Method and system for calculating transient power angle stability impact based on energy deviation rate

Abstract

The present invention discloses a method and system for calculating the impact of transient power angle stability based on energy deviation rate. The method includes determining the energy balance fulcrum of the power grid corresponding to each simulation time point after a preset disturbance; for each simulation time point, dividing the nodes in the network equation into two node sets with the energy balance fulcrum as the boundary, and determining the leading or lagging attribute of the node set based on the rotor inertia center angle of the synchronous generator associated with the nodes in the two node sets; determining the time period for calculating the impact of transient power angle stability based on the relative rotor inertia center angle between the leading node set and the lagging node set after the disturbance; calculating the energy deviation rate of the device connected to the grid through each node during the time period based on the grid-connected active power of the device through the connected nodes and the node set attributes at the simulation time point, and taking the sum of the energy deviation rates of the device connected to the grid through each connected node as the impact of the device on the transient power angle stability of the power grid after the disturbance. A corresponding system is also disclosed. The present invention provides an optimization direction for transient power angle stability control decision calculation.

Description

Technical Field

[0001] The present invention relates to a method and system for calculating transient power angle stability influence based on energy deviation rate, belonging to the technical field of power system stability analysis. Background Art

[0002] Calculating the impact of grid equipment on transient power angle stability is fundamental to optimizing transient power angle stability control decisions. The "Method for Identifying Component Participation Factors in Power System Transient Security and Stability Mode" (Patent No. ZL200910026801.4) uses the extended equal-area criterion to identify the dominant group of synchronous generators. Based on the proportional relationship between the kinetic energy of synchronous generators at critical moments during transient processes, this paper proposes a method for calculating the participation factors of synchronous generators influencing transient power angle stability. The "Method and System for Calculating the Impact of Non-synchronous Active Equipment on Transient Power Angle Stability" (Patent No. ZL 202011491158.5) also uses the extended equal-area criterion to identify the dominant group of synchronous generators. Based on the incremental energy injected into the grid at each connection point of non-synchronous active equipment during the initial period of the transient process, combined with oscillation center identification, this paper proposes a method for calculating the impact of non-synchronous active equipment (loads, renewable energy generation, DC systems, etc.) on transient power angle stability in the power grid. On the one hand, the existing technology lacks a method to calculate the impact of other equipment (AC lines, transformers, reactive equipment, etc.) on transient power angle stability. On the other hand, it does not consider that the impact of the energy injected into the grid by the equipment in each simulation step during the transient process on the transient power angle stability may be different. Summary of the Invention

[0003] The technical problem to be solved by the present invention is: there is a lack of a method for calculating the influence of multiple types of equipment such as AC lines, transformers, and reactive equipment on transient power angle stability.

[0004] A further technical problem to be solved by the present invention is to further consider the influence of the energy injected into the power grid by the equipment in each simulation step during the transient process on the transient power angle stability.

[0005] In order to solve the above technical problems, the technical solution adopted by the present invention is:

[0006] A method for calculating the transient power angle stability influence based on energy deviation rate comprises the following steps:

[0007] According to the preset grid operation state, determine the grid energy balance support corresponding to each simulation time point after the preset disturbance;

[0008] At each simulation time point, the nodes in the power grid network equation are divided into two node sets with the power grid energy balance support as the boundary, and two node sets corresponding to each simulation time point are obtained, wherein the nodes in the power grid network equation include the internal potential nodes of the synchronous generator;

[0009] At each simulation time point, by calculating the rotor inertia center angle of the synchronous generator associated with the nodes in the two node sets, the node set with the larger rotor inertia center angle is defined as the leading node set, and the node set with the smaller rotor inertia center angle is defined as the lagging node set, and the properties of the two node sets corresponding to each simulation time point are obtained;

[0010] Determine the time period for calculating the transient power angle stability influence according to the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point;

[0011] According to the grid-connected active power of the device through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the device connected to the grid through each connected node during the transient power angle stability impact calculation period is calculated respectively. The sum of the energy deviation rates of the device connected to the grid through each connected node is used as the impact of the device on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state.

[0012] According to the preset grid operation state, the grid energy balance fulcrum corresponding to each simulation time point after the preset disturbance is determined. The specific process is as follows:

[0013] Performing a time-domain simulation of a preset disturbance for a preset power grid operating state to obtain an admittance matrix for network equation calculation corresponding to each simulation time point, including the initial time point of the preset disturbance;

[0014] For the admittance matrix used for network equation calculation corresponding to each simulation time point, the power grid is equivalent to a 3-node 5-branch network that meets the set conditions through static network equivalence, where the nodes include two equivalent power nodes and one equivalent load node, and the network branches include three-node ground branches and branches between the two equivalent power nodes and the equivalent load node. The set condition is that the currents injected from the two equivalent power nodes to the equivalent load node are equal;

[0015] At each simulation time point, the position of the equivalent load node in the power grid before the equivalent process is determined by backtracking the static network equivalent process, and this position is used as the energy balance support of the power grid.

[0016] At each simulation time point, by calculating the central angle of inertia of the synchronous generator rotor associated with the nodes in the two node sets, the node set with the larger central angle of inertia of the rotor is defined as the leading node set, and the node set with the smaller central angle of inertia of the rotor is defined as the lagging node set. The properties of the two node sets corresponding to each simulation time point are obtained. The specific calculation formula of the central angle of inertia of the synchronous generator rotor associated with the nodes in the two node sets is:

[0017]

[0018]

[0019] Where A i 、B i are the two node sets corresponding to the i-th simulation time point, A i 、B i The rotor inertia center angle of the associated synchronous generator at the i-th simulation time point, M i.a , δ i.a A i The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node a of the synchronous generator at the i-th simulation time point, M i.b , δ i.b B i The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node b of the synchronous generator at the i-th simulation time point.

[0020] The node set attribute is referred to as "leading" or "lagging".

[0021] According to the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point, the time period for calculating the transient power angle stability influence is determined. The specific process is as follows;

[0022] If the A corresponding to the i-th simulation time point i If the leading node set is the set of nodes, then the relative rotor inertia center angle δ between the leading node set and the lagging node set corresponding to the i-th simulation time point is i Set to Otherwise, the relative rotor inertia center angle δ between the leading node set and the lagging node set corresponding to the i-th simulation time point is i Set to

[0023] Starting from the second simulation time point, the transient power angle stability impact calculation period is determined in the manner of increasing the simulation time points in sequence. Specifically, if δ i <δ s And δ i-1 <δ i , δ i >δ i+1 , then the period between the first simulation time point and the i-th simulation time point is used as the period for calculating the transient power angle stability influence. If δ i ≥δ s , then the period between the first simulation time point and the i-th simulation time point is used as the period for calculating the transient power angle stability influence, where δ s To set parameters;

[0024] In the process of determining the transient power angle stability influence calculation period, if the transient power angle stability influence calculation period is determined, the determination of the transient power angle stability influence calculation period is terminated.

[0025] Based on the grid-connected active power of the device through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the device connected to the grid through each connected node during the transient power angle stability impact calculation period is calculated respectively. The sum of the energy deviation rates of the device connected to the grid through each connected node is used as the impact of the device on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state. The specific calculation formula is:

[0026]

[0027]

[0028] Where B D The set of nodes that inject device D into the grid, D tas D is the influence of device D on the transient power angle stability of the grid after the preset disturbance in the preset grid operation state. j For device D via B D Energy deviation rate s of node j connected to the grid i.j Is positive or negative, if the ith simulation time point B D If the middle node j belongs to the preceding node set, the negative sign is taken, otherwise, s i.j Take the positive sign, I is the last simulation time point corresponding to the transient power angle stability influence calculation period, P i.j is the i-th simulation time point B D The grid-connected active power of node j, t i is the transient process moment corresponding to the i-th simulation time point.

[0029] The equipment mentioned herein refers to the primary equipment in the power grid except the synchronous generator.

[0030] In view of the situation that the power grid after the preset disturbance in the preset power grid operation state is composed of two or more asynchronous operation sub-grids, each asynchronous operation sub-grid is processed independently.

[0031] The transient power angle stability impact calculation system based on energy deviation rate includes:

[0032] Energy balance fulcrum determination module: for the preset grid operation state, determines the grid energy balance fulcrum corresponding to each simulation time point after the preset disturbance;

[0033] Node set generation module: For each simulation time point, the nodes in the power grid network equation, including the internal potential nodes of the synchronous generator, are divided into two node sets with the power grid energy balance support as the boundary, and the two node sets corresponding to each simulation time point are obtained;

[0034] Node set attribute determination module: For each simulation time point, by calculating the rotor inertia center angle of the synchronous generator associated with the nodes in the two node sets, the node set with the larger rotor inertia center angle is defined as the leading node set, and the node set with the smaller rotor inertia center angle is defined as the lagging node set, and the attributes of the two node sets corresponding to each simulation time point are obtained;

[0035] Impact calculation period determination module: determines the period for transient power angle stability impact calculation based on the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point;

[0036] Equipment impact calculation module: Based on the grid-connected active power of the equipment through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the equipment connected to the grid through each connected node during the period of transient power angle stability impact calculation is calculated respectively. The sum of the energy deviation rates of the equipment connected to the grid through each connected node is used as the impact of the equipment on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state.

[0037] A computer-readable storage medium storing one or more programs, wherein the one or more programs include instructions, which, when executed by a computing device, enable the computing device to perform a method for calculating a transient power angle stability influence based on an energy deviation rate.

[0038] A computing device includes one or more processors, one or more memories, and one or more programs, wherein the one or more programs are stored in the one or more memories and are configured to be executed by the one or more processors, and the one or more programs include instructions for executing a method for calculating the transient power angle stability influence based on an energy deviation rate.

[0039] The beneficial effects achieved by the present invention are as follows: the present invention realizes the quantitative evaluation of the impact of all primary equipment in the power grid except synchronous generators on transient power angle stability, provides a decision-making basis for the collaborative participation of large-scale multi-type equipment including new energy power generation stations, loads, DC systems, reactive compensation equipment, and even AC lines, transformers, etc. in transient power angle stability control, and can improve the computational efficiency of optimization decisions. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 Flowchart of the method of the present invention. DETAILED DESCRIPTION

[0041] The present invention will be further described below in conjunction with the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solutions of the present invention and are not intended to limit the scope of protection of the present invention.

[0042] like Figure 1 As shown, the method for calculating the transient power angle stability impact based on the energy deviation rate includes the following steps:

[0043] Step 1: for a preset power grid operating state, determine the power grid energy balance fulcrum corresponding to each simulation time point after a preset disturbance.

[0044] The specific process is as follows:

[0045] 11) performing a time domain simulation of a preset disturbance for a preset power grid operating state, and obtaining an admittance matrix for network equation calculation corresponding to each simulation time point including the initial time point of the preset disturbance;

[0046] 12) For the admittance matrix used for network equation calculation corresponding to each simulation time point, the power grid is equivalent to a 3-node 5-branch network that meets the set conditions through static network equivalent, where the nodes include two equivalent power nodes and one equivalent load node, and the network branches include three ground branches and branches between the two equivalent power nodes and the equivalent load node. The set condition is that the currents injected from the two equivalent power nodes to the equivalent load node are equal;

[0047] 13) For each simulation time point, the position of the equivalent load node in the power grid before the equivalent process is determined by backtracking the static network equivalent process, and this position is used as the energy balance support point of the power grid.

[0048] The grid energy balance fulcrum changes dynamically during the transient process, and each simulation time point has a corresponding grid energy balance fulcrum.

[0049] Step 2: For each simulation time point, the nodes in the power grid network equation including the internal potential nodes of the synchronous generator are divided into two node sets with the power grid energy balance support as the boundary, and the two node sets corresponding to each simulation time point are obtained.

[0050] Step 3: For each simulation time point, by calculating the central angle of inertia of the synchronous generator rotor associated with the nodes in the two node sets, the node set with the larger central angle of inertia of the rotor is defined as the leading node set, and the node set with the smaller central angle of inertia of the rotor is defined as the lagging node set, and the properties of the two node sets corresponding to each simulation time point are obtained.

[0051] The specific formula for calculating the central angle of inertia of the synchronous generator rotor associated with the nodes in the two node sets is:

[0052]

[0053]

[0054] Where A i 、B i are the two node sets corresponding to the i-th simulation time point, A i 、Bi The rotor inertia center angle of the associated synchronous generator at the i-th simulation time point, M i.a , δ i.a A i The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node a of the synchronous generator at the i-th simulation time point, M i.b , δ i.b B i The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node b of the synchronous generator at the i-th simulation time point.

[0055] The node set attribute is referred to as "leading" or "lagging".

[0056] Step 4: Determine the time period for calculating the transient power angle stability influence according to the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point.

[0057] The specific process is:

[0058] 41) If the A corresponding to the i-th simulation time point i If the leading node set is the set of nodes, then the relative rotor inertia center angle δ between the leading node set and the lagging node set corresponding to the i-th simulation time point is i Set to Otherwise, the relative rotor inertia center angle δ between the leading node set and the lagging node set corresponding to the i-th simulation time point is i Set to

[0059] 42) Starting from the second simulation time point, the transient power angle stability impact calculation period is determined in the manner of increasing the simulation time points. Specifically, if δ i <δ s And δ i-1 <δ i , δ i >δ i+1 , then the period between the first simulation time point and the i-th simulation time point is used as the period for calculating the transient power angle stability influence. If δ i ≥δ s , then the period between the first simulation time point and the i-th simulation time point is used as the period for calculating the transient power angle stability influence, where δ s To set the parameters, it is usually set to 180°;

[0060] In the process of determining the transient power angle stability influence calculation period, if the transient power angle stability influence calculation period is determined, the determination of the transient power angle stability influence calculation period is terminated.

[0061] Step 5: Based on the grid-connected active power of the device through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the device connected to the grid through each connected node during the transient power angle stability impact calculation period is calculated respectively, and the sum of the energy deviation rates of the device connected to the grid through each connected node is used as the impact of the device on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state.

[0062] The specific formula is;

[0063]

[0064]

[0065] Where B D The set of nodes that inject device D into the grid, D tas D is the influence of device D on the transient power angle stability of the grid after the preset disturbance in the preset grid operation state. j For device D via B D Energy deviation rate of node j connected to the grid, s i.j Is positive or negative, if the ith simulation time point B D If the middle node j belongs to the preceding node set, the negative sign is taken, otherwise, s i.j Take the positive sign, I is the last simulation time point corresponding to the transient power angle stability influence calculation period, P i.j is the i-th simulation time point B D The grid-connected active power of node j, t i is the transient process moment corresponding to the i-th simulation time point.

[0066] The equipment mentioned herein refers to the primary equipment in the power grid except the synchronous generator.

[0067] In view of the situation that the power grid after the preset disturbance in the preset power grid operation state is composed of two or more asynchronous operation sub-grids, each asynchronous operation sub-grid is processed independently.

[0068] If a device is injected into multiple asynchronous sub-grids, the impact of the device on the transient power angle stability includes its impact on the transient power angle stability of each asynchronous sub-grid.

[0069] Through the above method, if the obtained influence of the device on the transient power angle stability of the power grid is greater than 0, it means that the device is beneficial to the transient power angle stability, and the larger the value, the more beneficial it is. If the obtained influence of the device on the transient power angle stability of the power grid is less than 0, it means that the device is not conducive to the transient power angle stability, and the smaller the value, the more unfavorable it is. If the obtained influence of the device on the transient power angle stability of the power grid is equal to 0, it means that the device has no influence on the transient power angle stability of the power grid.

[0070] The theoretical basis of the above method is: based on the law of conservation of energy, the electric energy of the power grid is conserved in any period of time. If the original synchronous operation power grid is divided into two synchronous operation sub-grids with the grid energy balance fulcrum corresponding to each simulation time point in the transient process as the boundary, and the ground branch of the grid energy balance fulcrum at each simulation time point in the transient process is equivalent to two parallel branches with equal impedance, and respectively included in the two synchronous operation sub-grids, then the electric energy of the two synchronous operation sub-grids is conserved. Therefore, the energy deviation injected by the equipment into the synchronous operation sub-grid will inevitably affect the energy injected into the grid by the synchronous generator in the synchronous operation sub-grid, thereby affecting the transient power angle stability of the original synchronous operation power grid, and the energy deviation injected by the equipment into the two synchronous operation sub-grids has completely opposite effects on the transient power angle stability of the original synchronous operation power grid.

[0071] The above method realizes the quantitative evaluation of the impact of all primary equipment in the power grid except synchronous generators on transient power angle stability, provides a decision-making basis for the collaborative participation of large-scale multi-type equipment including renewable energy power generation stations, loads, DC systems, reactive compensation equipment, and even AC lines and transformers in transient power angle stability control, and can improve the computational efficiency of optimization decision-making.

[0072] The transient power angle stability impact calculation system based on energy deviation rate includes:

[0073] Energy balance fulcrum determination module: for the preset grid operation state, determines the grid energy balance fulcrum corresponding to each simulation time point after the preset disturbance;

[0074] Node set generation module: For each simulation time point, the nodes in the power grid network equation, including the internal potential nodes of the synchronous generator, are divided into two node sets with the power grid energy balance support as the boundary, and the two node sets corresponding to each simulation time point are obtained;

[0075] Node set attribute determination module: For each simulation time point, by calculating the rotor inertia center angle of the synchronous generator associated with the nodes in the two node sets, the node set with the larger rotor inertia center angle is defined as the leading node set, and the node set with the smaller rotor inertia center angle is defined as the lagging node set, and the attributes of the two node sets corresponding to each simulation time point are obtained;

[0076] Impact calculation period determination module: determines the period for transient power angle stability impact calculation based on the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point;

[0077] Equipment impact calculation module: Based on the grid-connected active power of the equipment through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the equipment connected to the grid through each connected node during the period of transient power angle stability impact calculation is calculated respectively. The sum of the energy deviation rates of the equipment connected to the grid through each connected node is used as the impact of the equipment on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state.

[0078] A computer-readable storage medium storing one or more programs, wherein the one or more programs include instructions, which, when executed by a computing device, enable the computing device to perform a method for calculating a transient power angle stability influence based on an energy deviation rate.

[0079] A computing device includes one or more processors, one or more memories, and one or more programs, wherein the one or more programs are stored in the one or more memories and are configured to be executed by the one or more processors, and the one or more programs include instructions for executing a method for calculating the transient power angle stability influence based on an energy deviation rate.

[0080] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0081] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the steps in the process. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0082] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 a process or multiple processes and / or boxes Figure 1The function specified in one or more boxes.

[0083] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operational steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing the instructions executed on the computer or other programmable device for implementing the process. Figure 1 a process or multiple processes and / or boxes Figure 1 The steps for the function specified in one or more boxes.

[0084] The above are merely embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention are included in the scope of the claims of the present invention to be approved.

Claims

1. A method for calculating the influence of transient power angle stability based on energy deviation rate, characterized in that: The following steps are involved: Based on the preset grid operation state, the grid energy balance fulcrum corresponding to each simulation time point after the preset disturbance is determined. The specific process is to perform a time domain simulation of the preset disturbance based on the preset grid operation state, and obtain the admittance matrix corresponding to each simulation time point including the initial time of the preset disturbance for network equation calculation; For the admittance matrix used for network equation calculation corresponding to each simulation time point, the power grid is equivalent to a 3-node 5-branch network that meets the set conditions through static network equivalence, where the nodes include two equivalent power nodes and one equivalent load node, and the network branches include three-node ground branches and branches between the two equivalent power nodes and the equivalent load node. The set condition is that the currents injected from the two equivalent power nodes to the equivalent load node are equal; At each simulation time point, the position of the equivalent load node in the power grid before the equivalent process is determined by backtracking the static network equivalent process, and this position is used as the energy balance support point of the power grid; At each simulation time point, the nodes in the power grid network equation are divided into two node sets with the power grid energy balance support as the boundary, and two node sets corresponding to each simulation time point are obtained; the nodes in the power grid network equation include the internal potential nodes of the synchronous generator; At each simulation time point, by calculating the rotor inertia center angle of the synchronous generator associated with the nodes in the two node sets, the node set with the larger rotor inertia center angle is defined as the leading node set, and the node set with the smaller rotor inertia center angle is defined as the lagging node set, and the properties of the two node sets corresponding to each simulation time point are obtained; According to the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point, the time period for calculating the transient power angle stability influence is determined. The specific process is as follows: if the i-th simulation time point corresponds to If the leading node set is the set of nodes, then the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to the i-th simulation time point is Set to Otherwise, the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to the i-th simulation time point is Set to ; Starting from the second simulation time point, the transient power angle stability impact calculation period is determined in the manner of increasing the simulation time points. Specifically: If and 、 , then the first simulation time point and the The period between the simulation time points is used as the period for calculating the transient power angle stability impact. , then the first simulation time point and the The period between the simulation time points is used as the period for calculating the transient power angle stability impact, where: To set parameters; in the process of determining the transient power angle stability influence calculation period, if the transient power angle stability influence calculation period is determined, then the determination of the transient power angle stability influence calculation period is terminated; 、 are the two node sets corresponding to the i-th simulation time point, 、 They are 、 The rotor inertia center angle of the associated synchronous generator at the ith simulation time point; According to the grid-connected active power of the device through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the device connected to the grid through each connected node during the transient power angle stability impact calculation period is calculated respectively. The sum of the energy deviation rates of the device connected to the grid through each connected node is used as the impact of the device on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state.

2. The method for calculating the transient power angle stability impact based on the energy deviation rate according to claim 1, characterized in that: At each simulation time point, by calculating the central angle of inertia of the synchronous generator rotor associated with the nodes in the two node sets, the node set with the larger central angle of inertia of the rotor is defined as the leading node set, and the node set with the smaller central angle of inertia of the rotor is defined as the lagging node set. The properties of the two node sets corresponding to each simulation time point are obtained. The specific calculation formula of the central angle of inertia of the synchronous generator rotor associated with the nodes in the two node sets is: ， ， Where, 、 They are The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node a of the synchronous generator at the i-th simulation time point, 、 They are The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node b of the synchronous generator at the i-th simulation time point; The node set attribute is referred to as "leading" or "lagging".

3. The method for calculating the transient power angle stability impact based on the energy deviation rate according to claim 1, characterized in that: Based on the grid-connected active power of the device through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the device connected to the grid through each connected node during the transient power angle stability impact calculation period is calculated respectively. The sum of the energy deviation rates of the device connected to the grid through each connected node is used as the impact of the device on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state. The specific calculation formula is: ， ， Where, The set of nodes that inject the grid for device D, is the influence of device D on the transient power angle stability of the grid after the preset disturbance in the preset grid operation state, For equipment D Energy deviation rate of grid-connected node j, Is positive or negative, if the ith simulation time point If the middle node j belongs to the preceding node set, the negative sign is taken, otherwise, Take the positive sign, I is the last simulation time point corresponding to the transient power angle stability influence calculation period, is the i-th simulation time point The grid-connected active power of node j is is the transient process moment corresponding to the i-th simulation time point.

4. The method for calculating the transient power angle stability impact based on the energy deviation rate according to claim 1, characterized in that: In view of the situation that the power grid after the preset disturbance in the preset power grid operation state is composed of two or more asynchronous operation sub-grids, each asynchronous operation sub-grid is processed independently.

5. A transient power angle stability impact calculation system based on energy deviation rate, characterized by: include, Energy balance fulcrum determination module: Based on the preset grid operation state, it determines the grid energy balance fulcrum corresponding to each simulation time point after the preset disturbance. The specific process is to perform a time domain simulation of the preset disturbance based on the preset grid operation state, and obtain the admittance matrix corresponding to each simulation time point including the initial moment of the preset disturbance for network equation calculation; For the admittance matrix used for network equation calculation corresponding to each simulation time point, the power grid is equivalent to a 3-node 5-branch network that meets the set conditions through static network equivalence, where the nodes include two equivalent power nodes and one equivalent load node, and the network branches include three-node ground branches and branches between the two equivalent power nodes and the equivalent load node. The set condition is that the currents injected from the two equivalent power nodes to the equivalent load node are equal; At each simulation time point, the position of the equivalent load node in the power grid before the equivalent process is determined by backtracking the static network equivalent process, and this position is used as the energy balance support point of the power grid; Node set generation module: For each simulation time point, the nodes in the power grid network equation, including the internal potential nodes of the synchronous generator, are divided into two node sets with the power grid energy balance support as the boundary, and the two node sets corresponding to each simulation time point are obtained; Node set attribute determination module: For each simulation time point, by calculating the rotor inertia center angle of the synchronous generator associated with the nodes in the two node sets, the node set with the larger rotor inertia center angle is defined as the leading node set, and the node set with the smaller rotor inertia center angle is defined as the lagging node set, and the attributes of the two node sets corresponding to each simulation time point are obtained; Influence calculation period determination module: According to the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point, the transient power angle stability influence calculation period is determined. The specific process is as follows: if the i-th simulation time point corresponds to If the leading node set is the set of nodes, then the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to the i-th simulation time point is Set to Otherwise, the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to the i-th simulation time point is Set to ; Starting from the second simulation time point, the transient power angle stability impact calculation period is determined in the manner of increasing the simulation time points. Specifically: If and 、 , then the first simulation time point and the The period between the simulation time points is used as the period for calculating the transient power angle stability impact. , then the first simulation time point and the The period between the simulation time points is used as the period for calculating the transient power angle stability impact, where: To set parameters; in the process of determining the transient power angle stability influence calculation period, if the transient power angle stability influence calculation period is determined, then the determination of the transient power angle stability influence calculation period is terminated; Equipment impact calculation module: Based on the grid-connected active power of the equipment through each connected node and the node set attributes at the simulation time point, the energy deviation rate of the equipment connected to the grid through each connected node during the period of transient power angle stability impact calculation is calculated respectively. The sum of the energy deviation rates of the equipment connected to the grid through each connected node is used as the impact of the equipment on the transient power angle stability of the grid after the preset disturbance under the preset grid operating state.

6. The transient power angle stability impact calculation system based on energy deviation rate according to claim 5, characterized in that: In the node set attribute determination module, the specific calculation formula for the central angle of inertia of the synchronous generator rotor associated with the nodes in the two node sets is: ， ， Where, 、 are the two node sets corresponding to the i-th simulation time point, 、 They are 、 The rotor inertia center angle of the synchronous generator at the i-th simulation time point is, 、 They are The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node a of the synchronous generator at the i-th simulation time point, 、 They are The moment of inertia and internal potential phase angle of the synchronous generator corresponding to the internal potential node b of the synchronous generator at the i-th simulation time point.

7. A computer-readable storage medium storing one or more programs, characterized in that: The one or more programs include instructions that, when executed by a computing device, cause the computing device to perform any one of the methods according to claims 1 to 4 .

8. A computing device, characterized in that: include, One or more processors, one or more memories, and one or more programs, wherein the one or more programs are stored in the one or more memories and configured to be executed by the one or more processors, the one or more programs comprising instructions for performing any of the methods according to claims 1 to 4.

Images