CN105470986A - Power system partitioning method - Google Patents

Power system partitioning method Download PDF

Info

Publication number
CN105470986A
CN105470986A CN201510955941.5A CN201510955941A CN105470986A CN 105470986 A CN105470986 A CN 105470986A CN 201510955941 A CN201510955941 A CN 201510955941A CN 105470986 A CN105470986 A CN 105470986A
Authority
CN
China
Prior art keywords
bus
isolated island
isolated
pmu
unobservable
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
Application number
CN201510955941.5A
Other languages
Chinese (zh)
Other versions
CN105470986B (en
Inventor
高峰
林昌年
魏文辉
王国平
金逸
费益军
马强
杨选怀
韩佳兵
陈捷
周书进
刘传良
慈文斌
史立勤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Shandong Electric Power Co Ltd
State Grid Jiangsu Electric Power Co Ltd
Beijing Kedong Electric Power Control System Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Shandong Electric Power Co Ltd
State Grid Jiangsu Electric Power Co Ltd
Beijing Kedong Electric Power Control System Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Shandong Electric Power Co Ltd, State Grid Jiangsu Electric Power Co Ltd, Beijing Kedong Electric Power Control System Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201510955941.5A priority Critical patent/CN105470986B/en
Publication of CN105470986A publication Critical patent/CN105470986A/en
Application granted granted Critical
Publication of CN105470986B publication Critical patent/CN105470986B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, 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 invention provides a power system partitioning method in a system recovery process. In order to increase the system recovery speed, the invention provides a partitioning method which considers various constraints such as a capacity constraint of a black start power supply, power balance of an island, and island observability based on distribution of PMU in a power grid in a partitioning process. According to the method, applying a wide area measurement system (WARMS) enables the island formed in the system recovery process satisfy the observability requirement. According to the method, after the network is preliminarily partitioned, buses which do not satisfy the observability requirement are processed according to a rule, thereby finally obtaining a power system partitioning scheme which satisfies the constraints. The invention provides the power network partitioning method which uses the optimum PMU distribution in the power network, thereby accelerating system recovery in a black start process. The method ensures the observability of the formed island and satisfies various operating constraints comprising a classical load constraint and an observability constraint in WARMS.

Description

A kind of electric power system partition method
Technical field
The present invention relates to a kind of electric power system partition method, especially relate to a kind of black starting-up partition method of base PMU best distribution.
Background technology
Power system blackstart refers to that whole network system is because of after fault stoppage in transit, by there is self-startup ability also referred to as the startup with black start-up ability unit in system, drive the unit of non self starting, and expand system power supply scope gradually, finally realize the recovery of whole system.
When black starting-up power supply point has been determined and system has multiple black starting-up power supply point, in order to improve the resume speed of electrical network, should according to the design feature of electrical network, large scale electric network is divided into two or more subsystem to recover independently, grid-connected after each subsystem stable operation a period of time, and then recover whole system.Take parallel recovery plan, greatly can shorten the recovery time of electrical network, if any one subsystem causes recovering unsuccessfully because of some unpredictable factor, also do not affect the recovering process of other subsystem, this provides for improved the speed of System recover.Therefore, Study system have a power failure on a large scale after recovery process in how subregion is carried out to whole network, namely how whole Iarge-scale system is divided into several subsystems, implement parallel recovery, have important meaning to the fast quick-recovery realizing whole electric power system.
In the middle of the research carried out at present, relate to the little of black starting-up partitioning problem.The general artificial appointment partitioned mode adopted by rule of thumb, lacked the foundation of science in the past.Conventional method has application ordered binary decision diagram to carry out model solution, the method carries out reasonable off-the-line according to established constraint on the basis of the unit and important load that set needs startup, determines that in backbone network, the connected state of circuit realizes subregion.Also have and use black starting-up power supplys all in dummy node equivalent system, call after traditional Bellman-ford algorithm and Prim algorithm construct and recover rack, then dummy node is reduced with the object reaching subregion.But the constraints that these researchs are considered is little, do not consider whether subnet meets the problem running constraint yet.Visible, the research in system subdivision is in the starting stage of a development, and a lot of detailed problem also needs deeply and carefully to inquire into.
Summary of the invention
Object of the present invention is exactly a kind of method in order to obtain electric power system subregion, ensures subregion observability, and meets institute's Constrained.
Object of the present invention can be achieved through the following technical solutions:
1 .the invention provides a kind of method of electric power system subregion, step is as follows:
1) first this method determines the optimum PMU distribution ensureing electric power networks observability;
2) the preliminary partition method of electric power networks is carried out;
3) partition scheme is modified be met constrained partition scheme.
2. step 1) described in problem for solving minimum PMU distribution problem:
If f is a 0-1 variable, be used for the distribution situation of PMU in expression system, i-th element definition of f is:
Optimum PMU distribution problem can represent with following mathematic(al) representation:
N represents n-th bus, A=[a ij] represent network adjacent matrix, be defined as follows:
3. after trying to achieve PMU best distribution, ask for the preliminary partition scheme of electric power system:
Do not contact each other to allow between all isolated islands, each transformer and coupled bus should be placed in same isolated island.Therefore transformer and the bus at its two ends, the generator be connected with these buses and load are all equivalent to the bus of a band equivalent generator be attached thereto and equivalent load.For New England 39 bus-bar system in accompanying drawing 1, bus 19,20,33,34 is by bus 17 equivalence in accompanying drawing 2, and bus 17 is used for replacing generator G4 and G5 in accompanying drawing 1 and the load on bus 20 with an equivalent generator and an equivalent load.When hereafter carrying out mission observability detection, the generator of equivalence, load and bus are reduced.Accompanying drawing 1 is carried out simplification and obtains accompanying drawing 2.
In order to carry out the System recover in black starting-up process, at least should comprise a black starting-up power supply in each isolated island, one is used for the load that recovers and ensure the PMU of mission observability.Determine possible maximum isolated island number m:
B iand l ibe defined as follows:
First two concepts are put forward:
Center bus: center bus is the center of each isolated island in this algorithm, it must have the combination of generator or PMU or load or three to hang on this bus.
Boundary line: the connecting line between isolated island, will open in the process of carrying out subregion.
To accompanying drawing 2 application of formula (5), assuming that black starting-up power supply is positioned at bus 2,17,20,22.Center bus number equals isolated island number.Due to so m=4.
The first step recursive matrix I of isolated island 1defining like that according to formula (8), is a upper triangular matrix.
If the i-th j in matrix element representation bus j belongs to isolated island i, be set to 1, otherwise be 0.I can be found out 1each provisional capital have a non-zero entry usually to represent the center bus of this isolated island, like this center bus is bus 2,17,20 and 22.
By I 1modify to carry out recursion next time, amendment specification is as follows:
I′ (k+1)=I k×A(9)
To I ' (k+1)apply following rule and obtain I (k+1):
Center bus rule: because each isolated island needs a center bus, then this bus position element zero setting in other isolated islands.
Historical rule: then next step can not be placed on other isolated islands if certain bus is placed in certain isolated island by back.So, show multiple 1 if a certain, except back being element zero setting except 1.
The minimum rule of bus: if certain bus is placed in how each isolated island and is first time appear in these isolated islands in a certain step, then bus will be placed in the isolated island of bus negligible amounts to accelerate resume speed.Idea bus quantity here refers to the bus quantity containing transformer bus in accompanying drawing 1.
Connect closely rule: if two isolated islands have identical bus number in previous rule, then this bus should be placed in the larger isolated island of the value of this bus position.Because representative has more bus to be connected with this bus by boundary line in this isolated island like this.
The shorter rule in boundary line: if in previous rule in all isolated islands this bus to be sitting in the value of position identical, then this bus should be placed in the shorter isolated island in the boundary line that is connected with this bus.
After above step terminates, the element being greater than 1 in matrix is put 1, so obtain next step matrix.Carry out this process always and know that all buses all belong to a not isolated island.
Accompanying drawing 3 represents the result of the preliminary subregion of 39 bus-bar system.
4. by step 3) try to achieve network primary partition after check the load constraint of isolated island:
After network is divided into different isolated island, each isolated island should meet and comprises WARMS and constrain in interior constraints, ensures completing smoothly of System recover.Constraint set is as follows:
1) load meets the minimum energy output of isolated island.
This constraint should represent by mathematical formulae as follows:
Wherein, J iisolated island i bus number, the minimum load of generator j, P liit is the workload demand of bus j.
If certain isolated island does not meet constraint (10), the boundary line of this this isolated island of inspection.If the non-central bus that adjacent isolated island exists the suitable load of band makes to meet constraint (10), and the adjacent isolated island removing this bus also meets constraint (10), then this bus moves on to the isolated island not meeting constraint.
To 39 node systems, suppose that the minimum load of every platform generator is 10% of its rated capacity, for the partition scheme of accompanying drawing 3, each isolated island meets constraint (10).
2) observability constraint.
In order to check the observability of isolated island, the element zero setting that the boundary line in formula (4) goes out is obtained matrix A i.The following observability representing isolated island median generatrix with variable ο:
ο=A I×f(11)
The bus that 0 element representation in ο is unobservable.If there is unobservable bus, by following formula definition matrix U:
U=A u×PIM(12)
Wherein, A ube the network connection matrix obtained by formula (4), u is bus number unobservable in all isolated islands, equals the quantity of all 0 elements in ο.
It is as follows that definition PMU installs matrix PIM:
Each line display unobservable bus in U, and have individual nonzero element to represent PMU is installed makes it possible to be that unobservable bus becomes Observable.
For making unobservable bus become observable, first define the observability degree of isolated island:
Wherein, ο ja jth element of ο, w jbe the observability degree of bus j, be defined as follows:
Wop ivalue higher showing has more bus to be observation in isolated island i.
Unobservable bus is become observable by following three rules of application:
1) unobservable bus Bu Shi center bus
Unobservable bus should be transferred to the bus that comprises and installed PMU and is connected with this bus by boundary line.If (10) do not meet after transfer, then calculate the wop of the unobservable bus of this isolated island ivalue.If be greater than 0.9, unobservable bus is retracted original isolated island, because observability degree accelerates resume speed enough greatly.If be less than or equal to 0.9, two isolated islands are merged.
2) unobservable bus is center bus but does not install PMU
The bus having installed PMU should be transferred in the isolated island comprising this unobservable bus.If (10) are broken after transfer, calculate the wop of unobservable bus place isolated island ivalue.If be greater than 0.9 the bus installing PMU is retracted original isolated island, if be less than or equal to, two isolated islands are merged.
3) unobservable bus and the bus Dou Shi center bus having installed PMU
If wop ibe less than or equal to 0.9 two isolated islands are merged.
Form the flow chart of complete partition scheme as shown in Figure 5.
Use formula (11) can show that bus 4,18,27 is unobservable to 39 node buses of accompanying drawing 3.Bus 4 Bu Shi center bus, the Observable that it can be become by the bus 14 having installed PMU.So bus 4 is moved on to the isolated island comprising bus 14.Bus 18 neither center bus, and bus 17 can make its Observable, so bus 18 should move on to bus 17 place isolated island.Bus 27 Bu Shi center bus also should move on to bus 17 place isolated island.The partition scheme of accompanying drawing 4 can be obtained through these adjustment.
Compared with prior art, the present invention has the following advantages:
The present invention proposes a kind of by utilizing the electric power networking partition method of optimum PMU distribution in electric power networks, accelerating the System recover in black starting-up process with this.This method ensure that the observability of formed isolated island, and meet the multiple observability constraint run in the load constraint and WARMS retraining and comprise classics.
Accompanying drawing explanation
Fig. 1 is New England 39 bus-bar system figure;
Fig. 2 is the simplified system diagram after New England 39 bus-bar system equivalent process transformer and load;
Fig. 3 is preliminary division result;
Fig. 4 is final division result;
Fig. 5 is subregion process flow diagram.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.The present embodiment is implemented premised on technical solution of the present invention, give detailed execution mode and concrete operating process, but protection scope of the present invention is not limited to following embodiment.
As shown in Figure 5, the present embodiment provides a kind of electric power system partition method, comprises step:
S1, solve minimum PMU distribution problem:
If f is a 0-1 variable, be used for the distribution situation of PMU in expression system, i-th element definition of f is:
Optimum PMU distribution problem can represent with following mathematic(al) representation:
N represents n-th bus, A=[a ij] represent network adjacent matrix, be defined as follows:
S2, try to achieve PMU best distribution after ask for the preliminary partition scheme of electric power system:
Do not contact each other to allow between all isolated islands, each transformer and coupled bus should be placed in same isolated island.Therefore transformer and the bus at its two ends, the generator be connected with these buses and load are all equivalent to the bus of a band equivalent generator be attached thereto and equivalent load.For New England 39 bus-bar system in accompanying drawing 1, bus 19,20,33,34 is by bus 17 equivalence in accompanying drawing 2, and bus 17 is used for replacing generator G4 and G5 in accompanying drawing 1 and the load on bus 20 with an equivalent generator and an equivalent load.When hereafter carrying out mission observability detection, the generator of equivalence, load and bus are reduced.Accompanying drawing 1 is carried out simplification and obtains accompanying drawing 2.
In order to carry out the System recover in black starting-up process, at least should comprise a black starting-up power supply in each isolated island, one is used for the load that recovers and ensure the PMU of mission observability.Determine possible maximum isolated island number m:
B iand l ibe defined as follows:
First two concepts are put forward:
Center bus: center bus is the center of each isolated island in this algorithm, it must have the combination of generator or PMU or load or three to hang on this bus.
Boundary line: the connecting line between isolated island, will open in the process of carrying out subregion.
To accompanying drawing 2 application of formula (5), assuming that black starting-up power supply is positioned at bus 2,17,20,22.Center bus number equals isolated island number.Due to so m=4.
The first step recursive matrix I of isolated island 1defining like that according to formula (8), is a upper triangular matrix.
If the i-th j in matrix element representation bus j belongs to isolated island i, be set to 1, otherwise be 0.I can be found out 1each provisional capital have a non-zero entry usually to represent the center bus of this isolated island, like this center bus is bus 2,17,20 and 22.
By I 1modify to carry out recursion next time, amendment specification is as follows:
I′ (k+1)=I k×A(9)
To I ' (k+1)apply following rule and obtain I (k+1):
Center bus rule: because each isolated island needs a center bus, then this bus position element zero setting in other isolated islands.
Historical rule: then next step can not be placed on other isolated islands if certain bus is placed in certain isolated island by back.So, show multiple 1 if a certain, except back being element zero setting except 1.
The minimum rule of bus: if certain bus is placed in how each isolated island and is first time appear in these isolated islands in a certain step, then bus will be placed in the isolated island of bus negligible amounts to accelerate resume speed.Idea bus quantity here refers to the bus quantity containing transformer bus in accompanying drawing 1.
Connect closely rule: if two isolated islands have identical bus number in previous rule, then this bus should be placed in the larger isolated island of the value of this bus position.Because representative has more bus to be connected with this bus by boundary line in this isolated island like this.
The shorter rule in boundary line: if in previous rule in all isolated islands this bus to be sitting in the value of position identical, then this bus should be placed in the shorter isolated island in the boundary line that is connected with this bus.
After above step terminates, the element being greater than 1 in matrix is put 1, so obtain next step matrix.Carry out this process always and know that all buses all belong to a not isolated island.
Accompanying drawing 3 represents the result of the preliminary subregion of 39 bus-bar system.
S3, try to achieve network primary partition after check isolated island load constraint:
After network is divided into different isolated island, each isolated island should meet and comprises WARMS and constrain in interior constraints, ensures completing smoothly of System recover.Constraint set is as follows:
1) load meets the minimum energy output of isolated island.
This constraint should represent by mathematical formulae as follows:
Wherein, J iisolated island i bus number, the minimum load of generator j, P liit is the workload demand of bus j.
If certain isolated island does not meet constraint (10), the boundary line of this this isolated island of inspection.If the non-central bus that adjacent isolated island exists the suitable load of band makes to meet constraint (10), and the adjacent isolated island removing this bus also meets constraint (10), then this bus moves on to the isolated island not meeting constraint.
To 39 node systems, suppose that the minimum load of every platform generator is 10% of its rated capacity, for the partition scheme of accompanying drawing 3, each isolated island meets constraint (10).
2) observability constraint.
In order to check the observability of isolated island, the element zero setting that the boundary line in formula (4) goes out is obtained matrix A i.The following observability representing isolated island median generatrix with variable ο:
ο=A I×f(11)
The bus that 0 element representation in ο is unobservable.If there is unobservable bus, by following formula definition matrix U:
U=A u×PIM(12)
Wherein, A ube the network connection matrix obtained by formula (4), u is bus number unobservable in all isolated islands, equals the quantity of all 0 elements in ο.
It is as follows that definition PMU installs matrix PIM:
Each line display unobservable bus in U, and have individual nonzero element to represent PMU is installed makes it possible to be that unobservable bus becomes Observable.
For making unobservable bus become observable, first define the observability degree of isolated island:
Wherein, ο ja jth element of ο, w jbe the observability degree of bus j, be defined as follows:
Wop ivalue higher showing has more bus to be observation in isolated island i.
Unobservable bus is become observable by following three rules of application:
1) unobservable bus Bu Shi center bus
Unobservable bus should be transferred to the bus that comprises and installed PMU and is connected with this bus by boundary line.If (10) do not meet after transfer, then calculate the wop of the unobservable bus of this isolated island ivalue.If be greater than 0.9, unobservable bus is retracted original isolated island, because observability degree accelerates resume speed enough greatly.If be less than or equal to 0.9, two isolated islands are merged.
2) unobservable bus is center bus but does not install PMU
The bus having installed PMU should be transferred in the isolated island comprising this unobservable bus.If (10) are broken after transfer, calculate the wop of unobservable bus place isolated island ivalue.If be greater than 0.9 the bus installing PMU is retracted original isolated island, if be less than or equal to, two isolated islands are merged.
3) unobservable bus and the bus Dou Shi center bus having installed PMU
If wop ibe less than or equal to 0.9 two isolated islands are merged.
Form the flow chart of complete partition scheme as shown in Figure 5.
Use formula (11) can show that bus 4,18,27 is unobservable to 39 node buses of accompanying drawing 3.Bus 4 Bu Shi center bus, the Observable that it can be become by the bus 14 having installed PMU.So bus 4 is moved on to the isolated island comprising bus 14.Bus 18 neither center bus, and bus 17 can make its Observable, so bus 18 should move on to bus 17 place isolated island.Bus 27 Bu Shi center bus also should move on to bus 17 place isolated island.The partition scheme of accompanying drawing 4 can be obtained through these adjustment.
Above-described specific embodiment is only and illustrates and of the present inventionly realize effect, not in order to limit the present invention.The amendment of any unsubstantiality done within all basic ideas in method proposed by the invention and framework, conversion and improvement, all should be included within protection scope of the present invention.

Claims (4)

1. a method for electric power system subregion, is characterized in that, comprises the steps:
1) the optimum PMU distribution ensureing electric power networks observability is determined;
2) the preliminary partition method of electric power networks is carried out;
3) partition scheme is modified be met constrained partition scheme.
2. the method for electric power system subregion according to claim 1, is characterized in that, described step 1) in, solve minimum PMU distribution problem:
If f is a 0-1 variable, be used for the distribution situation of PMU in expression system, i-th element definition of f is:
Optimum PMU distribution problem can represent with following mathematic(al) representation:
m i n Σ i = 1 n f i - - - ( 2 )
s . t . A × f ≥ 1 1 ... 1 1 × n T - - - ( 3 )
N represents n-th bus, A=[a ij] represent network adjacent matrix, be defined as follows:
3. the method for electric power system subregion according to claim 1, is characterized in that, asks for electric power system preliminary partition scheme concrete steps as follows after trying to achieve PMU best distribution:
Do not contact each other to allow between all isolated islands, each transformer and coupled bus should be placed in same isolated island; Therefore transformer and the bus at its two ends, the generator be connected with these buses and load are all equivalent to the bus of a band equivalent generator be attached thereto and equivalent load; In order to carry out the System recover in black starting-up process, at least should comprise a black starting-up power supply in each isolated island, one is used for the load that recovers and ensure the PMU of mission observability; Determine possible maximum isolated island number m:
m = m i n Σ i = 1 n b i , Σ i = 1 n l i Σ i = 1 n f i - - - ( 5 )
B iand l ibe defined as follows:
First two concepts are put forward:
Center bus: center bus is the center of each isolated island in this algorithm, it must have the combination of generator or PMU or load or three to hang on this bus;
Boundary line: the connecting line between isolated island, will open in the process of carrying out subregion;
The first step recursive matrix I of isolated island 1defining like that according to formula (8), is a upper triangular matrix;
I 1 = G = g 11 ... g 1 n . . . g m 1 ... g m n = b 1 b 2 ( 1 - g 11 ) ... b n ( 1 - g 11 ) ... ( 1 - g 1 ( n - 1 ) ) 0 b 2 ( 1 - g 12 ) b n ( 1 - g 1 n ) ( 1 - g 22 ) ... ( 1 - g 2 ( n - 1 ) ) 0 0 b n ( 1 - g 1 n ) ( 1 - g 2 n ) ( 1 - g 33 ) ... ( 1 - g 3 ( n - 1 ) ) . . . 0 0 b n ( 1 - g 1 n ) ... ( 1 - g ( n - 1 ) n ) - - - ( 8 )
If the i-th j in matrix element representation bus j belongs to isolated island i, be set to 1, otherwise be 0; I can be found out 1each provisional capital have a non-zero entry usually to represent the center bus of this isolated island;
By I 1modify to carry out recursion next time, amendment specification is as follows:
I′ (k+1)=I k×A(9)
To I ' (k+1)apply following rule and obtain I (k+1):
Center bus rule: because each isolated island needs a center bus, then this bus position element zero setting in other isolated islands;
Historical rule: then next step can not be placed on other isolated islands if certain bus is placed in certain isolated island by back; So, show multiple 1 if a certain, except back being element zero setting except 1;
The minimum rule of bus: if certain bus is placed in how each isolated island and is first time appear in these isolated islands in a certain step, then bus will be placed in the isolated island of bus negligible amounts to accelerate resume speed; Idea bus quantity here refers to the bus quantity containing transformer bus in accompanying drawing 1;
Connect closely rule: if two isolated islands have identical bus number in previous rule, then this bus should be placed in the larger isolated island of the value of this bus position; Because representative has more bus to be connected with this bus by boundary line in this isolated island like this;
The shorter rule in boundary line: if in previous rule in all isolated islands this bus to be sitting in the value of position identical, then this bus should be placed in the shorter isolated island in the boundary line that is connected with this bus;
After above step terminates, the element being greater than 1 in matrix is put 1, so obtain next step matrix; Carry out this process until all buses all belong to an isolated island always.
4. the method for electric power system subregion according to claim 1, is characterized in that, described step 3) try to achieve network primary partition after check isolated island load constraint:
After network is divided into different isolated island, each isolated island should meet and comprises WARMS and constrain in interior constraints, ensures completing smoothly of System recover; Constraint set is as follows:
1) load meets the minimum energy output of isolated island;
This constraint should represent by mathematical formulae as follows:
Σ j = 1 J i P G j min ≤ Σ j = 1 J i P L j , i = 1 , 2 ... m - - - ( 10 )
Wherein, J iisolated island i bus number, the minimum load of generator j, P liit is the workload demand of bus j;
If certain isolated island does not meet constraint (10), the boundary line of this this isolated island of inspection; If the non-central bus that adjacent isolated island exists the suitable load of band makes to meet constraint (10), and the adjacent isolated island removing this bus also meets constraint (10), then this bus moves on to the isolated island not meeting constraint;
To 39 node systems, suppose that the minimum load of every platform generator is 10% of its rated capacity, for the partition scheme of accompanying drawing 3, each isolated island meets constraint (10);
2) observability constraint;
In order to check the observability of isolated island, the element zero setting that the boundary line in formula (4) goes out is obtained matrix A i; The following observability representing isolated island median generatrix with variable ο:
ο=A I×f(11)
The bus that 0 element representation in ο is unobservable; If there is unobservable bus, by following formula definition matrix U:
U=A u×PIM(12)
Wherein, A ube the network connection matrix obtained by formula (4), u is bus number unobservable in all isolated islands, equals the quantity of all 0 elements in ο;
It is as follows that definition PMU installs matrix PIM:
P I M = f 1 f 2 ... f n n × n - - - ( 13 )
Each line display unobservable bus in U, and have individual nonzero element to represent PMU is installed makes it possible to be that unobservable bus becomes Observable;
For making unobservable bus become observable, first define the observability degree of isolated island:
wop i = Σ j = 1 J i w j × o j w j , i = 1 , ... m - - - ( 14 )
Wherein, ο ja jth element of ο, w jbe the observability degree of bus j, be defined as follows:
Wop ivalue higher showing has more bus to be observation in isolated island i;
Unobservable bus is become observable by following three rules of application:
1) unobservable bus Bu Shi center bus
Unobservable bus should be transferred to the bus that comprises and installed PMU and is connected with this bus by boundary line; If (10) do not meet after transfer, then calculate the wop of the unobservable bus of this isolated island ivalue; If be greater than 0.9, unobservable bus is retracted original isolated island, because observability degree accelerates resume speed enough greatly; If be less than or equal to 0.9, two isolated islands are merged;
2) unobservable bus is center bus but does not install PMU
The bus having installed PMU should be transferred in the isolated island comprising this unobservable bus; If (10) are broken after transfer, calculate the wop of unobservable bus place isolated island ivalue; If be greater than 0.9 the bus installing PMU is retracted original isolated island, if be less than or equal to, two isolated islands are merged;
3) unobservable bus and the bus Dou Shi center bus having installed PMU
If wop ibe less than or equal to 0.9 two isolated islands are merged.
CN201510955941.5A 2015-12-17 2015-12-17 Power system partitioning method Expired - Fee Related CN105470986B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510955941.5A CN105470986B (en) 2015-12-17 2015-12-17 Power system partitioning method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510955941.5A CN105470986B (en) 2015-12-17 2015-12-17 Power system partitioning method

Publications (2)

Publication Number Publication Date
CN105470986A true CN105470986A (en) 2016-04-06
CN105470986B CN105470986B (en) 2020-06-19

Family

ID=55608431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510955941.5A Expired - Fee Related CN105470986B (en) 2015-12-17 2015-12-17 Power system partitioning method

Country Status (1)

Country Link
CN (1) CN105470986B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106410789A (en) * 2016-10-13 2017-02-15 东南大学 Power system black-start portioning method for isolated network
CN108734599A (en) * 2017-04-18 2018-11-02 中国科学院深圳先进技术研究院 Power grid island effect detecting system and method
CN108808670A (en) * 2018-07-03 2018-11-13 中国电力科学研究院有限公司 It is a kind of have a power failure on a large scale after power system recovery partition scheme optimization method and device
CN113949099A (en) * 2021-10-28 2022-01-18 国网上海市电力公司 Real-time controlled island detection and recovery method based on state estimation

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102509168A (en) * 2011-10-06 2012-06-20 武陟县电业总公司 Improved integer programming method-based optimum PMU (Power Management Unit) placement method
CN102611102A (en) * 2012-03-13 2012-07-25 河海大学 Power system partition modeling method based on actual measurement by aid of PMU (phasor measurement units)
CN103248043A (en) * 2013-04-28 2013-08-14 西南交通大学 Power system multi-zone distributed state estimation method based on synchronous phase angle measurement device
CN104092212A (en) * 2014-07-24 2014-10-08 河海大学 Electric system multi-domain distributed state estimation method based on PMU measurement
CN104362613A (en) * 2014-09-28 2015-02-18 国家电网公司 Power system observability judgment method by considering transformer tap position estimation
US20150355655A1 (en) * 2014-06-06 2015-12-10 Shanghai Jiao Tong University Method for optimizing the flexible constraints of an electric power system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102509168A (en) * 2011-10-06 2012-06-20 武陟县电业总公司 Improved integer programming method-based optimum PMU (Power Management Unit) placement method
CN102611102A (en) * 2012-03-13 2012-07-25 河海大学 Power system partition modeling method based on actual measurement by aid of PMU (phasor measurement units)
CN103248043A (en) * 2013-04-28 2013-08-14 西南交通大学 Power system multi-zone distributed state estimation method based on synchronous phase angle measurement device
US20150355655A1 (en) * 2014-06-06 2015-12-10 Shanghai Jiao Tong University Method for optimizing the flexible constraints of an electric power system
CN104092212A (en) * 2014-07-24 2014-10-08 河海大学 Electric system multi-domain distributed state estimation method based on PMU measurement
CN104362613A (en) * 2014-09-28 2015-02-18 国家电网公司 Power system observability judgment method by considering transformer tap position estimation

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106410789A (en) * 2016-10-13 2017-02-15 东南大学 Power system black-start portioning method for isolated network
CN106410789B (en) * 2016-10-13 2018-12-28 东南大学 A kind of power system blackstart partition method for isolated network
CN108734599A (en) * 2017-04-18 2018-11-02 中国科学院深圳先进技术研究院 Power grid island effect detecting system and method
CN108808670A (en) * 2018-07-03 2018-11-13 中国电力科学研究院有限公司 It is a kind of have a power failure on a large scale after power system recovery partition scheme optimization method and device
CN113949099A (en) * 2021-10-28 2022-01-18 国网上海市电力公司 Real-time controlled island detection and recovery method based on state estimation
CN113949099B (en) * 2021-10-28 2024-05-10 国网上海市电力公司 Real-time controlled island detection and recovery method based on state estimation

Also Published As

Publication number Publication date
CN105470986B (en) 2020-06-19

Similar Documents

Publication Publication Date Title
Sarmadi et al. A sectionalizing method in power system restoration based on WAMS
CN105470986A (en) Power system partitioning method
CN101860037B (en) Determination method of network side harmonic current of high-voltage direct current power transmission system
CN109742779B (en) Distributed energy storage scheduling method and device
CN105870949B (en) A kind of micro-capacitance sensor energy-storage units optimal control method based on distributed gradient algorithm
CN104537258A (en) Cone optimization modeling method for allowing distributed stored energy to participate in running adjustment of active power distribution network
CN107666155A (en) System Stochastic Stability Analysis method of providing multiple forms of energy to complement each other based on Markov model
CN108493930A (en) The load restoration two-phase optimization method of meter and wind power integration
CN107834602A (en) A kind of micro- source half-bridge current transformer tandem type micro-grid system
CN110797874A (en) State estimation method for alternating current-direct current hybrid power distribution network containing power electronic transformer
CN107294130A (en) A kind of wind storage generating black starting-up system and its energy storage configuration method
CN113708366A (en) Power distribution network optimal scheduling method and system considering new energy and electric automobile
CN106505575A (en) A kind of Line Flow economic load dispatching method based on Granule Computing
CN107332277B (en) Active power distribution network island operation method considering source load storage operation characteristics
CN115333110A (en) Power distribution network-microgrid group collaborative distributed optimization scheduling method and system based on ADMM
CN106159955A (en) Based on the power system distributed optimal power flow method punishing Duality Decomposition continuously
CN105406517B (en) Economic Dispatch method based on finite time average homogeneity algorithm
Liu et al. Consensus-based energy management of multi-microgrid: An improved SoC-based power coordinated control method
CN104090496A (en) Smart grid control operation continuous analog simulation method
Liu et al. Start-up sequence of generators in power system restoration avoiding the backtracking algorithm
CN109031946A (en) A kind of control method and device of mixed energy storage system
CN110417040B (en) Control method of energy management system for large-scale VRB (virtual router bus) optimized operation
CN115313438B (en) AC/DC power transmission network and energy storage collaborative planning method and medium
CN108336731A (en) A kind of computational methods of power distribution network distributed generation resource allowed capacity
CN101964526B (en) Electric network recovery and reconstruction method for underpowered network

Legal Events

Date Code Title Description
C06 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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200619

Termination date: 20201217

CF01 Termination of patent right due to non-payment of annual fee