CN110350596A - It is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin - Google Patents
It is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin Download PDFInfo
- Publication number
- CN110350596A CN110350596A CN201910542412.0A CN201910542412A CN110350596A CN 110350596 A CN110350596 A CN 110350596A CN 201910542412 A CN201910542412 A CN 201910542412A CN 110350596 A CN110350596 A CN 110350596A
- Authority
- CN
- China
- Prior art keywords
- synchronous generator
- transient stability
- bus nodes
- receiving end
- 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.)
- Pending
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005457 optimization Methods 0.000 title claims abstract description 35
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 230000035945 sensitivity Effects 0.000 claims abstract description 8
- 230000001360 synchronised effect Effects 0.000 claims description 67
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000005070 sampling Methods 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 3
- 210000004209 hair Anatomy 0.000 claims description 3
- 230000003466 anti-cipated effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/40—Synchronising a generator for connection to a network or to another generator
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- 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)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin, the receiving end power grid of this method operating status known for one, calculate the sensitivity of the transient stability margin and transient stability margin of each generator to unit active power output, the mathematical model that Generator Transient Stability constraint is proposed to unit output optimization as constraint condition, obtains the optimal case of each generator output.The method of the present invention has comprehensively considered the various operation constraint conditions during optimization unit output, consider power output limit value, the limit value of branch power and the limit value of node voltage amplitude of unit, therefore to the adaptable of the various methods of operation of system, and it is easy to engineering practical operation, it can guarantee that the transient stability characteristic of generator meets anticipated conditions by optimization unit output.
Description
Technical field
The invention belongs to technical field of power systems, and in particular to a kind of for improving receiving end Power Network Transient Stability nargin
Unit output optimization method.
Background technique
The security reliability of electric system is the primary goal of operation of power networks, to the analysis of power system transient stability with
Control is always one of research emphasis in the industry.China is vast in territory, and the energy is unevenly distributed, the part under the overall background of Power System Interconnection
Disturbance may cause large-scale transient stability problem.Meanwhile HVDC Transmission Technology is because its dramatic benefit is in Power System Interconnection, long distance
It is widely applied from large capacity transmission etc..Alternating current-direct current mixed connection receiving end Power grid structure is increasingly complicated, and system occurs
The probability of failure is soaring；The commutation failure even direct current that the Transient Instability of receiving end power grid frequently can lead to DC transmission system closes
Lock, causes chain reaction.Therefore the transient stability analysis of receiving end power grid is one of critical issue urgently to be resolved.
The appraisal procedure of transient stability analysis of power system mainly has energy function method and generator rotor angle difference criterion method, wherein energy
Whether function method is stable to judge system by construction energy function, but is not suitable for the large scale electric network of current complexity；Generator rotor angle
Poor criterion method by judge the opposite generator rotor angle between generator amature whether more than one threshold value come assessment system whether transient state lose
Surely.
Xue Feng, Ding Chun, Xue Yusheng et al. are in document " algorithm and its engineer application that numerical integration terminates automatically, power train
System automation, 2001 (20): proposed in 9~13. " based on extended equal area criterion (EEAC) stability Quantity Analysis Theory with
And the automatic termination algorithm of numerical integration, integral each time step by multimachine trajectory map to two-dimentional phase plane, and then utilize
The stability of equal area criterion analysis system, but this method the shortcomings that there are model applicabilities is not strong, low precision.
Su Fu, Yang Songhao, Wang Huaiyuan, Zhang Baohui et al. are in document " electric power system transient stability time-domain-simulation Quick stop
Algorithm research, Proceedings of the CSEE, 2017,37 (15): 4372~4378+4577 " in propose transient stability margin it is general
It reads, judges the transient stability of system rapidly by the transient stability margin size of phase path concavity and convexity, avoid for a long time
Time-domain-simulation, improve the assessment efficiency of electric power system transient stability；Transient stability margin tool based on phase path concavity and convexity
There is the characteristics of accurate, convenience of calculation, and can be applied to the power grid with strong nonlinearity characteristic.
In addition, China's power system staff proposes " three lines of defence " for resisting failure, it is that may be present to power grid
Transient safe and stable problem cuts the control measure such as machine, cutting load or off-the-line by control strategy implementation, keeps the transient state of system steady
It is fixed；But " three lines of defence " is carved need to lose certain loads and generator if necessary there are operation difficulty.
Flexible DC transmission technology is also widely used in the transient stability characteristic in power grid with regulating system at present, but soft
Straight technological investment is huge, while may bring new stable problem.The transient stability of system and the variation of its trend mode are close
Cut phase is closed, and under different generator output combinations, its transient stability can not yet when system is by a certain identical disturbance
Together；Therefore a kind of unit output optimization side that receiving end grid stability can be only improved by change generator output combination is needed
Method, while needing to find a Transient Stability Criterion accurate, easy to power system transient stability.
Summary of the invention
In view of above-mentioned, the present invention provides a kind of for improving the unit output optimization side of receiving end Power Network Transient Stability nargin
Method can guarantee that power system transient stability characteristic is met the requirements by optimization unit output.
It is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin, include the following steps:
(1) the receiving end power grid of initial operating state known for one, according to synchronous generator in receiving end power grid respective
The initial transient that power-angle curve and angular speed offset curve under most serious short trouble calculate each synchronous generator is stablized
Nargin；
(2) by changing the active power output of synchronous generator, calculate each synchronous generator transient stability margin has relatively
The sensitivity of function power output；
(3) it establishes following unit output mathematical model and it is optimized, it is active out to obtain each synchronous generator
The optimal case of power；
Wherein: SGFor the synchronous generator set in receiving end power grid, i is synchronous generator set SGIn any synchronization
Generator, PGiFor the active power output of synchronous generator i, c2i、c1i、c0iRespectively synchronous generator i consumption is characteristic secondary
Term coefficient, Monomial coefficient, constant term.
Further, the concrete methods of realizing of the step (1) are as follows: for any synchronous generator in receiving end power grid
By bus nodes each in traverse scanning power grid the power-angle curve of the synchronous generator and angle speed under three phase short circuit fault occur for machine
Offset curve is spent, the i.e. maximum one group of power-angle curve of the generator rotor angle amplitude of oscillation and angular speed offset are taken under wherein most serious fault condition
Curve is bent about the phase path of the synchronous generator to draw in power angle of synchronous generator-angular speed offset phase plane
Line, and calculate by following formula the initial transient stability margin η of the synchronous generator0；
η0=min | ln|
Wherein: Δ ωnFor the angular speed offset of n-th of sampled point corresponding in phase path curve, δnFor in phase path curve
The generator rotor angle of corresponding n-th of sampled point, knFor the slope of n-th of sampled point corresponding in phase path curve, lnIt is phase path curve with function
Angle δnChange rate, n is natural number greater than 0 and 1≤n≤N, N be number of sampling points and N=T/ Δ T, T are to sample duration, Δ
T is the sampling interval.
Further, the concrete methods of realizing of the step (2) are as follows: for any synchronous generator in receiving end power grid
Machine changes its active power output size under its initial operating state, and recalculates the synchronous generator under current active power output
The transient stability margin of machine, and then it is calculated by the following formula spirit of the synchronous generator transient stability margin with respect to active power output
Sensitivity s；
Wherein: η is the transient stability margin of the synchronous generator under current active power output, η0For the synchronous generator
Initial transient stability margin, Δ PGFor given active power output variable quantity.
Further, the unit output mathematical model in the step (3) need to meet following economize on electricity power-balance constraint；
Wherein: SBFor the bus nodes set in receiving end power grid, k and j are bus nodes set SBIn bus nodes and
K ≠ j, VkAnd VjThe respectively voltage magnitude of bus nodes k and j, PHkBy on bus nodes k even synchronous generator it is active go out
Power, the P if without synchronous generator connectionHk=0；QHkIdle power output by connecting synchronous generator on bus nodes k, if without same
Walk generator connection then QHk=0；PDkBurden with power amount by connecting load on bus nodes k, the P if zero load connectionDk=0；
QDkLoad or burden without work amount by connecting load on bus nodes k, the Q if zero load connectionDk=0；θkjBetween bus nodes k and j
Phase angle difference, GkjAnd BkjRespectively xth in receiving end network system node admittance matrixkRow yjThe real and imaginary parts of column element,
xkFor bus nodes k in system node admittance matrix corresponding row serial number, yjIt is bus nodes j in system node admittance matrix
In corresponding column serial number.
Further, the unit output mathematical model in the step (3) need to meet following unit active power output bound
Constraint；
PGimin≤PGi≤PGimax
Wherein: PGiminAnd PGimaxThe respectively lower limit value and upper limit value of synchronous generator i active power output.
Further, the unit output mathematical model in the step (3) need to meet following bus nodes voltage bound
Constraint；
Vkmin≤Vk≤Vkmax(k∈SB)
Wherein: SBFor the bus nodes set in receiving end power grid, k is bus nodes set SBIn any bus nodes, Vk
For the voltage magnitude of bus nodes k, VkminAnd VkmaxThe respectively Lower Limit Amplitude and upper limit magnitude of bus nodes k voltage.
Further, the unit output mathematical model in the step (3) need to meet following Line Flow constraint；
Wherein: SlFor the line set in receiving end power grid, l is line set SlIn any bar route, PlFor route l's
Effective power flow, p and q are respectively the bus nodes that the both ends route l connect, VpAnd VqThe respectively voltage magnitude of bus nodes p and q,
θpqFor the phase angle difference between bus nodes p and q, GpqAnd BpqRespectively xth in receiving end network system node admittance matrixpRow yq
The real and imaginary parts of column element, xpFor bus nodes p in system node admittance matrix corresponding row serial number, yqFor bus nodes q
Corresponding column serial number, P in system node admittance matrixlmaxFor the effective power flow upper limit value of route l.
Further, it is steady need to meet following unit active power output transient state for the unit output mathematical model in the step (3)
Conclude a contract or treaty beam；
ηi0+si(PGi-PGi0)≥ηiset
Wherein: siSensitivity for synchronous generator i transient stability margin with respect to active power output, ηi0For synchronous generator i
Initial transient stability margin, PGi0For active power output of the synchronous generator i under initial operating state, ηisetFor synchronous generator
The transient stability margin desired value of machine i.
The method of the present invention has comprehensively considered the various operation constraint conditions during optimization unit output, it is contemplated that unit
Power output limit value, the limit value of branch power and the limit value of node voltage amplitude, therefore to the adaptability of the various methods of operation of system
By force, and it is easy to engineering practical operation, can guarantees that the transient stability characteristic of generator meets expected item by optimization unit output
Part.
Detailed description of the invention
Fig. 1 is the attachment structure schematic diagram of 10 machine system of IEEE39 node.
Fig. 2 is the initial power-angle curve schematic diagram of generator 35.
Fig. 3 is that the initial concavity and convexity index of generator 35 changes schematic diagram.
Fig. 4 is the power-angle curve schematic diagram after the optimization power output of generator 35.
Fig. 5 is that the concavity and convexity index after the optimization power output of generator 35 changes schematic diagram.
Specific embodiment
In order to more specifically describe the present invention, with reference to the accompanying drawing and specific embodiment is to technical solution of the present invention
It is described in detail.
This example verifies the feasibility of unit output optimization method of the present invention using 10 machine system of IEEE39 node.Fig. 1
In give the wiring diagram of 39 node, 10 machine system, generator 39 is ideal source in Fig. 1, to equivalent infinite bulk power grid,
Therefore generator 39 is not considered when considering transient stability margin constraint.For above-mentioned network system, using the present invention for mentioning
The unit output optimization method of high receiving end Power Network Transient Stability nargin, detailed process is as follows:
(1) power-angle curve according to synchronous generator each in power grid under respective most serious short trouble and angular speed deviate
Amount curve calculates the transient stability margin of each generator.By occurring under three phase short circuit fault everywhere in scanning power grid, respectively
Power-angle curve and angular speed offset curve of each synchronous generator respectively under most serious fault condition are obtained, in synchronous generator
Generator rotor angle-angular speed (δi—Δωi) draw using time t as the curve (phase path) of parameter, with following formula meter in phase plane
Calculate the initial transient stability margin η of each generatori:
ηi=| li|min
Wherein: Δ ωiFor the angular speed offset of i-th synchronous generator in power grid, δiFor i-th in power grid synchronous hair
The generator rotor angle of motor, kiFor i-th power angle of synchronous generator-angular speed (δ in power gridi—Δωi) it with time t is ginseng in phase plane
The slope of the curve (phase path) of variable, liFor the phase path concavity and convexity index of i-th synchronous generator in power grid, it is defined as phase
Track slope is with generator rotor angle δiChange rate.
By taking generator 35 as an example, (i.e. No. 22 bus three-phase shortcircuits) obtain the generator rotor angle song of generator 35 under most serious failure
Line is as shown in Fig. 2, concavity and convexity index changes as shown in figure 3, being η by the initial transient stability margin that Fig. 3 can calculate generator 3535
=0.1281.The initial active power output of all generators and transient stability margin are as shown in table 1 in power grid.
Table 1
(2) to calculate each unit transient state according to the different transient stability margins under each synchronous generator difference active power output steady
Nargin is determined to the sensitivity s of unit active power outputi.Specifically, an active power output changes delta is given to each unit respectively
PGi, calculate its transient stability margin η under the power output under respective most serious short troublei1, calculated respectively using following formula
Sensitivity s of the unit transient stability margin to unit active power outputi:
Wherein: ηi1It is i-th synchronous generator in power grid in power output changes delta PGiUnder it is temporary under most serious short trouble
State stability margin, ηiFor the initial transient stability margin of i-th synchronous generator in power grid.All hairs in IEEE39 node power grid
Motor transient stability margin is as shown in table 1 to the sensitivity of unit active power output.
(3) mathematical model for proposing unit output optimization, obtains the optimal case of each generator output；Unit output optimization
Model are as follows:
Wherein: PGiFor the active power output after i-th generator optimization in power grid, SGFor generator collections all in power grid,
c2i、c1i、c0iFor the consumption characteristic curve parameter of i-th generator in power grid.
Equality constraint:
This formula is economize on electricity power balance equation, wherein SBFor node sets all in system, QGiFor institute's running fire in node i
The idle power output of motor, PDi、QDiFor the active and reactive load in node i, Vi、θiFor the voltage magnitude and phase angle of node i, θij
Phase angle difference between node i, j, Gij、BijFor the real part and imaginary part of node admittance matrix the i-th row jth column element.
Inequality constraints condition:
PGimin≤PGi≤PGimax(i∈SG)
Vimin≤Vi≤Vimax(i∈SB)
ηi0+si(PGi-PGi0)≥ηiset(i∈SG)
The above inequality constraints is followed successively by generated power power output or more bound constrained, node voltage bound constrained, route up and down
The Transient Stability Constraints of trend constraint, unit active power output.Wherein PGimin、PGimaxIt is minimum for the active power output of i-th generator
Limit value and threshold limit value, Vimin、VimaxFor the minimum limitation of node i voltage magnitude and threshold limit value, SlFor all routes of power grid
Set, PlFor the effective power flow of route l, two end node of route is i, j, PlmaxFor Line Flow limit value, ηisetRespectively to be sent out in power grid
The transient stability margin desired value of motor, PGi0For the active power output under i-th synchronous generator initial operating state in power grid.
Take ηiset=1 is the Power Network Transient Stability nargin desired value, obtains all power generations in IEEE39 node system after solving the model
The optimization active power output of machine is as shown in table 2.
Table 2
(4) it calculates each synchronous generator to contribute the transient stability margin under optimal case, verifies whether to meet transient state steady
Determine nargin desired value.The optimization power generating value that each unit is obtained after solving optimization model, under respective most serious short circuit fault condition
Transient stability margin value is calculated again, verifies whether to meet Transient Stability Constraints condition:
ηi0+si(PGi-PGi0)≥ηiset(i∈SG)
In IEEE39 node system, solve all units optimization power output after in method described in step (1) again
It is as shown in table 2 to calculate transient stability margin of all units under respective most serious short trouble.
With generator 35 come verify optimization after transient stability margin whether meet the requirements, (i.e. 22 under most serious failure
Number bus three-phase shortcircuit) power-angle curve after the optimization power output of generator 35 is obtained as shown in figure 4, phase path concavity and convexity index is as schemed
Shown in 5.The transient stability margin η of the optimization of generator 35 power output as shown in Figure 535=1.128, compared to initial transient stability margin
0.1281 is greatly improved；As shown in Figure 4 under identical disturbance, the generator rotor angle swing curve of generator 35 is compared to initial fortune
Row state (Fig. 2) fluctuation it is smaller, therefore the transient characterisitics of generator 35 really arrived improvement.
It is not as shown in Table 2 that the transient stability margin of each unit all increases, because the transient state of certain units is steady
Determine that nargin is sufficiently large, but it can be seen that the transient stability margin of each unit all meets desired value ηiset=1, thus
Demonstrate effectiveness of the invention.In practice scene can each unit is specified according to the actual situation it is different
Transient stability margin desired value flexibly uses unit output optimization method of the present invention.
The above-mentioned description to embodiment is for that can understand and apply the invention convenient for those skilled in the art.
Person skilled in the art obviously easily can make various modifications to above-described embodiment, and described herein general
Principle is applied in other embodiments without having to go through creative labor.Therefore, the present invention is not limited to the above embodiments, ability
Field technique personnel announcement according to the present invention, the improvement made for the present invention and modification all should be in protection scope of the present invention
Within.
Claims (8)
1. it is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin, include the following steps:
(1) the receiving end power grid of initial operating state known for one, according to synchronous generator in receiving end power grid respectively most tight
Power-angle curve and angular speed offset curve under weight short trouble calculate the initial transient stability margin of each synchronous generator；
(2) by change synchronous generator active power output, calculate each synchronous generator transient stability margin relatively it is active go out
The sensitivity of power；
(3) it establishes following unit output mathematical model and it is optimized, obtain each synchronous generator active power output
Optimal case；
Wherein: SGFor the synchronous generator set in receiving end power grid, i is synchronous generator set SGIn any synchronous generator
Machine, PGiFor the active power output of synchronous generator i, c2i、c1i、c0iThe respectively characteristic secondary term system of synchronous generator i consumption
Number, Monomial coefficient, constant term.
2. unit output optimization method according to claim 1, it is characterised in that: the concrete methods of realizing of the step (1)
Are as follows: for any synchronous generator in receiving end power grid, three-phase shortcircuit occurs by bus nodes each in traverse scanning power grid
The power-angle curve of the synchronous generator and angular speed offset curve under failure take under wherein most serious fault condition i.e. generator rotor angle pendulum
Maximum one group of power-angle curve is with angular speed offset curve to equal in power angle of synchronous generator-angular speed offset
The phase path curve about the synchronous generator is drawn on face, and the initial temporary of the synchronous generator is calculated by following formula
State stability margin η0；
Wherein: Δ ωnFor the angular speed offset of n-th of sampled point corresponding in phase path curve, δnIt is corresponding in phase path curve
The generator rotor angle of n-th of sampled point, knFor the slope of n-th of sampled point corresponding in phase path curve, lnIt is phase path curve with generator rotor angle δn
Change rate, n is natural number greater than 0 and 1≤n≤N, N be number of sampling points and N=T/ Δ T, T are to sample duration, and Δ T is
Sampling interval.
3. unit output optimization method according to claim 1, it is characterised in that: the concrete methods of realizing of the step (2)
Are as follows: for any synchronous generator in receiving end power grid, changes its active power output size under its initial operating state, lay equal stress on
The new transient stability margin for calculating the synchronous generator under current active power output, and then it is calculated by the following formula synchronization hair
Sensitivity s of the motor transient stability margin with respect to active power output；
Wherein: η is the transient stability margin of the synchronous generator under current active power output, η0For the initial of the synchronous generator
Transient stability margin, Δ PGFor given active power output variable quantity.
4. unit output optimization method according to claim 1, it is characterised in that: the unit output number in the step (3)
Following economize on electricity power-balance constraint need to be met by learning model；
Wherein: SBFor the bus nodes set in receiving end power grid, k and j are bus nodes set SBIn bus nodes and k ≠
J, VkAnd VjThe respectively voltage magnitude of bus nodes k and j, PHkActive power output by connecting synchronous generator on bus nodes k,
The P if the connection of no synchronous generatorHk=0；QHkIdle power output by connecting synchronous generator on bus nodes k, if being sent out without synchronous
Motor connects then QHk=0；PDkBurden with power amount by connecting load on bus nodes k, the P if zero load connectionDk=0；QDkFor
The load or burden without work amount of connected load on bus nodes k, the Q if zero load connectionDk=0；θkjFor the phase between bus nodes k and j
Angular difference, GkjAnd BkjRespectively xth in receiving end network system node admittance matrixkRow yjThe real and imaginary parts of column element, xkFor
Bus nodes k corresponding row serial number, y in system node admittance matrixjIt is right in system node admittance matrix for bus nodes j
The column serial number answered.
5. unit output optimization method according to claim 1, it is characterised in that: the unit output number in the step (3)
Following unit active power output bound constraint need to be met by learning model；
PGimin≤PGi≤PGimax
Wherein: PGiminAnd PGimaxThe respectively lower limit value and upper limit value of synchronous generator i active power output.
6. unit output optimization method according to claim 1, it is characterised in that: the unit output number in the step (3)
Following bus nodes voltage bound constraint need to be met by learning model；
Vkmin≤Vk≤Vkmax(k∈SB)
Wherein: SBFor the bus nodes set in receiving end power grid, k is bus nodes set SBIn any bus nodes, VkFor mother
The voltage magnitude of line node k, VkminAnd VkmaxThe respectively Lower Limit Amplitude and upper limit magnitude of bus nodes k voltage.
7. unit output optimization method according to claim 1, it is characterised in that: the unit output number in the step (3)
Following Line Flow constraint need to be met by learning model；
Wherein: SlFor the line set in receiving end power grid, l is line set SlIn any bar route, PlFor the active of route l
Trend, p and q are respectively the bus nodes that the both ends route l connect, VpAnd VqThe respectively voltage magnitude of bus nodes p and q, θpq
For the phase angle difference between bus nodes p and q, GpqAnd BpqRespectively xth in receiving end network system node admittance matrixpRow yqColumn
The real and imaginary parts of element, xpFor bus nodes p in system node admittance matrix corresponding row serial number, yqExist for bus nodes q
Corresponding column serial number, P in system node admittance matrixlmaxFor the effective power flow upper limit value of route l.
8. unit output optimization method according to claim 1, it is characterised in that: the unit output number in the step (3)
Following unit active power output Transient Stability Constraints need to be met by learning model；
ηi0+si(PGi-PGi0)≥ηiset
Wherein: siSensitivity for synchronous generator i transient stability margin with respect to active power output, ηi0For the first of synchronous generator i
Beginning transient stability margin, PGi0For active power output of the synchronous generator i under initial operating state, ηisetFor synchronous generator i's
Transient stability margin desired value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910542412.0A CN110350596A (en) | 2019-06-21 | 2019-06-21 | It is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910542412.0A CN110350596A (en) | 2019-06-21 | 2019-06-21 | It is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110350596A true CN110350596A (en) | 2019-10-18 |
Family
ID=68182718
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910542412.0A Pending CN110350596A (en) | 2019-06-21 | 2019-06-21 | It is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110350596A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110571809A (en) * | 2019-10-28 | 2019-12-13 | 华北电力大学 | Commutation sequence control method and system based on energy function |
| CN111786386A (en) * | 2020-07-15 | 2020-10-16 | 华北电力大学 | Control system and method for preventing direct current blocking based on transient energy method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107689638A (en) * | 2017-09-28 | 2018-02-13 | 中国电力科学研究院 | A kind of control method for coordinating of electrical power system transient containing wind-powered electricity generation based on phase path analysis |
| CN108183484A (en) * | 2017-12-14 | 2018-06-19 | 河海大学 | A kind of Semi linearization optimal load flow model of decoupling based on thermal starting environment |
-
2019
- 2019-06-21 CN CN201910542412.0A patent/CN110350596A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107689638A (en) * | 2017-09-28 | 2018-02-13 | 中国电力科学研究院 | A kind of control method for coordinating of electrical power system transient containing wind-powered electricity generation based on phase path analysis |
| CN108183484A (en) * | 2017-12-14 | 2018-06-19 | 河海大学 | A kind of Semi linearization optimal load flow model of decoupling based on thermal starting environment |
Non-Patent Citations (1)
| Title |
|---|
| 杨松浩 等: "计及暂态稳定约束的断面功率调整方法研究", 《中国电机工程学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110571809A (en) * | 2019-10-28 | 2019-12-13 | 华北电力大学 | Commutation sequence control method and system based on energy function |
| CN110571809B (en) * | 2019-10-28 | 2020-12-01 | 华北电力大学 | Commutation sequence control method and system based on energy function |
| CN111786386A (en) * | 2020-07-15 | 2020-10-16 | 华北电力大学 | Control system and method for preventing direct current blocking based on transient energy method |
| CN111786386B (en) * | 2020-07-15 | 2021-09-07 | 华北电力大学 | Control system and method for preventing direct current blocking based on transient energy method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102879668B (en) | Asymmetric fault analysis method for power distribution network including inverted distribution type power supply | |
| CN110797874B (en) | State estimation method for AC/DC hybrid power distribution network containing power electronic transformer | |
| CN103760400B (en) | Voltage stabilization on-line monitoring method based on wide area measurement system | |
| Ameur et al. | Steady state and dynamic analysis of renewable energy integration into the grid using PSS/E software | |
| CN105244901B (en) | A kind of nonlinear decentralized control method of HVDC transmission system | |
| CN109038641A (en) | Calculation method and system for rapidly improving power capacity of high-voltage direct-current transmission system | |
| CN110350596A (en) | It is a kind of for improving the unit output optimization method of receiving end Power Network Transient Stability nargin | |
| CN108233408B (en) | A kind of MMC-MTDC system self-adaption droop control method | |
| CN109546677A (en) | A kind of scale offshore wind farm flexibility transmitting system safety control strategy method for solving | |
| CN112865167A (en) | Method and system for determining transient stability margin of alternating current-direct current hybrid system | |
| CN104617576A (en) | Multi-DC-feed (multi-direct current-feed) AC (alternate current) grid fault calculating method taking DC control characteristics into consideration | |
| CN108667027A (en) | Ac and dc systems power flow transfer search containing flexible direct current and quantitative analysis method | |
| CN103413043B (en) | A kind of method for solving of electric system hyperspace section thermal stability limit boundary | |
| CN1804644B (en) | Method for capturing out-of-step interface based on reactive power | |
| Yu et al. | Dynamic performance of embedded HVDC with frequency control strategy | |
| CN106451421A (en) | Modeling and ring-opening method for electromagnetic ring network-containing alternating current-direct current power grid | |
| CN110880780B (en) | Optimal emergency power control method of multi-terminal direct current system | |
| Maya et al. | A generalised three phase power flow algorithm incorporating the uncertainty of Photo Voltaic (PV) source for unbalanced distribution network | |
| CN111952995B (en) | Direct-current capacitor voltage balance control method under single-phase fault of power grid | |
| CN113270871B (en) | Flexible interconnection device capacity configuration optimization method based on power distribution network N-1 safety assessment | |
| Spallarossa et al. | Assessment of disturbance propagation between AC grids through VSC HVDC links using reduced Great Britain model | |
| CN107959276A (en) | A kind of adaptive reclosing method of single loop line | |
| Djehaf et al. | Modeling of a multi-level converter based VSC HVDC supplying a dead load | |
| Li et al. | Interval Electrical Betweenness Method for Power Grid Vulnerability Assessment Considering Wind Power | |
| CN105119281A (en) | UPFC control method with rapid terminal sliding mode variable structure |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191018 |