CN111181153B - Method and system for determining Thevenin equivalent parameters based on branch response information - Google Patents
Method and system for determining Thevenin equivalent parameters based on branch response information Download PDFInfo
- Publication number
- CN111181153B CN111181153B CN201911141762.2A CN201911141762A CN111181153B CN 111181153 B CN111181153 B CN 111181153B CN 201911141762 A CN201911141762 A CN 201911141762A CN 111181153 B CN111181153 B CN 111181153B
- Authority
- CN
- China
- Prior art keywords
- power system
- thevenin
- internal potential
- reactance
- iteration
- 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.)
- Active
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
-
- 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
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a method and a system for determining Thevenin equivalent parameters based on branch response information, and belongs to the technical field of voltage stability online analysis and monitoring of safety and stability of an electric power system. The method comprises the following steps: selecting any one of a plurality of buses in the power system as a boundary, and taking a part in the buses as a target to be detected; acquiring a voltage amplitude, active power and reactive power of a target bus node to be detected; acquiring thevenin internal potential amplitude, thevenin internal potential phase angle and thevenin reactance of the power system with iteration rounds of t and t + 1; determining whether the absolute values of the difference values of the Thevenin internal potential amplitude, the Thevenin internal potential phase angle and the Thevenin reactance of the power system with the iteration round number t and the t +1 round are smaller than a preset value; and determining the Thevenin internal potential amplitude, the Thevenin internal potential phase angle and the Thevenin reactance of the power system with the iteration round number of t +1 rounds as Thevenin equivalent parameters of the power system. The invention ensures the voltage stability of the power system.
Description
Technical Field
The invention relates to the technical field of voltage stability online analysis and monitoring of safety and stability of a power system, in particular to a method and a system for determining Thevenin equivalent parameters based on branch response information.
Background
Most provincial power grids in the eastern region of China are receiving-end power grids. In order to ensure safe and reliable operation of a receiving-end power grid without voltage collapse, online voltage stability monitoring and analysis are required. Thevenin equivalent of an electric power system is a basic technology for realizing online voltage stability monitoring and analysis.
Thevenin equivalence equates an arbitrarily selected sub-network in a power system to the form of a series connection of an internal potential and a reactance. Thevenin equivalence can effectively simplify a complex power system, and can also accurately indicate the voltage stability and stability margin of the power system by using equivalent internal potential and internal impedance obtained by the equivalent internal potential and the equivalent internal impedance, so that effective perception of a dispatcher on the safety situation of a power grid is promoted, and voltage instability of the system is prevented. The Thevenin equivalent parameters are closely related to the power grid structure and the system operation mode and change along with the power grid structure and the system operation mode, and the accuracy of the Thevenin equivalent parameters directly determines the accuracy of recognizing the safety characteristics of the power grid by a dispatcher, so that the accurate determination of the Thevenin equivalent parameters under various working conditions is related to the correct judgment of the voltage stability of the power grid, and the Thevenin equivalent parameters are extremely important work.
The traditional Thevenin equivalent parameter determination method is suitable for the following scenes: during a plurality of adjacent time steps of continuous sampling, the internal potential amplitude, the internal potential phase angle and the Thevenin reactance of a selected certain sub-network are kept unchanged; at the same time, the active, reactive, and voltage magnitudes of the external networks outside the sub-network need to change significantly between these adjacent time steps.
In fact, the actual operation condition of the power system often differs greatly from the above scenario. A typical scenario is: in a plurality of time steps of Thevenin equivalence, which are carried out on a certain sub-network of an electric power system, in order to adjust the voltage level of the system to keep the voltage level in a reasonable range, the switching operation of reactive equipment (inductor or capacitor) often exists in an external network except the sub-network. Under the condition, the Thevenin equivalent internal potential of the sub-network can be changed obviously, and the traditional method has larger error and is not applicable any more.
The traditional method can be used only in the scene that the magnitude of thevenin internal potential, the phase angle of the internal potential and thevenin reactance are kept unchanged in a plurality of adjacent time steps. Under the working condition that the potential amplitude in thevenin changes, thevenin equivalent parameters obtained by the traditional method have poor precision and are difficult to support and accurately judge the voltage stability.
Disclosure of Invention
Aiming at the problems, the invention provides a method for determining thevenin equivalent parameters based on branch response information, which comprises the following steps:
selecting any one of a plurality of buses in the power system as a boundary, and taking the part in the buses as a target to be measured;
according to the preset time steps, acquiring the voltage amplitude of a target bus node to be detected, the active power of a target bus line to be detected and the reactive power of the target bus line to be detected;
according to the amplitude to be subjected to voltage, the active power and the reactive power, iteration is carried out on the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin, the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin are obtained, the number of iteration rounds t and t +1 rounds are obtained, and t is larger than or equal to 1;
obtaining the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin in the iteration rounds t and t +1, and determining whether the absolute value of the difference value of the internal potential amplitudes, the absolute value of the internal potential phase angle difference and the absolute value of the difference value of the Thevenin reactance of the power system Thevenin in the iteration rounds t and t +1 are smaller than a preset value or not;
and when the absolute value of the difference value of the internal potential amplitude value of the power system Thevenin, the absolute value of the internal potential phase angle difference value and the absolute value of the difference value of the reactance of the iteration round number t and the iteration round number t +1 are all smaller than preset values, determining the internal potential amplitude value of the power system Thevenin, the internal potential phase angle and the reactance of the iteration round number t +1 as equivalent parameters of the power system Thevenin.
Optionally, the preset time step is greater than or equal to 2.
Optionally, the preset value is 0.001.
Optionally, the thevenin internal potential amplitude, the thevenin internal potential phase angle and the power system thevenin reactance of the iteration round number t and the t +1 round are obtained, specifically:
iteration is carried out by enabling the Dyvenin internal potential amplitude of the power system in 1 iteration round number to be 1, the Dyvenin internal potential phase angle of the power system in 1 iteration round number to be 0 and the Dyvenin reactance of the power system in 1 iteration round number to be 0, and the Dyvenin internal potential amplitude of the power system, the Dyvenin internal potential phase angle of the power system and the Dyvenin reactance of the power system in t iteration round number are obtained;
according to the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t, the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t +1 are obtained, and the formula is as follows:
wherein eta is 0.001-0.01,
The number of iteration rounds is t +1 round, namely the magnitude of the internal potential of the Thevenin of the power system is delta th (t + 1) iteration round number t +1 round of power system Thevenin internal potential phase angle X th (t + 1) iterative turns t +1 turns of Thevenin reactance of power system,
The internal potential amplitude value delta of the Thevenin of the power system is t rounds of iteration th (t) iterative round number t round power system Thevenin internal potential phase angle X th (t) iterative round number t round Thevenin reactance, U of power system k Voltage amplitude, P, of a bus node of a target to be measured k For the active power, Q, of the target busbar line to be tested k For the reactive power of the target bus line to be measured, K is a preset time step, K is any one time step, f k As the sum of the active residuals e k Is a reactive residual.
The invention also provides a system for determining Thevenin equivalent parameters based on branch response information, which comprises the following steps:
the target determining module selects any one of the buses in the power system as a boundary, and takes the part in the buses as a target to be detected;
the parameter acquisition module is used for acquiring the voltage amplitude of a node of a target bus to be detected, the active power of a target bus line to be detected and the reactive power of the target bus line to be detected according to a preset time step;
the iteration module is used for iterating the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin according to the amplitude to be voltage, the active power and the reactive power to obtain the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin with iteration turns t and t +1 turns, wherein t is greater than or equal to 1;
the comparison module is used for obtaining the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin in the iteration rounds t and t +1, and determining whether the absolute value of the difference value of the internal potential amplitude values of the power system Thevenin, the absolute value of the difference value of the internal potential phase angles of the power system Thevenin and the absolute value of the difference value of the Thevenin reactance of the power system are smaller than a preset value or not;
and the judgment module is used for determining that the Thevenin internal potential amplitude, the internal potential phase angle and the reactance of the power system in the iteration round number t +1 are equivalent parameters of the power system when the absolute value of the difference value of the Thevenin internal potential amplitude, the absolute value of the internal potential phase angle difference and the absolute value of the difference value of the reactance of the power system in the iteration round number t +1 are smaller than preset values.
Optionally, the preset time step is greater than or equal to 2.
Optionally, the preset value is 0.001.
Optionally, the iteration module is configured to:
iteration is carried out by enabling the Dyvenin internal potential amplitude of the power system in 1 iteration round number to be 1, the Dyvenin internal potential phase angle of the power system in 1 iteration round number to be 0 and the Dyvenin reactance of the power system in 1 iteration round number to be 0, and the Dyvenin internal potential amplitude of the power system, the Dyvenin internal potential phase angle of the power system and the Dyvenin reactance of the power system in t iteration round number are obtained;
according to the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t round, acquiring the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t +1 round, wherein the formula is as follows:
wherein eta is 0.001-0.01,
The number of iteration rounds is t +1 round, namely the magnitude and delta of the Thevenin internal potential of the power system th (t + 1) iteration round number t +1 round of power system Thevenin internal potential phase angle X th (t + 1) iteration turns t +1 turns of Thevenin reactance of power system,
The internal potential amplitude value delta of the Thevenin of the power system is t rounds of iteration th (t) iterative round number t round power system Thevenin internal potential phase angle X th (t) iterative round number t round Thevenin reactance, U of power system k Voltage amplitude, P, of a bus node of a target to be measured k For the active power, Q, of the target busbar line to be tested k For the reactive power of the target bus line to be measured, K is a preset time step, K is any one time step, f k As the sum of the active residuals e k Is a reactive residual.
The method can obtain accurate Thevenin equivalent parameters, thereby supporting accurate voltage stability judgment, improving the effective perception of dispatching personnel on the safety situation of the power grid, preventing the voltage instability of the power system and ensuring the voltage stability of the power system.
According to the method, thevenin equivalent parameters are calculated only by using the bus and the local measurement data on the outgoing line of the bus, a system model and parameters are not depended on, and the Thevenin equivalent parameters obtained by the method are more accurate compared with the traditional method for the working condition of the change of the potential amplitude in Thevenin.
Drawings
FIG. 1 is a flow chart of a method of determining Thevenin equivalent parameters based on branch response information according to the present invention;
FIG. 2 is a Thevenin equivalence schematic diagram of the power system of the method for determining Thevenin equivalence parameters based on branch response information;
FIG. 3 is a graph showing an active power curve of a target busbar line to be tested and a reactive power curve of the target busbar line to be tested according to the method for determining Thevenin equivalent parameters based on branch response information;
FIG. 4 is a voltage amplitude curve diagram of a target bus node to be measured by the method for determining Thevenin equivalent parameters based on branch response information;
FIG. 5 is a Thevenin internal potential amplitude curve diagram of a power system of the method for determining Thevenin equivalent parameters based on branch response information;
FIG. 6 is a Thevenin internal potential phase angle graph of a power system of the method for determining Thevenin equivalent parameters based on branch response information of the invention;
FIG. 7 is a Thevenin reactance graph of a force system of a method for determining Thevenin equivalent parameters based on branch response information according to the present invention;
FIG. 8 is a diagram of a system for determining Thevenin equivalent parameters based on branch response information in accordance with the present invention.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the embodiments described herein, which are provided for complete and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, the same units/elements are denoted by the same reference numerals.
Unless otherwise defined, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Further, it will be understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
The invention provides a method for determining Thevenin equivalent parameters based on branch response information, which comprises the following steps of:
the thevenin equivalent diagram of the power system is shown in fig. 2, any one of a plurality of buses in the power system is selected as a boundary, the part inside the bus is used as a target to be detected, the target to be detected is a subsystem, the boundary outside the subsystem is used as an external system, and U is adopted k Voltage amplitude, P, of a bus node of a target to be measured k Active power, Q, for the target busbar line to be tested k The reactive power of the bus line of the target to be measured,
For the magnitude, delta, of the Thevenin internal potential of the power system th Thevenin internal potential phase angle and X of power system th Is a Thevenin reactance of a wheel power system;
according to the preset time step, the preset time step is more than or equal to 2;
the voltage amplitude of the target bus node to be detected and the active power of the target bus line to be detected are obtained, and the voltage amplitude of the target bus node to be detected and the active power curve of the target bus line to be detected are shown in fig. 3;
obtaining the reactive power of a target busbar line to be tested, wherein a reactive power curve chart of the target busbar line to be tested is shown in fig. 4;
iterating the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system according to the voltage amplitude of the target bus node to be detected, the active power of the target bus line to be detected and the reactive power of the target bus line to be detected, and acquiring the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system for t +1 rounds of iteration rounds, wherein t is greater than or equal to 1;
the method comprises the following steps of obtaining thevenin internal potential amplitude, thevenin internal potential phase angle and power system thevenin reactance of iteration rounds t and t +1, and specifically comprises the following steps:
iteration is carried out by enabling the Dyvenin internal potential amplitude of the power system in 1 iteration round number to be 1, the Dyvenin internal potential phase angle of the power system in 1 iteration round number to be 0 and the Dyvenin reactance of the power system in 1 iteration round number to be 0, and the Dyvenin internal potential amplitude of the power system, the Dyvenin internal potential phase angle of the power system and the Dyvenin reactance of the power system in t iteration round number are obtained;
according to the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t round, acquiring the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t +1 round, wherein the formula is as follows:
wherein eta is 0.001-0.01,
The number of iteration rounds is t +1 round, namely the magnitude of the internal potential of the Thevenin of the power system is delta th (t + 1) iteration round number t +1 round of power system Thevenin internal potential phase angle X th (t + 1) iterative turns t +1 turns of Thevenin reactance of power system,
The internal potential amplitude value delta of the Thevenin of the power system is t rounds of iteration th (t) iterative round number t round power system Thevenin internal potential phase angle X th (t) iterative turns t turns of Thevenin reactance, U of power system k Voltage amplitude, P, of a bus node of a target to be measured k For the active power, Q, of the target busbar line to be tested k For the reactive power of the target bus line to be measured, K is a preset time step, K is any one time step, f k Is the sum of active residual errors e k Is a reactive residual.
Acquiring the internal potential amplitude value of the power system Thevenin, the internal potential phase angle of the power system Thevenin and the Thevenin reactance of the power system according to the iteration turns t and t +1, and determining whether the absolute value of the difference value of the internal potential amplitude values of the power system Thevenin, the absolute value of the difference value of the internal potential phase angles of the power system Thevenin and the absolute value of the difference value of the Thevenin reactance of the power system Thevenin of the iteration turns t and t +1 are smaller than a preset value or not;
when the absolute value of the difference value of the internal potential amplitude value of the power system Thevenin, the absolute value of the phase angle difference value of the internal potential of the power system Thevenin and the absolute value of the difference value of the Thevenin reactance of the power system, which are the iteration turns t and t +1, are all smaller than preset values, the preset values are 0.001, and the internal potential amplitude value of the power system Thevenin, the internal potential phase angle of the power system Thevenin and the Thevenin reactance of the power system, which are the iteration turns t +1, are determined as equivalent parameters of the Thevenin of the power system;
graphs of the magnitude of the Thevenin internal potential of the power system, the phase angle of the Thevenin internal potential of the power system, and the Thevenin reactance of the power system are shown in FIGS. 5, 6, and 7, respectively.
The present invention further provides a system 200 for determining thevenin equivalent parameters based on branch response information, as shown in fig. 8, comprising:
the target determining module 201 selects any one of a plurality of buses in the power system as a boundary, and takes a part in the buses as a target to be detected;
the parameter obtaining module 202 is used for obtaining the voltage amplitude of the node of the target bus to be detected, the active power of the target bus line to be detected and the reactive power of the target bus line to be detected according to the preset time step, wherein the preset time step is greater than or equal to 2;
the iteration module 203 is used for iterating the Thevenin internal potential amplitude value of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system according to the voltage amplitude value of the target bus node to be detected, the active power of the target bus line to be detected and the reactive power of the target bus line to be detected, and acquiring the Thevenin internal potential amplitude value of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system in t +1 iteration rounds, wherein t is greater than or equal to 1;
the iteration module 203 is configured to:
iteration is carried out by enabling the Dyinan internal potential amplitude of the power system with 1 iteration round number, the Dyinan internal potential phase angle of the power system with 1 iteration round number to be 1, and the Dyinan reactance of the power system with 1 iteration round number to be 0, and the Dyinan internal potential amplitude of the power system, the Dyinan internal potential phase angle of the power system and the Dyinan reactance of the power system with t iteration round numbers are obtained;
according to the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t, the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t +1 are obtained, and the formula is as follows:
wherein eta is 0.001-0.01,
The number of iteration rounds is t +1 round, namely the magnitude and delta of the Thevenin internal potential of the power system th (t + 1) iteration round number t +1 round of power system Thevenin internal potential phase angle X th (t + 1) iteration turns t +1 turns of Thevenin reactance of power system,
The internal potential amplitude value delta of the Thevenin of the power system is t rounds of iteration th (t) iterative round number t round power system Thevenin internal potential phase angle X th (t) iterative turns t turns of Thevenin reactance, U of power system k Voltage amplitude, P, of a target bus node to be measured k For the active power, Q, of the target busbar line to be tested k For the reactive power of the target bus line to be measured, K is a preset time step, K is any one time step, f k As the sum of the active residuals e k Is a reactive residual.
The comparison module 204 is used for obtaining the internal potential amplitude value of the power system thevenin, the internal potential phase angle of the power system thevenin and the thevenin reactance of the power system according to iteration rounds t and t +1, and determining whether the absolute value of the difference value of the internal potential amplitude values of the power system thevenin, the absolute value of the difference value of the internal potential phase angle of the power system thevenin and the absolute value of the difference value of the thevenin reactance of the power system are smaller than a preset value or not;
the judgment module 205 determines that when the absolute value of the difference value of the internal potential amplitude of the power system thevenin, the absolute value of the phase angle difference value of the internal potential of the power system thevenin and the absolute value of the difference value of the thevenin reactance of the power system are all smaller than preset values, the preset values are 0.001, and determines that the internal potential amplitude of the power system thevenin, the internal potential phase angle of the power system thevenin and the thevenin reactance of the power system thevenin in the iteration round t +1 are equivalent parameters of the power system thevenin.
The method can obtain accurate Thevenin equivalent parameters, thereby supporting accurate voltage stability judgment, improving the effective perception of dispatching personnel on the safety situation of the power grid, preventing the voltage instability of the system and ensuring the voltage stability of the power system.
According to the method, thevenin equivalent parameters are calculated only by using local measurement data on the bus and the outgoing line thereof, a system model and parameters are not depended on, and the Thevenin equivalent parameters obtained by the method are more accurate compared with the traditional method for the working condition of the change of the internal potential amplitude of Thevenin.
As will be appreciated by one skilled in the art, 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 so forth) 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 flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams 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.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.
Claims (6)
1. A method of determining thevenin equivalent parameters based on branch response information, the method comprising:
selecting any one of a plurality of buses in the power system as a boundary, and taking the part in the buses as a target to be measured;
according to the preset time step, acquiring the voltage amplitude of a target bus node to be detected, the active power of a target bus line to be detected and the reactive power of the target bus line to be detected;
iterating the internal potential amplitude, the internal potential phase angle and the reactance of the power system thevenin according to the voltage amplitude, the active power and the reactive power of the target bus node to be detected, and acquiring the internal potential amplitude, the internal potential phase angle and the reactance of the power system thevenin with iteration rounds t and t +1 rounds, wherein t is greater than or equal to 1;
the method for acquiring the Thevenin internal potential amplitude, the internal potential phase angle and the reactance of the power system with the iteration round number t and the t +1 round specifically comprises the following steps:
iteration is carried out by enabling the Dyvenin internal potential amplitude of the power system in 1 iteration round number to be 1, the Dyvenin internal potential phase angle of the power system in 1 iteration round number to be 0 and the Dyvenin reactance of the power system in 1 iteration round number to be 0, and the Dyvenin internal potential amplitude of the power system, the Dyvenin internal potential phase angle of the power system and the Dyvenin reactance of the power system in t iteration round number are obtained;
according to the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t, the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t +1 are obtained, and the formula is as follows:
wherein eta is 0.001-0.01,
The number of iteration rounds is t +1 round, namely the magnitude and delta of the Thevenin internal potential of the power system th (t + 1) iterative turns t +1 turns of the internal potential phase angle of thevenin in the power system,X th (t + 1) iterative turns t +1 turns of Thevenin reactance of power system,
The internal potential amplitude value delta of the Thevenin of the power system is t rounds of iteration th (t) iterative round number t round power system Thevenin internal potential phase angle X th (t) iterative turns t turns of Thevenin reactance, U of power system k Voltage amplitude, P, of a bus node of a target to be measured k For the active power, Q, of the target busbar line to be tested k For the reactive power of the target bus line to be measured, K is a preset time step, K is any one time step, f k As the sum of the active residuals e k Is a reactive residual error;
obtaining the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin in the iteration rounds t and t +1, and determining whether the absolute value of the difference value of the internal potential amplitudes, the absolute value of the internal potential phase angle difference and the absolute value of the difference value of the Thevenin reactance of the power system Thevenin in the iteration rounds t and t +1 are smaller than a preset value or not;
and when the absolute value of the difference value of the internal potential amplitude value of the power system Thevenin, the absolute value of the internal potential phase angle difference value and the absolute value of the difference value of the reactance of the iteration round number t and the iteration round number t +1 are all smaller than preset values, determining the internal potential amplitude value of the power system Thevenin, the internal potential phase angle and the reactance of the iteration round number t +1 as equivalent parameters of the power system Thevenin.
2. The method of claim 1, the preset time step being greater than or equal to 2.
3. The method of claim 1, wherein the preset value is 0.001.
4. A system for determining thevenin equivalent parameters based on branch response information, the system comprising:
the target determining module selects any one of a plurality of buses in the power system as a boundary, and takes the part in the buses as a target to be detected;
the parameter acquisition module is used for acquiring the voltage amplitude of a node of a target bus to be detected, the active power of a target bus line to be detected and the reactive power of the target bus line to be detected according to a preset time step;
the iteration module is used for iterating the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin according to the voltage amplitude, the active power and the reactive power of the target bus node to be detected to obtain the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin for t and t +1 iteration rounds, wherein t is greater than or equal to 1;
the method for acquiring the Thevenin internal potential amplitude, the internal potential phase angle and the reactance of the power system with the iteration round number t and the t +1 round specifically comprises the following steps:
iteration is carried out by enabling the Dyvenin internal potential amplitude of the power system in 1 iteration round number to be 1, the Dyvenin internal potential phase angle of the power system in 1 iteration round number to be 0 and the Dyvenin reactance of the power system in 1 iteration round number to be 0, and the Dyvenin internal potential amplitude of the power system, the Dyvenin internal potential phase angle of the power system and the Dyvenin reactance of the power system in t iteration round number are obtained;
according to the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t, the Thevenin internal potential amplitude of the power system, the Thevenin internal potential phase angle of the power system and the Thevenin reactance of the power system of the iteration round number t +1 are obtained, and the formula is as follows:
wherein eta is 0.001-0.01,
The number of iteration rounds is t +1 round, namely the magnitude and delta of the Thevenin internal potential of the power system th (t + 1) iteration round number t +1 round of power system Thevenin internal potential phase angle X th (t + 1) iterative turns t +1 turns of Thevenin reactance of power system,
The internal potential amplitude value delta of the Thevenin of the power system is t rounds of iteration th (t) iterative round number t round power system Thevenin internal potential phase angle X th (t) iterative turns t turns of Thevenin reactance, U of power system k Voltage amplitude, P, of a bus node of a target to be measured k For the active power, Q, of the target busbar line to be tested k For the reactive power of the target bus line to be measured, K is a preset time step, K is any one time step, f k As the sum of the active residuals e k Is a reactive residual error;
the comparison module is used for obtaining the internal potential amplitude, the internal potential phase angle and the reactance of the power system Thevenin in the iteration rounds t and t +1, and determining whether the absolute value of the difference value of the internal potential amplitude values of the power system Thevenin, the absolute value of the difference value of the internal potential phase angles of the power system Thevenin and the absolute value of the difference value of the Thevenin reactance of the power system are smaller than a preset value or not;
and the judgment module is used for determining that the Thevenin internal potential amplitude, the internal potential phase angle and the reactance of the power system in the iteration round number t +1 are equivalent parameters of the power system when the absolute value of the difference value of the Thevenin internal potential amplitude, the absolute value of the internal potential phase angle difference and the absolute value of the difference value of the reactance of the power system in the iteration round number t +1 are smaller than preset values.
5. The system of claim 4, the preset time step being greater than or equal to 2.
6. The system of claim 4, wherein the preset value is 0.001.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911141762.2A CN111181153B (en) | 2019-11-20 | 2019-11-20 | Method and system for determining Thevenin equivalent parameters based on branch response information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911141762.2A CN111181153B (en) | 2019-11-20 | 2019-11-20 | Method and system for determining Thevenin equivalent parameters based on branch response information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111181153A CN111181153A (en) | 2020-05-19 |
| CN111181153B true CN111181153B (en) | 2022-10-11 |
Family
ID=70651889
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911141762.2A Active CN111181153B (en) | 2019-11-20 | 2019-11-20 | Method and system for determining Thevenin equivalent parameters based on branch response information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111181153B (en) |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102522742B (en) * | 2011-10-26 | 2014-02-26 | 重庆大学 | Estimation method of external network Thevenin equivalent parameters based on single point measurement information |
| CN102420427B (en) * | 2011-12-02 | 2014-05-14 | 余杭供电局 | Method for calculating voltage stability margin of regional power grid by considering outer network equivalence |
| CN103872679A (en) * | 2014-03-07 | 2014-06-18 | 国家电网公司 | Identification method for power grid Thevenin equivalent model parameter under condition of weak signals |
| CN104778367B (en) * | 2015-04-20 | 2017-08-08 | 山东大学 | Wide area Thevenin's equivalence parameter on-line calculation method based on a single state section |
| GB2567460B (en) * | 2017-10-12 | 2020-02-12 | Univ Cape Town | Method and apparatus for obtaining Thévenin equivalent parameters of a power network |
| CN107958120B (en) * | 2017-11-30 | 2021-10-08 | 上海电力学院 | System Thevenin equivalent parameter calculation method based on power series expansion |
| CN108336730B (en) * | 2018-03-14 | 2021-07-13 | 哈尔滨工业大学 | Thevenin equivalent parameter identification method based on reduced order self-sensitivity |
-
2019
- 2019-11-20 CN CN201911141762.2A patent/CN111181153B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN111181153A (en) | 2020-05-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110579682B (en) | Fault recording data transient homologous comparison method and device | |
| CN106505557B (en) | Remote measurement error identification method and device | |
| CN107016489A (en) | A kind of electric power system robust state estimation method and device | |
| CN109633366B (en) | Power distribution network single-phase earth fault phase selection method and phase selection device | |
| CN106709147B (en) | BPA-based automatic batch calculation method for limit clearing time of short circuit fault | |
| CN109284933B (en) | Electronic transformer state evaluation system and method based on mathematical statistics | |
| CN111130074A (en) | Direct-current transmission line ultra-high-speed protection system and method based on voltage first extreme value time | |
| CN109829246B (en) | Line parameter identification method based on parameter comprehensive suspicion | |
| CN114065118B (en) | Power system robust state estimation method based on exponential function | |
| CN113484804B (en) | Wrong wiring identification method applied to feeder terminal device | |
| CN111181153B (en) | Method and system for determining Thevenin equivalent parameters based on branch response information | |
| CN110417009A (en) | Power distribution network hybrid robust state estimation method based on data of different sampling periods | |
| CN111612232B (en) | Power distribution network line re-jump probability prediction optimization method and device based on gradient descent | |
| CN112180314A (en) | Anti-interference self-correction GIS electronic transformer acquisition method and device | |
| CN110190617B (en) | Evaluation method, system, device and storage medium for multi-feed-in direct current power system | |
| CN104849684A (en) | Oscilloscope, correcting unit, and method of automatically correcting horizontal center of oscilloscope | |
| CN114386510A (en) | Method and system for identifying measurement errors of power system | |
| CN115879040A (en) | Capacitor failure prediction method, capacitor failure prediction device, electronic device, and storage medium | |
| CN113051741B (en) | Method and device for solving reliability parameters of high-voltage direct-current transmission system element | |
| CN113036716B (en) | Generator winding fault protection method and device | |
| CN108767821B (en) | Switching method and system for phase selection switch-on | |
| CN112421616A (en) | Method and system for identifying device parameters in main network of power grid | |
| CN112231981A (en) | Method for establishing large-scale electromagnetic transient simulation example | |
| CN114123164B (en) | Method and device for calculating short-circuit capacity of power system and terminal equipment | |
| CN110661261A (en) | PMU data processing method and device for low-frequency oscillation online analysis |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |