CN116186979A - Transient power angle stability influence degree calculation method and system based on energy deviation rate - Google Patents
Transient power angle stability influence degree calculation method and system based on energy deviation rate Download PDFInfo
- Publication number
- CN116186979A CN116186979A CN202211568869.7A CN202211568869A CN116186979A CN 116186979 A CN116186979 A CN 116186979A CN 202211568869 A CN202211568869 A CN 202211568869A CN 116186979 A CN116186979 A CN 116186979A
- Authority
- CN
- China
- Prior art keywords
- node
- time point
- simulation time
- node set
- power grid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Power Engineering (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a transient power angle stability influence degree calculation method and a transient power angle stability influence degree calculation system based on energy deviation rate, wherein the method comprises the steps of determining power grid energy balance supporting points corresponding to simulation time points after preset disturbance; for each simulation time point, dividing nodes in a network equation into two node sets by taking an energy balance fulcrum as a boundary, and determining the leading or lagging attribute of the node sets according to the inertia central angle of the synchronous generator rotor associated with the nodes in the two node sets; determining a time period for calculating the stability influence degree of the transient power angle according to the relative rotor inertia center angle between the front node set and the hysteresis node set after disturbance; according to the grid-connected active power of the simulation time point equipment through the connected nodes and the node set attribute, calculating the energy deviation rate of the equipment which is connected through each node in the period, and taking the sum of the energy deviation rates of the equipment which is connected through each node as the influence of the equipment on the stability of the transient power angle of the power grid after the disturbance. Corresponding systems are also disclosed. The invention provides an optimization direction for transient power angle stability control decision calculation.
Description
Technical Field
The invention relates to a transient power angle stability influence degree calculation method and system based on energy deviation rate, and belongs to the technical field of power system stability analysis.
Background
The influence degree calculation of the power grid equipment on the transient power angle stability is a basis for realizing the transient power angle stability control optimization decision. The method for identifying element participation factors in a transient safety and stability mode of an electric power system (patent number: ZL 200910026801.4) is based on synchronous generator dominant group identification of area criteria such as expansion, a participation factor calculation method for influencing transient power angle stability of a synchronous generator is provided according to a proportional relation between the kinetic energy of a synchronous generator at key moments in a transient process, a transient power angle stability influence degree calculation method and a system of asynchronous active equipment (patent number: ZL 202011491158.5) are also based on synchronous generator dominant group identification of area criteria such as expansion, and a transient power angle stability influence degree calculation method for asynchronous active equipment (load, new energy power generation, direct current system and the like) except the synchronous generator in a power grid is provided according to increment of energy injected into a power grid in an initial period of the asynchronous active equipment in the transient process and combined with oscillation center identification. On one hand, the prior art lacks a method for calculating the influence degree of other equipment (alternating current lines, transformers, reactive power equipment and the like) on the transient state power angle stability, and on the other hand, the influence of the energy injected into the power grid by the equipment in each simulation step length in the transient state process on the transient state power angle stability is not considered to be different.
Disclosure of Invention
The technical problems to be solved by the invention are as follows: the method for calculating the influence degree of various equipment such as an alternating current circuit, a transformer, reactive equipment and the like on the transient power angle stability is lacked.
The invention further aims to solve the technical problems that: further consider the influence of the energy of the power grid injected by the equipment in each simulation step in the transient process on the transient power angle stability.
In order to solve the technical problems, the invention adopts the following technical scheme:
a transient power angle stability influence degree calculation method based on energy deviation rate comprises the following steps:
determining power grid energy balance supporting points corresponding to simulation time points after preset disturbance according to a preset power grid running state;
aiming at each simulation time point, dividing nodes in a power grid network equation into two node sets by taking a power grid energy balance fulcrum as a boundary, and obtaining two node sets corresponding to each simulation time point, wherein the nodes in the power grid network equation comprise potential nodes in a synchronous generator;
aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator which are related to the nodes in the two node sets, defining a node set with a large central angle of the rotor inertia as a leading node set, defining a node set with a small central angle of the rotor inertia as a lagging node set, and obtaining the attributes of the two node sets corresponding to each simulation time point;
according to the relative rotor inertia center angles between the leading node set and the lagging node set corresponding to each simulation time point, determining the time period for calculating the transient power angle stability influence degree;
according to the grid-connected active power of the equipment at the simulation time point and the node set attribute of each connected node, respectively calculating the energy deviation rate of the equipment which is connected through each connected node in the transient power angle stability influence calculation period, and taking the sum of the energy deviation rates of the equipment which is connected through each connected node as the influence of the equipment on the transient power angle stability of the power grid after the disturbance is preset in the preset power grid running state.
Aiming at the running state of a preset power grid, determining power grid energy balance supporting points corresponding to each simulation time point after the preset disturbance, wherein the specific process is that,
performing time domain simulation of preset disturbance on a preset power grid running state to obtain an admittance matrix for network equation calculation corresponding to each simulation time point including a preset disturbance initial moment;
aiming at admittance matrixes for network equation calculation corresponding to simulation time points, respectively equating a power grid into a 3-node 5-branch network meeting set conditions through static network equivalence, wherein the nodes comprise two equivalent power supply nodes and 1 equivalent load node, the network branch comprises a ground branch of 3 nodes and a branch between the two equivalent power supply nodes and the equivalent load node, and the set conditions are that currents injected into the equivalent load node by the two equivalent power supply nodes are equal;
and for each simulation time point, determining the position of the equivalent load node in the power grid before the equivalent through backtracking of the equivalent process of the static network, and taking the position as an energy balance pivot of the power grid.
Aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator associated with the nodes in the two node sets, defining the node set with the large central angle of the rotor inertia as a leading node set, defining the node set with the small central angle of the rotor inertia as a lagging node set, and obtaining the attribute of the two node sets corresponding to each simulation time point, wherein the specific calculation formula of the central angles of the rotor inertia of the synchronous generator associated with the nodes in the two node sets is as follows:
wherein A is i 、B i Respectively two node sets corresponding to the ith simulation time point,
respectively A i 、B i Ith simulation time point rotor inertia center angle of medium-association synchronous generator, M i.a 、δ i.a Respectively A i Moment of inertia and internal potential phase angle, M, of the ith simulation time point of the synchronous generator corresponding to internal potential node a of the medium synchronous generator i.b 、δ i.b Respectively B i And the moment of inertia and the internal potential phase angle of the ith simulation time point of the synchronous generator corresponding to the internal potential node b of the medium synchronous generator.
The node set attribute is referred to as "lead" or "lag".
According to the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point, determining the time period for transient power angle stability influence degree calculation, wherein the specific process is as follows;
if A corresponding to the ith simulation time point i For the leading node set, the relative rotor inertia center angle delta between the leading node set and the lagging node set corresponding to the ith simulation time point i Is arranged as
Otherwise, the relative rotor inertia center angle delta between the leading node set and the lagging node set corresponding to the ith simulation time point i Set to->
Starting from the 2 nd simulation time point, determining a transient power angle stability influence degree calculation period in a mode that the simulation time points are sequentially increased, wherein the method specifically comprises the following steps: if delta i <δ s And delta i-1 <δ i 、δ i >δ i+1 Taking the time period between the 1 st simulation time point and the i-th simulation time point as the time period for calculating the transient power angle stability influence degree, if delta i ≥δ s Taking the time period between the 1 st simulation time point and the i-th simulation time point as the time period for calculating the transient power angle stability influence degree, wherein delta s Setting 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, the determination of the transient power angle stability influence calculation period is terminated.
According to the grid-connected active power of the equipment at the simulation time point and the node set attribute of each connected node, respectively calculating the energy deviation rate of the equipment connected through each connected node in the transient power angle stability influence calculation period, taking the sum of the energy deviation rates of the equipment connected through each connected node as the influence of the equipment on the transient power angle stability of the power grid after the equipment is subjected to preset disturbance in the preset power grid running state, wherein a specific calculation formula is as follows;
wherein B is D Injecting a device D into a node set of the power grid, D tas The influence of the equipment D on the stability of the transient power angle of the power grid after the disturbance is preset under the preset power grid running state is shown as D j For apparatus D via B D Energy deviation rate s of grid connection of middle node j i.j For positive or negative sign, if the ith simulation time point B D The middle node j belongs to the leading node set, and takes a negative sign, otherwise, s i.j Taking positive sign, wherein I is the last simulation time point corresponding to the transient power angle stability influence calculation period, and P i.j For the ith simulation time point B D Grid-connected active power of middle node j, t i The transient process time corresponding to the ith simulation time point.
Wherein the device refers to a primary device in the power grid except for the synchronous generator.
Aiming at the situation that the power grid is composed of two or more asynchronous operation sub-power grids after the preset disturbance in the operation state of the preset power grid, each asynchronous operation sub-power grid is independently processed.
The transient power angle stability influence degree calculation system based on the energy deviation rate comprises,
the energy balance pivot point determining module: determining power grid energy balance supporting points corresponding to simulation time points after preset disturbance according to a preset power grid running state;
the node set generation module: aiming at each simulation time point, dividing nodes in a power grid network equation containing potential nodes in the synchronous generator into two node sets by taking a power grid energy balance fulcrum as a boundary, and obtaining two node sets corresponding to each simulation time point;
the node set attribute determining module: aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator which are related to the nodes in the two node sets, defining a node set with a large central angle of the rotor inertia as a leading node set, defining a node set with a small central angle of the rotor inertia as a lagging node set, and obtaining the attributes of the two node sets corresponding to each simulation time point;
influence degree calculation period determination module: according to the relative rotor inertia center angles between the leading node set and the lagging node set corresponding to each simulation time point, determining the time period for calculating the transient power angle stability influence degree;
and the equipment influence degree calculating module is used for: according to the grid-connected active power of the equipment at the simulation time point and the node set attribute of each connected node, respectively calculating the energy deviation rate of the equipment connected through each connected node in the period of transient power angle stability influence calculation, and taking the sum of the energy deviation rates of the equipment connected through each connected node as the influence of the equipment on the transient power angle stability of the power grid after the disturbance in the preset power grid running state.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform a transient power angle stability influence calculation method based on an energy deviation rate.
A computing device comprising one or more processors, one or more memories, and one or more programs, wherein 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 a transient power angle stability influence calculation method based on an energy deviation rate.
The invention has the beneficial effects that: the invention realizes the quantitative evaluation of the transient power angle stability influenced by all primary equipment except the synchronous generator in the power grid, provides decision basis for large-scale multi-class equipment including new energy power generation stations, loads, direct current systems, reactive compensation equipment, even alternating current lines, transformers and the like to cooperatively participate in the transient power angle stability control, and can improve the calculation efficiency of optimization decisions.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for more clearly illustrating the technical aspects of the present invention, and are not intended to limit the scope of the present invention.
As shown in fig. 1, the transient power angle stability influence degree calculation method based on the energy deviation rate comprises the following steps:
step 1, determining power grid energy balance supporting points corresponding to simulation time points after preset disturbance according to a preset power grid running state.
The specific process is as follows:
11 Performing time domain simulation of preset disturbance on a preset power grid running state to obtain an admittance matrix for network equation calculation corresponding to each simulation time point including a preset disturbance initial moment;
12 Aiming at admittance matrixes for network equation calculation corresponding to simulation time points, respectively equating a power grid into a 3-node 5-branch network meeting set conditions through static network equivalence, wherein the nodes comprise two equivalent power supply nodes and 1 equivalent load node, the network branch comprises a ground branch of 3 nodes and a branch between the two equivalent power supply nodes and the equivalent load node, and the set conditions are that currents injected into the equivalent load node by the two equivalent power supply nodes are equal;
13 For each simulation time point, determining the position of the equivalent load node in the power grid before the equivalent through backtracking of the equivalent process of the static network, and taking the position as an energy balance pivot of the power grid.
The power grid energy balance pivot point is dynamically changed in the transient process, and each simulation time point is provided with a corresponding power grid energy balance pivot point.
And 2, dividing nodes in a power grid network equation containing potential nodes in the synchronous generator into two node sets by taking a power grid energy balance fulcrum as a boundary aiming at each simulation time point, and obtaining the two node sets corresponding to each simulation time point.
And 3, aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator associated with the nodes in the two node sets, defining the node set with the large central angle of the rotor inertia as a leading node set, defining the node set with the small central angle of the rotor inertia as a lagging node set, and obtaining the attributes of the two node sets corresponding to each simulation time point.
The specific formula for calculating the inertia center angle of the synchronous generator rotor associated with the two nodes is as follows:
wherein A is i 、B i Respectively two node sets corresponding to the ith simulation time point,
respectively A i 、B i Ith simulation time point rotor inertia center angle of medium-association synchronous generator, M i.a 、δ i.a Respectively A i Moment of inertia and internal potential phase angle, M, of the ith simulation time point of the synchronous generator corresponding to internal potential node a of the medium synchronous generator i.b 、δ i.b Respectively B i The same corresponding to the internal potential node b of the medium synchronous generatorThe i-th simulation of the step generator simulates the moment of inertia and the internal potential phase angle.
The node set attribute is referred to as "lead" or "lag".
And 4, determining a time period for calculating the stability influence degree of the transient power angle 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 specific process is as follows:
41 If A corresponding to the ith simulation time point i For the leading node set, the relative rotor inertia center angle delta between the leading node set and the lagging node set corresponding to the ith simulation time point i Is arranged as
Otherwise, the relative rotor inertia center angle delta between the leading node set and the lagging node set corresponding to the ith simulation time point i Set to->
42 From the 2 nd simulation time point, determining a transient power angle stability influence degree calculation period in a mode that the simulation time points are sequentially increased, specifically: if delta i <δ s And delta i-1 <δ i 、δ i >δ i+1 Taking the time period between the 1 st simulation time point and the i-th simulation time point as the time period for calculating the transient power angle stability influence degree, if delta i ≥δ s Taking the time period between the 1 st simulation time point and the i-th simulation time point as the time period for calculating the transient power angle stability influence degree, wherein delta s For setting parameters, it is usually set to 180 °;
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.
And 5, respectively calculating the energy deviation rate of the equipment connected with each node in the transient power angle stability influence calculation period according to the grid-connected active power and the node set attribute of the equipment connected with each node at the simulation time point, and taking the sum of the energy deviation rates of the equipment connected with each node in the grid connection as the influence of the equipment on the transient power angle stability of the power grid after the disturbance is preset in the preset power grid running state.
The specific formula is as follows;
wherein B is D Injecting a device D into a node set of the power grid, D tas The influence of the equipment D on the stability of the transient power angle of the power grid after the disturbance is preset under the preset power grid running state is shown as D j For apparatus D via B D Energy deviation rate s of grid connection of middle node j i.j For positive or negative sign, if the ith simulation time point B D The middle node j belongs to the leading node set, and takes a negative sign, otherwise, s i.j Taking positive sign, wherein I is the last simulation time point corresponding to the transient power angle stability influence calculation period, and P i.j For the ith simulation time point B D Grid-connected active power of middle node j, t i The transient process time corresponding to the ith simulation time point.
Wherein the device refers to a primary device in the power grid except for the synchronous generator.
Aiming at the situation that the power grid is composed of two or more asynchronous operation sub-power grids after the preset disturbance in the operation state of the preset power grid, each asynchronous operation sub-power grid is independently processed.
If the device is injected into a plurality of asynchronously operating sub-grids, the influence degree of the device on the transient power angle stability comprises the influence degree of the device on the transient power angle stability of each asynchronously operating sub-grid.
By the method, if the obtained influence degree of the equipment on the transient power angle stability of the power grid is larger than 0, the equipment is favorable for the transient power angle stability, the larger the numerical value is, the more favorable, if the obtained influence degree of the equipment on the transient power angle stability of the power grid is smaller than 0, the equipment is unfavorable for the transient power angle stability, the smaller the numerical value is, the more unfavorable, and if the obtained influence degree of the equipment on the transient power angle stability of the power grid is equal to 0, the equipment is not influenced for the transient power angle stability of the power grid.
The theoretical basis of the method is as follows: the grid is conservation of electrical energy over any period of time based on the law of conservation of energy. If the original synchronous operation power grid is divided into two synchronous operation sub-power grids by taking the power grid energy balance pivot corresponding to each simulation time point in the transient process as a boundary, and the ground support of the power grid energy balance pivot of each simulation time point in the transient process is equivalent to two parallel branches with equal impedance, and the parallel branches are respectively incorporated into the two synchronous operation sub-power grids, the electric energy of the two synchronous operation sub-power grids is conserved. Therefore, the energy deviation of the device injected into the synchronous operation sub-power grid can necessarily influence the energy of the synchronous generator in the synchronous operation sub-power grid injected into the power grid, so that the transient power angle stability of the original synchronous operation power grid is influenced, and the influence of the energy deviation of the device injected into the two synchronous operation sub-power grids on the transient power angle stability of the original synchronous operation power grid is completely opposite.
The method realizes the quantitative evaluation of the transient power angle stability influenced by all primary equipment except the synchronous generator in the power grid, provides decision basis for the cooperative participation of large-scale multi-class equipment including new energy power generation stations, loads, direct current systems, reactive compensation equipment, even alternating current circuits, transformers and the like in the transient power angle stability control, and can improve the calculation efficiency of optimization decisions.
The transient power angle stability influence degree calculation system based on the energy deviation rate comprises,
the energy balance pivot point determining module: determining power grid energy balance supporting points corresponding to simulation time points after preset disturbance according to a preset power grid running state;
the node set generation module: aiming at each simulation time point, dividing nodes in a power grid network equation containing potential nodes in the synchronous generator into two node sets by taking a power grid energy balance fulcrum as a boundary, and obtaining two node sets corresponding to each simulation time point;
the node set attribute determining module: aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator which are related to the nodes in the two node sets, defining a node set with a large central angle of the rotor inertia as a leading node set, defining a node set with a small central angle of the rotor inertia as a lagging node set, and obtaining the attributes of the two node sets corresponding to each simulation time point;
influence degree calculation period determination module: according to the relative rotor inertia center angles between the leading node set and the lagging node set corresponding to each simulation time point, determining the time period for calculating the transient power angle stability influence degree;
and the equipment influence degree calculating module is used for: according to the grid-connected active power of the equipment at the simulation time point and the node set attribute of each connected node, respectively calculating the energy deviation rate of the equipment connected through each connected node in the period of transient power angle stability influence calculation, and taking the sum of the energy deviation rates of the equipment connected through each connected node as the influence of the equipment on the transient power angle stability of the power grid after the disturbance in the preset power grid running state.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform a transient power angle stability influence calculation method based on an energy deviation rate.
A computing device comprising one or more processors, one or more memories, and one or more programs, wherein 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 a transient power angle stability influence calculation method based on an energy deviation rate.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing is illustrative of the present invention and is not to be construed as limiting thereof, but rather as providing for the use of additional embodiments and advantages of all such modifications, equivalents, improvements and similar to the present invention are intended to be included within the scope of the present invention as defined by the appended claims.
Claims (10)
1. The transient power angle stability influence degree calculation method based on the energy deviation rate is characterized by comprising the following steps of:
determining power grid energy balance supporting points corresponding to simulation time points after preset disturbance according to a preset power grid running state;
aiming at each simulation time point, dividing nodes in a power grid network equation into two node sets by taking a power grid energy balance fulcrum as a boundary, and obtaining two node sets corresponding to each simulation time point; nodes in the grid network equation comprise potential nodes in the synchronous generator;
aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator which are related to the nodes in the two node sets, defining a node set with a large central angle of the rotor inertia as a leading node set, defining a node set with a small central angle of the rotor inertia as a lagging node set, and obtaining the attributes of the two node sets corresponding to each simulation time point;
according to the relative rotor inertia center angles between the leading node set and the lagging node set corresponding to each simulation time point, determining the time period for calculating the transient power angle stability influence degree;
according to the grid-connected active power of the equipment at the simulation time point and the node set attribute of each connected node, respectively calculating the energy deviation rate of the equipment which is connected through each connected node in the transient power angle stability influence calculation period, and taking the sum of the energy deviation rates of the equipment which is connected through each connected node as the influence of the equipment on the transient power angle stability of the power grid after the disturbance is preset in the preset power grid running state.
2. The method for calculating the transient power angle stability influence degree based on the energy deviation rate according to claim 1, wherein the method comprises the following steps of: aiming at the running state of a preset power grid, determining power grid energy balance supporting points corresponding to each simulation time point after the preset disturbance, wherein the specific process is that,
performing time domain simulation of preset disturbance on a preset power grid running state to obtain an admittance matrix for network equation calculation corresponding to each simulation time point including a preset disturbance initial moment;
aiming at admittance matrixes for network equation calculation corresponding to simulation time points, respectively equating a power grid into a 3-node 5-branch network meeting set conditions through static network equivalence, wherein the nodes comprise two equivalent power supply nodes and 1 equivalent load node, the network branch comprises a ground branch of 3 nodes and a branch between the two equivalent power supply nodes and the equivalent load node, and the set conditions are that currents injected into the equivalent load node by the two equivalent power supply nodes are equal;
and for each simulation time point, determining the position of the equivalent load node in the power grid before the equivalent through backtracking of the equivalent process of the static network, and taking the position as an energy balance pivot of the power grid.
3. The method for calculating the transient power angle stability influence degree based on the energy deviation rate according to claim 1, wherein the method comprises the following steps of: aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator associated with the nodes in the two node sets, defining the node set with the large central angle of the rotor inertia as a leading node set, defining the node set with the small central angle of the rotor inertia as a lagging node set, and obtaining the attribute of the two node sets corresponding to each simulation time point, wherein the specific calculation formula of the central angles of the rotor inertia of the synchronous generator associated with the nodes in the two node sets is as follows:
wherein A is i 、B i Respectively two node sets corresponding to the ith simulation time point,
respectively A i 、B i Ith simulation time point rotor inertia center angle of medium-association synchronous generator, M i.a 、δ i.a Respectively A i Internal potential node a of medium synchronous generatorMoment of inertia and internal potential phase angle, M, of the ith simulation time point of the corresponding synchronous generator i.b 、δ i.b Respectively B i The moment of inertia and the internal potential phase angle of the ith simulation time point of the synchronous generator corresponding to the internal potential node b of the medium synchronous generator;
the node set attribute is referred to as "lead" or "lag".
4. A method for calculating the transient power angle stability influence degree based on the energy deviation rate according to claims 1 and 3, wherein: according to the relative rotor inertia center angle between the leading node set and the lagging node set corresponding to each simulation time point, determining the time period for transient power angle stability influence degree calculation, wherein the specific process is as follows;
if A corresponding to the ith simulation time point i For the leading node set, the relative rotor inertia center angle delta between the leading node set and the lagging node set corresponding to the ith simulation time point i Is arranged as
Otherwise, the relative rotor inertia center angle delta between the leading node set and the lagging node set corresponding to the ith simulation time point i Set to->
Starting from the 2 nd simulation time point, determining a transient power angle stability influence degree calculation period in a mode that the simulation time points are sequentially increased, wherein the method specifically comprises the following steps: if delta i <δ s And delta i-1 <δ i 、δ i >δ i+1 Taking the time period between the 1 st simulation time point and the i-th simulation time point as the time period for calculating the transient power angle stability influence degree, if delta i ≥δ s Taking the time period between the 1 st simulation time point and the i-th simulation time point as the time period for calculating the transient power angle stability influence degree, wherein delta s Setting 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, the determination of the transient power angle stability influence calculation period is terminated.
5. The method for calculating the transient power angle stability influence degree based on the energy deviation rate according to claim 1, wherein the method comprises the following steps of: according to the grid-connected active power of the equipment at the simulation time point and the node set attribute of each connected node, respectively calculating the energy deviation rate of the equipment connected through each connected node in the transient power angle stability influence calculation period, taking the sum of the energy deviation rates of the equipment connected through each connected node as the influence of the equipment on the transient power angle stability of the power grid after the equipment is subjected to preset disturbance in the preset power grid running state, wherein a specific calculation formula is as follows;
wherein B is D Injecting a device D into a node set of the power grid, D tas The influence of the equipment D on the stability of the transient power angle of the power grid after the disturbance is preset under the preset power grid running state is shown as D j For apparatus D via B D Energy deviation rate s of grid connection of middle node j i.j For positive or negative sign, if the ith simulation time point B D The middle node j belongs to the leading node set, and takes a negative sign, otherwise, s i.j Taking positive sign, wherein I is the last simulation time point corresponding to the transient power angle stability influence calculation period, and P i.j For the ith simulationTime point B D Grid-connected active power of middle node j, t i The transient process time corresponding to the ith simulation time point.
6. The method for calculating the transient power angle stability influence degree based on the energy deviation rate according to claim 1, wherein the method comprises the following steps of: aiming at the situation that the power grid is composed of two or more asynchronous operation sub-power grids after the preset disturbance in the operation state of the preset power grid, each asynchronous operation sub-power grid is independently processed.
7. Transient state power angle stability influence degree computing system based on energy deviation rate, its characterized in that: comprising the steps of (a) a step of,
the energy balance pivot point determining module: determining power grid energy balance supporting points corresponding to simulation time points after preset disturbance according to a preset power grid running state;
the node set generation module: aiming at each simulation time point, dividing nodes in a power grid network equation containing potential nodes in the synchronous generator into two node sets by taking a power grid energy balance fulcrum as a boundary, and obtaining two node sets corresponding to each simulation time point;
the node set attribute determining module: aiming at each simulation time point, calculating the central angles of the rotor inertia of the synchronous generator which are related to the nodes in the two node sets, defining a node set with a large central angle of the rotor inertia as a leading node set, defining a node set with a small central angle of the rotor inertia as a lagging node set, and obtaining the attributes of the two node sets corresponding to each simulation time point;
influence degree calculation period determination module: according to the relative rotor inertia center angles between the leading node set and the lagging node set corresponding to each simulation time point, determining the time period for calculating the transient power angle stability influence degree;
and the equipment influence degree calculating module is used for: according to the grid-connected active power of the equipment at the simulation time point and the node set attribute of each connected node, respectively calculating the energy deviation rate of the equipment connected through each connected node in the period of transient power angle stability influence calculation, and taking the sum of the energy deviation rates of the equipment connected through each connected node as the influence of the equipment on the transient power angle stability of the power grid after the disturbance in the preset power grid running state.
8. The energy deviation rate-based transient power angle stability influence calculation system of claim 7, wherein: in the node set attribute determining module, a specific calculation formula of the synchronous generator rotor inertia center angle related to the nodes in the two node sets is as follows:
wherein A is i 、B i Respectively two node sets corresponding to the ith simulation time point,
respectively A i 、B i Ith simulation time point rotor inertia center angle of medium-association synchronous generator, M i.a 、δ i.a Respectively A i Moment of inertia and internal potential phase angle, M, of the ith simulation time point of the synchronous generator corresponding to internal potential node a of the medium synchronous generator i.b 、δ i.b Respectively B i And the moment of inertia and the internal potential phase angle of the ith simulation time point of the synchronous generator corresponding to the internal potential node b of the medium synchronous generator.
9. A computer readable storage medium storing one or more programs, characterized by: the one or more programs include instructions, which when executed by a computing device, cause the computing device to perform any of the methods of claims 1-6.
10. A computing device, characterized by: comprising the steps of (a) a step of,
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 of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211568869.7A CN116186979A (en) | 2022-12-08 | 2022-12-08 | Transient power angle stability influence degree calculation method and system based on energy deviation rate |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211568869.7A CN116186979A (en) | 2022-12-08 | 2022-12-08 | Transient power angle stability influence degree calculation method and system based on energy deviation rate |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116186979A true CN116186979A (en) | 2023-05-30 |
Family
ID=86437302
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211568869.7A Pending CN116186979A (en) | 2022-12-08 | 2022-12-08 | Transient power angle stability influence degree calculation method and system based on energy deviation rate |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116186979A (en) |
-
2022
- 2022-12-08 CN CN202211568869.7A patent/CN116186979A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103825269B (en) | Rapid probabilistic load flow calculation method considering static power frequency characteristics of electric power system | |
| CN111797510A (en) | Method and system for calculating short circuit ratio of new energy station | |
| CN112260326B (en) | Equivalent short circuit ratio calculation method and system for new energy cluster | |
| CN109783895A (en) | A kind of electric system hybrid simulation method and system containing distributed generation resource | |
| Lu et al. | Decentralised nonlinear H∞ excitation control based on regulation linearisation | |
| CN104836260B (en) | Calculation method for allowed capacity of DGs in active power distribution network under voltage constraint | |
| CN110970905B (en) | Multi-power-supply reactive power coordination method and system for optimizing voltage control capability of wind power plant | |
| CN115481856A (en) | Comprehensive energy system multi-scale scheduling method and system considering comprehensive demand response | |
| CN106712032A (en) | Optimal power flow model construction method considering active power voltage regulation capacity of wind turbine generator set | |
| CN112072692A (en) | Impedance equivalence method and device for new energy power generation station | |
| CN105896547B (en) | A kind of bulk power grid hierarchical voltage control method under wind power integration | |
| CN109390959B (en) | Storage battery energy storage control method based on virtual synchronous machine technology | |
| CN110380423A (en) | A kind of wind farm group Reactive-power control object priority determines method and device | |
| CN116186979A (en) | Transient power angle stability influence degree calculation method and system based on energy deviation rate | |
| CN111884234B (en) | Emergency control method and device for precise cutting machine, electronic equipment and storage medium | |
| CN114649826A (en) | Black-start partitioning method and system considering wind power plant | |
| CN113346553A (en) | Evaluation method and device for output ratio of renewable energy power station | |
| CN109698516A (en) | The maximum capacity computing system and method for renewable energy access power distribution network | |
| CN109873454B (en) | Method and device for evaluating group voltage supporting capability in receiving end power grid | |
| CN116127742A (en) | Method and system for calculating influence degree of synchronous generator on transient power angle stability of power grid | |
| Li et al. | Newton-Raphson Method Based Power Flow Analysis and Dynamic Security Assessment [J] | |
| Ding et al. | An Integer Variable Reduction and Neighborhood Search Cutting Plane Method for Security Constrained Unit Commitment | |
| CN117200266B (en) | Distribution line photovoltaic collection point energy storage capacity configuration method, device and equipment | |
| Xu et al. | Intelligent Commutation Strategy Based on Improved Discrete Particle Swarm Optimization | |
| CN111092447B (en) | Optimal energy storage grid connection point configuration method and device under system small interference stability boundary |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |