CN105552899B - A kind of method of power system restoration power after calculating is had a power failure on a large scale - Google Patents
A kind of method of power system restoration power after calculating is had a power failure on a large scale Download PDFInfo
- Publication number
- CN105552899B CN105552899B CN201610038041.9A CN201610038041A CN105552899B CN 105552899 B CN105552899 B CN 105552899B CN 201610038041 A CN201610038041 A CN 201610038041A CN 105552899 B CN105552899 B CN 105552899B
- Authority
- CN
- China
- Prior art keywords
- unit
- power
- node
- beeline
- black starting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000005540 biological transmission Effects 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 18
- 230000005611 electricity Effects 0.000 claims description 14
- 238000011084 recovery Methods 0.000 claims description 12
- 230000001131 transforming effect Effects 0.000 claims description 8
- 238000003745 diagnosis Methods 0.000 claims description 7
- 238000011161 development Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 230000007257 malfunction Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 230000008450 motivation Effects 0.000 claims 1
- 238000011002 quantification Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 abstract 1
- 230000008569 process Effects 0.000 description 3
- 238000000205 computational method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 238000012546 transfer Methods 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
A kind of method of power system restoration power after having a power failure on a large scale the invention discloses calculating, including：1) particular state of system after identification has a power failure, obtains data used in calculating；2) available black starting-up unit in system is filtered out, calculates the sum of capacity of all black starting-up units and corresponding perunit value；3) each unit to be restored in system is calculated to the average beeline and load bus of black starting-up unit to the average beeline of each machine group node, and the average distance both calculated；4) using the capacity of black starting-up unit and perunit value and step 3) in average distance business as power grid after having a power failure on a large scale restoring force.Advantageous effect of the present invention：Provide a kind of means of quantification to assess and analyzing restoring force of the power grid after having a power failure on a large scale, and can weak spot of the power grid easy to find in terms of restoring force is improved, proposition is correspondingly improved measure.
Description
Technical field
The meter of power system restoration power after having a power failure on a large scale the present invention relates to a kind of power grid index calculating method more particularly to a kind of calculating
Calculation method.
Background technology
With the fast development of economic society, the scale of electric system is increasing, and with extra-high voltage alternating current-direct current electricity
The large-scale renewable energy source base access power grid of the construction of net and large-scale wind power field and photovoltaic plant etc., the dynamic of electric system
Behavior becomes increasingly complicated, and to consider environmental problem and economic factor during operation of power networks, gets over the operating point of power grid
Come closer to its safety margins point, which greatly increases operations of power networks and the complexity safeguarded.Therefore at system local fault
It manages not at that time, it is most likely that trigger system occurrence of large-area power outage.Such as in August in 2003 14 days, North America history has occurred
The power outage of upper most serious, this power outage have fed through to many areas of America ＆ Canada, lose load 61.8GW,
It is 50,000,000 to influence population, and occurrence cause is dealt with improperly for a 345kV transmission line short-circuit fault, causes trend
A wide range of transfer, causes a plurality of transmission line of electricity cascading trip due to overload, ultimately results in generation large-scale blackout；In November, 2006
4, large-scale blackout occurred for European interconnected network, and the Kuomintang-Communists such as moral, method, meaning lose about 16,000,000 kW of load, 1500 general-purpose families by
To influence；Large-scale blackout occurs on November 10th, 2009, Brazil, Paraguay power grid, causes three times 750kV line trippings, and two
Item ± 600kV DC lines are latched, and Brazilian power network loses about 17,000,000 kW power supplys, and power failure range involves Brazilian 12 states and neighbouring country
The most area of Paraguay, 5000 to 60,000,000 people are affected；On July 30th, 2012 and 31 days, north India and east
Area recurs large area blackout twice, covers the territory of more than half, directly affects the life of people more than 600,000,000.State
The operating experience of inside and outside electric system shows that although the extensive application of new technology and new equipment in electric system can improve
The stability and reliability of system operation, but it is still unavoidable from the generation having a power failure on a large scale.
The interdependent degree that modern society supplies firm power is higher and higher, once the power outage of occurrence of large-area, it will
It can be to the influence and huge economic loss that social production life is brought.To the data analysis table of domestic and international repeatedly large-scale blackout
It is bright, the loss that power outage is brought and power off time exponent function relation.The system blackout time is longer, caused by adversely affect it is tighter
Weight.After generation power outage, under the guidance of corresponding recovery scheme, carry out corresponding recovery operation as early as possible, be conducive to accelerate system
The recovering process of system shortens the power off time of system, so as to reduce the influence and loss that system blackout accident is brought.System is recovered
The duration of process not only with whether recovery scheme complete and the factors such as dispatcher's qualification are related, also with power grid by
The fast restoration capabilities of itself after external disturbance, the i.e. restoring force of power grid are related.It is big that power system restoration power mainly considers that power grid occurs
In the case of the extreme accident such as area power failure, system utilizes the ability of the fast quick-recovery of various resources, this is conventional electric power systems organization
Never it is related to.Restoring force is the key character of intelligent grid, and the research for carrying out power system restoration power is intelligent grid development
Inexorable trend.The many factors such as the grid structure and power supply of power system restoration power and system, sharing of load are related, extensive to improve power grid
Multiple power can instruct planning for the expansion of power grids for target, and also existing power grid can be analyzed, and proposition, which is correspondingly improved, to be arranged
It applies.
At this stage, the research in relation to power system restoration power is at the early-stage, and Assessment theory and technique study are still up for into one
The systemic achievement of step, Study on Relative Factors lack the instrument of quantitative analysis also in the qualitative analysis stage.
The content of the invention
The purpose of the present invention is exactly that a kind of side of power system restoration power after calculating is had a power failure on a large scale is provided to solve above-mentioned problem
Method.Consider after large-scale blackout occurs, the capacity and distribution of available black starting-up unit and unit to be restored in power grid,
The factors such as the distribution of grid structure and load website propose a kind of quantitative calculation method of power system restoration power, to instruct power grid
Planning and power system restoration power, which are promoted, provides foundation.
To achieve the above object, concrete scheme of the invention is as follows：
The method of power system restoration power, comprises the following steps after a kind of calculating is had a power failure on a large scale：
(1) carry out availability diagnosis to all hair power transmission and transforming equipment in system, black starting-up unit and treated in acquisition system
Recover in the quantity of unit and its corresponding unit capacity, system the connection relation of transmission line of electricity and its corresponding circuit can be used to join
Number；
(2) according to the data of acquisition, the corresponding weighted network connection matrix of system is established；
(3) available black starting-up unit in system is filtered out, calculates the sum of the capacity of all available black starting-up units
And corresponding perunit value；
(4) each unit to be restored in system is calculated to the average beeline and load bus of black starting-up unit to each
The average beeline of machine group node, and calculate above-mentioned two and be averaged the average distance of beeline；
(5) according to the perunit value of the sum of the capacity of all available black starting-up units asked in step (3) and step (4)
The average distance of middle calculating asks for the restoring force of power grid after having a power failure on a large scale, and reference frame is provided for the development plan of power grid.
In the step (1), the method that availability diagnosis is carried out to hair power transmission and transforming equipment all in system is：
Operating status residing for hair power transmission and transforming equipment all in system is identified and judged, is filtered out in failure
Not available equipment in the equipment of state or inspecting state, i.e. recovery process.
In the step (2), when establishing the weighted network connection matrix M corresponding to system, M is a square formation, and line number
The quantity of available power plant and power transformation tiny node in system is equal to columns；
When there are available transmission line of electricity l between system interior joint i and node jijWhen, if circuit lijLine reactance
Perunit value is xij, corresponding element is M in connection matrix Mij, and have Mij=Mji=xij；If conversely, system interior joint i
When available transmission line of electricity being not present between node j, Mij=Mji=∞.
In the step (3), in system each unit to be restored to the average beeline of black starting-up unit according to each treating
The minimum value for recovering unit to the beeline of each black starting-up unit determines.
Further, the beeline Dis of unit i to black starting-up unit j to be restored is being calculatedijWhen, by black starting-up unit j
Place node is as source point s, and node is as destination node t where unit i to be restored, by the section in weighted network connection matrix M
Point is divided into two classes, is denoted as set S and set U respectively：Set S is the destination node for the shortest path for having sought out source point s
Set, set U is the destination node set having not yet been calculated out to the shortest path of source point s.
Further, the beeline Dis of unit i to black starting-up unit j to be restored is calculated using critical path method (CPM)ij, tool
Body method is as follows：
1) during initial calculation, source point s has been contained only in set S, all nodes in addition to source point s are included in set U, for
Any point x in set U, D [x] are the distance of source point s to point x, and as circuit is right in connection matrix M between source point s and node x
Answer the numerical value M of elementsx；
2) the point x with minimum D [x] is selected from set U, then moves on to point x in set S from set U；
3) for any node y in set U residue nodes, if there is D [y]>D[x]+Mxy, then by source point s to point y
Distance D [y] be arranged to D [x]+Mxy, wherein, MxyCircuit corresponding element in connection matrix M between node x and node y
Numerical value；
4) the point y with minimum D [y] is selected from remaining node, then moves on to point y in set S from set U；
5) step 3) and step 4) are repeated, until all nodes are all shifted in set S in set U；
6) beeline of source point s to destination node t is D [t].
In the step (4), the average beeline of load bus to each machine group node is according to unit i to be restored to all
The minimum value of the beeline of machine group node determines.
Further, calculated load node to each machine group node average beeline when, load bus is divided into two classes：
One kind is to be included in each unit to be restored into the corresponding restoration path of minimum value of the beeline of each black starting-up unit
Load bus；Load bus that is another kind of, being not included in above-mentioned path；
First kind load bus is already contained in restoration path, therefore is not considered, and only considers the second type load node, and
And the node that will be included in the corresponding restoration path of minimum value of each unit to be restored to the beeline of each black starting-up unit
It is clustered.
Further, unit i to be restored is calculated to the beeline of all machine group nodes using critical path method (CPM).
In the step (5), the restoring force of power grid is after having a power failure on a large scale：The sum of the capacity of all available black starting-up units
Perunit value and the quotient of the average distance calculated in step (4).
Beneficial effects of the present invention：
Firstth, the method for power system restoration power after calculating proposed by the present invention is had a power failure on a large scale, it is convenient to calculate different power grids and exist
The size of the restoring force after power outage occurs, and is compared to each other.
Secondth, the method for power system restoration power after calculating proposed by the present invention is had a power failure on a large scale, is occurring greatly to assess and analyzing power grid
Restoring force after power failure provides a kind of means of quantification, and can weakness of the power grid easy to find in terms of restoring force is improved
Point, proposition are correspondingly improved measure.
3rd, the quantitative calculation method of power system restoration power after calculating proposed by the present invention is had a power failure on a large scale so that with power system restoration
Power, which is up to objectives guide Electric Power Network Planning, becomes more directly perceived and effective.
4th, the present invention has considered black starting-up unit in power grid and has held in the restoring force after calculating power grid and having a power failure on a large scale
The factors such as amount and position so that the result of calculation of power system restoration power is more objective, and layouting for black starting-up unit is excellent available for instructing
Change and capacity optimizes.
5th, the position relationship of unit to be restored and black starting-up unit is considered in the present invention, can effectively reflect black open
Unit starting process to be restored during dynamic so that power system restoration power can effectively embody actual recovery of the power grid after breaking down
Ability.
6th, load bus is classified, is more in line in the restoring force after calculating power grid and having a power failure on a large scale by the present invention
Actual recovery process after power outage occurs for power grid.
Description of the drawings
Fig. 1 is the method flow diagram of power system restoration power after present invention calculating is had a power failure on a large scale；
Fig. 2 is IEEE30 node system structure charts of the embodiment of the present invention.
Specific embodiment
The present invention is described in detail below in conjunction with the accompanying drawings：
A kind of method of power system restoration power after calculating is had a power failure on a large scale, as shown in Figure 1, comprising the following steps：
(1) carry out availability diagnosis to all hair power transmission and transforming equipment in system, black starting-up unit and treated in acquisition system
Recover in the quantity of unit and its corresponding unit capacity, system the connection relation of transmission line of electricity and its corresponding circuit can be used to join
Number；
(2) according to the data of acquisition, the corresponding weighted network connection matrix of system is established；
When establishing the weighted network connection matrix M corresponding to system, M is a square formation, and line number is equal to columns and is
The quantity of available website in system.When there are available interconnection l between system interior joint i and node jijWhen, and circuit lij's
The perunit value of line reactance is xij, then the element M in connection matrix Mij=Mji=xij；If conversely, system interior joint i and node j
Between be not present available interconnection when, Mij=Mji=∞.For the system containing n website, corresponding connection matrix M
For following form.
In formula, for arbitrary element MijAnd Mji, there are Mij=Mji。
(3) available black starting-up unit in system is filtered out, calculates the sum of the capacity of all available black starting-up units
And corresponding perunit value；
Recovery process after system occurrence of large-area power outage includes black starting-up stage, rack reconstruct and load restoration rank
Section.In the black starting-up stage, available black starting-up unit can be started first in system.After black starting-up unit successfully starts up,
The unit for not possessing self-startup ability for remaining provides required startup power, makes its restarting grid-connected.Rank is reconstructed in rack
Section, the progressively important transmission line of electricity in recovery system and important load website, establish a strong rack, be follow-up
Load restoration lays the first stone, and in the load restoration stage, on the basis of early period recovers, carries out large-scale load restoration, as early as possible
Recover the power supply again to dead electricity load, shorten the power off time of system.Successfully starting up for black starting-up unit is entirely to recover
The beginning of journey, the capacity of black starting-up unit represent the rapid ability recovered after system occurrence of large-area power outage,
It can be as an important indicator of reflection power system restoration power.It filters out in system after all available black starting-up units, owns
The sum of the capacity of black starting-up unit PblackCalculation formula is as follows
In formula, n is the number of available black starting-up unit in system；PiFor the capacity of i-th black starting-up unit.Pblack's
The calculation formula of perunit value is as follows
In formula, SBFor the power base value chosen in step S1.
(4) each unit to be restored in system is calculated to the average beeline and load bus of black starting-up unit to each
The average beeline of machine group node, and calculate above-mentioned two and be averaged the average distance of beeline；
Each unit to be restored is to the average beeline of black starting-up unitCalculation formula is as follows：
LSi=min (Disij), j=1,2 ..., n
In formula, m be unit number to be restored, n be black starting-up unit number, DisijFor unit i to be restored to black starting-up machine
The beeline of group j, LSiFor unit i to be restored to the minimum value of the beeline of each black starting-up unit.
Calculating the beeline Dis of unit i to black starting-up unit j to be restoredijWhen, by node where black starting-up unit j
As source point s, node is calculated as destination node t using shortest path first where unit i to be restored.For weighting
Node in M is divided into two classes by network connection matrix M, shortest path first, is denoted as set S and set U respectively.Set S is
The set of the point of shortest path through seeking out source point s, and the collection that set U is the shortest path having not yet been calculated out to source point s
It closes.The calculation procedure of shortest path first is：
1st, during initial calculation, source point s has been contained only in set S, all nodes in addition to source point s are included in set U, for
Any point x in set U, D [x] are the distance of source point s to point x, and as circuit is right in connection matrix M between source point s and node x
Answer the numerical value M of elementsx；
2nd, the point x with minimum D [x] is selected from set U, then moves on to point x in set S from set U；
3rd, for any node y in set U residue nodes, if there is D [y]>D[x]+Mxy, then by source point s to point y
Distance D [y] be arranged to D [x]+Mxy；
4th, step 2 and step 3 are repeated, until all nodes are all shifted in set S in set U；
5th, the beeline of source point s to destination node t are D [t].
In beeline of the calculated load node to all machine group nodes, load bus is divided into two classes：One kind is
Load section included in the corresponding restoration path of minimum value of each unit to be restored to the beeline of each black starting-up unit
Point；Load bus that is another kind of, being not included in above-mentioned path.For in first kind load bus included restoration path,
Therefore do not consider, only consider the second type load node herein, and by each unit to be restored to the most short of each black starting-up unit
The node included in the corresponding restoration path of minimum value of distance is clustered.Load bus is most short to being averaged for each machine group node
DistanceCalculation formula is
Lldi=min (Disij), j=1,2 ..., nut
In formula, nld is the number of the second type load node；Nut is the number of all machine group nodes；DisijFor load bus
I is to the beeline of machine group node j；LldiFor unit i to be restored to the minimum value of the beeline of all machine group nodes.It is counting
Calculate DisijWhen, equally using shortest path first.
WithAverage valueFor
(5) according to the perunit value of the sum of the capacity of all available black starting-up units asked in step (3) and step (4)
The average distance of middle calculating asks for the restoring force of power grid after having a power failure on a large scale, and reference frame is provided for the development plan of power grid.
Power system restoration power R after having a power failure on a large scaleesComputational methods are
Wherein,The perunit value of the sum of capacity for all available black starting-up units,ForWithAverage value.
As a kind of embodiment, the present invention has carried out simulation calculation for IEEE30 node systems, illustrates to calculate power failure
The method flow of power system restoration power after accident.IEEE30 node system structure charts are as shown in Figure 2.Power grid after having a power failure on a large scale according to calculating
The method of restoring force calculates the power system restoration power after IEEE30 node systems generation large-scale blackout and is as follows：
S1：Particular state residing for system is identified, data used in calculating is obtained, and establishes the corresponding weighting of system
Network connection matrix.
Pair the state of system is identified and refers to carry out availability diagnosis to all hair power transmission and transforming equipment in system, i.e.,
Equipment state in which is identified and judges, filters out the equipment in malfunction or inspecting state (in recovery process
Not available equipment)；Calculating required data has：The quantity and unit of black starting-up unit and unit to be restored are held in system
It measures, the connection relation of transmission line of electricity and corresponding line parameter circuit value can be used in system.
Calculate circuit lijLine reactance perunit value xijWhen, a reference value SBIt is taken as 100MVA, UBTake each voltage class
Normal voltage.
In IEEE30 node systems, unit on node 1 is black starting-up unit, capacity 50MW.Node 2,13,22,23
It is unit to be restored with the unit on 27, unit capacity is respectively 60MW, 25MW, 45MW, 50MW and 55MW.IEEE30 nodes system
The voltage class of system is 220kV.
In this step, after power outage occurring to IEEE30 node systems, the availability of various kinds of equipment carries out in power grid
Diagnosis, it was therefore concluded that as follows：Various kinds of equipment is available in recovery process in IEEE30 node systems.And it gets black in system
Start the quantity and unit capacity of unit and unit to be restored, the connection relation of transmission line of electricity and corresponding circuit ginseng in system
Number, and calculate the perunit value of each line reactance.According to the system topological of IEEE30 nodes and each line electricity calculated
Anti- perunit value establishes the corresponding connection matrix M of IEEE30 node systems.
S2：Available black starting-up unit in system is selected, calculates the sum of capacity of all black starting-up units and right
The perunit value answered；
In the black starting-up stage, available black starting-up unit can be started first in system.Black starting-up unit successfully starts up
Afterwards, the unit for not possessing self-startup ability for remaining provides required startup power, makes its restarting grid-connected.Black starting-up machine
The capacity of group represents the rapid ability recovered after system occurrence of large-area power outage, can be used as reflection power system restoration
One important indicator of power.
There was only the unit on node 1 for black starting-up unit in this step, in IEEE30 node systems, easily calculate system
The sum of the capacity of middle black starting-up unit PblackFor 50MW, corresponding perunit value is
S3：Each unit to be restored in system is calculated to the average beeline and load bus of black starting-up unit to each
The average beeline of machine group node, and calculate the average distance of the two.
According to the corresponding connection matrix M of IEEE30 node systems, node 2 where calculating unit to be restored first, 13,22,
23 and 27 arrive the beeline of black starting-up unit place node 1.The calculation procedure of shortest path first is：
1. during initial calculation, node 1 is only contained in set S, all nodes in addition to node 1 are included in set U, for
Any point x in set U, D [x] are the distance that node 1 arrives point x；
2. the point x with minimum D [x] is node 2 in set U, node 2 is moved on to from set U in set S；
3. for any node y in set U residue nodes, if there is D [y]>D[x]+Mxy, then node 1 is arrived into point y
Distance D [y] be arranged to D [x]+Mxy；
4. step 2 and step 3 are repeated, until all nodes are all shifted in set S in set U；
5. the beeline of node 2,13,22,23 and 27 to node 1 where black starting-up unit is D [2], D [13], D
[22], D [23] and D [27].
Unit to be restored on node 2,13,22,23 and 27 is to the average beeline of black starting-up unitCalculation formula
It is as follows：
LSi=Dis [i]
The beeline path of the black starting-up unit on unit to be restored to node 1 on node 2 is：1-2, beeline
For 0.0575.
The beeline path of the black starting-up unit on unit to be restored to node 1 on node 13 is：1-3,3-4,4-
12,12-13, beeline 0.6191.
The beeline path of the black starting-up unit on unit to be restored to node 1 on node 22 is：1-2,2-6,6-
9,9-10,10-21,21-22, beeline 0.6503.
The beeline path of the black starting-up unit on unit to be restored to node 1 on node 23 is：1-3,3-4,4-
12,12-15,15-23, beeline 0.8115.
The beeline path of the black starting-up unit on unit to be restored to node 1 on node 27 is：1-2,2-6,6-
28,28-27, beeline 0.6897.
Unit to be restored on node 2,13,22,23 and 27 is to the average beeline of black starting-up unitFor 0.5656.
Load bus can be divided into two classes：The first kind is that the shortest path of node 1 is arrived included in node 2,13,22,23 and 27
In node, be comprising node altogether：3、4、6、9、10、15、21、28.
Second class is remaining load bus, is comprising node altogether：5、7、8、11、14、16、17、18、19、20、24、25、
26、29、30.In the beeline of the second class node of calculating to each unit, by circuit 1-2,1-3,3-4,4-12,12-13,2-
The node clustering that 6,6-9,9-10,10-21,21-22,12-15,15-23,2-6,6-28,28-27 and circuit are connected is one
A point 1, the second class node to the beeline of each machine group node, as arrive cluster after point 1 beeline.Using most short
The beeline that routing algorithm calculates the second class stage to each machine group node is as shown in table 1.
1 second class node of table is to the beeline of each machine group node
| Node number | Beeline | Node number | Beeline |
| 5 | 0.1980 | 19 | 0.2770 |
| 7 | 0.0820 | 20 | 0.2090 |
| 8 | 0.0420 | 24 | 0.1790 |
| 11 | 0.2080 | 25 | 0.2087 |
| 14 | 0.1997 | 26 | 0.5887 |
| 16 | 0.1987 | 29 | 0.4153 |
| 17 | 0.0845 | 30 | 0.6027 |
| 18 | 0.2185 |
Load bus is to the average beeline of each machine group nodeFor 0.2475.ThereforeWithAverage valueFor
S4：The business that will ask for average distance in the perunit value and S3 of the sum of the capacity of black starting-up unit calculates greatly as quantitative
The method of power system restoration power after power failure.
Power system restoration power R after having a power failure on a large scaleesComputational methods are
Power system restoration power of the IEEE30 node systems after having a power failure on a large scale is 1.23.
Above-mentioned, although the foregoing specific embodiments of the present invention is described with reference to the accompanying drawings, not protects model to the present invention
The limitation enclosed, those skilled in the art should understand that, based on the technical solutions of the present invention, those skilled in the art are not
Need to make the creative labor the various modifications or changes that can be made still within protection scope of the present invention.
Claims (3)
1. a kind of method of power system restoration power after calculating is had a power failure on a large scale, it is characterized in that, comprise the following steps：
(1) availability diagnosis is carried out to all hair power transmission and transforming equipment in system, black starting-up unit and to be restored in acquisition system
The connection relation of transmission line of electricity and its corresponding line parameter circuit value can be used in the quantity of unit and its corresponding unit capacity, system；
(2) according to the data of acquisition, the corresponding weighted network connection matrix M of system is established；
When establishing the weighted network connection matrix M corresponding to system, weighted network connection matrix M be a square formation, and line number with
Columns is equal to the quantity of available power plant and power transformation tiny node in system；
When there are available transmission line of electricity l between system interior joint i and node jijWhen, if circuit lijLine reactance perunit value
For xij, corresponding element is M in weighted network connection matrix Mij, and have Mij=Mji=xij；If conversely, system interior joint
When available transmission line of electricity being not present between i and node j, Mij=Mji=∞；
(3) filter out available black starting-up unit in system, calculate the sum of capacity of all available black starting-up units and
Corresponding perunit value；
(4) each unit to be restored in system is calculated to the average beeline and load bus of black starting-up unit to each unit
The average beeline of node, and calculate above-mentioned two and be averaged the average distance of beeline；
The average beeline of each unit to be restored to black starting-up unit black is opened according to each unit to be restored to each in system
The minimum value of the beeline of motivation group determines；
Calculating the beeline Dis of unit i to black starting-up unit j to be restoredijWhen, using node where black starting-up unit j as
Node in weighted network connection matrix M is divided into two classes, divided by source point s, node where unit i to be restored as destination node t
Set S and set U are not denoted as：Set S is the set of the destination node for the shortest path for having sought out source point s, and set U is
The destination node set to the shortest path of source point s has not yet been calculated out；
The beeline Dis of unit i to black starting-up unit j to be restored is calculated using critical path method (CPM)ij, specific method is as follows：
1) during initial calculation, source point s has been contained only in set S, all nodes in addition to source point s have been included in set U, for set
Any point x in U, D [x] are the distance of source point s to point x, and circuit is in weighted network connection matrix M as between source point s and node x
The numerical value M of middle corresponding elementsx；
2) the point x with minimum D [x] is selected from set U, then moves on to point x in set S from set U；
3) for any node y in set U residue nodes, if there is D [y]>D[x]+Mxy, then by source point s to point y away from
D [x]+M is arranged to from D [y]xy, wherein, MxyCircuit is corresponding in weighted network connection matrix M between node x and node y
The numerical value of element；
4) the point y with minimum D [y] is selected from remaining node, then moves on to point y in set S from set U；
5) step 3) and step 4) are repeated, until all nodes are all shifted in set S in set U；
6) beeline of source point s to destination node t is D [t]；
Load bus to each machine group node average beeline according to unit i to be restored to the beeline of all machine group nodes
Minimum value determine；
Calculated load node to each machine group node average beeline when, load bus is divided into two classes：One kind is to have included
Load bus in the corresponding restoration path of minimum value of each unit to be restored to the beeline of each black starting-up unit；Separately
It is a kind of then be not included in the load bus in above-mentioned path；
First kind load bus is already contained in restoration path, therefore is not considered, and only considers the second type load node, and will
The node that each unit to be restored is included into the corresponding restoration path of minimum value of the beeline of each black starting-up unit carries out
Cluster；
Unit i to be restored is calculated to the beeline of all machine group nodes using critical path method (CPM)；
(5) fallen into a trap according to the perunit value of the sum of the capacity of all available black starting-up units asked in step (3) with step (4)
The average distance of calculation asks for the restoring force of power grid after having a power failure on a large scale, and reference frame is provided for the development plan of power grid.
The method of power system restoration power after 2. a kind of calculating as described in claim 1 is had a power failure on a large scale, it is characterized in that, the step (1)
In, the method that availability diagnosis is carried out to hair power transmission and transforming equipment all in system is：
Operating status residing for hair power transmission and transforming equipment all in system is identified and judged, is filtered out in malfunction
Or the equipment of inspecting state, i.e., not available equipment in recovery process.
The method of power system restoration power after 3. a kind of calculating as described in claim 1 is had a power failure on a large scale, it is characterized in that, the step (5)
In, the restoring force of power grid is after having a power failure on a large scale：The perunit value of the sum of the capacity of all available black starting-up units is fallen into a trap with step (4)
The quotient of the average distance of calculation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610038041.9A CN105552899B (en) | 2016-01-20 | 2016-01-20 | A kind of method of power system restoration power after calculating is had a power failure on a large scale |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610038041.9A CN105552899B (en) | 2016-01-20 | 2016-01-20 | A kind of method of power system restoration power after calculating is had a power failure on a large scale |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105552899A CN105552899A (en) | 2016-05-04 |
| CN105552899B true CN105552899B (en) | 2018-05-25 |
Family
ID=55831913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610038041.9A Expired - Fee Related CN105552899B (en) | 2016-01-20 | 2016-01-20 | A kind of method of power system restoration power after calculating is had a power failure on a large scale |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105552899B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108631306B (en) * | 2018-05-21 | 2020-03-13 | 西安交通大学 | Method for evaluating recovery capability of power system after disaster |
| CN108512226B (en) * | 2018-05-21 | 2020-06-19 | 西安交通大学 | Method for evaluating resilience of power system under disaster |
| CN110676869B (en) * | 2018-07-02 | 2022-10-11 | 国网山东省电力公司潍坊供电公司 | Method for determining distribution points of newly-increased black start unit |
| CN108879785B (en) * | 2018-08-10 | 2021-09-07 | 国网山东省电力公司潍坊供电公司 | Power grid resilience calculation method and system based on load recovery process |
| CN109193751B (en) * | 2018-08-10 | 2022-02-22 | 国网山东省电力公司潍坊供电公司 | Power grid resilience calculation method and system based on black start and load recovery process |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102904256A (en) * | 2012-10-23 | 2013-01-30 | 广东电网公司电力调度控制中心 | Method and system for rapid self-healing of power grid |
| CN102904251A (en) * | 2012-10-23 | 2013-01-30 | 广东电网公司电力调度控制中心 | Method for system for configuring thermal power generating units to carry out self-healing of power grid |
| CN103746366A (en) * | 2012-11-29 | 2014-04-23 | 浙江省电力公司 | Electric power system recovery method based on node importance degree |
| CN103792924A (en) * | 2014-02-14 | 2014-05-14 | 暨南大学 | Black-start method for expansion of electric power system with micro-grids |
| CN104424510A (en) * | 2013-09-06 | 2015-03-18 | 华北电力大学（保定） | Power transmission grid frame recovering method based on plurality of intelligent bodies and cluster simulation of thermal power generating unit |
| WO2015090375A1 (en) * | 2013-12-18 | 2015-06-25 | Abb Technology Ltd | Microgrid black-start |
-
2016
- 2016-01-20 CN CN201610038041.9A patent/CN105552899B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102904256A (en) * | 2012-10-23 | 2013-01-30 | 广东电网公司电力调度控制中心 | Method and system for rapid self-healing of power grid |
| CN102904251A (en) * | 2012-10-23 | 2013-01-30 | 广东电网公司电力调度控制中心 | Method for system for configuring thermal power generating units to carry out self-healing of power grid |
| CN103746366A (en) * | 2012-11-29 | 2014-04-23 | 浙江省电力公司 | Electric power system recovery method based on node importance degree |
| CN104424510A (en) * | 2013-09-06 | 2015-03-18 | 华北电力大学（保定） | Power transmission grid frame recovering method based on plurality of intelligent bodies and cluster simulation of thermal power generating unit |
| WO2015090375A1 (en) * | 2013-12-18 | 2015-06-25 | Abb Technology Ltd | Microgrid black-start |
| CN103792924A (en) * | 2014-02-14 | 2014-05-14 | 暨南大学 | Black-start method for expansion of electric power system with micro-grids |
Non-Patent Citations (1)
| Title |
|---|
| 《考虑机组恢复成功率的黑启动阶段待恢复机组优选》;朱海南等;《电力***自动化》;20131125;第37卷(第22期);第28-34页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105552899A (en) | 2016-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105552899B (en) | A kind of method of power system restoration power after calculating is had a power failure on a large scale | |
| CN105656040B (en) | A kind of power system restoration power computational methods for considering black starting-up process | |
| CN103296677B (en) | A kind of online bulk power grid recovers aid decision-making system | |
| CN108631306B (en) | Method for evaluating recovery capability of power system after disaster | |
| CN101917001B (en) | Fault sequencing method for on-line static security analysis of power system | |
| CN108226713B (en) | Concentric relaxation depression domain analysis method for voltage sag | |
| CN109066659B (en) | Microgrid island operation reliability assessment method and terminal equipment | |
| CN109698500B (en) | Power distribution network power supply reliability improving method based on intelligent soft switch | |
| CN104716646B (en) | A kind of node Coupling Degrees method based on Injection Current | |
| CN104750878A (en) | Mixed searching strategy-based topology fault diagnosis method | |
| CN100550564C (en) | Electric power system transient stability and thermally-stabilised coordination emergency control method | |
| CN107622360A (en) | A kind of critical circuits recognition methods for considering subjective and objective factor | |
| CN110907754A (en) | Fault line severity evaluation method based on PSD-BPA | |
| CN106972517A (en) | Reliability of UHVDC transmission system computational methods based on bipolar symmetrical feature | |
| CN105656036A (en) | Probability static safety analysis method considering flow-and-sensitivity consistency equivalence | |
| CN113328437A (en) | CPS (control performance Standard) topology construction method and fault recovery method for intelligent power distribution network | |
| CN110412417B (en) | Micro-grid data fault diagnosis method based on intelligent power monitoring instrument | |
| CN109193751B (en) | Power grid resilience calculation method and system based on black start and load recovery process | |
| CN113098012B (en) | Regional autonomy capability assessment method for DG-containing power distribution network | |
| CN103296678B (en) | A kind of online bulk power grid recovers aid decision-making method | |
| CN111276967B (en) | Elasticity capability evaluation method and device of power system | |
| CN116014790A (en) | Elasticity evaluation method for power distribution network containing distributed energy sources | |
| Bretas et al. | Bi-level linear programming model for automatic load shedding: A distributed wide-area measurement system-based solution | |
| CN110189038B (en) | Power grid risk analysis method and system, computer device and storage medium | |
| Lang et al. | Optimal allocation of voltage sag monitors in multi voltage grades distribution network based on improved MRA method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180525 |