CN106026083B - A kind of Power Network Transient Stability analysis method based on MATLAB - Google Patents
A kind of Power Network Transient Stability analysis method based on MATLAB Download PDFInfo
- Publication number
- CN106026083B CN106026083B CN201610479004.1A CN201610479004A CN106026083B CN 106026083 B CN106026083 B CN 106026083B CN 201610479004 A CN201610479004 A CN 201610479004A CN 106026083 B CN106026083 B CN 106026083B
- Authority
- CN
- China
- Prior art keywords
- equation
- generator
- analysis
- matlab
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention belongs to Simulating technique in Electric Power System field more particularly to a kind of Power Network Transient Stability analysis methods based on MATLAB.Using following step: determining electricity net safety stable analysis and calculation model；Power balance and trend mode select；Tidal current analysis；Power system transient stability analysis；The method has given full play to the features such as MATLAB programming efficiency height, good reliability, strong real-time, fining emulation can be carried out to the control method of electric system complexity, fast implement the customized of model, scalability is strong, it is clear, just with implementation, can satisfy the needs of electric system under the new situation.
Description
Technical field
The invention belongs to Simulating technique in Electric Power System field more particularly to a kind of Power Network Transient Stabilities based on MATLAB point
Analysis method.
Background technique
Electric system is a complicated dynamical system, at every moment among continuous disturbance, disturb generation when
Between, place, type, seriousness have randomness, once stability problem occurs in the system dynamic course after disturbance occurs,
In several seconds serious consequence may occur for system, cause great economic loss and social influence.Therefore, it is advised in electric system
It draws, design, operation etc. need to carry out a large amount of transient stability analysis in work.
But with the continuous expansion of power grid scale, the equation order of transient stability problem is higher and higher, and new element is continuous
Power grid is accessed, more stringent requirements are proposed to electric system simulation calculating.Some power system simulation softwares are in simulation calculation speed
Degree, new equipment model etc. have been unable to meet the needs of current electric grid development.
With the quick raising that electric system simulation fining requires, correlative study is weighed by university and scientific research institution
Depending on, but the research at present all concentrates on above theoretical and small test, and there are no the electric system for putting into actual motion.
Summary of the invention
Technical problem to be solved by the present invention lies in a kind of Power Network Transient Stability analysis method is provided, simulation calculation is improved
Efficiency and reliability can be realized the customized of complicated control methods and model, and scalability is strong, can satisfy electric system
Needs under the new situation.
The invention is realized in this way
A kind of Power Network Transient Stability analysis method based on MATLAB, using following step:
Determine electricity net safety stable analysis and calculation model；
Power balance and trend mode select；
Tidal current analysis；
Power system transient stability analysis；
It is characterized by: power system transient stability analysis includes the steps that are as follows:
1) according to initial value y before trend method of operation calculation perturbation(0): the system-wide mathematical model of multilayer output feedback network includes
The differential equation and electric power networks equation of the electric system in relation to element dynamic characteristic, the related element dynamic characteristic of electric system are described
The differential equation and the general type of electric power networks equation be
In formula: x, y are respectively state variable and algebraic variable；Algebraic variable y indicates the operating parameter of electric power networks, for electricity
The variable that power network equation solves, state variable x are the variable of differential equation,Indicate the differential form of state variable x, g
(x, y) indicates electric power networks equation；
2) by initial value y before disturbing(0)Find out the initial value x in system in relation to element state variable(0)；
3) algebraic equation of differential equation and electric power networks in relation to element dynamic characteristic is formed；
4) t=0 at the time of setting system operation；
5) check there is failure-free operation in t moment system；
6) if so, then modifying the related equation in the differential equation or algebraic equation in formula (1)；
7) electric power networks equation g (x, y)=0 is solved, acquires the operating parameter y of system after disturbance(t+0), state variable x(t)No
It can be mutated；
8) (x is used(t), y(t+0)) or (x(t), y(t)) substitute into formula 1. in the differential equation and electric power networks equation solution；
9) according to the differential equation and electric power networks equation joint solving result and judgment basis, judge whether system is steady
It is fixed；
10) t=t+ Δ t is enabled, calculates the operating status of t+ time Δt system, time interval Δ t is known as step-length, and transient state is steady
Devise a stratagem takes 0.05s in calculating, if meeting error: output generator power-angle curve, node voltage curve；
11) determine whether the time reaches previously given maximum moment tmax；
12) calculated result is exported.
Further, in step 5, if undisturbed, operating procedure 8.
Further, in step 9, such as system unstability has been distinguished, has then terminated operation, return step 12 exports calculated result.
Further, in step 11, as the time does not reach previously given maximum moment tmax, then: return step 5, inspection
Whether there is or not failed operation outputs for etching system when looking into this.
Further, it is included in simulink in tidal current analysis and establishes area power grid main electrical scheme Load flow calculation simulation model
Figure, reads Eexcel table and TXT text file content in MATLAB, and input power grid various parameters carry out Load flow calculation, use
MATLAB works out Newton-Raphson approach or PQ decomposition method flow calculation program, carries out preliminary Load flow calculation.
Further, comprising:
The static stability analysis under each trend mode is completed with MATLAB establishment " N-1 " checking routine, according to transmission system
Static security standard " N-l " criterion carries out line energizing flow ability and voltage level analysis.
Further, comprising:
Static stability analysis is carried out to design level year engineering location, static stability analysis is static using transmission system
Safety standard " N-l " criterion, on the basis of selected trend mode, fault-free open-circuit line, transformer single element one by one,
Load flow calculation is carried out again, obtains the trend distribution after N-1 is cut-off.
Compared with prior art, the present invention beneficial effect is:
1) the Power Network Transient Stability analysis method based on MATLAB that the present invention provides a kind of, the method give full play to
The features such as MATLAB programming efficiency height, good reliability, strong real-time, can carry out the control method of electric system complexity fine
To change emulation, fast implements the customized of model, scalability is strong, and it is clear, just with implementation, it can satisfy the new shape of electric system
Needs under gesture.
2) before the operation of newly-built project of transmitting and converting electricity, security and stability analysis is carried out to it using this method, it is possible to find engineering is thrown
Postpartum there are the problem of, adopt an effective measure in time, ensure safe and stable operation of power system.
Detailed description of the invention
Fig. 1 is the Power Network Transient Stability analysis method implementation flow chart the present invention is based on MATLAB；
Fig. 2 is Power Network Transient Stability analysis flow chart diagram provided in an embodiment of the present invention；
Fig. 3 (a) is the Eq and P when connecting between discussion generator and infinitely great system through complex network by taking salient-pole machine as an exampleEWith
The relationship of its dependent variable；
Fig. 3 (b) is the equivalent circuit between a salient pole generator and infinitely great system, and generator can only be withAnd xqAs
Its electromotive force and equivalent reactance；
Fig. 3 (c) is Y11…Y22For node admittance matrix element；
Fig. 3 (d) is y10、y12、y20For the admittance in equivalent π circuit；
Fig. 4 is IEEE9 bus test system figure；
Fig. 5 is generator 21 in 9 meshed network routes generation three phase short circuit fault, and the time that failure occurs cuts when being 1 second
The power-angle curve obtained except three phase short circuit fault；
Fig. 6 is generator 31 in 9 meshed network routes generation three phase short circuit fault, and the time that failure occurs is 1 second,
Three phase short circuit fault is cut off at 1.5 seconds, the power-angle curve obtained；
Fig. 7 is that 21 time of failure of generator is 1 second, and three phase short circuit fault occurs for route, and the mute time is 2.5 seconds
Power-angle curve；
Fig. 8 is that 31 time of failure of generator is 1 second, and three phase short circuit fault occurs for route, and the mute time is 2.5 seconds
Power-angle curve.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
Step 1 is to determine electricity net safety stable analysis and calculation model in Fig. 1, builds regional network framework for the later period and does standard
It is standby.System main models include: generator (containing excitation system, speed-regulating system model), load, route, transformer, direct current, nothing
The models such as function compensation.In MATLAB/Simulink, generator model mainly includes simplified Synchronous Machine Models, substantially same
Step machine model, standard synchronisation motor model, can satisfy the requirement of Transient-State Analysis of Power System different accuracy.Model of power transmission system
There are " π " equivalent module and distribution parameter equivalence module.Load model has static load model, dynamic load model and asynchronous electric
Machine model can simulate actual load characteristic to greatest extent.
Simpowersystems includes various component models needed for circuit simulation in Simulink, as power module,
The modules such as tandem circuit module, electric power electronic module, motor module, line-connecting machine module, detection module and secondary power module
Library.It include various primary element models in each module library, user can according to need oneself for carrying out electric system complex model
The fining emulation of circuit and electric system is realized in definition.
Model is chosen, and to consider element model, power system capacity used in real system, the factors such as simulation accuracy, with power generation
For machine utility model is chosen, utility model is chosen according to generator type: for the non-salient pole turbine generator with solid rotor
Machine should choose 6 ranks, 4 ranks and 2 rank models；5 ranks, 3 ranks and 2 rank models should be chosen for the hydrogenerator with field spider,
Wherein Eq ', Eq 〞, Ed 〞, 5 ranks of Ed ' potential change and 6 rank models consider time transient state and transient process (detailed model), examine
Consider 3 rank models of Eq ' potential change and consider Eq ', 4 rank models of Ed ' potential change only consider transient process, consider Eq ',
2 rank models of Ed ' constant potential and 2 rank model approximations of consideration E ' constant potential consider the model of adjustment of field excitation.In engineering reality
In trampling, general heavy-duty generator group is changed detailed using the secondary transient state (Eq 〞, Ed 〞) of Damper Winding the considerations of field measurement
Model, the generating set of intermediate size using consideration transient potential (Eq ', Ed ') variation model E, adopt by the generating set of low capacity
With transient potential (Eq ', Ed ') constant model.
In building rack, each parameter of input system element, needs that operation can be simplified by Excel table when reduction, mention
High efficiency.
Step 2 is that power balance and trend mode select in Fig. 1, determines designs quality year, regional rack, to area
Load is predicted: because area may have planning in the project of transmitting and converting electricity built, it is thus determined that designs quality year, ability is really
Determine the network architecture after engineering is gone into operation, and then regional load growth could be predicted.The method of load prediction mainly has bullet
Property coefficient method, unit consumption method, regression analysis etc..It can use different methods to be predicted, comprehensive ratio carried out to prediction result
Compared with, obtain basic, normal, high three kinds of schemes, according to different situations use different load prediction schemes.
Determine installation planning, spare capacity etc.: according to the load level of system prediction, necessary spare capacity and station-service
Installed capacity needed for electricity, network loss capacity determine system is horizontal.The generating equipment capacity that system needs should be that system is comprehensive most
The sum of big load and capacity needed for station service in the comprehensive spare capacity of system and system and network loss, i.e. equipment level=prediction is negative
Lotus Pmax+ spare capacity+station service, network loss.Spare capacity is divided into spare load, emergency duty by spare effect and overhauls standby
With.Load is spare generally to take the 2%~5% of system peak load；Emergency duty generally take system peak load 5%~
10%, and cannot be less than the capacity of a maximum unit in system；Overhaul reserve viewing system yearly load curve low ebb area it is big
Small, taking repair reserve of capacity is the 8%~15% of the comprehensive peak load of system；
Power balance is carried out to area, works out power balance table: according to area installation planning and load prediction, carrying out area
220kV layers of electrical reticulation design forcasted years and following grid power balance, work out power balance table.
System wiring is set for the rough sledding that actual capabilities in system operation occur according to the purpose for calculating analysis
And the method for operation, selection is likely to occur respectively from the following three method of operation (normal mode, maintenance mode, particular form)
To system safety and stability adverse conditions, calculating analysis is carried out.The selection of trend mode is made every effort to sufficiently expose after engineering is gone into operation
Safety and stability problem existing for uniting.
Step 3 is static system safety analysis in Fig. 1.Before carrying out static security analysis, route in system will be clear that, become
Depressor parameter, model and the thermostabilization limit, as whether overloading the foundation of judgement.It is provided according to related regulation, 220kV substation
220kV bus, when normal operating mode, voltage tolerance be system nominal voltage -3%~+7% (213.4kV~
235.4kV), emergency operating mode when for system nominal voltage -5%~+10% (209kV~242kV), as judging bus
Whether voltage is in tolerance range.
Static security analysis uses transmission system static security standard " N-l " criterion, i.e., in the trend mode base studied
On plinth, check that the system after disconnecting the single element (including route, transformer etc.) in electric system one by one carries out Load flow calculation,
After judging that N-I is cut-off, apparatus of load current value, system busbar voltage value, whether the load that route is carried be safely operated
Range etc.It for DC transmission system, is then latched with transient stability simulated program analog DC line monopole, observes direct current
Route whether there is safety problem after transient state and dynamic process.Generally route, transformer are carried out in engineering practice " single
Element fault " and " the element simultaneous faults of paired running " etc. are verified, and judge each in power grid after each element fault set
Whether standby to overload, whether system busbar voltage is out-of-limit.
Step 4 is power system transient stability analysis in Fig. 1, and time-domain-simulation method is that current development is most mature, most popular
Transient stability analysis method is widely used in electric system research, planning, operation, design etc..Multilayer output feedback network total system
Mathematical model include description the differential equation and electric power networks equation of the electric system in relation to element dynamic characteristic.Its general type
It can be written as
In formula: x, y are respectively state variable and algebraic variable；Algebraic variable y indicates that the operating parameter of electric power networks is electricity
The variable that power network equation solves, such as node admittance matrix, system node voltage, each generator node Injection Current.Temporary
The transient process of electric power networks is not considered in state calculating generally, therefore, y can be mutated in calculating process.State variable x is micro-
Divide the variable of equation solution, mainly including electromotive force in rotor angle δ, rotational speed omega, generator etc..
It is excessively not cumbersome to make to discuss, it is assumed that after short circuit, the excitation voltage of exciter is stood under the effect of automatic field forcing device
It is raised to maximum value uffmax.There are four due to the effect of encouraging by force, after short circuit, the generator differential equation to be solved is total to: first is that exciter
The differential equation, i.e.,
Second is that the excitation winding differential equation, i.e.,
Other two is the equation of rotor motion of generator:
Wherein Tff、Td′o、TJFor time constant, Eq' it is transient internal voltage, EqeFor the electromotive force of synchronous generator, EqemFor
Corresponding uffmaxNo-load electromotive force, Eq is no-load electromotive force, and δ is generator amature angle, ω0For system benchmark angular frequency, ω
For system angular frequency, PTFor the mechanical output of generator prime machine output, PEFor the electromagnetic power of generator output.
By taking salient-pole machine as an example, discuss between generator and infinitely great system through Eq when complex network connection (shown in such as Fig. 3 (a))
And PEWith the relationship of its dependent variable.Fig. 3 (b) is the equivalent circuit between a salient pole generator and infinitely great system, and generator is only
It can be withAnd xqAs its electromotive force and equivalent reactance.Y in Fig. 3 (c)11...Y22For node admittance matrix element, in Fig. 3 (d)
y10、y12、y20For the admittance in equivalent π circuit.The two has following relationship Y11=y10+y12；Y22=y20+y21；Y12=Y21=-y12
=-y21。
According to EQDefinition, EQFor generator-motor gesture, then EQAnd Eq、Eq' relationship be
I in formuladFor dynamo currentIn the component of d axis.If hypothesis q axis is in vertical axis direction, i.e., Hanging down
Axis direction, then IdForComponent on a horizontal axis, i.e.,
Substitution formula (5)
E in subtractive (7)QE can be obtainedqWith Eq', the relationship of δ.
The power equation of generator are as follows:
If it is non-salient pole machine, due to EQ=Eq, therefore have
Non-salient pole machine power equation is rewritten as
In formula, YijTransadmittance (G between generator-motor gesture node i and jij+jBij), | Yij| it is YijModulus value；δij
For the voltage of node iWith the voltage of node jAngle between phasor, i.e. δi-δj；δiAnd δjFor electromotive forceWithOpposite Mr. Yu
The angle of one reference phasor；βijFor
Step 1 and step 2 are to carry out initial value calculating.Initial value calculating includes dynamic to all generators and its excitation system, original
The calculating of machine and its speed-regulating system and load correlated variables initial value.
Before multilayer output feedback network, first has to carry out Load flow calculation according to the trend method of operation of selection, determine disturbance
Preceding operation state finds out the operating parameter y of electric power networks(0), and thus find out the state parameter in electric system in relation to element
Initial value x(0)。x(0)It mainly include Eqe(i.e. Eq|0|)、Eq′|0|、δ|0|、ω|0|.Their namely Eqe、Eq', the initial value of δ and ω.
These state variables are not mutated before and after disturbance.
Step 3 is to form the algebraic equation of differential equation and electric power networks in relation to element dynamic characteristic.From step 4
Start to enter the recurrence calculation for disturbing operation states of electric power system change procedure later, sets t=0 first, it is assumed that transient stability
Calculating have progressed to t moment, variable x(t), y(t)It, can E after calculation perturbation using formula (7), formula (8) for known quantityQ、EqWith
PE, these three variables are mutation before and after disturbance.Then differential equation (2), (3), (4) right-hand vector are it is known that t=0 can be acquired
When Eqe、Eq', the change rate of δ and ω.
Step 5 first checks for having failure-free operation in t moment system, if so, then enter step the 6 modification differential equations or
Related equation in algebraic equation, and then calculate the operating parameter y of electric power networks after failure(t+0), state variable x(t)It will not dash forward
Become, then uses (x(t), y(t+0)) substitute into the differential equation and algebraic equation solving.Otherwise directly by (x(t), y(t)) substitute into solution.
Step 8 carries out the solution of the differential equation and algebraic equation, and MATLAB provides many algorithms for ODE
Numerical solution, can be using Runge-Kutta method, improved Euler method, trapezoidal integration etc..The differential equation and algebra side can be used
Journey alternating or the method for simultaneous solution calculate the x of t+ time Δt(t+Δt)And y(t+Δt).Specifically it is achieved in that solution differential
Equation (2), (3), (4) calculate the E after Δ t according to the recurrence formula of selected method of value solvingqe、Eq', δ and ω
Estimated value recycles formula (6), (7) to calculate the E after Δ tQ、EqAnd PEEstimated value.It is similar to acquire correspondence using result above
The status variation rate of estimated value when Δ t.The E after Δ t is calculated further according to recurrence formulaqe、Eq', the corrected value of δ and ω.In this way,
After extrapolating t+ time Δt system running state, 10 are entered step, judges whether system is stable according to criterion appropriate, such as
Distinguish system unstability, then terminates operation.Otherwise, continue the operating status of subsequent time system, until previously given maximum
Moment tmax.The change curve of δ (and its dependent variable)~t so can be obtained.The method of the above analysis automatic adjustment excitation system is complete
It is suitable for automatic speed regulation system entirely, only needs additional notes PTThe differential equation of variation, then simultaneous solution.Step 5 in Fig. 1
For control measure, different control measure are taken according to the different faults of appearance, the relatively low measure taken of node voltage is installation
Reactive power compensator, increase power factor or cutting load etc..The measure that transformer, circuit overload are taken is cutting load.It is synchronous
The measure that stable problem is taken is to cut machine.
Embodiment: matlab program is write by the data that Fig. 4 is provided, Fig. 4 is IEEE9 bus test system figure.9 node nets
Three phase short circuit fault occurs for winding thread road, and the time that failure occurs is 1 second, cuts off three phase short circuit fault at 1.5 seconds, the function obtained
Angular curve is as shown in Figure 5, Figure 6.Using generator 1 as reference, the function with other two generators 21 and generator 31 is obtained
Angular difference.It can be seen from the curve obtained after breaking down, a series of concussions occur for system.After failure removal, system shake
The trend for presenting and constantly decaying is swung, final system is attributed to the state before failure occurs substantially.Therefore, it is possible to determine that go out system transient modelling
It is stable.
Time of failure is 1 second, and three phase short circuit fault occurs for route, and the power-angle curve that the mute time is 2.5 seconds is as schemed
7, shown in Fig. 8.
The power-angle curve as obtained in figure does not cut off failure it is found that after system jam in time, leads to generator
Power-angle curve diverging, can not return to original or new steady s tate.Therefore, it is possible to determine that Transient Instability occurs for system.
Claims (7)
1. a kind of Power Network Transient Stability analysis method based on MATLAB, using following step:
Determine electricity net safety stable analysis and calculation model；
Power balance and trend mode select；
Tidal current analysis；
Power system transient stability analysis；
It is characterized by: power system transient stability analysis includes the steps that are as follows:
1) according to initial value y before trend method of operation calculation perturbation(0): the system-wide mathematical model of multilayer output feedback network includes description
The differential equation and electric power networks equation of the electric system in relation to element dynamic characteristic, electric system are micro- in relation to element dynamic characteristic
The general type of point equation and electric power networks equation is
In formula: x, y are respectively state variable and algebraic variable；Algebraic variable y indicates the operating parameter of electric power networks, is power network
The variable of network equation solution, state variable x are the variable of differential equation,The differential form of expression state variable x, g (x,
Y) electric power networks equation is indicated；
2) by initial value y before disturbing(0)Find out the initial value x in system in relation to element state variable(0)；
3) algebraic equation of differential equation and electric power networks in relation to element dynamic characteristic is formed；
4) t=0 at the time of setting system operation；
5) check there is failure-free operation in t moment system；
6) if so, then modifying the related equation in the differential equation or algebraic equation in formula (1)；
7) electric power networks equation g (x, y)=0 is solved, acquires the operating parameter y of system after disturbance(t+0), state variable x(t)It will not dash forward
Become；
8) (x is used(t), y(t+0)) or (x(t), y(t)) substitute into formula (1) in the differential equation and electric power networks equation solution；
It include: the solution differential equation
According to the recurrence formula of selected method of value solving, the E after Δ t is calculatedqe、E′q, δ and ω estimated value, then benefit
Use formula
Calculate the E after Δ tQ、EqAnd PEEstimated value；Using result above, similar to the state of estimated value when acquiring corresponding Δ t
Change rate；The E after Δ t is calculated further according to recurrence formulaqe、E′q, δ and ω corrected value；
Wherein Tff、T′do、TJFor time constant, E 'qFor transient internal voltage, EqeFor the electromotive force of synchronous generator, EqemFor correspondence
uffmaxNo-load electromotive force, uffmaxFor exciter excitation voltage under the effect of automatic field forcing device the maximum value that is raised to immediately,
EqFor no-load electromotive force, δ is generator amature angle, ω0For system benchmark angular frequency, ω is system angular frequency, PTFor generator
The mechanical output of prime mover output, PEFor the electromagnetic power of generator output, EQFor generator-motor gesture, IdForIn horizontal axis
On component, Y11...Y22For node admittance matrix element；
9) according to the differential equation and electric power networks equation joint solving result and judgment basis, judge whether system is stable；
10) t=t+ Δ t is enabled, calculates the operating status of t+ time Δt system, time interval Δ t is known as step-length, transient stability meter
0.05s is taken in calculation, if meeting error: output generator power-angle curve, node voltage curve；
11) determine whether the time reaches previously given maximum moment tmax；
12) calculated result is exported.
2. according to the method for claim 1, which is characterized in that in step 5, if undisturbed, operating procedure 8.
3. according to the method for claim 1, which is characterized in that in step 9, such as distinguishes system unstability, then terminated operation,
Return step 12 exports calculated result.
4. according to the method for claim 1, which is characterized in that in step 11, as the time do not reach it is previously given most
Big moment tmax, then: return step 5, whether there is or not failed operation outputs for etching system when checking this.
5. according to the method for claim 1, which is characterized in that be included in simulink in tidal current analysis and establish area electricity
Host's wiring Load flow calculation simulation model figure reads Eexcel table and TXT text file content in MATLAB, inputs power grid
Various parameters carry out Load flow calculation, work out Newton-Raphson approach or PQ decomposition method flow calculation program with MATLAB, carry out just
Walk Load flow calculation.
6. according to the method for claim 5 characterized by comprising
The static stability analysis under each trend mode is completed with MATLAB establishment " N-1 " checking routine, according to transmission system static state
Safety standard " N-l " criterion carries out line energizing flow ability and voltage level analysis.
7. according to the method for claim 4 characterized by comprising
Static stability analysis is carried out to design level year engineering location, static stability analysis uses transmission system static security
Standard " N-l " criterion, on the basis of selected trend mode, fault-free open-circuit line, transformer single element one by one, then into
Row Load flow calculation obtains the trend distribution after N-1 is cut-off.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610479004.1A CN106026083B (en) | 2016-06-27 | 2016-06-27 | A kind of Power Network Transient Stability analysis method based on MATLAB |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610479004.1A CN106026083B (en) | 2016-06-27 | 2016-06-27 | A kind of Power Network Transient Stability analysis method based on MATLAB |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106026083A CN106026083A (en) | 2016-10-12 |
| CN106026083B true CN106026083B (en) | 2019-02-19 |
Family
ID=57084553
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610479004.1A Expired - Fee Related CN106026083B (en) | 2016-06-27 | 2016-06-27 | A kind of Power Network Transient Stability analysis method based on MATLAB |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106026083B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107202971B (en) * | 2017-04-24 | 2019-07-23 | 国网浙江省电力公司电力科学研究院 | Unconventional low voltage electric network electric energy measuring equipment operation characteristic simulation testing device |
| CN108365629B (en) * | 2018-03-16 | 2020-03-24 | 邓宏伟 | Electric power system electromechanical transient simulation method based on direct algorithm |
| CN110135031A (en) * | 2019-04-30 | 2019-08-16 | 东南大学 | Electric power system transient stability calculation method based on half implicit runge kutta method |
| CN113409643A (en) * | 2021-06-28 | 2021-09-17 | 江苏核电有限公司 | Transformer substation secondary circuit simulation training method and system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104240151A (en) * | 2014-09-11 | 2014-12-24 | 广西大学 | Transient stability optimal correcting and control system and method for power system |
| CN104732033A (en) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | Simulation analysis method of electromechanical transient process of electrical power system |
| CN104852394A (en) * | 2015-05-07 | 2015-08-19 | 陈鸽 | Power load balancing system applied to low-voltage power distribution network |
-
2016
- 2016-06-27 CN CN201610479004.1A patent/CN106026083B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104240151A (en) * | 2014-09-11 | 2014-12-24 | 广西大学 | Transient stability optimal correcting and control system and method for power system |
| CN104732033A (en) * | 2015-03-31 | 2015-06-24 | 河南行知专利服务有限公司 | Simulation analysis method of electromechanical transient process of electrical power system |
| CN104852394A (en) * | 2015-05-07 | 2015-08-19 | 陈鸽 | Power load balancing system applied to low-voltage power distribution network |
Non-Patent Citations (3)
| Title |
|---|
| 用于在线暂态稳定分析的电网实时运行状态获取方法研究;陈永稳;《电工技术》;20110310;第36-43页 |
| 风电场接入电网的安全稳定分析;张红光;《中国电力》;20070505;第105-110页 |
| 风电系统动态稳定分析与计算;王晓东;《中国优秀硕士学位论文全文数据库(电子期刊) 工程科技Ⅱ辑》;20040215;C042-118 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106026083A (en) | 2016-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Du et al. | Modeling of grid-forming and grid-following inverters for dynamic simulation of large-scale distribution systems | |
| Karimi et al. | A new centralized adaptive underfrequency load shedding controller for microgrids based on a distribution state estimator | |
| Kalcon et al. | Small-signal stability analysis of multi-terminal VSC-based DC transmission systems | |
| Van et al. | Virtual synchronous generator: An element of future grids | |
| Kundur | Power system stability | |
| Soultanis et al. | A stability algorithm for the dynamic analysis of inverter dominated unbalanced LV microgrids | |
| Phurailatpam et al. | Estimation of non-synchronous inertia in AC microgrids | |
| CN106026083B (en) | A kind of Power Network Transient Stability analysis method based on MATLAB | |
| Bangash et al. | Investigation of energy storage batteries in stability enforcement of low inertia active distribution network | |
| Collados-Rodriguez et al. | Stability and operation limits of power systems with high penetration of power electronics | |
| Kirakosyan et al. | Selective frequency support approach for MTDC systems integrating wind generation | |
| Vournas et al. | Application of QSS modeling to stabilizer design for interarea oscillations | |
| Elizondo et al. | Model reduction, validation, and calibration of wind power plants for dynamic studies | |
| Bharati et al. | Role of distributed inverter-based-resources in bulk grid primary frequency response through HELICS based SMTD co-simulation | |
| Mohammed et al. | Grid-forming power inverters: control and applications | |
| Erlich et al. | Selective damping of inter area oscillations using phasor measurement unit (PMU) signals | |
| Vergara et al. | Feasibility and performance assessment of commercial PV inverters operating with droop control for providing voltage support services | |
| Miranda et al. | Impact of energy storage on island grid dynamics: A case study of Bonaire | |
| Osika | Stability of micro-grids and inverter-dominated grids with high share of decentralised sources | |
| Dicorato et al. | Voltage compensation for wind integration in power systems | |
| Panwar et al. | Significance of dynamic and transient analysis in the design and operation of hybrid energy systems | |
| Teigenes | Provision of Primary Frequency Control from Electric Vehicles in the Nordic Power System | |
| Jacobsen et al. | Synchronous islanding control in the presence of hydro-electric generation | |
| Eyisi et al. | Load Sharing Strategy Incorporating Power Limits in Islanded Inverter-Based Microgrids | |
| Toigo | Investigation of stability aspects in a Pan-European interconnected grid with different wind power penetration |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190219 Termination date: 20200627 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |