CN102545216B - Projection method for generator node voltage during electric power system transient stability simulation process - Google Patents
Projection method for generator node voltage during electric power system transient stability simulation process Download PDFInfo
- Publication number
- CN102545216B CN102545216B CN 201210017819 CN201210017819A CN102545216B CN 102545216 B CN102545216 B CN 102545216B CN 201210017819 CN201210017819 CN 201210017819 CN 201210017819 A CN201210017819 A CN 201210017819A CN 102545216 B CN102545216 B CN 102545216B
- Authority
- CN
- China
- Prior art keywords
- generator
- node voltage
- generator node
- constantly
- electric power
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000001052 transient effect Effects 0.000 title claims abstract description 41
- 238000004088 simulation Methods 0.000 title abstract description 12
- 230000010354 integration Effects 0.000 claims description 28
- 238000002347 injection Methods 0.000 claims description 11
- 239000007924 injection Substances 0.000 claims description 11
- 238000005530 etching Methods 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 230000002265 prevention Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013213 extrapolation Methods 0.000 description 1
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a projection method for a generator node voltage during an electric power system transient stability simulation process. Compared with the prior art, the generator node voltage at a next integral step time can be more accurately projected through the time generator node voltage and the correlative state variable. The basic concept adopted by the invention is that when the generator node voltage is required to be projected, the whole electric power system is separately regarded as an electric generating set and an external system represented by equivalent admittance which are connected; the system equivalent admittance outside each generator is approximately regarded as to maintain invariability in an integral step; and a time state variable estimation value is projected through an explicit integral, and the predominant value of each generator node voltage can be obtained according to an approximate system equivalent impedance corresponding to each generator andthe state variable projection value. The generator node voltage projected by the method provided by the invention serves as the initial estimation value of each integral step simulation, and the iteration frequency during the solution process can be remarkably reduced.
Description
Technical field
The invention belongs to Automation of Electric Systems, particularly be used for the method for prediction generator node voltage in a kind of electric power system transient stability simulation process.
Background technology
Transient stability analysis of power system is one of content basic, the most most crucial during power system analysis is calculated.Online dynamic security analysis, the advanced technologies such as safety and stability emergency control, prevention and control, intelligent scheduling progressively require to apply in electric power system.The precondition that realizes these advanced technologies is can be fast, accurately, the transient stability emulation reliably large-scale electrical power system carried out.
Numerical integration is accurate, the most reliable method of electric power system transient stability simulation calculation.The disadvantage of numerical integrating is that amount of calculation is large, although computer speed has had raising at full speed, for large-scale electrical power system, still is difficult to satisfy the requirement that online dynamic security analysis, prevention and control, emergency control etc. are calculated computing time.
In existing electric power system transient stability numerical integration method, in order to guarantee the numerical stability of algorithm, usually adopt implicit integration algorithm.Implicit integration algorithm is in basis
State variable constantly
With the operation variable
Value is asked for
State variable constantly
With the operation variable
During value, need first given
The constantly discreet value of state variable and operation variable
,
, pass through again iterative.When considering the generator saliency, usually also need iterative when finding the solution network equation, namely need to provide the initial discreet value of generator and load bus voltage.If can provide good initial estimate, can reduce iterations, thereby significantly reduce the amount of calculation of transient stability emulation.
At present, usually use explicit numerical integration algorithm for state variable and just can provide preferably discreet value, and for the active node voltage in the system, particularly the generator node voltage then is difficult to dope more accurately discreet value.Initial estimate for each integration step generator node voltage in the transient stability emulation is provided by two kinds of methods usually: a kind of is with the preceding integration calculated value in step, as this integration initial estimate in step, namely gets
Another kind is the calculated value of getting former integration steps, obtains by the curve extrapolation, namely gets
Matched curve
Usually obtained by least square method or generalized extended algorithm,
Usually get 1~3.The method of above two kinds of initial estimates, in each integration step state variable and operation variable change hour, effect still can.But, when variable change is larger, will significantly increase iterations, sometimes not even convergence.The second estimation algorithm particularly, sometimes even poorer than first method effect.
Along with the continuous expansion of electric power system scale, the fault that may occur also becomes increasingly complex, and causes that thus system loses the loss of stablizing and causing also more serious.This just can carry out dynamic security analysis, implement safety and stability prevention and control and emergency control electric power system rapidly in the urgent need to us, and everything all is to be calculated as the basis with Fast numerical integration transient stability.In the transient stability computational process, the good predict of each integration step initial value is one of key technology that realizes the fast transient stability Calculation.
Summary of the invention
The present invention seeks in order to solve in the electric power system transient stability simulation process, the initial discreet value of generator node voltage is inaccurate, so that the iterative number of times increases, amount of calculation is large, computational speed can not satisfy the shortcoming of the online calculation requirement of electric power system, invented in a kind of electrical power system transient computational process, can dope more exactly the method for the initial discreet value of each integration step generator node voltage, reduce so that the iterations in each integration step of power system transient simulation when calculating can have significantly, and then significantly improve simulation velocity.
The present invention seeks to be achieved through the following technical solutions: a kind of generator node voltage Forecasting Methodology for electric power system transient stability emulation is characterized in that comprising the steps:
Step 1: according to
The time etching system in each generator node voltage and Injection Current, be calculated as follows the corresponding whole external system Equivalent admittance of each generator:
In the formula,
Be respectively
Constantly the
Real part, the imaginary part of the real part of platform generator node voltage, imaginary part and node Injection Current;
Expression
Constantly with the
The whole external system Equivalent admittance that generator joins;
Step 2: use explicit integration scheme, by
Moment state variable
With the operation variable
Value is calculated
The initial estimate of state variable constantly
Obtain the initial estimate of each generator's power and angle and each electromotive force
The generator built-in potential
Different and different according to the selected model of generator, their corresponding value is referring to table 1;
Step 3: be calculated as follows
Each generator of moment generator node admittance corresponding with the state variable initial estimate:
In the formula,
Expression is resistance in the generator,
,
,
,
What represent is that the admittance relevant with the generator built-in potential exists
Initial discreet value constantly,
,
,
,
What represent is that the admittance relevant with the generator node voltage exists
Initial discreet value constantly,
,
Represent respectively each generator
Axle and
The reactance of axle, along with the used model of generator is different and different, their corresponding value can see Table 1.
Table 1 motor model and the parameter table of comparisons
In the table 1
What represent is generator transient state reactance after-potential,
,
,
,
Represent respectively generator
Axle,
The transient internal voltage of axle and time transient potential,
What represent is generator
The axle synchronous reactance,
,
,
,
Represent respectively generator
Axle,
Transient state reactance and the subtranient reactance of axle.
In the formula
Represent respectively
The platform generator exists
Real part, the imaginary part of moment node voltage initial estimate; The inventive method is Single-step Prediction, namely only uses
The voltage of generator node and relevant generator state variable just can dope next integration step more exactly constantly
Generator node voltage constantly.
The inventive method when prediction generator node voltage just according to the generator node voltage before the prediction constantly, but integrated application generator node voltage and relevant generator state variables before the prediction constantly.
The inventive method is when each generator node voltage of prediction, regard one by one whole electric power system as a generating set and an external system that represents with Equivalent admittance links, namely one by one whole complication system is become the single system that the system equivalent admittance forms outside a power supply and the generator.
The inventive method is in prediction during each generator node voltage, approximately thinks that the Equivalent admittance with the whole external system of each generator electric remains unchanged in the step at an integration.
The present invention only needs according to following simple forecast formula
Just can dope more exactly the generator node voltage.
The invention has the beneficial effects as follows: a generator node voltage Forecasting Methodology foundation
State variable constantly and operation variate-value just can be obtained more rapidly
The constantly discreet value of generator node voltage can reduce the iterations of network algebra equation solution in the transient stability emulation effectively, improves the speed of power system transient simulation.
Description of drawings
Fig. 1 is flow chart of the present invention;
Fig. 2 is the equivalent schematic diagram of the present invention;
Fig. 3 is the equivalent schematic diagram of unit corresponding to system.
Embodiment
The invention will be further described below in conjunction with accompanying drawing.
A kind of generator node voltage Forecasting Methodology for electric power system transient stability emulation that the present invention proposes may further comprise the steps:
Step 1: according to
The time etching system in each generator node voltage and Injection Current, be calculated as follows the corresponding whole external system Equivalent admittance of each generator:
In the formula,
Be respectively
Constantly the
Real part, the imaginary part of the real part of platform generator node voltage, imaginary part and node Injection Current;
Expression
Constantly with the
The whole external system Equivalent admittance that generator joins;
Step 2: use explicit integration scheme, by
Moment state variable
With the operation variable
Value is calculated
The initial estimate of state variable constantly
Obtain the initial estimate of each generator's power and angle and each electromotive force
The generator built-in potential
Different and different according to the selected model of generator, their corresponding value is referring to table 1;
Step 3: be calculated as follows
Each generator of moment generator node admittance corresponding with the state variable initial estimate:
In the formula,
Expression is resistance in the generator,
,
,
,
What represent is that the admittance relevant with the generator built-in potential exists
Initial discreet value constantly,
,
,
,
What represent is that the admittance relevant with the generator node voltage exists
Initial discreet value constantly,
,
Represent respectively each generator
Axle and
The reactance of axle, along with the used model of generator is different and different, their corresponding value can see Table 1.
Table 1 motor model and the parameter table of comparisons
In the table 1
What represent is generator transient state reactance after-potential,
,
,
,
Represent respectively generator
Axle,
The transient internal voltage of axle and time transient potential,
What represent is generator
The axle synchronous reactance,
,
,
,
Represent respectively generator
Axle,
Transient state reactance and the subtranient reactance of axle.
In the formula
Represent respectively
The platform generator exists
Real part, the imaginary part of moment node voltage initial estimate;
The inventive method is Single-step Prediction, namely only uses
The voltage of generator node and relevant generator state variable just can dope next integration step more exactly constantly
Generator node voltage constantly.
The inventive method when prediction generator node voltage just according to the generator node voltage before the prediction constantly, but integrated application generator node voltage and relevant generator state variables before the prediction constantly.
The inventive method is when each generator node voltage of prediction, regard one by one whole electric power system as a generating set and an external system that represents with Equivalent admittance links, namely one by one whole complication system is become the single system that the system equivalent admittance forms outside a power supply and the generator.
The inventive method is in prediction during each generator node voltage, approximately thinks that the Equivalent admittance with the whole external system of each generator electric remains unchanged in the step at an integration.
The present invention only needs according to following simple forecast formula
Just can dope more exactly the generator node voltage.
Below introduce in detail the detailed process of the inventive method.
As shown in Figure 2 whole electric power system is regarded as a generating set one by one and an external system that represents with Equivalent admittance links, i.e. whole electric power system can be expressed as respectively equivalent circuit shown in Figure 3.Then
The moment all can be tried to achieve with voltage and the Injection Current of each generator node by (1) formula with the Equivalent admittance of the whole external system of each generator electric:
In the formula,
Be respectively
Constantly the
Real part, the imaginary part of the real part of platform generator node voltage, imaginary part and node Injection Current;
Expression
Constantly with the
The whole external system Equivalent admittance that generator joins;
Use explicit integration scheme, by
Moment state variable
With the operation variable
Value is calculated
The initial estimate of state variable constantly
Obtain the initial estimate of each generator's power and angle and each electromotive force
The generator built-in potential
According to the difference of the selected model of generator and difference.
The constantly discreet value of the Injection Current of generator node can be shown at this state variable of constantly estimating and the voltmeter of estimating with generator:
In the formula (2)
Moment generator
Internal impedance under the coordinate transforms to
Admittance under the coordinate
,
But formula (3), formula (4) are asked for:
Simultaneously
The constantly discreet value of each generator node Injection Current also estimating with the external system Equivalent admittance of available this moment generator node voltage is tried to achieve:
Suppose that each generator exists
Constantly with
The Equivalent admittance approximately equal of valve system is waited in the outside constantly, namely has:
Then by (2)~(6) Shi Kede
The voltage of estimating of every generator node is constantly:
Each generator exists in the transient stability simulation process like this
The constantly voltage discreet value of each generator node, can try to achieve as follows:
1.According to
The time etching system in each generator node voltage and Injection Current, calculate the corresponding whole external system Equivalent admittance of each generator by (1) formula
2.Use explicit integration scheme, by
Moment state variable
With the operation variable
Value is calculated
The initial estimate of state variable constantly
Obtain the initial estimate of each generator's power and angle and each electromotive force
By calculating suc as formula (3), formula (4)
Each generator of moment generator node admittance corresponding with the state variable initial estimate:
,
Below be an embodiment of the inventive method, carry out emulation experiment with East China 529 machines 3918 systems and make embodiment, further specify as follows: all generators all adopt in the test macro
Constant model, load adopts the constant impedance model, and iteration error is
The generator node voltage Forecasting Methodology that this paper is proposed is used for implicit expression trapezoidal integration transient stability simulation calculation, and with the calculated value that adopts the preceding integration step as this integration initial estimate in step, namely get
Method compare.
Suppose near the SNQ node of the SNQ--SSJ branch road of East China 529 machines, 3918 node systems three phase short circuit fault occurs, continue 0.3s fault time after, with the faulty line excision, simulation time is 3 seconds.
The amount of calculation of table 1. pair two kinds of algorithms compares.As can be seen from Table 1, the number of times that adopts the present invention to predict that the generator node voltage can make transient emulation find the solution network equation reduces 52%~117%, and the CPU time saves 30%~44%, has reduced significantly the amount of calculation of transient stability emulation.
Claims (5)
1. be used for the generator node voltage Forecasting Methodology of electric power system transient stability emulation, it is characterized in that comprising the steps:
Step 1: each generator node voltage and Injection Current in the etching system during according to t are calculated as follows the corresponding whole external system Equivalent admittance of each generator:
In the formula, V
XGi(t), V
YGi(t), I
XGi(t), I
YGi(t) be respectively t constantly real part, the imaginary part of real part, imaginary part and the node Injection Current of i platform generator node voltage; G
Gi(t), B
Gi(t) represent the whole external system Equivalent admittance that the t moment and i platform generator join;
Step 2: use explicit integration scheme, by the value of t moment state variable x (t) and operation variable y (t), calculate the initial estimate x of t+h state variable constantly
[0](t+h); Obtain the initial estimate of each generator's power and angle
Initial estimate with built-in potential
The generator built-in potential
Different and different according to the selected model of generator, their corresponding value is referring to table 1;
Step 3: be calculated as follows constantly each generator generator node admittance corresponding with the state variable initial estimate of t+h:
In the formula, R
AiExpression is resistance in the generator,
Expression be the admittance relevant with the generator built-in potential in the t+h initial discreet value in the moment,
Expression be the admittance relevant with the generator node voltage in the t+h initial discreet value in the moment,
Represent respectively the reactance of each generator d axle and q axle, along with the used model of generator is different and different, their corresponding value can see Table 1;
Table 1 motor model and the parameter table of comparisons
What E' represented in the table 1 is generator transient state reactance after-potential, E'
d, E'
q, E "
d, E "
qThe transient internal voltage and the inferior transient potential that represent respectively generator d axle, q axle, X
qWhat represent is generator q axle synchronous reactance, X'
d, X'
q, X "
d, X "
qThe transient state reactance and the subtranient reactance that represent respectively generator d axle, q axle;
Step 4: be calculated as follows the constantly initial estimate of each generator node voltage of t+h
2. the generator node voltage Forecasting Methodology for electric power system transient stability emulation according to claim 1, it is characterized in that: the method is Single-step Prediction, and voltage and the relevant state variable of generator of generator node just can dope next integration step t+h generator node voltage constantly more exactly constantly namely only to use t.
3. the generator node voltage Forecasting Methodology for electric power system transient stability emulation according to claim 1, it is characterized in that: the method when prediction generator node voltage just according to the generator node voltage before the prediction constantly, but integrated application generator node voltage and relevant generator state variables before the prediction constantly.
4. the generator node voltage Forecasting Methodology for electric power system transient stability emulation according to claim 1, it is characterized in that: the method is when each generator node voltage of prediction, regarding one by one whole electric power system as a generating set and an external system that represents with Equivalent admittance links, namely is the single system that the system equivalent admittance forms outside a power supply and the generator one by one with whole complication system equivalence.
5. the generator node voltage Forecasting Methodology for electric power system transient stability emulation according to claim 1, it is characterized in that: the method is in prediction during each generator node voltage, thinks that the Equivalent admittance with the whole external system of each generator electric remains unchanged in the step at an integration.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210017819 CN102545216B (en) | 2012-01-19 | 2012-01-19 | Projection method for generator node voltage during electric power system transient stability simulation process |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210017819 CN102545216B (en) | 2012-01-19 | 2012-01-19 | Projection method for generator node voltage during electric power system transient stability simulation process |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102545216A CN102545216A (en) | 2012-07-04 |
CN102545216B true CN102545216B (en) | 2013-10-23 |
Family
ID=46351432
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201210017819 Expired - Fee Related CN102545216B (en) | 2012-01-19 | 2012-01-19 | Projection method for generator node voltage during electric power system transient stability simulation process |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102545216B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105375475B (en) * | 2015-11-26 | 2017-11-24 | 浙江大学 | A kind of power system fast transient stability simulation method based on Real-time Dynamic Equivalence |
CN107025335B (en) * | 2017-03-13 | 2019-12-10 | 清华大学 | simulation calculation method and simulation system based on state variable discretization |
CN112290550A (en) * | 2020-12-24 | 2021-01-29 | 中国电力科学研究院有限公司 | Emergency control method and system based on branch transient state power transmission capacity |
CN115659849A (en) * | 2022-12-09 | 2023-01-31 | 南方电网数字电网研究院有限公司 | Transient simulation method of power system based on data model hybrid drive |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4602288B2 (en) * | 2005-12-01 | 2010-12-22 | 中部電力株式会社 | Method and apparatus for determining transient stability of power system |
CN102157936A (en) * | 2011-03-31 | 2011-08-17 | 清华大学 | Coordination control method |
-
2012
- 2012-01-19 CN CN 201210017819 patent/CN102545216B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4602288B2 (en) * | 2005-12-01 | 2010-12-22 | 中部電力株式会社 | Method and apparatus for determining transient stability of power system |
CN102157936A (en) * | 2011-03-31 | 2011-08-17 | 清华大学 | Coordination control method |
Non-Patent Citations (5)
Title |
---|
JP特许4602288B2 2010.12.22 |
基于输电线路单端量测信息的电力***暂态稳定性预估;罗远翔等;《东北电力技术》;20100831(第8期);全文 * |
电力***实时等值及电压稳定性分析;邱晓燕等;《电网技术》;20040430;第28卷(第7期);全文 * |
罗远翔等.基于输电线路单端量测信息的电力***暂态稳定性预估.《东北电力技术》.2010,(第8期),全文. |
邱晓燕等.电力***实时等值及电压稳定性分析.《电网技术》.2004,第28卷(第7期),全文. |
Also Published As
Publication number | Publication date |
---|---|
CN102545216A (en) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102609575B (en) | Power system transient stability simulating method based on implicit numerical integration | |
CN102403720B (en) | Hyper-real-time setting method for superposed time sequences based on transient voltage safety margin | |
CN102175922B (en) | Phasor measurement unit (PMU) measurement data-based power line parameter identification and estimation method | |
CN102420427B (en) | Method for calculating voltage stability margin of regional power grid by considering outer network equivalence | |
CN104866665A (en) | Hybrid simulation method including power electronic equipment based on interface equivalence and interaction | |
CN102545216B (en) | Projection method for generator node voltage during electric power system transient stability simulation process | |
CN106374488A (en) | Fractional order terminal sliding mode-based AFNN control method of active power filter | |
CN108711846B (en) | AC/DC system long-term voltage stabilization model prediction control method | |
CN103606921A (en) | Transient voltage stabilization judgment method based on real-time generalized Thevenin equivalence | |
CN102545263B (en) | Power system transient stability simulation method based on explicit numerical integration | |
CN105184027A (en) | Power load modeling method based on interactive multi-model algorithm | |
Ma et al. | Application of dynamic equivalencing in large-scale power systems | |
CN105512502A (en) | Weight function least square state estimation method based on residual normalization | |
CN102508086A (en) | Parameter estimation method for electric power circuit in T-type connection | |
Kim et al. | Parameter estimation for the composite load model | |
CN105375475B (en) | A kind of power system fast transient stability simulation method based on Real-time Dynamic Equivalence | |
CN109245122B (en) | Power system regional dynamic equivalence method based on wide-area measurement information | |
Bila | Power system dynamic state estimation and load modeling | |
CN103825270B (en) | A kind of power distribution network three-phase state estimates the processing method of Jacobian matrix constant | |
CN105305392A (en) | Symmetrical component method for short circuit calculation of voltage-controlled type IIDG included power distribution network | |
CN103606952B (en) | A kind of machine of cutting control measure quantization method based on system acceleration energy | |
CN103944158B (en) | A kind of online Poewr control method based on small interference stability risk supervision and system | |
CN104779613B (en) | Test-based equivalent modeling method for electric element comprising converter | |
CN102136105A (en) | Phase measurement unit-based power grid information graph parameter estimation method | |
CN102798751B (en) | A kind of Novel voltage stability detection method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131023 Termination date: 20210119 |
|
CF01 | Termination of patent right due to non-payment of annual fee |