CN104680323B - Electric system energy stabilization domain builds system and method - Google Patents
Electric system energy stabilization domain builds system and method Download PDFInfo
- Publication number
- CN104680323B CN104680323B CN201510092163.1A CN201510092163A CN104680323B CN 104680323 B CN104680323 B CN 104680323B CN 201510092163 A CN201510092163 A CN 201510092163A CN 104680323 B CN104680323 B CN 104680323B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- energy
- stabilization domain
- generator
- 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
Abstract
A kind of electric system energy stabilization domain builds system and method.The system comprises data acquisition module, system dynamical equation generation module, damping capacities to solve module, energy stabilization domain structure module and result output module；Wherein, data acquisition module for gathering network architecture parameters, generator frequency, generator rotor angle in system；System dynamical equation generation module is used to obtain the dynamical equation of system；Damping capacity solves module and is used to ask for the damping capacity that generator armature reaction accumulates in a cycle oscillatory process in system；Energy stabilization domain structure module reflects the electric system energy stabilization domain of system dynamic stability and transient stability for building；As a result output module is used for output power system capacity stable region.System and method is built by the electric system energy stabilization domain of the present invention, can reflect the static stability and dynamic stability of system, there is higher accuracy.
Description
Technical field
The present invention relates to technical field of power systems more particularly to stability of power system analysis systems and method.
Background technology
With the fast development of modern power systems, grid equipment variation, information diversification, the method for operation changeableization etc.
The deficiencies of feature is further apparent, and deterministic type security conceives the information dualization exposed, model determinization has been difficult to meet people
For the more safe and reliable requirement of electric energy, under these circumstances, the analysis means in " domain " are increasingly becoming a kind of electric system
The important measures of safety and stability research.
According to the division that IEEE analyzes power system security, the stability of power system research for being currently based on domain is main
It is divided into Steady State Security Region (Steady-State Security Region) and Dynamic Security Region (Dynamic Security
Region).Steady State Security Region is defined in a set in node injecting power space, for any one note in the set
Incoming vector, it is necessary at the same meet trend can solve with two constraintss of small disturbance stability, therefore, Steady State Security Region is that trend is feasible
Domain and the intersection of Small Signal Stability Region.And Dynamic Security Region is in node injecting power space, system is still protected after making accident
Hold the set for being possible to operating point of transient stability.From above-mentioned definition can be seen that Steady State Security Region can analysis system it is steady
Security during state, but the considerations of lack to disturbance factor, and Dynamic Security Region can ensure the first pendulum stability of system, but lack
The weary analysis to system subsequent dynamic stability.As it can be seen that single Steady State Security Region or Dynamic Security Region can not analyze failure
How the dynamic stability of system afterwards, build that both both to take into account advantage and make up insufficient comprehensive safety domain be to solve the problems, such as this
Key.
The content of the invention
In view of this, the purpose of the present invention is intended to overcome deficiency of the prior art, avoid single Steady State Security Region or
Dynamic Security Region method can not analyze the problem of dynamic stability of post-fault system, play the work that structure takes into account the two advantage
With.
System is carried out linearisation expansion by the present invention first with the wide area signal data of reading at non-equilibrium, is obtained
System dynamical equation.Secondly, time-varying damping component in generator armature reaction process is defined as dynamic antivibration, it is basic herein
On, it takes into full account the time-varying factor in oscillatory process, asks for generator armature in system and react in a cycle oscillatory process
The damping capacity of accumulation.Then, in parameter-energy hyperspace, according to the electricity for considering disturbance size and system mode simultaneously
Force system dynamic stability criterion draws out the electric system energy that can intuitively reflect system dynamic stability and transient stability
Stable region.It is higher that simulation result based on four machines, two district system shows that the system capacity stable region constructed by the present invention has
Accuracy.Using the energy stabilization domain that the present invention is formed for system call is instructed to run, learn system potential risk, formulate
Control program ensures that system safe and stable operation is of great significance.
In order to realize this purpose, the technical solution that the present invention takes is as follows.
A kind of electric system energy stabilization domain builds system, the system comprises sequentially connected data acquisition module, is
Uniting, dynamical equation generation module, damping capacity solve module, energy stabilization domain builds module and result output module；
Wherein, data acquisition module for gathering network architecture parameters, generator frequency, generator rotor angle in system, and will acquisition
Data sending is to system dynamical equation generation module；
System dynamical equation generation module is used at non-equilibrium point be unfolded to obtain system by system using the data of reading
Dynamical equation；
Damping capacity solves module and is accumulated for asking for generator armature reaction in system in a cycle oscillatory process
Damping capacity；
Energy stabilization domain structure module is used for the power system dynamic stability criterion according to disturbance size and system mode, structure
Build the electric system energy stabilization domain of reflection system dynamic stability and transient stability；
As a result output module is used for output power system capacity stable region.
Wherein, the system dynamical equation of system dynamical equation generation module acquisition is：
Wherein, c=Ak -1f(xk,yk),
xkAnd ykThe quantity of state and algebraic quantity of system, Δ x at respectively point kkFor xkKnots modification.
In addition, the damping capacity that damping capacity solution module is asked for is：
Wherein,
Δδi、ΔωiRespectively i-th generator's power and angle variable quantity and rotation speed change amount,
KekAnd DekRespectively tkMoment generator synchronization factor and damped coefficient.
A kind of electric system energy stabilization domain construction method, the method includes the steps：
A, time lag system network architecture parameters are gathered, generator frequency, generator rotor angle in system；
B, system is unfolded at non-equilibrium point using the data of acquisition to obtain the dynamical equation of system；
C, the damping capacity that generator armature reaction accumulates in a cycle oscillatory process in system is asked for；
D, according to disturbance size and the power system dynamic stability criterion of system mode, structure reflection system dynamic stability
With the electric system energy stabilization domain of transient stability；
E, the electric system energy stabilization domain is exported.
The system dynamical equation obtained in stepb is：
Wherein, c=Ak -1f(xk,yk),
xkAnd ykThe quantity of state and algebraic quantity of system, Δ x at respectively point kkFor xkKnots modification..
In addition the damping capacity obtained in step C is：
Wherein,
Δδi、ΔωiRespectively i-th generator's power and angle and rotation speed change amount,
KekAnd DekRespectively tkMoment generator synchronization factor and damped coefficient.
Size and the power system dynamic stability criterion of system mode are disturbed in step D to be included：
System ensures the damping in equalization point for just under current operating condition；
After disturbance in first cycle, damping capacity and generator windings, machinery that generator armature reaction generates hinder
The sum of damping capacity that Buddhist nun generates is more than zero；
Transient state energy accumulation is no more than the max-thresholds that system can be born caused by failure.
Determine that electric system energy stabilization domain includes in step D：
D1, energy stabilization domain bottom surface D, the energy stabilization are determined using Small Signal Stability Region method for solving in parameter space
Domain bottom surface D meets system and ensures that the damping in equalization point is positive condition under current operating condition；
D2, any operating point in the bottom surface D of energy stabilization domain is obtained to step D1, which is determined by direct method
Transient Instability energy threshold E1, the Transient Instability energy threshold E1 meets transient state energy accumulation and is no more than system and can bear
Max-thresholds condition；
D3, the disturbance for applying different-energy size in the range of 0-E1 to operating point in step D2, ask for disturbing system
Damping is 0 energy value after dynamic, remembers the energy value for E2, and E2 meets after disturbance in first cycle, generator armature reaction
The sum of damping capacity that the damping capacity of generation is generated with generator windings, mechanical damping is more than zero condition；
All operating points on D4, traversal energy stabilization domain bottom surface D, repeat step D2 and D3, determine that all operating points correspond to
E1 and E2, the top surface S1 being made of all E1 is related to power system transient stability, if disturbance energy be higher than the top surface, system will
Transient Instability, if disturbance energy is less than the top surface, system will keep transient stability, the top surface S2 institutes structure being made of all E2
Into top surface it is related with the dynamic stability during system oscillation, if disturbance energy, higher than S2 and less than S1, system will be gradual
Divergence instability, if disturbance energy is less than S2, system oscillation process will gradually decay；
D5, the energy stabilization domain bottom surface D institutes encircling space determined in the top surface S2 and step D1 determined in step D4 is utilized
As energy stabilization domain.
In addition in step D1, energy stabilization domain bottom surface D is determined using point by point method or protection mapping method.
System and method is built by using the electric system energy stabilization domain of the present invention, can be unfolded at non-equilibrium point
The dynamical equation for obtaining system solves damping capacity, on this basis, it is steady to construct a kind of new comprehensive safety domain-energy
Localization.Electric system energy stabilization domain constructed by the present invention can reflect the static stability and dynamic stability of system,
With higher accuracy.Using the electric system energy stabilization domain that the present invention is formed for instructing system call operation, obtaining
System potential risk is known, control program is formulated, ensures that system safe and stable operation has advantageous effect.
Description of the drawings
Fig. 1 builds system structure diagram for electric system energy stabilization domain of the present invention.
Fig. 2 is the signal in the energy stabilization domain that electric system energy stabilization domain of the present invention is built constructed by system and method
Figure.
Fig. 3 is the two district system structure chart of IEEE4 machines applied in one embodiment of the present invention.
Fig. 4 is the two district system parameter space small disturbance stability of IEEE4 machines applied in one embodiment of the present invention
Domain schematic diagram.
Fig. 5 is identified energy stabilization domain in one embodiment of the present invention.
Fig. 6 is two district system energy stabilization domain section of the IEEE4 machines signal applied in one embodiment of the present invention
Figure.
Fig. 7 be section I section of two district system energy stabilization domain of IEEE4 machines applied in one embodiment of the present invention with
The time-domain-simulation figure of point in IV sections.
Fig. 8 is section II sections of two district system energy stabilization domain of IEEE4 machines applied in one embodiment of the present invention
The result in time domain figure of interior point.
Fig. 9 is section III sections of two district system energy stabilization domain of IEEE4 machines applied in one embodiment of the present invention
The result in time domain figure of interior point.
Specific embodiment
Below in conjunction with the accompanying drawings, elaborate to the present invention.
The detailed example embodiment of following discloses.However, concrete structure disclosed herein and function detail merely for the sake of
The purpose of example embodiment is described.
It should be appreciated, however, that the present invention is not limited to disclosed particular exemplary embodiment, but covering falls into disclosure model
Enclose interior all modifications, equivalent and alternative.In the description to whole attached drawings, identical reference numeral represents identical member
Part.
It will also be appreciated that term "and/or" includes the arbitrary of one or more relevant list items as used in this
With all combinations.It will further be appreciated that when component or unit are referred to as " connecting " or during " coupled " to another component or unit, it
It can be directly connected or coupled to other component or unit or there may also be intermediate member or units.In addition, for describing
Between component or unit other words of relation should understand in the same fashion (for example, " between " to " directly between ",
" adjacent " is to " direct neighbor " etc.).
In order to introduce technical scheme, illustrate the principle of the present invention first.
The model of electric system is represented by with differential-algebraic equation group：
In formula, x is quantity of state, and y is algebraic quantity.
Values of the x and y in each integration step starting point is obtained by analysis.X in making the integration step at certain point k
=xk+ Δ x, y=yk, differential equation f (x, y) is substituted into, can be obtained：
By formula (2) equation right side in (xk,yk) at carry out Taylor expansion：
When Δ x is sufficiently small, formula (3) can be neglected high-order term H.O.T. and obtain：
In formula：
Formula (4) is made into following deformation：
Make c=Ak -1f(xk,yk), it obtains：
Formula (6) is the system dynamical equation that system is unfolded at non-equilibrium point.
In addition, in one machine infinity bus system, the t moment generator amature equation of motion is as follows：
Wherein, Pm、PeRespectively act on the mechanical output and electromagnetic power on rotor, PmIt is determined by prime mover output, and
PeCome from the interaction for the electromagnetic field that stator and rotor current generates in magnetic gap.DMBe meter and D windings, Q windings in dynamic process
In damping action and mechanical damping in rotor motion after, the permanent damped coefficient of introducing.
In [tk,tk+1] in time step, calculated using second order Runge-Kutta method and inject generator by generator interior nodes
The oscillation energy of branch is：
Two parts are included by the energy of network injection generator branch：One is the variation of generator transient state energy, i.e.,Secondly it is the energy of Generator Damping consumption, i.e.,
With the electromagnetic power Δ P of the generator of Incremental Equation expressionekIt can be expressed as by Heffron-Phillips models：
ΔPek=KekΔδ+DekΔω (9)
Wherein, KekAnd DekRespectively tkMoment generator synchronization factor and damped coefficient.The present invention is by each moment generator
Armature-reaction is defined as the damping action size that generator amature moves the dynamic antivibration of generator, it is with electromagnetic power increment
In linear representation Δ ω it is related, in one machine infinity bus system, DekAs dynamic antivibration.
Due to DMPerseverance is positive value, and positive damping capacity is only sent in oscillatory process, to studying because of negative damping energy production
Unstable oscillation process do not contact directly, for ease of analysis, generator makes D using classical second-order modelM=0, at this point, being
The dynamical equation at non-equilibrium point k of uniting is：
In formula：
In formula (10), because generator armature reaction is only related with Δ ω to the damping action of generator amature movement, extraction
Part in the formula with Δ ω same-phases can obtain：
-TJΔ ω=Kekvδk+Dek(Δω+rωk) (12)
The damping energy that No. i-th generator armature reaction accumulates in a cycle oscillatory process in system is asked for by formula (12)
Measure EDiFor：
Knowable to analysis mode (13), the damping capacity E of certain generator after a failure in multi-computer systemDiExcept its in by system
The influence of his generator speed, it is also related with the generator rotor angle of other generators.After disturbance during system oscillation, as disturbance is big
Small difference, generator parameter KekAnd DekPositive and negative values ratio in a cycle damps energy also by difference so as to cause the machine
The positive negative variance of amount.If being vibrated at first in cycle, the damping capacity that generator generates is negative, in the movement of follow-up cycle
In, system is by gradual divergence instability；It is 0 that if the damping capacity that generator generates in cycle is vibrated at first, in follow-up cycle
Movement in, system will occur self-sustained oscillation phenomenon；If the damping capacity that generator generates in first cycle is just, it is
System will gradually tend towards stability in the movement of follow-up cycle.
Therefore, whether stablize under a certain disturbance or failure for judging a system, it is impossible to only dotted by system balancing
State judges, it is necessary to be codetermined with reference to disturbance size.In order to intuitively analyze different disturbance sizes comprehensively to system stability
It influences, the present invention will provide in parameter-energy space structure electric system energy stabilization domain for system stability judgement after disturbance
One quick and accurate instrument.
Energy stabilization domain is defined in all of parameter-energy space and can guarantee that system keeps dynamic stability after disturbance
Operating point-disturbance energy to set.Electric system being capable of following 3 conditions of safe and stable operation needs satisfaction under disturbance：
Condition 1：System must assure that the damping in equalization point is just, to ensure that system is being born under current operating condition
With certain stability margin without being to bear that Unstable oscillation phenomenon occurs because of damping during certain disturbance.
Condition 2：After disturbance in first cycle, the damping capacity and generator windings of generator armature reaction generation,
The sum of damping capacity that mechanical damping generates is more than 0, to ensure power system dynamic stability after disturbance.
Condition 3：Transient state energy accumulation is no more than the max-thresholds that system can be born caused by failure, to ensure system not
Transient Instability problem can occur.
Therefore, Fig. 1 builds system structure diagram for electric system energy stabilization domain of the present invention.It is as shown in Figure 1, of the invention
Electric system energy stabilization domain structure system includes sequentially connected data acquisition module in embodiment, system dynamical equation is given birth to
Module, energy stabilization domain structure module and result output module are solved into module, damping capacity；Wherein, data acquisition module is used
System dynamical equation generation mould is sent in generator frequency, generator rotor angle in acquisition network architecture parameters, system, and by gathered data
Block；System dynamical equation generation module is used to that system to be unfolded at non-equilibrium point using the data of reading to obtain the dynamic of system
Equation；Damping capacity solves module and is used to ask for the resistance that generator armature reaction accumulates in a cycle oscillatory process in system
Buddhist nun's energy；Energy stabilization domain structure module is used for the power system dynamic stability criterion according to disturbance size and system mode, structure
Build the electric system energy stabilization domain of reflection system dynamic stability and transient stability；As a result output module is used for output power
System capacity stable region.
Particularly, in the specific embodiment of the invention, the system dynamical equation of system dynamical equation generation module acquisition
For：
Wherein, c=Ak -1f(xk,yk),
xkAnd ykThe quantity of state and algebraic quantity of system at respectively point k.
In addition, in the embodiment of the invention, damping capacity solves the damping capacity that module is asked for and is：
Wherein,
Δδi、ΔωiRespectively i-th generator's power and angle variable quantity and rotation speed change amount,
KekAnd DekRespectively tkMoment generator synchronization factor and damped coefficient.
Present invention additionally comprises a kind of electric system energy stabilization domain construction method, the method includes the steps：
A, time lag system network architecture parameters are gathered, generator frequency, generator rotor angle in system；
B, system is unfolded at non-equilibrium point using the data of acquisition to obtain the dynamical equation of system；
C, the damping capacity that generator armature reaction accumulates in a cycle oscillatory process in system is asked for；
D, according to disturbance size and the power system dynamic stability criterion of system mode, structure reflection system dynamic stability
With the electric system energy stabilization domain of transient stability；
E, the electric system energy stabilization domain is exported.
Correspondingly, the system dynamical equation obtained in step B is：
Wherein, c=Ak -1f(xk,yk),
xkAnd ykThe quantity of state and algebraic quantity of system at respectively point k.
In addition, the damping capacity obtained in step C is：
Wherein,
Δδi、ΔωiRespectively i-th generator's power and angle and rotation speed change amount,
KekAnd DekRespectively tkMoment generator synchronization factor and damped coefficient.
It is corresponding with analyzing above, it is characterised in that size is disturbed in step D and the Electrical Power System Dynamic of system mode is steady
Determining criterion includes：
System ensures the damping in equalization point for just under current operating condition；
After disturbance in first cycle, damping capacity and generator windings, machinery that generator armature reaction generates hinder
The sum of damping capacity that Buddhist nun generates is more than zero；
Transient state energy accumulation is no more than the max-thresholds that system can be born caused by failure.
Correspondingly, determine that electric system energy stabilization domain includes in step D：
D1, energy stabilization domain bottom surface D, the energy stabilization are determined using Small Signal Stability Region method for solving in parameter space
Domain bottom surface D meets system and ensures that the damping in equalization point is positive condition under current operating condition；
D2, any operating point in the bottom surface D of energy stabilization domain is obtained to step D1, which is determined by direct method
Transient Instability energy threshold E1, the Transient Instability energy threshold E1 meets transient state energy accumulation and is no more than system and can bear
Max-thresholds condition；
D3, the disturbance for applying different-energy size in the range of 0-E1 to operating point in step D2, ask for disturbing system
Damping is 0 energy value after dynamic, remembers the energy value for E2, and E2 meets after disturbance in first cycle, generator armature reaction
The sum of damping capacity that the damping capacity of generation is generated with generator windings, mechanical damping is more than zero condition；
All operating points on D4, traversal energy stabilization domain bottom surface D, repeat step D2-D3, determine that all operating points are corresponding
E1 and E2.Fig. 2 is the schematic diagram in the energy stabilization domain that electric system energy stabilization domain of the present invention is built constructed by system and method,
As shown in Fig. 2, the top surface S1 being made of all E1 is related to power system transient stability, if disturbance energy is higher than the top surface, system
By Transient Instability, if disturbance energy is less than the top surface, system will keep transient stability, the top surface S2 institutes being made of all E2
The top surface of composition is related with the dynamic stability during system oscillation, if disturbance energy, higher than S2 and less than S1, system will be by
Gradually divergence instability, if disturbance energy is less than S2, system oscillation process will gradually decay；
D5, the energy stabilization domain bottom surface D institutes encircling space determined in the top surface S2 and step D1 determined in step D4 is utilized
As energy stabilization domain.
Particularly, in a detailed embodiment, in step D1, energy is determined using point by point method or protection mapping method
Stable region bottom surface D.
Illustrate the technique effect of the present invention followed by specific application example.
Two district system of IEEE4 machines as shown in Figure 3 is built first, using it as analysis foundation.
In system, generator uses six rank transient Models.Using selecting model analysis (Selective mode
Analysis, SMA) method is to high order system progress depression of order, and during depression of order, reservation is related to low-frequency oscillation and is easy to survey
The quantity of state of amount.Excitation system uses high-speed excitation, the load under benchmark model using 50% constant-impedance and 50% perseverance electricity
Flow mixed model.
No. 1 generator energy stable region is established according to electric system energy stabilization domain construction method of the present invention.
Step 1：Using on busbar 4 and 14 band burden with power P4 and P14 as parametric variable, asked for using mapping method is protected
Border of the equalization point damping more than 0 is as stable region bottom surface on the parameter space for system, as shown in Figure 4.
Step 2：To the operating point required by step 1 in stable region bottom surface, the transient state for determining the operating point by direct method is lost
Steady energy threshold, it is E1 to remember the energy.By taking (- 0.25,0.25) is put in stable region bottom surface required by step 1 as an example, direct method is utilized
This Transient Instability energy is obtained as 0.60MW*S, as shown in Figure 5.
Step 3：Apply the disturbance of different-energy size in the range of 0-E1 to the operating point, according to formula (13) searching make be
The value that the damping capacity after disturbance of uniting is 0, it is E2 to remember the energy.With (- 0.25,0.25) point in stable region bottom surface required by step 1
Exemplified by, the different disturbance energies between 0-0.60MW*S are applied to system, ask for sending out in the first cycle after disturbing using formula (13)
Motor damping energy size if damping capacity is less than 0, reduces initial disturbance energy, if damping capacity is more than 0, increases
Add initial disturbance energy, until hunting out the initial disturbance energy value that damping capacity is 0, acquired results 0.53MW*S.
Step 4：It to all operating points on required stable region bottom surface in step 1, is traveled through with 0.05 for step-length, repeats to walk
Rapid 2-3 asks for all operating points corresponding E1 and E2 respectively, as shown in the figure of Fig. 6 middle sections.
Step 5：Take E2 form energy top surface and step 1 in the space that is surrounded of stable region bottom surface as energy stabilization
Domain.
In order to further verify the correctness in constructed energy stabilization domain, analysis verification is carried out to the point inside and outside stable region.
Using P14=-0.25 sections as analysis object, as shown in Figure 6.Wherein, it is 4 sections by -1.5 to 1 points by P4, is respectively
I sections ([- 1.5, -0.75]), II sections ([- 0.75, -0.05]), III sections ([- 0.05,0.5]), IV sections ([0.5,1]), and
III sections are drawn top surface shown in E2, E1 respectively.In Section II section, top surface shown in E2, E1 overlaps.
To A in Fig. 6 I sectionsIA in point (- 0.25, -1,0.1) and sectionIV(- 0.25,0.6,0.1) point is emulated respectively,
The results are shown in Figure 7.Reflect in Fig. 7 be generator 1 and 3 generator rotor angle difference of generator changes with time situation.It can be with from figure
Find out, system diverging oscillation due to underdamping, at this point, system is operated in unsafe condition, thus be excluded that in energy stabilization domain
Outside.
To A in Fig. 6 II sectionsIIPoint (- 0.25, -0.4,0.3) and BIIPoint (- 0.25, -0.4,0.4) is emulated respectively,
The results are shown in Figure 8.What is reflected in figure is that the generator rotor angle difference of generator 1 and generator 3 changes with time situation.It can be with from figure
Find out, when system is in A in II sectionsIIDuring dotted state, after disturbing, damping capacity is just generator 1 and hair in oscillatory process
The generator rotor angle difference of motor 3 gradually tends towards stability, as shown in Fig. 8 (a).And when system is located at B in III sectionsIIDuring dotted state, because of oscillation energy
Amount is more than Transient Instability threshold value, and system is by direct Transient Instability after disturbance occurs, as shown in Fig. 8 (b).
To A in Fig. 6 III sectionsIIIPoint (- 0.25,0.25,0.48), BIIIPoint (- 0.25,0.25,0.0.53), CIIIPoint (-
0.25,0.25,0.58)、DIIIPoint (- 0.25,0.25,0.63) is emulated respectively, and the results are shown in Figure 9.What is reflected in figure is
Generator 1 and 3 generator rotor angle difference of generator change with time situation.
As seen from Figure 9, when system is in A in III sectionsIIIDuring dotted state, after disturbing, damped in oscillatory process
Energy is just, the generator rotor angle difference of generator 1 and generator 3 gradually tends towards stability, as shown in Fig. 9 (a).And positioned at B in III sectionsIIIPoint
After disturbing, damping capacity is just 0 in oscillatory process, and the generator rotor angle between generator 1 and generator 3 occurs constant amplitude and persistently shakes
It swings, as shown in Fig. 9 (b).When system running state is in C in III sectionsIIIDuring point, after disturbing, energy is damped in oscillatory process
It measures and is no more than Transient Instability threshold value for negative and oscillation energy, system is because damping capacity is negative ever-diverging oscillations, although the short time
Interior system is unlikely to unstability, but the process continues and triggers system transient modelling unstability by final, as shown in Fig. 9 (c).When system is transported
Row state is located at D in III sectionsIIIDuring point, because oscillation energy is more than Transient Instability threshold value, system is by direct transient state after disturbance occurs
Unstability, as shown in Fig. 9 (d).
It follows that there is higher accuracy using the system capacity stable region constructed by the present invention.
It should be noted that the above embodiment is only the preferable embodiment of the present invention, it is impossible to is understood as to this
The limitation of invention distance protection scope, under the premise of without departing from present inventive concept, to any minor variations for being done of the present invention with
Modification belongs to the distance protection scope of the present invention.
Claims (5)
1. a kind of electric system energy stabilization domain builds system, the system comprises sequentially connected data acquisition module, systems
Dynamical equation generation module, damping capacity solve module, energy stabilization domain structure module and result output module；
Wherein, data acquisition module for gathering network architecture parameters, generator frequency, generator rotor angle in system, and by gathered data
It is sent to system dynamical equation generation module；
System dynamical equation generation module is used to that system to be unfolded at non-equilibrium point using the data of reading to obtain the dynamic of system
State equation；
Damping capacity solves module and is used to ask for the resistance that generator armature reaction accumulates in a cycle oscillatory process in system
Buddhist nun's energy；
Energy stabilization domain structure module is used for the power system dynamic stability criterion according to disturbance size and system mode, and structure is anti-
Reflect the electric system energy stabilization domain of system dynamic stability and transient stability；
As a result output module is used for output power system capacity stable region；
The system dynamical equation that the system dynamical equation generation module obtains is：
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;x</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein, c=Ak -1f(xk,yk),
<mrow>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mi>k</mi>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
xkAnd ykThe quantity of state and algebraic quantity of system, Δ x at respectively point kkFor xkKnots modification；
The damping capacity solves the damping capacity that module is asked for：
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>T</mi>
</munderover>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&Delta;&delta;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>e</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>&delta;</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>e</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>&Delta;&omega;</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>&omega;</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein,
Δδi、ΔωiRespectively i-th generator's power and angle variable quantity and rotation speed change amount,
KekAnd DekRespectively tkMoment generator synchronization factor and damped coefficient.
2. a kind of electric system energy stabilization domain construction method, the method includes the steps：
A, time lag system network architecture parameters are gathered, generator frequency, generator rotor angle in system；
B, system is unfolded at non-equilibrium point using the data of acquisition to obtain the dynamical equation of system；
C, the damping capacity that generator armature reaction accumulates in a cycle oscillatory process in system is asked for；
D, according to disturbance size and system mode power system dynamic stability criterion, structure reflection system dynamic stability and temporarily
The electric system energy stabilization domain of state stability；
E, the electric system energy stabilization domain is exported；
Wherein, the system dynamical equation obtained in step B is：
<mrow>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;x</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein, c=Ak -1f(xk,yk),
<mrow>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>f</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mi>k</mi>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
xkAnd ykThe quantity of state and algebraic quantity of system, Δ x at respectively point kkFor xkKnots modification；
The damping capacity obtained in step C is：
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>D</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>T</mi>
</munderover>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msub>
<mi>&Delta;&delta;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>e</mi>
<mi>k</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>&delta;</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>e</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>&Delta;&omega;</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>&omega;</mi>
<mi>k</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
Wherein,
Δδi、ΔωiRespectively i-th generator's power and angle and rotation speed change amount,
KekAnd DekRespectively tkMoment generator synchronization factor and damped coefficient.
3. according to the electric system energy stabilization domain construction method described in claim 2, it is characterised in that disturbed in step D big
Small and system mode power system dynamic stability criterion includes：
System ensures the damping in equalization point for just under current operating condition；
After disturbance in first cycle, damping capacity and generator windings, mechanical damping that generator armature reaction generates are produced
The sum of raw damping capacity is more than zero；
Transient state energy accumulation is no more than the max-thresholds that system can be born caused by failure.
4. according to the electric system energy stabilization domain construction method described in claim 3, it is characterised in that determine electricity in step D
Force system energy stabilization domain includes：
D1, energy stabilization domain bottom surface D, energy stabilization domain bottom are determined using Small Signal Stability Region method for solving in parameter space
Face D meets system and ensures that the damping in equalization point is positive condition under current operating condition；
D2, any operating point in the bottom surface D of energy stabilization domain is obtained to step D1, the temporary of the operating point is determined by direct method
State unstability energy threshold E1, the Transient Instability energy threshold E1 meets transient state energy accumulation can bear most no more than system
The condition of big threshold value；
D3, the disturbance for applying different-energy size in the range of 0-E1 to operating point in step D2, ask for making system after disturbance
The energy value for 0 is damped, remembers the energy value for E2, E2 meets after disturbance in first cycle, and generator armature reaction generates
Damping capacity and generator windings, the sum of the damping capacity that generates of mechanical damping be more than zero condition；
All operating points on D4, traversal energy stabilization domain bottom surface D, repeat step D2 and D3, determine the corresponding E1 of all operating points
And E2, the top surface S1 being made of all E1 is related to power system transient stability, if disturbance energy is higher than the top surface, system is by transient state
Unstability, if disturbance energy is less than the top surface, system will keep transient stability, what the top surface S2 being made of all E2 was formed
Top surface is related with the dynamic stability during system oscillation, if disturbance energy, higher than S2 and less than S1, system will gradually dissipate
Unstability, if disturbance energy is less than S2, system oscillation process will gradually decay；
D5, by the use of the energy stabilization domain bottom surface D institute's encircling spaces determined in the top surface S2 and step D1 determined in step D4 as
Energy stabilization domain.
5. according to the electric system energy stabilization domain construction method described in claim 4, it is characterised in that in step D1, utilize
Point by point method or protection mapping method determine energy stabilization domain bottom surface D.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510092163.1A CN104680323B (en) | 2015-03-02 | 2015-03-02 | Electric system energy stabilization domain builds system and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510092163.1A CN104680323B (en) | 2015-03-02 | 2015-03-02 | Electric system energy stabilization domain builds system and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104680323A CN104680323A (en) | 2015-06-03 |
| CN104680323B true CN104680323B (en) | 2018-05-25 |
Family
ID=53315331
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510092163.1A Active CN104680323B (en) | 2015-03-02 | 2015-03-02 | Electric system energy stabilization domain builds system and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104680323B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104933636B (en) * | 2015-06-12 | 2019-02-05 | 华北电力大学 | Multi-state method for analyzing stability of power system and device |
| CN104950676B (en) * | 2015-06-12 | 2017-11-28 | 华北电力大学 | Time-varying power system self-adaptation control method and device |
| CN104915892B (en) * | 2015-06-12 | 2018-12-18 | 华北电力大学 | Electric system stochastic Time-Delay method for analyzing stability |
| CN107742892B (en) * | 2017-10-30 | 2019-07-12 | 国家电网公司 | A kind of energy storage damping control method inhibiting the oscillation of New-energy power system broadband |
| CN110311414B (en) * | 2019-07-30 | 2020-11-20 | 华北电力大学 | Doubly-fed wind turbine control parameter optimization method and system based on energy stability domain |
| CN111769575B (en) * | 2020-07-15 | 2021-09-28 | 华北电力大学 | Fan parameter optimization oscillation suppression system and method based on modal stability domain |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226795A (en) * | 2013-05-06 | 2013-07-31 | 中国南方电网有限责任公司 | Method for quantitatively evaluating damping of generator |
| CN103473478A (en) * | 2013-09-30 | 2013-12-25 | 电子科技大学 | Energy function-based assessment method for transient stability of grid |
| CN103812123A (en) * | 2014-02-28 | 2014-05-21 | 华北电力大学 | Time-lag stability control system considering jump characteristic of power system and method thereof |
-
2015
- 2015-03-02 CN CN201510092163.1A patent/CN104680323B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226795A (en) * | 2013-05-06 | 2013-07-31 | 中国南方电网有限责任公司 | Method for quantitatively evaluating damping of generator |
| CN103473478A (en) * | 2013-09-30 | 2013-12-25 | 电子科技大学 | Energy function-based assessment method for transient stability of grid |
| CN103812123A (en) * | 2014-02-28 | 2014-05-21 | 华北电力大学 | Time-lag stability control system considering jump characteristic of power system and method thereof |
Non-Patent Citations (3)
| Title |
|---|
| Power system small-signal stability region calculation method based on the guardian map theory;Jing Ma，Shangxing Wang，Zengping Wang，James S. Thorp;《IET Generation, Transmission & Distribution》;20140831;第8卷(第8期);第1479-1488页 * |
| 基于能量函数法的电网暂态稳定性分析;汪小明;《电网技术》;20110831;第35卷(第8期);第114-118页 * |
| 电力系统结构保留模型及上界函数法暂态稳定域构建的研究;刘峰;《中国优秀博硕士学位论文全文数据库 (博士) 工程科技Ⅱ辑》;20050915(第05期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104680323A (en) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104680323B (en) | Electric system energy stabilization domain builds system and method | |
| Yang et al. | Data-driven estimation of inertia for multiarea interconnected power systems using dynamic mode decomposition | |
| Chavan et al. | Identification and predictive analysis of a multi-area WECC power system model using synchrophasors | |
| Shah et al. | Direct quadrate (dq) modeling of 3-phase induction motor using matlab/simulink | |
| Liu et al. | Impedance modeling and controllers shaping effect analysis of PMSG wind turbines | |
| Bharati et al. | SMTD co-simulation framework with HELICS for future-grid analysis and synthetic measurement-data generation | |
| Huang et al. | Faster than real-time dynamic simulation for large-size power system with detailed dynamic models using high-performance computing platform | |
| Ma et al. | Model order reduction analysis of DFIG integration on the power system small‐signal stability considering the virtual inertia control | |
| Baek et al. | Nonlinear parameter optimization of FACTS controller via real-time digital simulator | |
| Abdulrahman et al. | Simulink-based programs for power system dynamic analysis | |
| CN104915892B (en) | Electric system stochastic Time-Delay method for analyzing stability | |
| Vournas et al. | Application of QSS modeling to stabilizer design for interarea oscillations | |
| CN108988320A (en) | Electrical Power System Dynamic element responds characteristic is to Enhancement of Transient Voltage Stability impact analysis method | |
| Kyesswa et al. | A Matlab-based dynamic simulation module for power system transients analysis in the eASiMOV framework | |
| Saborio Romano | Small-signal modelling and stability analysis of a traditional generation unit and a virtual synchronous machine in grid-connected operation | |
| Marjanian et al. | Designing of a TEF‐based neuro‐fuzzy controller for IPFC to improve the transient stability | |
| Feng et al. | Dynamic equivalencing of distribution network with high penetration of distributed generation | |
| Abdulrahman et al. | Simulink-based program for simulating multi-machine power systems | |
| CN110309625A (en) | A kind of energy stabilization domain of double-fed fan motor grid-connected system determines method and system | |
| Haghighat et al. | Mitigating communication delay impact on microgrid stability using a compensator based on smith predictor | |
| Concepcion et al. | On extended-term dynamic simulations with high penetrations of photovoltaic generation | |
| Argüello et al. | Dynamic modeling in OpenDSS: An implementation sequence for object Pascal | |
| Liao et al. | Measurement-based dynamic equivalent modeling for small and medium hydropower generator group | |
| Jakobsen et al. | Development of a test system for identification of turbine dynamics using the dc power flow | |
| Bravo‐López et al. | Simplified dynamic models for three‐phase microgrids |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |