CN109256769A - A kind of uncertainty transient stability evaluation in power system method - Google Patents
A kind of uncertainty transient stability evaluation in power system method Download PDFInfo
- Publication number
- CN109256769A CN109256769A CN201811135263.8A CN201811135263A CN109256769A CN 109256769 A CN109256769 A CN 109256769A CN 201811135263 A CN201811135263 A CN 201811135263A CN 109256769 A CN109256769 A CN 109256769A
- Authority
- CN
- China
- Prior art keywords
- state
- transient
- power system
- generator
- equation
- 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.)
- Granted
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000011156 evaluation Methods 0.000 title claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 16
- 230000004044 response Effects 0.000 claims abstract description 10
- 238000009795 derivation Methods 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims 1
- 238000000547 structure data Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 206010034719 Personality change Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005183 dynamical system Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention relates to a kind of uncertain transient stability evaluation in power system methods, belong to power system stability field.Method of the invention is: the basic structure data of input electric power system first, obtains the electromechanical transient equation of description generator, and further obtain the state equation of its vector form；Determine the disturbance quantity λ of system；The uncertain transient stability evaluation in power system model for introducing uncertain variable is established, relevant response coefficient is calculated using the method for section Taylor expansion, obtains the bound of quantity of state x in transient state process of electric power system.Show that method proposed by the present invention can effectively assess uncertain electric power system transient stability in one machine infinity bus system emulation.
Description
Technical field
The present invention relates to a kind of uncertain transient stability evaluation in power system methods, belong to power system stability field.
Background technique
Modern power systems development is large-scale regional internet power grid, and this development can bring huge economy,
Also transient stability problem can be made more complicated simultaneously.Transient state unstability is still the biggest threat that modern power systems face
One of.The effective prediction of transient stability in real time and emergency flight control are most important.Traditional Transient Stability Control strategy is mainly adopted
With the control method of " table of generating strategy offline, real-time matching ".And there are some deviations for system model and parameter, to influence to count
Calculate the precision of result.Currently, after common transient stability analysis method has time-domain-simulation, time-domain-simulation method to can be used for device design
The inspection of phase stability is regarded as the standard of other method for analyzing stability inspection, but there are computationally intensive, imitative for simulation method
The true time is long and the disadvantages of can not emulating all operating statuses, therefore is unable to get the quantitative information of system degree of stability.Separately
A kind of main Transient Stability Analysis method is the direct method based on modern differentiable dynamical system, and this method, which can be furtherd investigate, is
The mechanism of system transient stability, has high learning value.Direct method is without integrating post-fault system, and by by fault clearance
The system transient modelling energy at moment is compared with transition energy, the transient stability of direct decision-making system.And with intelligent algorithm
As a result will for the temporary steady analysis method of guidance if being that there are errors between online data and preset data the shortcomings that this method
It is undesirable.
With the extensive access of power electronics, new energy etc., the operation of electric system is provided with very big uncertainty, because
This traditional Power System Analysis method be not applicable in.For the uncertain input amount often occurred in electric system and disturbance
The influence that variable generates power system transient stability, system need to capture probabilistic input quantity, analyze it
Influence to electric system this paper presents a kind of uncertain electrical power system transient appraisal procedure based on Taylor series expansion.
Summary of the invention
To solve existing deficiency, the present invention provides a kind of uncertain transient stability evaluation in power system method.
The technical scheme is that the basic structure data of input electric power system first, obtain the machine of description generator
Electric transient state equation, and further obtain the state equation of its mathematic(al) representation；Establish the uncertainty based on section Taylor series
Transient stability evaluation in power system model；Establish the disturbance quantity of system variable；Relevant response coefficient is calculated, electric system is obtained
Characteristic quantity bound in transient process.Show that method proposed by the present invention can be assessed effectively in one machine infinity bus system emulation
Uncertain electric power system transient stability.
A kind of uncertainty transient stability evaluation in power system method, includes the following steps:
Step 1 establishes the electromechanical transient equation of generator, and further obtains the electromechanical transient state side of its vector form
Journey；
Step 2 determines the disturbance quantity λ of system；
Step 3 establishes the electromechanical transient state equation for introducing disturbance quantity λ；
Step 4 calculates relevant response coefficient using the method for section Taylor expansion, obtains in transient state process of electric power system
The bound of quantity of state x: min (max) x (t, λ) (1)
Wherein x is the quantity of state of system；T is the time；λ is the disturbance quantity of quantity of state x.
Detailed process is as follows for the step 1:
Establish the electromechanical transient equation of following generator:
Wherein, δiIt is the generator rotor angle of i-th generator, ωiIt is the angular speed of i-th generator；ω0It is opposite angular speed；Di
It is the damped coefficient of i-th generator；PmiAnd PeiIt is the mechanical output and electromagnetic power of i-th generator, M respectivelyiIt is i-th
The mechanical output rotary inertia of generator；
The electromechanical transient state equation of its vector form is obtained according to above formula:
Wherein, x=[δ, ω]T；δ=[δ1,...δi,...,δm]；ω=[ω1,...ωi,...,ωm], m is total hair
Number of motors.
Detailed process is as follows for the step 3:
The disturbance quantity of input system:
λ=[λ1,...λj,...,λn]T (4)
Wherein, t changes quantity of state x at any time, λ1... λj... λnFor the disturbance quantity of the quantity of state at corresponding each moment, n
For the number of disturbance quantity,
Establish the electromechanical transient state equation introduced under disturbance quantity:
Detailed process is as follows for the step 4:
Establish the state equation of the first derivative containing probabilistic quantity of state x under the Taylor series expansion of section:
Establish the state equation containing probabilistic quantity of state x under the Taylor series expansion of section:
The further derivation of above formula is obtained:
Comprehensive above formula obtains:
Respective items are equal in this formula obtains:
The bound that t changes at any time during the quantity of state x dynamic response of system is obtained accordingly, wherein xcIt is system
Quantity of state when normal operation, λcIt is corresponding xcDisturbance quantity.
The beneficial effects of the present invention are:
1) this paper presents a kind of novel uncertain electric power system transient stabilities based on section Taylor series expansion
Appraisal procedure, this method can solve the problems, such as to be introduced by the intermittent renewable energy and measurement error.
2) calculation method of the present invention is simple, and precision of prediction is high, can be used for system transient modelling unstability early warning, further increases unstability
The cost that the timeliness of differentiation, as far as possible reduction subsequent control measure are paid.
Detailed description of the invention
Fig. 1 one machine infinity bus system and its basic structure data；
Fig. 2 introduces the generator's power and angle upper lower limit value after uncertainty.
Specific embodiment
The present invention will be further explained below with reference to the attached drawings and specific examples.
Embodiment 1: a kind of uncertainty transient stability evaluation in power system method, the specific steps of the method for discrimination are such as
Under:
Step 1 establishes the electromechanical transient equation of generator, and further obtains the electromechanical transient state side of its vector form
Journey；
Step 2 determines the disturbance quantity λ of system；
Step 3 establishes the electromechanical transient state equation for introducing disturbance quantity λ；
Step 4 calculates relevant response coefficient using the method for section Taylor expansion, obtains in transient state process of electric power system
The bound of quantity of state x: min (max) x (t, λ) (1)
Wherein x is the quantity of state of system；T is the time；λ is the disturbance quantity of quantity of state x.
Detailed process is as follows for the step 1:
Establish the electromechanical transient equation of following generator:
Wherein, δiIt is the generator rotor angle of i-th generator, ωiIt is the angular speed of i-th generator；ω0It is opposite angular speed；Di
It is the damped coefficient of i-th generator；PmiAnd PeiIt is the mechanical output and electromagnetic power of i-th generator, M respectivelyiIt is i-th
The mechanical output rotary inertia of generator；
The electromechanical transient state equation of its vector form is obtained according to above formula:
Wherein, x=[δ, ω]T；δ=[δ1,...δi,...,δm]；ω=[ω1,...ωi,...,ωm], m is total hair
Number of motors.
Detailed process is as follows for the step 3:
The disturbance quantity of input system:
λ=[λ1,...λj,...,λn]T (4)
Wherein, t changes quantity of state x at any time, λ1... λj... λnFor the disturbance quantity of the quantity of state at corresponding each moment, n
For the number of disturbance quantity,
Establish the electromechanical transient state equation introduced under disturbance quantity:
Detailed process is as follows for the step 4:
Establish the state equation of the first derivative containing probabilistic quantity of state x under the Taylor series expansion of section:
Establish the state equation containing probabilistic quantity of state x under the Taylor series expansion of section:
The further derivation of above formula is obtained:
Comprehensive above formula obtains:
Respective items are equal in this formula obtains:
The bound that t changes at any time during the quantity of state x dynamic response of system is obtained accordingly, wherein xcIt is system
Quantity of state when normal operation, λcIt is corresponding xcDisturbance quantity.
Embodiment 2: as shown in Figs. 1-2, proposed in this paper based on phase path and the knowledge of the stability of transient state energy in order to verify
Other method is verified using MATLAB software in one machine infinity bus system.Wherein the basic parameter of system is as follows:
A reference value is selected, SB=250MVA, UB (220)=209kV is chosen.
E'=1.5164
δ0=27.75 °
Wherein SB is the reference power of system, and UB is the reference voltage of system,For the equivalent reactance of generator unit stator, x2
For negative sequence neactance, xL1、xL2For the reactance per unit value of route 1 and route 2, xT1、xT2The respectively reactance of transformer 1 and transformer 2
Per unit value, TJFor time constant,For transient potential, δ0For generator rotor angle value.
Step S1, the state equation of the mathematic(al) representation of generator transient process is described
Step S2, the disturbance parameter that system is arranged is the fluctuation of mechanical output generation ± 5%.
Step S3, the state equation of the mathematic(al) representation of the generator transient process introduced under uncertain variables is established；
Step S4, obtain System describe system character changes over time bound, as shown in Figure 2.
Fig. 2 show when mechanical output occur ± 5% fluctuation when one machine infinity bus system angle stability curve and
The up-and-down boundary of dynamic response, it can be seen that power-angle curve when system stable operation dynamic response up-and-down boundary curve it
Interior, when mechanical output changes, minor change occurs for generator rotor angle versus time curve, and the operation of electric system is by crowd
The influence of more uncertain factors, for these uncertainties, the transient stability evaluation in power system based on section Taylor expansion
Method can accurately judge the transient stability of electric system.
Above in conjunction with attached drawing, the embodiment of the present invention is explained in detail, but the present invention is not limited to above-mentioned
Embodiment within the knowledge of a person skilled in the art can also be before not departing from present inventive concept
It puts and makes a variety of changes.
Claims (4)
1. a kind of uncertainty transient stability evaluation in power system method, characterized by the following steps:
Step 1 establishes the electromechanical transient equation of generator, and further obtains the electromechanical transient state equation of its vector form；
Step 2 determines the disturbance quantity λ of system；
Step 3 establishes the electromechanical transient state equation for introducing disturbance quantity λ；
Step 4 calculates relevant response coefficient using the method for section Taylor expansion, obtains state in transient state process of electric power system
Measure the bound of x: min (max) x (t, λ), x are the quantity of states of system；T is the time；λ is the disturbance quantity of quantity of state x.
2. uncertainty transient stability evaluation in power system method according to claim 1, it is characterised in that: the step
1 detailed process is as follows:
Establish the electromechanical transient equation of following generator:
Wherein, δiIt is the generator rotor angle of i-th generator, ωiIt is the angular speed of i-th generator；ω0It is opposite angular speed；DiIt is i-th
The damped coefficient of platform generator；PmiAnd PeiIt is the mechanical output and electromagnetic power of i-th generator, M respectivelyiIt is i-th power generation
The mechanical output rotary inertia of machine；
The electromechanical transient state equation of its vector form is obtained according to above formula:
Wherein, x=[δ, ω]T；δ=[δ1,...δi,...,δm]；ω=[ω1,...ωi,...,ωm], m is total generator
Quantity.
3. uncertainty transient stability evaluation in power system method according to claim 1, is characterized in that: the step 3
Detailed process is as follows:
The disturbance quantity of input system:
λ=[λ1,...λj,...,λn]T
Wherein, t changes quantity of state x at any time, λ1... λj... λnFor the disturbance quantity of the quantity of state at corresponding each moment, n is to disturb
The number of momentum,
Establish the electromechanical transient state equation introduced under disturbance quantity:
4. uncertainty transient stability evaluation in power system method according to claim 1, is characterized in that: the step 4
Detailed process is as follows:
Establish the state equation of the first derivative containing probabilistic quantity of state x under the Taylor series expansion of section:
Establish the state equation containing probabilistic quantity of state x under the Taylor series expansion of section:
The further derivation of above formula is obtained:
Comprehensive above formula obtains:
Respective items are equal in this formula obtains:
The bound that t changes at any time during the quantity of state x dynamic response of system is obtained accordingly, wherein xcIt is that system is normally transported
Quantity of state when row, λcIt is corresponding xcDisturbance quantity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811135263.8A CN109256769B (en) | 2018-09-28 | 2018-09-28 | Transient stability evaluation method for uncertain power system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811135263.8A CN109256769B (en) | 2018-09-28 | 2018-09-28 | Transient stability evaluation method for uncertain power system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109256769A true CN109256769A (en) | 2019-01-22 |
| CN109256769B CN109256769B (en) | 2021-08-20 |
Family
ID=65047995
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811135263.8A Active CN109256769B (en) | 2018-09-28 | 2018-09-28 | Transient stability evaluation method for uncertain power system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109256769B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009120190A1 (en) * | 2008-03-26 | 2009-10-01 | The Tokyo Electric Power Company, Incorporated | Stable equilibrium point (sep) calculation apparatus of power system |
| CN102915471A (en) * | 2012-09-26 | 2013-02-06 | 中国电力科学研究院 | Wind power disturbance online safety precaution method |
| CN104158191A (en) * | 2014-09-09 | 2014-11-19 | 武汉大学 | Dispersed coordination control method for stable running of multi-machine power system |
| CN104217090A (en) * | 2013-06-04 | 2014-12-17 | 刘光晔 | Method for analyzing and calculating voltage stabilization critical point of electric power system through Taylor series |
| CN105305425A (en) * | 2015-10-20 | 2016-02-03 | 国家电网公司 | UPFC buffeting control method of additional interference observer |
-
2018
- 2018-09-28 CN CN201811135263.8A patent/CN109256769B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2009120190A1 (en) * | 2008-03-26 | 2009-10-01 | The Tokyo Electric Power Company, Incorporated | Stable equilibrium point (sep) calculation apparatus of power system |
| CN102915471A (en) * | 2012-09-26 | 2013-02-06 | 中国电力科学研究院 | Wind power disturbance online safety precaution method |
| CN104217090A (en) * | 2013-06-04 | 2014-12-17 | 刘光晔 | Method for analyzing and calculating voltage stabilization critical point of electric power system through Taylor series |
| CN104158191A (en) * | 2014-09-09 | 2014-11-19 | 武汉大学 | Dispersed coordination control method for stable running of multi-machine power system |
| CN105305425A (en) * | 2015-10-20 | 2016-02-03 | 国家电网公司 | UPFC buffeting control method of additional interference observer |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109256769B (en) | 2021-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gomez et al. | Support vector machine-based algorithm for post-fault transient stability status prediction using synchronized measurements | |
| Resende et al. | Application of dynamic equivalence techniques to derive aggregated models of active distribution network cells and microgrids | |
| CN107132772A (en) | Real-time simulation system and method for AC/DC power grid | |
| Zhang et al. | Synchrophasor measurement-based wind plant inertia estimation | |
| CN104992017B (en) | Meter and the broad sense load modeling method and device of wind power random fluctuation | |
| Shabani et al. | Real-time transient stability detection in the power system with high penetration of DFIG-based wind farms using transient energy function | |
| CN103033716B (en) | Calculation method of proportion of each lode component in power grid comprehensive load model | |
| CN111130135B (en) | Power system inertia calculation method suitable for high-proportion new energy access | |
| Meliopoulos et al. | PMU data characterization and application to stability monitoring | |
| Bai et al. | Automatic modeling and optimization for the digital twin of a regional multi-energy system | |
| CN108054768A (en) | Transient stability evaluation in power system method based on principal component analysis | |
| He et al. | Power system frequency situation prediction method based on transfer learning | |
| CN108988320A (en) | Electrical Power System Dynamic element responds characteristic is to Enhancement of Transient Voltage Stability impact analysis method | |
| La Scala et al. | A qualitative approach to the transient stability analysis [of power systems] | |
| CN104779613B (en) | Test-based equivalent modeling method for electric element comprising converter | |
| CN109256769A (en) | A kind of uncertainty transient stability evaluation in power system method | |
| CN103279641A (en) | Method for conducting multi-modal dynamic equivalence on complex electrical power systems except for generators | |
| Cokkinides et al. | Visualization and characterization of stability swings via GPS-synchronized data | |
| Le-Thanh et al. | Test bench for self-healing functionalities applied on distribution network with distributed generators | |
| CN113113908A (en) | Time domain analysis method and system suitable for frequency response of modern large power grid | |
| Xinhui et al. | A data-driven slow dynamic characteristic extraction and state estimation method for large power grid | |
| Abdelaziz et al. | Phasor measurement units for out of step detection | |
| CN108964122A (en) | A kind of verification method and device of Wind turbines electrical model frequency response characteristic | |
| CN109088413A (en) | A kind of direct current receiving end power grid dynamic reactive Site Selection method and apparatus | |
| Senyuk et al. | A PMU-based algorithm of synchronous generator stability prediction during a disturbance |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |