CN103259265B - Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix - Google Patents

Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix Download PDF

Info

Publication number
CN103259265B
CN103259265B CN201310223708.9A CN201310223708A CN103259265B CN 103259265 B CN103259265 B CN 103259265B CN 201310223708 A CN201310223708 A CN 201310223708A CN 103259265 B CN103259265 B CN 103259265B
Authority
CN
China
Prior art keywords
matrix
distribution network
vector
reachability matrix
island
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
Application number
CN201310223708.9A
Other languages
Chinese (zh)
Other versions
CN103259265A (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.)
Shandong University
Original Assignee
Shandong University
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 Shandong University filed Critical Shandong University
Priority to CN201310223708.9A priority Critical patent/CN103259265B/en
Publication of CN103259265A publication Critical patent/CN103259265A/en
Application granted granted Critical
Publication of CN103259265B publication Critical patent/CN103259265B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a power network partitioning and detached island identification method of an intelligent power distribution network based on a reachability matrix. The method has good power network partitioning capability, and can quickly and accurately detect the quantity and a range of detached islands, and a specific parallel-in switch. The method comprises the following steps that 1, a branch-node diagram is established according to the intelligent power distribution network and used for describing a structure of the intelligent power distribution network; 2, an adjacent matrix A is established and used for describing connection relations between various branches and nodes in the branch-node diagram; 3, each idempotent matrix of the adjacent matrix A is calculated, and a partitioning situation of the network is obtained; and 4, a switch property vector K is defined, a block vector Pi for describing network partitioning information is obtained according to the partitioning situation and value situations of elements in a reachability matrix P; corresponding elements in the vector K and the vector Pi are subjected to logic and calculation; vectors for describing properties of various areas of the network are obtained; whether the detached islands exist is judged; and the quantity of the detached islands, and power supplies corresponding to each detached island are judged.

Description

Based on intelligent distribution network sub area division and the island identification method of reachability matrix
Technical field
The present invention relates to a kind of sub area division and island detection method, especially a kind of intelligent distribution network sub area division based on reachability matrix and island identification method.
Background technology
Power distribution network node is numerous, various DG throws and moves back frequently, cause the operational mode of intelligent distribution network flexible and changeable, especially when some switch in electrical network disconnects, electrical network can be divided into some distinct, separate subnets, if the load in subnet is powered by DG or microgrid, so this subnet just runs in the mode of isolated island.The appearance of isolated island brings considerable influence to the scheduling of intelligent distribution network, automated system, relaying protection system and automatic reclosing etc. on the one hand, and the safety and stability of isolated island self-operating also exists problems on the other hand.
Current island detection method can be divided into passive type to detect and active detection two class generally.Passive type detection method refers to that the change of the power output by directly measuring distributed generation system inverter or the change of PCC place voltage or frequency judge isolated island situation.Passive detection method principle is simple, is easy to realize, can not the quality of power supply of influential system and stability.But the threshold value of these class methods is all more difficult to be determined, and whether they mate with the power of load the power output of inverter and have strict requirement, there is certain check frequency.Active detection method is that when distributed generation system and electrical network are normally incorporated into the power networks, microvariations are not enough to have an impact to system by exporting the little disturbance of generation one to the inverter of distributed generation system; When isolated island produces, the microvariations initiatively produced will cause the change that island state is large.The check frequency of active detecting method has had certain improvement, but owing to introducing disturbance, have certain influence to the transient response of system, in multiple stage grid-connected inverters situation, the consistency of disturbance is also ubiquitous problem.
Current islanding detect and discrimination method study mainly based on and the local information at site place, if using overall for power distribution network as research object, the position of the switch state information in network is utilized in conjunction with network topology structure, electric network composition situation of change can be obtained fast and accurately, and identification can be carried out to isolated island reliably.
Have a kind of based on normalization spectral clustering and the two benches Active Splitting method retraining spectral clustering at present.The method utilizes normalization spectrum to birds of the same feather flock together to identify Coherent Generator Group in the first phase, and utilizes these Coherent Generator Groups to be embodied as constraint; Second stage utilization constraint spectral clustering finds satisfied constraint in pairs, has the off-the-line section of minimum effective power flow impact.When system needs to be split as multiple isolated island, recursive bisection is used to realize.The method relates to generating and impacts two factors with mediation effective power flow, and algorithm reliability has much room for improvement.
In addition based on the parallel network active alone island detection system that M sequence is modulated and variance judges and method.This detection method is the pseudo-random signal M sequence of mechanics of communication, modulates inverter modulating wave according to the regularity of distribution of M sequence, and with the characteristic making output also have M sequence, thus output signal also has the feature of random signal.A variance threshold values is set, when detecting that variance exceedes threshold value and can be judged as isolated island during detection.This algorithm is simple, cost is low, and little to the harmonic effects exporting electric energy, islanding detect blind area is little, but amount of calculation is high, cannot meet the demands completely for large scale system.
Summary of the invention
The object of the invention is for overcoming above-mentioned the deficiencies in the prior art, a kind of intelligent distribution network sub area division based on reachability matrix and island identification method are provided, the method is not by network topology structure variable effect, clear principle, highly sensitive, there is good sub area division ability, the quantity of isolated island, the scope of isolated island and concrete grid-connected switch can be detected fast and accurately.
For achieving the above object, the present invention adopts following technical proposals:
Based on intelligent distribution network sub area division and the island identification method of reachability matrix, comprise the following steps:
The first step: be corresponding node by each switch equivalent in intelligent distribution network, the line equivalent between switch is branch road, sets up a road – node diagram and describes intelligent power distribution web frame;
Second step: set up adjacency matrix A to describe the annexation of each branch road and node in Zhi Lu – node diagram, if certain branch road is directly connected with certain node, then corresponding in matrix A element value is got " 1 ", otherwise gets " 0 ";
3rd step: each the power matirx calculating adjacency matrix A, carries out logic summation operation to corresponding element in each power matirx, obtains reachability matrix P, analyzes the value rule of element in reachability matrix P, obtains the partitioning scenario of network;
4th step: be power supply, circuit or load according to what be connected with switch, definition switching behavior vector K, according to the value condition of element in partitioning scenario and reachability matrix P, obtains the block vector P describing net partition information i, to vectorial K and P iin corresponding element carry out logical AND calculating, obtain the vector describing network each region character, and then judge whether there is isolated island, which power supply corresponding in the quantity of isolated island and each isolated island.
The define method of each element of adjacency matrix A in described second step is:
Reachability matrix computational methods in the 3rd described step are: each the power matirx first obtaining adjacency matrix A: A, A 2, A 3..., A n, wherein n is the nodes of Zhi Lu – node diagram, and then reachability matrix P is obtained in calculating, and computational methods are: P=A+A 2+ A 3+ ... + A n, addition involved is here logic add (∨), i.e. 0 ∨ 0=0,0 ∨ 1=1,1 ∨ 0=1,1 ∨ 1=1.
Switching behavior vector K=[k is defined in the 4th described step 1k 2k n], element k in K ivalue is defined as follows:
If a road – node diagram is divided into some regions, block vector P can be used in each region 1, P 2p mrepresent, will judge whether to comprise in each region the character that power supply and institute comprise power supply further, need with switching behavior vector K to the vectorial P of block irevise:
p i'=k i·p i
Analyze vectorial P i' in element value, just can judge the character of each isolated area fast, thus determine whether there is isolated island, the quantity of isolated island and the scope of isolated island.
The present invention adopts reachability matrix to carry out sub area division and the isolated island identification of intelligent distribution network, and tool has the following advantages compared with existing analogous algorithms:
1. clear principle, information needed is few, is easy to realize;
2. matrix degree of rarefication is high, and amount of calculation is little, is convenient to adopt centralized structure;
3. can define by flexible expansion switchtype as requested, meet the sub area division requirement of different object;
4. can detect the quantity of isolated island, the scope of isolated island and concrete grid-connected switch fast and accurately.
Accompanying drawing explanation
Fig. 1 is the structural representation of intelligent distribution network;
Fig. 2 is node-branch road schematic diagram that intelligent distribution network simplifies;
Embodiment
Below in conjunction with drawings and Examples, the present invention is further described.
1. based on the intelligent distribution network subregion of reachability matrix
In Fig. 1, S represents backbone switch and branch switch, corresponds to node in fig. 2; K is interconnection switch; D represents offside line switching.
When switch S 2 disconnects, adjacency matrix corresponding to the system shown in Fig. 2 is:
A = 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 0 0 0 1 0 1 0 0 0 1 1 0 0 0 1 0 0 0 0 0 1 1 0 0 0 0 0 1 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0
Can calculate each power of adjacency matrix A, and then the reachability matrix P trying to achieve its correspondence is:
P = 1 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 1 0 0 0 0 1 1 0 0 0 1 0 0 0 0 1 1 0 0 0 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1 0 0 1 1 1
The value analyzing element in reachability matrix P is visible, 1st, 6,7 row (column) element values are consistent, and the 3rd, 4,5,8,9,10 row (column) element values are consistent, and the 2nd row (column) element is complete " 0 ", therefore known switch S 2 disconnects, now Zhi Lu – node diagram is divided into 2 relatively independent regions.First region comprises the branch road of switch S 1, S6 and S7 and correspondence, and corresponding element value is " 1 ", shows can reach between S1, S6 and S7; Two Areas comprises switch S 3, the branch road of S4, S5, S8, S9 and S10 and correspondence, and can reach each other.Because S2 disconnects, therefore not any with other node is communicated with.The division result obtained by calculating reachability matrix P and the connected relation of each node, completely the same with the practical structures of the Lu of Zhi shown in Fig. 2 – node diagram.
2. Grid type judges and isolated island identification
After the switch S 2 in Fig. 2 disconnects, electrical network is divided into two relatively independent regions, and corresponding block vector is designated as P1 and P2 respectively, and based on the result of calculation of reachability matrix P, P1 and P2 obtained is respectively:
P 1=[1 0 0 0 0 1 1 0 0 0]
P 2=[0 0 1 1 1 0 0 1 1 1]
According to the switchtype vector K that system configuration and running status obtain be:
K=[1 0 0 0 0 -1 -2 -2 -2 -2]
Utilize vectorial K to carry out correction to vectorial P1 and P2 respectively to obtain:
P 1'=[1 0 0 0 0 -1 -2 0 0 0]
P 2'=[0 0 0 0 0 0 0 -2 -2 -2]
At P 1' in, owing to comprising the element that value is " 1 ", therefore can judge in this region, have switch (i.e. switch S 1) to be connected with system power supply, belong to region of being incorporated into the power networks; At P 2' in, except the element that value is " 0 ", only comprise the element of value for "-2 ", therefore can judge that this region is islet operation region, and be respectively S8 with the switch that DG is directly connected, S9 and S10.

Claims (4)

1., based on intelligent distribution network sub area division and the island identification method of reachability matrix, it is characterized in that, comprise the following steps:
The first step: be corresponding node by each switch equivalent in intelligent distribution network, the line equivalent between switch is branch road, sets up a road – node diagram and describes intelligent power distribution web frame;
Second step: set up adjacency matrix A to describe the annexation of each branch road and node in Zhi Lu – node diagram, if certain branch road is directly connected with certain node, then corresponding in matrix A element value is got " 1 ", otherwise gets " 0 ";
3rd step: each the power matirx calculating adjacency matrix A, carries out logic summation operation to corresponding element in each power matirx, obtains reachability matrix P, analyzes the value rule of element in reachability matrix P, obtains the partitioning scenario of network;
4th step: be power supply, circuit or load according to what be connected with switch, definition switching behavior vector K, according to the value condition of element in partitioning scenario and reachability matrix P, obtains the block vector P describing net partition information i, to vectorial K and P iin corresponding element carry out logical AND calculating, obtain the vector describing network each region character, and then judge whether there is isolated island, the quantity of isolated island and each isolated island which power supply corresponding.
2., as claimed in claim 1 based on intelligent distribution network sub area division and the island identification method of reachability matrix, it is characterized in that, the define method of each element of adjacency matrix A in described second step is:
3., as claimed in claim 1 based on intelligent distribution network sub area division and the island identification method of reachability matrix, it is characterized in that, the reachability matrix computational methods in the 3rd described step are: each the power matirx first obtaining adjacency matrix A: A, A 2, A 3..., A n, wherein n is the nodes of Zhi Lu – node diagram, and then reachability matrix P is obtained in calculating, and computational methods are: P=A+A 2+ A 3+ ... + A n, addition involved is here logic add (∨), i.e. 0 ∨ 0=0,0 ∨ 1=1,1 ∨ 0=1,1 ∨ 1=1.
4. as claimed in claim 1 based on intelligent distribution network sub area division and the island identification method of reachability matrix, it is characterized in that, in the 4th described step, define switching behavior vector K=[k 1k 2k n], element k in K ivalue is defined as follows:
If a road – node diagram is divided into some regions, block vector P can be used in each region 1, P 2p mrepresent, will judge whether to comprise in each region the character that power supply and institute comprise power supply further, need with switching behavior vector K to the vectorial P of block irevise:
p i′=k i·p i
Analyze P ielement value in ' vector, just can judge the character of each isolated area fast, thus determine whether there is isolated island, the quantity of isolated island and the scope of isolated island.
CN201310223708.9A 2013-06-06 2013-06-06 Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix Expired - Fee Related CN103259265B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310223708.9A CN103259265B (en) 2013-06-06 2013-06-06 Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310223708.9A CN103259265B (en) 2013-06-06 2013-06-06 Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix

Publications (2)

Publication Number Publication Date
CN103259265A CN103259265A (en) 2013-08-21
CN103259265B true CN103259265B (en) 2015-02-18

Family

ID=48963039

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310223708.9A Expired - Fee Related CN103259265B (en) 2013-06-06 2013-06-06 Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix

Country Status (1)

Country Link
CN (1) CN103259265B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104124686B (en) * 2014-07-28 2016-04-27 国家电网公司 Based on the network optimization operation method of the blocking theory of distribution
CN106571622A (en) * 2015-10-09 2017-04-19 云南电网有限责任公司 Method for evaluating reliability of distribution network based on traversal algorithm and component hierarchy
CN105470934B (en) * 2015-12-21 2018-07-24 华中科技大学 A kind of power distribution network automatic reclosing control method suitable for distributed generation resource access
CN106199288B (en) * 2016-08-22 2019-01-25 中国电力科学研究院 Active distribution network electric power isolated island recognition methods and device based on network model
CN108090616A (en) * 2017-12-22 2018-05-29 东南大学 A kind of electric system Active Splitting optimal section searching method
CN108536764B (en) * 2018-03-21 2020-12-29 山东昭元信息科技有限公司 Topology and island analysis method and device based on temporary topology analysis table
CN109560546A (en) * 2018-12-10 2019-04-02 四川大学 A kind of quick partition method of urban high voltage distribution network
CN112904239B (en) * 2019-11-19 2023-02-17 中国商用飞机有限责任公司 Cable conduction network determining method and device, storage medium and equipment

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3825172B2 (en) * 1998-04-06 2006-09-20 関西電力株式会社 Power flow calculator for distribution system control system.
CN102611190A (en) * 2012-03-13 2012-07-25 华北电力大学 Combination method for automatic throw-in device of emergency power supply
CN102904334B (en) * 2012-11-05 2015-07-08 江苏省电力公司镇江供电公司 Matrix locking method of spare power automatic switching device
CN103093097A (en) * 2013-01-15 2013-05-08 天津大学 Electrical power system fragile section identification method based on normalized-cut

Also Published As

Publication number Publication date
CN103259265A (en) 2013-08-21

Similar Documents

Publication Publication Date Title
CN103259265B (en) Power network partitioning and detached island identification method of intelligent power distribution network based on reachability matrix
Bolognani et al. Identification of power distribution network topology via voltage correlation analysis
Liao et al. Urban MV and LV distribution grid topology estimation via group lasso
CN109274095B (en) Mutual information-based low-voltage distribution network user topology estimation method and system
Gurusinghe et al. Post-disturbance transient stability status prediction using synchrophasor measurements
Ustun et al. A microgrid protection system with central protection unit and extensive communication
Shams et al. Active power imbalance detection, size and location estimation using limited PMU measurements
CN107104416B (en) The DC line relay protecting method of direct-flow distribution system is pressed during a kind of multiterminal are flexible
CN108022002B (en) Neutral point grounding mode decision method and system based on tensor machine support
CN103138227B (en) Power distribution network fast power restoration method containing distributed power connected grid
CN102521478B (en) Power network topology dynamic coloring algorithm
Simonov Dynamic partitioning of DC microgrid in resilient clusters using event-driven approach
Rahnamay-Naeini et al. On the role of power-grid and communication-system interdependencies on cascading failures
Gutierrez et al. Vulnerability analysis of power grids using modified centrality measures
CN103093097A (en) Electrical power system fragile section identification method based on normalized-cut
CN104901308A (en) Power system critical circuit identification method
Duong et al. Online voltage stability monitoring based on PMU measurements and system topology
Garau et al. ICT reliability modelling in co-simulation of smart distribution networks
CN105046011A (en) Rapid analysis method used for state of electric device and based on topological computation of distributed power grid
Sayeekumar et al. Graph theory and its applications in power systems-a review
Deka et al. Structure learning and statistical estimation in distribution networks-Part II
CN104578057B (en) A kind of consider source flow path containing wind-powered electricity generation line voltage Pilot bus recognition methods
CN108964282B (en) Method and device for deploying power quality monitoring device and computer readable medium
Anderson et al. A minimum cover algorithm for PMU placement in power system networks under line observability constraints
CN105356457A (en) Standby electrical network accident recovery capacity adequacy evaluation 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
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150218

Termination date: 20200606

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