Specific embodiment
As shown in figure 1, a kind of radial topological constraints determination methods of many electrical power distribution nets, including step:
S100:Obtain the power supply point in many electrical power distribution nets, the connection between branch road, node and branch road and node to close
System, and the adjacency matrix for representing the direct connection relational between each node is formed, adjacency matrix includes not considering on off state
Element adjacency matrix and the electric energy adjacency matrix for considering on off state.
All power point in many electrical power distribution nets, all branches and whole nodes are obtained, wherein, power supply point is with spy
The node of different mark, in addition, also obtaining the annexation between branch road and node, then forms expression based on these data for obtaining
The adjacency matrix of the direct connection relational between each node, adjacency matrix includes not considering the element adjacency matrix of on off state
With the electric energy adjacency matrix for considering on off state.Non-essential, recording element adjacency matrix is A0, recording electric energy adjacency matrix is
A.More specifically, according to element situation between each node in power distribution network, N is formedb×NbRank (NbFor nodes) element adjacency matrix
A0, represented with element 1 and there is element between two nodes.In element adjacency matrix A0On the basis of, by the switch in separated position
Corresponding element is entered as 0, and remaining element keeps constant, obtains electric energy adjacency matrix A.
S200:Interconnection switch number is calculated according to element adjacency matrix and electric energy adjacency matrix.
The element adjacency matrix obtained by step S100 is calculated with electric energy adjacency matrix, obtains interconnection switch number.It is non-
It is necessary, interconnection switch number can be calculated according to the difference element of element adjacency matrix and electric energy adjacency matrix.Specially:Contact
Switch number Nts=(1/2) * matrix (A0- A) in all elements sum.
S300:Judge whether interconnection switch number meets judgement equation, judge that equation subtracts for interconnection switch number equal to circuitry number
Go difference between nodes and power supply number.
Based on interconnection switch number, circuitry number, nodes and power supply number that step S100 and step S200 are obtained, connection is judged
Whether network switch number meets judgement equation, wherein, equation is judged for interconnection switch number=circuitry number-(nodes-power supply number).
S400:If being unsatisfactory for, judge that distribution net work structure is unsatisfactory for radial topological constraints.
When interconnection switch number is unsatisfactory for judging equation, judge that distribution net work structure is unsatisfactory for radial topological constraints.
The radial topological constraints determination methods of many electrical power distribution nets of the present invention, obtain power supply point in many electrical power distribution nets,
Annexation between branch road, node and branch road and node, and element adjacency matrix and electric energy adjacency matrix are formed, calculate connection
Network switchs number, when interconnection switch number is unsatisfactory for judging equation, judges that distribution net work structure is unsatisfactory for radial topological constraints, its
In, equation is judged for interconnection switch number=circuitry number-(nodes-power supply number).In whole process, branch road in many electrical power distribution nets
Number, nodes and power supply number are judged, without the need for being reconstructed by switching division and the backtracking of node power point, therefore nothing
Need the process of big data quantity and calculate, can efficiently to being unsatisfactory for the judgement of radial topological constraints in many electrical power distribution nets.
As shown in Fig. 2 also including after step S300 in one embodiment wherein:
S520:When interconnection switch number meets judges equation, from each power supply point, according to electric energy adjacency matrix to away from
The branch road in power supply point direction carries out node searching.
After step S300 judges, when interconnection switch number meets judges equation, in addition it is also necessary to further located
Reason, now from each power supply point, node searching is carried out according to electric energy adjacency matrix to the branch road away from power supply point direction.It is non-must
Want, from each power supply point, directions of the Zhi Luxiang represented by each nonzero element in electric energy adjacency matrix away from power supply point
Search, and record power supply point source, you can the connected relation between each node and each power supply point is obtained, power supply point record is built
Matrix.
S540:Record it is searched to node and the searched node for arriving power supply point situation.
From each power supply point, node searching is carried out to the branch road away from power supply point direction according to electric energy adjacency matrix
During, by it is searched to node and the power supply point situation of the node record.It is non-essential, structure can be adopted
Build power supply point matrix-style to record the power supply point situation of the searched node for arriving and the node.
S560:Judge whether it is searched to node and each power supply point between there is power supply relation and search each section
The power supply point number of point is 1.
Can be based on step S540 obtain it is searched to node and the searched node for arriving power supply point feelings
Condition, arranges relation of powering between egress and node corresponding power point, further judges whether the node that searches and each electricity
Whether there is power supply relation between source point.To find that power supply has the respective element note of power supply relation to node in search procedure
Be that 0 expression does not search power supply relation for 1, it is logical go over search by record all nodes in the power distribution network for searching
Power supply point situation.
S580:If so, judge that distribution net work structure meets radial topological constraints.
When it is searched to node and each power supply point between there is power supply relation when, judge distribution net work structure meet radiation
Shape topological constraints.Non-essential, in actual applications, step 580 is specially:After judging all branch roads of traversal, each node is given
Whether the power supply points of power supply are 1, if so, then provide the conclusion of " distribution net work structure meets radial ".In addition, when each section
When the power supply points of point power supply are not to be 1, the conclusion of " distribution net work structure is unsatisfactory for radial " is given.
As shown in Fig. 2 including before step S520 in an example wherein:
S510:Power supply point source is recorded, power supply point record matrix is built.
Before step S520 is carried out, from each power supply point, represented by each nonzero element in electric energy adjacency matrix A
Zhi Luxiang is searched for away from the direction of power supply point, and records power supply point source, you can obtained between each node and each power supply point
Connected relation, is denoted as power supply point record matrix S.Specifically, power supply point record matrix S is K × NbRank matrix, K is power supply point,
NbFor node.Further, power supply point record matrix can also be initialized, K power supply point in establishing electrical network
Node serial number is respectively m1, m2, m3..., mK, then element s in SijFor 1 when represent i-th power supply point, i.e. miIndividual node is to jth
Individual node is powered;It is not difficult to find out, when S is initialized, should be by si, mi(i=1,2,3 ..., K) is set to 1, represents each power supply point extremely
Less by power itself.
Wherein in one embodiment, from each power supply point, according to electric energy adjacency matrix to away from power supply point direction
The step of branch road carries out node searching includes:
Step one:Remaining adjacency matrix in record search procedure, searches the node composition searched in search procedure
Rope knot vector, wherein, remaining adjacency matrix is electric energy adjacency matrix when initial, and search node vector includes each power supply point.
Remaining adjacency matrix in search procedure is designated as AS, A when initialS=A;The section searched in search procedure
Point composition of vector is designated as d, and when initial, d is only made up of power supply point node.
Step 2:Each node from the vector of search node, to away from the remaining adjacent square of power supply point direction search
Remaining branch road in battle array.
Each node from d, to away from power supply point direction A is searched forSIn remaining branch road.
Record it is searched to node and it is searched to node power supply point situation the step of include:
Step one:Power supply point of the node for searching in power supply point is remembered relevant position in power supply point record matrix
Under record.
Step 2:The branch road searched for is deleted from the remaining adjacency matrix of search, to update remaining adjacency matrix, will
It is searched to node be added to search node vector in.
Power supply point of the node for searching in power supply point is recorded relevant position in matrix S, and searching for
The branch road crossed is from matrix ASMiddle deletion, by it is searched to node be added in d.
In order to further explain in detail the radial topological constraints determination methods of many electrical power distribution nets of the invention and system, under
Face will be using being described in more details in a specific embodiment.In this specific embodiment, many electrical power distribution nets of the invention
Radial topological constraints determination methods include step:
1st, power supply point, the annexation between branch road, node and branch road and node in many electrical power distribution nets, shape are obtained
Into element adjacency matrix A0With electric energy adjacency matrix A.According to element situation between each node in power distribution network, N is formedb×NbRank (NbFor
Nodes) element adjacency matrix A0, represented with element 1 and there is element between two nodes.In element adjacency matrix A0On the basis of,
The corresponding element of switch in separated position is entered as into 0, remaining element keeps constant, obtains electric energy adjacency matrix A.
2nd, by element adjacency matrix A0Interconnection switch number N is calculated with the difference element of electric energy adjacency matrix Ats;By matrix
(A0- A) in all elements summation, then divided by 2, you can obtain Nts。
3rd, judge whether " interconnection switch number=circuitry number-(nodes-power supply number) " are set up, if not, be then given
The conclusion of " distribution net work structure is unsatisfactory for radial ", and terminate radial distribution networks and judge flow process;If so, then continue.
4th, power supply point record matrix S (K × N are initializedbRank).The node serial number difference of K power supply point in establishing electrical network
For m1, m2, m3..., mK, then element s in SijFor 1 when represent i-th power supply point, i.e. miIndividual node is powered to j-th node;
It is not difficult to find out, when S is initialized, should be by si, mi(i=1,2,3 ..., K) is set to 1, represents that each power supply point is at least supplied by itself
Electricity.
5th, the remaining adjacency matrix in search procedure is designated as AS, initialize AS=A;What is searched in search procedure
Node composition of vector is designated as d, and when initial, d is only made up of power supply point node.
6th, each node from d, to away from power supply point direction A is searched forSIn remaining branch road, the node for searching
Record relevant position in matrix S in the power supply point of power supply point, and the branch road searched for from matrix ASMiddle deletion,
By it is searched to node be added in d.
If the 7, matrix ASElement be all zero or unchanged, then stop search;Otherwise repeat step 6.
8th, after judging all branch roads of traversal, 1 whether is to the power supply points that each node is powered, is if so, then given and " matches somebody with somebody
Electric network composition meets radial " conclusion;Otherwise, the conclusion of " distribution net work structure is unsatisfactory for radial ", program determination are given.
As shown in figure 3, a kind of radial topological constraints of many electrical power distribution nets judge system, including:
Analysis module 100, for obtaining many electrical power distribution nets in power supply point, branch road, node and branch road and node it
Between annexation, and form the adjacency matrix for representing the direct connection relational between each node, adjacency matrix includes not considering
The element adjacency matrix of on off state and the electric energy adjacency matrix for considering on off state.
Computing module 200, for calculating interconnection switch number according to element adjacency matrix and electric energy adjacency matrix.
Judge module 300, for judging whether interconnection switch number meets judgement equation, judges that equation is interconnection switch number etc.
Difference between nodes and power supply number is deducted in circuitry number.
First determination module 400, for when interconnection switch number is unsatisfactory for judging equation, judging that distribution net work structure is unsatisfactory for
Radial topological constraints.
The radial topological constraints of many electrical power distribution nets of the present invention judge system, and analysis module 100 obtains many electrical power distribution nets
In power supply point, the annexation between branch road, node and branch road and node, and it is adjacent to form element adjacency matrix and electric energy
Matrix, computing module 200 calculates interconnection switch number, and judge module 300 judges whether interconnection switch number meets judgement equation, its
In, equation is judged for interconnection switch number=circuitry number-(nodes-power supply number), when interconnection switch number is unsatisfactory for judging equation
When, the first determination module 400 judges that distribution net work structure is unsatisfactory for radial topological constraints,.In whole process, many electrical power distribution nets
Middle circuitry number, nodes and power supply number are judged, without the need for being reconstructed by switching division and the backtracking of node power point,
Therefore without the need for the process and calculating of big data quantity, efficiently can sentence to being unsatisfactory for radial topological constraints in many electrical power distribution nets
It is fixed.
Wherein in one embodiment, difference unit of the computing module 200 according to element adjacency matrix and electric energy adjacency matrix
Element calculates interconnection switch number.
As shown in figure 4, wherein in one embodiment, the radial topological constraints of many electrical power distribution nets judge that system is also wrapped
Include:
Search module 520, it is adjacent according to electric energy from each power supply point for when interconnection switch number satisfaction judges equation
Connect matrix carries out node searching to the branch road away from power supply point direction.
Logging modle 540, for record it is searched to node and the searched node for arriving power supply point situation.
Power supply relationship module 560, for judge whether it is searched to node and each power supply point between be exist power supply close
The power supply point number for being and searching each node is 1.
There is power supply relation and search each in the second determination module 580, the power supply point for powering when each node
When the power supply point number of node is 1, judge that distribution net work structure meets radial topological constraints.
As shown in figure 4, wherein in one embodiment, the radial topological constraints of many electrical power distribution nets judge that system is also wrapped
Include:
Power supply point records matrix module 510, for recording power supply point source, builds power supply point record matrix.
Wherein in one embodiment, search module 520 includes the first search unit and the second search unit, logging modle
540 include recording unit and updating block;
First search unit, for recording search procedure in remaining adjacency matrix, the section searched in search procedure
Search node is vectorial for point composition, wherein, remaining adjacency matrix is electric energy adjacency matrix when initial, and search node vector includes
Each power supply point.
Second search unit, for each node from the vector of search node, to away from the search of power supply point direction
Remaining branch road in remaining adjacency matrix.
Recording unit, for by the node for searching power supply point power supply point power supply point record matrix in phase
Position is answered to record.
Updating block, it is remaining adjacent to update for the branch road searched for be deleted from the remaining adjacency matrix of search
Matrix, by it is searched to node be added to search node vector in.
In order to further explain in detail the radial topological constraints determination methods of many electrical power distribution nets of the invention with system
Technical scheme and its effect brought, will be described in detail below using instantiation and with reference to Fig. 5 to Figure 10.
For the distribution network topology embodiment shown in Fig. 5, " Bus " represents node, and " Line " represents branch road, and solid line is normal
The block switch branch road for closing, dotted line is the normally opened interconnection for being furnished with interconnection switch, and Bus 1 and Bus 10 are power supply nodes, feedback
Line F1 is connected by contact branch road 6 with F2 with branch road 7, there is contact branch road 5 inside feeder line F1, there is interconnection 9 inside feeder line F2;
Knowable to known conditions, node 1 and 10 be power supply point, power supply number K=2, nodes Nb=10, circuitry number Nbr=
12, element adjacency matrix is:
With reference to A0And the state of each branch switch, electric energy adjacency matrix can be obtained:
So as to,
By (A0- A) N can be obtainedts=4;Nts=Nbr-(Nb- K) set up, continue;
Source node vector d=[1,10] has been searched in initializationT(element number is designated as N in dd, d and NdIn search procedure not
It is disconnected to update), the remaining adjacency matrix A of initializationS=A, and initialize power supply point record matrix:
Start the power conditions for searching for each node from power supply point, until powered branch lives through search.
By taking network structure shown in accompanying drawing 5 as an example, search procedure is:
Original state d=[1,10]T,
→ d=[1,2,9,10]T,
→ d=[1,2,3,5,6,7,9,10]T,
→ d=[1,2,3,4,5,6,7,8,9,10]T,
So far, flag=0, AS=0, search terminates.
It can be seen that, the power supply points of each node are 1, and the Distributing network structure meets radial topological constraints.
It is exactly above the whole process of the tree topology constraint qualification of Fig. 5 embodiments.
In order to embody the applicability of the present invention, radial topology is carried out about to the different running method of power distribution network shown in Fig. 5
Beam judgement (under the various methods of operation, element adjacency matrix A0Keep constant, K=2, Nb=10, Nbr=12 is constant, and electric energy
Adjacency matrix A is different):
Fig. 6 is on the basis of Fig. 5, the switch of branch road 5 to be closed the corresponding method of operation;Fig. 7 is on the basis of Fig. 5
On, the switch of branch road 6 is closed the corresponding method of operation;Nts=3, Nts≠Nbr-(Nb- K), can immediately judge both operations
Mode is unsatisfactory for radial topological constraints;This is mutual with power supply node 1 and 10 in branch road 2-3-5-4 cyclizations, Fig. 7 in Fig. 6
Reachable intuitive judgment is also what is be consistent;
Fig. 8 is on the basis of Fig. 1, the switch of branch road 6 is closed, the switch of branch road 4 separates the corresponding method of operation;Nts
=4, Nts=Nbr-(Nb- K), continue executing with judgement flow process;Until ASWhen=0, have:
D=[1,2,3,4,5,6,7,8,9,10]T,
It can be seen that, the power supply points of each node are 1, Fig. 4 Distributing network structures and meet radial topological constraints;
Fig. 9 is on the basis of Fig. 5, the switch of branch road 6 is closed, the switch of branch road 2 separates the corresponding method of operation;Nts
=4, Nts=Nbr-(Nb- K), continue executing with judgement flow process;Until ASWhen no longer changing, have:
D=[1,2,5,6,7,8,9,10]T,
It is 0 to the power supply number of node 3 and the power supply of node 4, and is 2, Fig. 9 distributions to the power supply number that node 5 and node 6 are powered
Structure is unsatisfactory for radial topological constraints;This isolated with Fig. 9 interior joints 3 and 4 and power supply node 1 and 10 is mutually reachable intuitively to be sentenced
Disconnected is also what is be consistent;
Figure 10 is on the basis of Fig. 5, the switch of branch road 5 and 6 is closed, the switch of branch road 2 and 4 separates corresponding operation
Mode;Nts=4, Nts=Nbr-(Nb- K), continue executing with judgement flow process;Until ASWhen=0, have:
D=[1,2,3,4,5,6,7,8,9,10]T,
It can be seen that, the power supply points of each node are 1, and Distributing network structure meets radial topological constraints.
Embodiment described above only expresses the several embodiments of the present invention, and its description is more concrete and detailed, but and
Can not therefore be construed as limiting the scope of the patent.It should be pointed out that for one of ordinary skill in the art comes
Say, without departing from the inventive concept of the premise, some deformations and improvement can also be made, these belong to the protection of the present invention
Scope.Therefore, the protection domain of patent of the present invention should be defined by claims.