CN110222966B - Distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method - Google Patents
Distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method Download PDFInfo
- Publication number
- CN110222966B CN110222966B CN201910453337.0A CN201910453337A CN110222966B CN 110222966 B CN110222966 B CN 110222966B CN 201910453337 A CN201910453337 A CN 201910453337A CN 110222966 B CN110222966 B CN 110222966B
- Authority
- CN
- China
- Prior art keywords
- node
- sub
- nodes
- distribution network
- state estimation
- 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.)
- Active
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 35
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 34
- 238000000638 solvent extraction Methods 0.000 title claims abstract description 33
- 238000005192 partition Methods 0.000 claims abstract description 62
- 239000011159 matrix material Substances 0.000 claims abstract description 39
- 230000002068 genetic effect Effects 0.000 claims abstract description 13
- 210000000349 chromosome Anatomy 0.000 claims description 9
- 230000035772 mutation Effects 0.000 claims description 3
- 108090000623 proteins and genes Proteins 0.000 claims description 3
- 238000004891 communication Methods 0.000 abstract description 6
- 238000004364 calculation method Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 238000010223 real-time analysis Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/067—Enterprise or organisation modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/22—Flexible AC transmission systems [FACTS] or power factor or reactive power compensating or correcting units
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Power Engineering (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
A distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method comprises the following steps: for a selected power distribution system, acquiring a topological connection relation of a power distribution network, constructing an adjacent matrix A, and setting the number M of sub-regions; establishing a synchronous phasor measurement configuration partition model facing to the distributed state estimation of the power distribution network; solving a synchronous phasor measurement configuration partition model facing the distribution network distributed state estimation through a genetic algorithm; and outputting the obtained synchronous phasor measurement configuration partition scheme oriented to the distribution network distributed state estimation. The invention ensures connectivity of the subareas by using a method of selecting the central node first and then gradually expanding outwards, takes the minimum difference of the number of the nodes among the subareas as a target function and solves the problem through a genetic algorithm, thereby ensuring the minimum difference of the number of the nodes among the subareas and ensuring the similar state estimation time of different areas, thereby reducing the waiting time of adjacent area communication and ensuring the real-time property of distributed state estimation.
Description
Technical Field
The invention relates to a power distribution network partitioning method. In particular to a distribution network distributed state estimation oriented synchronous phasor measurement configuration partitioning method.
Background
With the enlargement of the scale of a power distribution system, a large number of distributed power sources are accessed and users participate widely, the communication burden of the traditional centralized state estimation is gradually increased, the dimension of a state variable is high, the overall solving time is long, and the requirement of real-time analysis and control of the power distribution network is difficult to meet. In order to reduce communication burden and increase computation speed, distributed state estimation is becoming an effective solution to the problem. The premise of state estimation convergence is that the system is considerable, and the existing power distribution network cannot meet the integral observability due to the fact that the data acquisition points are multiple, the distribution is wide, the coverage of monitoring points is not comprehensive, and the configuration is unreasonable. In addition, most of the measurement used for state estimation comes from a data acquisition and monitoring control System (SCADA) and an advanced measurement system (AMI), and not only is the data itself low in precision, poor in synchronism and long in acquisition period, but also the measurement and state variables are in a nonlinear relation, and iterative solution needs to be performed through algorithms such as a weighted least square method, and even if distributed state estimation is adopted, the accuracy and the real-time performance of system state solution are difficult to guarantee.
Due to the introduction of the synchronous phasor measurement device, the operation monitoring level of the power distribution network is greatly improved. Compared with the traditional measuring device, the synchronous phasor measuring device can synchronously acquire the measurement information of the node voltage phasor, the branch current phasor and the system frequency, is high in uploading speed, effectively solves the problems of poor measurement data quality, low synchronism and long acquisition period of the traditional power distribution system, improves the calculation speed and accuracy of applications such as model parameter verification, state estimation, system protection and operation control, and is an important link for the technical development of the intelligent power distribution network. Particularly, in the application of state estimation, the voltage and current phasor measurement acquired by the synchronous phasor measurement device is in a linear relation with the state variable of the system, so that the calculation time is greatly reduced, and the real-time analysis and control of the power distribution system are guaranteed.
If all the measurement data adopted in the distributed state estimation come from the synchronous phasor measurement device, the rapidity of state variable solving can be ensured, the communication burden can be reduced, and the accuracy and the instantaneity of the estimation result can be ensured due to strong data synchronism and high uploading frequency. However, before configuring the synchronous phasor measurement apparatus, the partitioning of the power distribution network needs to be completed. The existing power distribution network distributed state estimation partitioning method mainly comprises the following steps: and manually completing the partition of the power distribution network based on the topological structure and the geographical position partition, the existing measurement partition and the distributed state estimation partition principle. Partitioning is difficult to ensure that the number of nodes in each partition is similar based on a topological structure and a geographic position; the method is based on the fact that the partition cannot be completed when the measurement configuration of the power distribution network is insufficient or configuration information cannot be obtained by the existing measurement partition, and the partition effect is related to the setting of relevant parameters, so that the method is not suitable for the situation that synchronous phasor measurement is configured after the partition is performed; the optimal partitioning scheme is difficult to obtain after the power distribution network partitioning is completed manually. Therefore, finding an optimization method which is not based on existing measurement and meets the principle of distributed state estimation to complete power distribution network partitioning, enabling the number of nodes between each region to be similar and enabling adjacent regions to have only one overlapped node is a prerequisite condition for configuring a synchronous phasor measurement device to enable each sub-network to be completely observable so as to achieve distributed state estimation.
Disclosure of Invention
The invention aims to solve the technical problem of providing a synchronous phasor measurement configuration partitioning method for distributed state estimation of a power distribution network, which can reduce the waiting time of adjacent area communication and ensure the real-time property of the distributed state estimation.
The technical scheme adopted by the invention is as follows: a distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method comprises the following steps:
1) For a selected power distribution system, acquiring a topological connection relation of a power distribution network, constructing an adjacent matrix A, and setting the number M of sub-regions;
2) Establishing a synchronous phasor measurement configuration partition model facing to the distributed state estimation of the power distribution network;
3) Solving the synchronous phasor measurement configuration partition model facing the power distribution network distributed state estimation in the step 2) through a genetic algorithm;
4) And outputting the obtained synchronous phasor measurement configuration partition scheme oriented to the distribution network distributed state estimation.
The number M of the subregions in the step 1) is as follows:
wherein N represents the total node number of the power distribution system to be partitioned,representing number->And rounding down.
The distribution network distributed state estimation-oriented synchronous phasor measurement configuration partition model in the step 2) takes the minimum difference of the number of nodes among the sub-partitions as a target function, and the mathematical expression is as follows:
wherein N is i Indicates the number of nodes included in the sub-region i.
The step 3) comprises the following steps:
(3.1) setting the length of the chromosome of an individual in a genetic algorithm to be equal to the number M of sub-regions, setting a value set of genes on the chromosome to be a set omega formed by all nodes of a power distribution system, wherein omega = {1, 2.., N }, N represents the total number of the nodes of the power distribution system to be partitioned, randomly generating an initial population, setting an evolutionary algebra d =1, and setting a cross probability p C Genetic probability p M And a maximum evolution algebra D;
(3.2) obtaining an initial center node through chromosomes of individuals in a group, partitioning the power distribution network by using a center expansion partitioning method, and further calculating the fitness of the individuals;
(3.3) acting a selection operator on the population, and reserving the individual with the highest fitness;
(3.4) acting the crossover operator and the mutation operator on the selected population to generate a next generation, wherein the evolution algebra d = d +1;
(3.5) if the evolution algebra D is less than D, entering the step (3.2), otherwise, entering the step (3.6);
and (3.6) outputting the individual with the highest fitness as an optimal solution.
The method for center-expanded partition in step (3.2) comprises:
(3.2.1) setting all the nodes not to be partitioned, randomly acquiring M nodes as initial central nodes, and marking the initial central nodes as partitioned nodes;
(3.2.2) the number of expansion times k =0, updating the sub-region set corresponding to the initial central node according to the initial central node and constructing a partition matrix E k The number of expansion times k = k +1; wherein the partition matrix E k Is a matrix of dimension M × N, for l and i respectively represent a node number and a sub-region number, and if the node l is in the sub-region i, E k Ith row and ith column element E k,i,l =1, otherwise E k,i,l Is 0, set of subregions Γ l The middle element is the number of the subregion to which the node l belongs, if the node l is in the subregion i, i belongs to the gamma l ;
(3.2.3) setting a partition matrix E 'to be adjusted after the k expansion' k =E k-1 A and A are adjacent matrixes, and a partition matrix E 'to be adjusted' k All the non-zero elements in the sub-area are set to be 1, i =1, j =2, and i and j are the numbers of the sub-areas;
(3.2.4) if j > M, i = i +1, j = i +1, and entering the (3.2.5) th step, and otherwise, entering the (3.2.6) th step;
(3.2.5) if i = M, entering the (3.2.8) step, otherwise, entering the (3.2.6) step;
(3.2.6) obtaining a set Λ of overlapping nodes of the sub-region i and the sub-region j i,j If the node set is overlapped i,j If the current is an empty set, j = j +1, and the step (3.2.4) is carried out, otherwise, the step (3.2.7) is carried out;
(3.2.7) if overlapping node set Λ i,j In the method, only two nodes are provided, and the two nodes are partitioned nodes, the two nodes are respectively used as the shared boundary nodes of the sub-regions i and j, and the comparison is maximumSelecting a boundary node by taking the minimum difference as a target according to the difference between the number of the area nodes and the minimum number of the area nodes, and adjusting a to-be-adjusted partition matrix E 'by taking a node in the pre-expansion sub-area i as a boundary if the differences are the same' k Updating a sub-region set corresponding to the boundary node, wherein j = j +1, and entering the (3.2.4) th step;
(3.2.8) set the set of all partitioned nodes to Ω 1 To aIf i ∈ Γ l Then the partition matrix E 'is to be adjusted' k Line i of l element E' k,i,l =1, otherwise E' k,i,l =0;
(3.2.9) set the set of all non-partitioned overlapping nodes to Ω 2 To forIf the overlapped subareas have boundaries, dividing the node l into the subarea with the minimum node number, if the subarea with the minimum node number has a plurality of subareas, dividing the node l into the subarea with the minimum serial number, and adjusting the subarea matrix E 'to be adjusted' k Updating the sub-region set corresponding to the node l;
(3.2.10) assuming that i =1, j =2, i and j are both sub-region numbers;
(3.2.11) if j > M, i = i +1, j = i +1, and entering the (3.2.12) th step, and otherwise entering the (3.2.13) th step;
(3.2.12) if i = M, entering the (3.2.15) step, otherwise, entering the (3.2.13) step;
(3.2.13) obtaining the set Λ of overlapping nodes of the sub-regions i and j i,j If the node set is overlapped i,j If the current is an empty set, j = j +1, and the step (3.2.11) is entered, otherwise, the step (3.2.14) is entered;
(3.2.14) if overlapping node set Λ i,j The node is not partitioned, and the sub-region i and the sub-region j are overlapped for the first time, the node is used as a boundary node shared by the sub-region i and the sub-region j, and a partition matrix E 'to be adjusted is adjusted' k Updating the pairs of edge nodesThe corresponding subregion set, j = j +1, enters the step (3.2.11);
(3.2.15) setting the adjusted partition matrix E k =E′ k If E is k =E k-1 Entering the step (3.2.16), otherwise, E k The node contained in the node list is marked as a partitioned node, the expansion times k = k +1, and the step (3.2.3) is carried out;
(3.2.16) according to partition matrix E k And outputting a partition result.
The synchronous phasor measurement configuration partitioning method for the distributed state estimation of the power distribution network ensures connectivity of sub-regions by using a method of selecting a central node first and then gradually expanding outwards, takes the minimum difference of the number of nodes among the sub-regions as a target function and solves the target function through a genetic algorithm, ensures the minimum difference of the number of nodes among the regions, ensures that state estimation time of different regions is close, reduces the waiting time of communication of adjacent regions and ensures the real-time property of the distributed state estimation.
Drawings
FIG. 1 is a flow chart of a method for partitioning a synchronous phasor measurement configuration for distributed state estimation of a power distribution network according to the present invention;
FIG. 2 is an IEEE 33 node calculation diagram;
FIG. 3 is a graph of the PG & E69 node calculation;
FIG. 4 is an IEEE 33 node example partitioning result;
FIG. 5 is the PG & E69 node example partitioning results.
Detailed Description
The following describes the partition method for the measurement configuration of the synchronous phasor for the distributed state estimation of the power distribution network in detail with reference to the embodiments and the drawings.
As shown in fig. 1, the method for partitioning a synchronized phasor measurement configuration for distributed state estimation of a power distribution network according to the present invention includes the following steps:
1) For a selected power distribution system, acquiring a topological connection relation of a power distribution network, constructing an adjacent matrix A, and setting the number M of sub-regions, wherein the number M of the sub-regions is as follows:
wherein N represents the total node number of the power distribution system to be partitioned,means number->And rounding down.
2) Establishing a synchronous phasor measurement configuration partition model facing the power distribution network distributed state estimation, wherein the synchronous phasor measurement configuration partition model facing the power distribution network distributed state estimation takes the minimum node number difference among sub-regions as a target function, and the mathematical expression is as follows:
wherein N is i Indicates the number of nodes included in the sub-region i.
3) Solving the synchronous phasor measurement configuration partition model facing the power distribution network distributed state estimation in the step 2) through a genetic algorithm; the method comprises the following steps:
(3.1) setting the chromosome length of an individual in a genetic algorithm to be equal to the number M of subregions, setting a value set of genes on the chromosome to be a set omega formed by all nodes of a power distribution system, wherein omega = {1, 2.., N }, and N represents the total number of the nodes of the power distribution system to be partitioned, randomly generating an initial population, setting an evolution algebra d =1, and setting a cross probability p C Genetic probability p M And a maximum evolution algebra D;
(3.2) obtaining an initial center node through chromosomes of individuals in a group, partitioning the power distribution network by using a center expansion partitioning method, and further calculating the fitness of the individuals; the center expansion partition method comprises the following steps:
(3.2.1) setting all the nodes not to be partitioned, randomly acquiring M nodes as initial central nodes, and marking the initial central nodes as partitioned nodes;
(3.2.2) the number of expansion times k =0, updating the sub-region set corresponding to the initial central node according to the initial central node and constructing a partition matrix E k The number of expansion times k = k +1; wherein the partition matrix E k Is a matrix of dimension M × N, for l and i respectively represent a node number and a sub-region number, and if the node l is in the sub-region i, E k Row i and column l element E of k,i,l =1, otherwise E k,i,l Is 0, the set of sub-regions Γ l The middle element is the number of the subregion to which the node l belongs, if the node l is in the subregion i, i belongs to the gamma l ;
(3.2.3) setting a partition matrix E 'to be adjusted after the k expansion' k =E k-1 A and A are adjacent matrixes, and a partition matrix E 'to be adjusted' k All the non-zero elements in the list are set to be 1, i =1, j =2, and i and j are all subarea numbers;
(3.2.4) if j > M, i = i +1, j = i +1, and the step (3.2.5) is entered, otherwise, the step (3.2.6) is entered;
(3.2.5) if i = M, entering the (3.2.8) th step, otherwise, entering the (3.2.6) th step;
(3.2.6) obtaining the overlapping node set Lambda of the sub-regions i and j i,j If the node set is overlapped i,j If the current is an empty set, j = j +1, and the step (3.2.4) is carried out, otherwise, the step (3.2.7) is carried out;
(3.2.7) if overlapping node set Λ i,j If there are only two nodes in the partition matrix E ' and the two nodes are all partitioned nodes, the two nodes are respectively used as shared boundary nodes of the sub-area i and the sub-area j, the difference between the maximum area node number and the minimum area node number is compared, the boundary node is selected by taking the minimum difference as a target, and if the differences are the same, the node in the sub-area i before expansion is taken as a boundary, and the partition matrix E ' to be adjusted is adjusted ' k Updating a sub-region set corresponding to the boundary node, wherein j = j +1, and entering the (3.2.4) step;
(3.2.8) set the set of all partitioned nodes to Ω 1 To aIf i ∈ Γ l Then the partition matrix E 'is to be adjusted' k Line ith element E' k,i,l =1, otherwise E' k,i,l =0;
(3.2.9) set the set of all non-partitioned overlapping nodes to Ω 2 To forIf the overlapped subareas have boundaries, dividing the node l into the subarea with the minimum node number, if the subarea with the minimum node number has a plurality of subareas, dividing the node l into the subarea with the minimum serial number, and adjusting the subarea matrix E 'to be adjusted' k Updating the sub-region set corresponding to the node l;
(3.2.10) assuming that i =1, j =2, i and j are both sub-region numbers;
(3.2.11) if j > M, i = i +1, j = i +1, and entering the (3.2.12) th step, and otherwise entering the (3.2.13) th step;
(3.2.12) if i = M, entering the (3.2.15) step, otherwise, entering the (3.2.13) step;
(3.2.13) obtaining a set Λ of overlapping nodes of the sub-region i and the sub-region j i,j If the node sets are overlapped, Λ i,j If the current is an empty set, j = j +1, and the step (3.2.11) is entered, otherwise, the step (3.2.14) is entered;
(3.2.14) if overlapping node set Λ i,j The node is not partitioned, and the sub-area i and the sub-area j are overlapped for the first time, then the node is used as a boundary node shared by the sub-area i and the sub-area j, and a partition matrix E 'to be adjusted is adjusted' k Updating a subregion set corresponding to the boundary node, wherein j = j +1, and entering the step (3.2.11);
(3.2.15) setting the adjusted partition matrix E k =E′ k If E is k =E k-1 Then step (3.2.16) is entered, otherwise E is entered k Node marking contained inFor the partitioned node, the expansion times k = k +1, and the step (3.2.3) is entered;
(3.2.16) according to partition matrix E k Outputting a partition result;
(3.3) acting a selection operator on the population, and reserving the individuals with the highest fitness;
(3.4) acting the crossover operator and the mutation operator on the selected population to generate a next generation, wherein the evolution algebra d = d +1;
(3.5) if the evolution algebra D is less than D, entering the step (3.2), otherwise, entering the step (3.6);
and (3.6) outputting the individual with the highest fitness as an optimal solution.
Specific examples are given below:
employing IEEE 33 node algorithm and PG&The method provided by the invention is verified by an example of the E69 node, and the example topological connection relation is shown in FIGS. 2 and 3. Let the cross probability p C =0.3, genetic probability p M =0.1, the number of initial population individuals is 5000, and the maximum evolution algebra D =100, after the partitioning process according to the present invention is performed, the partitioning result of the IEEE 33 node example is shown in fig. 4, the specific partitioning scheme is shown in table 1, and PG is&The partitioning result of the E69 node algorithm is shown in fig. 5, and the specific partitioning scheme is shown in table 2. Therefore, the partitioning result obtained by solving by using the method provided by the invention can ensure that the number of nodes between the sub-regions has smaller difference, and the validity of the method is verified.
TABLE 1 IEEE 33 node example partition scheme obtained by solving by the method of the invention
Region numbering | Center node numbering | Node numbering within a region | Region(s)Number of |
1 | 1 | 1,2,3,4,5,19,20,21,22,23,24,25 | 12 |
2 | 29 | 5,6,7,26,27,28,29,30,31,32,33 | 11 |
3 | 12 | 7,8,9,10,11,12,13,14,15,16,17,18 | 12 |
Table 2 PG and E69 node sample partitioning scheme
Claims (3)
1. A distribution network distributed state estimation oriented synchronous phasor measurement configuration partitioning method is characterized by comprising the following steps:
1) For a selected power distribution system, acquiring a topological connection relation of a power distribution network, constructing an adjacent matrix A, and setting the number M of sub-regions;
2) Establishing a synchronous phasor measurement configuration partition model facing to the distributed state estimation of the power distribution network;
3) Solving the synchronous phasor measurement configuration partition model facing the power distribution network distributed state estimation in the step 2) through a genetic algorithm; the method comprises the following steps:
(3.1) setting genetic AlgorithmThe length of a chromosome of a middle individual is equal to the number M of subregions, the value set of genes on the chromosome is a set omega formed by all nodes of a power distribution system, omega = {1, 2.., N }, N represents the total number of the nodes of the power distribution system to be partitioned, an initial population is generated randomly, an evolution algebra d =1, and a cross probability p is set C Genetic probability p M And a maximum evolution algebra D;
(3.2) obtaining an initial center node through chromosomes of individuals in a group, partitioning the power distribution network by using a center expansion partitioning method, and further calculating the fitness of the individuals; the central extension partition method comprises the following steps:
(3.2.1) setting all the nodes not to be partitioned, randomly acquiring M nodes as initial central nodes, and marking the initial central nodes as partitioned nodes;
(3.2.2) the number of expansion times k =0, updating the sub-region set corresponding to the initial central node according to the initial central node and constructing a partition matrix E k The number of expansion times k = k +1; wherein the partition matrix E k Is a matrix of dimension M × N, for l and i respectively represent a node number and a sub-region number, and if the node l is in the sub-region i, E k Row i and column l element E of k,i,l =1, otherwise E k,i,l Is 0, the set of sub-regions Γ l The middle element is the number of the subregion to which the node l belongs, if the node l is in the subregion i, i belongs to the gamma l ;
(3.2.3) setting a partition matrix E 'to be adjusted after the k expansion' k =E k-1 A and A are adjacent matrixes, and a partition matrix E 'to be adjusted' k All the non-zero elements in the sub-area are set to be 1, i =1, j =2, and i and j are the numbers of the sub-areas;
(3.2.4) if j > M, i = i +1, j = i +1, and entering the (3.2.5) th step, and otherwise, entering the (3.2.6) th step;
(3.2.5) if i = M, entering the (3.2.8) step, otherwise, entering the (3.2.6) step;
(3.2.6) obtaining a set Λ of overlapping nodes of the sub-region i and the sub-region j i,j If the node set is overlapped i,j If the current is an empty set, j = j +1, and the step (3.2.4) is entered, otherwise, the step (3.2.7) is entered;
(3.2.7) if overlapping node set Λ i,j If there are only two nodes in the partition matrix E ' and the two nodes are all partitioned nodes, the two nodes are respectively used as shared boundary nodes of the sub-area i and the sub-area j, the difference between the maximum area node number and the minimum area node number is compared, the boundary node is selected by taking the minimum difference as a target, and if the differences are the same, the node in the sub-area i before expansion is taken as a boundary, and the partition matrix E ' to be adjusted is adjusted ' k Updating a sub-region set corresponding to the boundary node, wherein j = j +1, and entering the (3.2.4) th step;
(3.2.8) let the set of all partitioned nodes be Ω 1, forIf i ∈ Γ l Then the partition matrix E 'is to be adjusted' k Line i of l element E' k,i,l =1, otherwise E' k,i,l =0;
(3.2.9) set the set of all non-partitioned overlapping nodes to Ω 2 To aIf the overlapped subareas have boundaries, dividing the node l into the subarea with the minimum node number, if the subarea with the minimum node number has a plurality of subareas, dividing the node l into the subarea with the minimum serial number, and adjusting the subarea matrix E 'to be adjusted' k Updating the sub-region set corresponding to the node l;
(3.2.10) assuming that i =1, j =2, i and j are both sub-region numbers;
(3.2.11) if j > M, i = i +1, j = i +1, and the step (3.2.12) is entered, otherwise, the step (3.2.13) is entered;
(3.2.12) if i = M, entering the (3.2.15) otherwise entering the (3.2.13) step;
(3.2.13) obtaining the sub-regionOverlap node set Λ of domain i and sub-region j i,j If the node set is overlapped i,j If the current is an empty set, j = j +1, and the step (3.2.11) is entered, otherwise, the step (3.2.14) is entered;
(3.2.14) if overlapping node set Λ i,j The node is not partitioned, and the sub-area i and the sub-area j are overlapped for the first time, then the node is used as a boundary node shared by the sub-area i and the sub-area j, and a partition matrix E 'to be adjusted is adjusted' k Updating a subregion set corresponding to the boundary node, wherein j = j +1, and entering the step (3.2.11);
(3.2.15) setting the adjusted partition matrix E k =E′ k If E is k =E k-1 Then step (3.2.16) is entered, otherwise E is entered k The nodes contained in the data are marked as partitioned nodes, the expansion times k = k +1, and the step (3.2.3) is carried out;
(3.2.16) according to partition matrix E k Outputting a partition result;
(3.3) acting a selection operator on the population, and reserving the individual with the highest fitness;
(3.4) acting the crossover operator and the mutation operator on the selected population to generate a next generation, wherein the evolution algebra d = d +1;
(3.5) if the evolution algebra D is less than D, entering the step (3.2), otherwise, entering the step (3.6);
(3.6) outputting the individual with the highest fitness as an optimal solution;
4) And outputting the obtained synchronous phasor measurement configuration partition scheme oriented to the distribution network distributed state estimation.
2. The distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method according to claim 1, wherein the number M of the sub-regions in the step 1) is as follows:
3. The distribution network distributed state estimation-oriented synchronous phasor measurement configuration partition method according to claim 1, wherein the distribution network distributed state estimation-oriented synchronous phasor measurement configuration partition model in step 2) is a target function with minimum node number difference between sub-areas, and a mathematical expression is as follows:
wherein, N i Indicating the number of nodes contained in sub-region i.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910453337.0A CN110222966B (en) | 2019-05-28 | 2019-05-28 | Distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910453337.0A CN110222966B (en) | 2019-05-28 | 2019-05-28 | Distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110222966A CN110222966A (en) | 2019-09-10 |
CN110222966B true CN110222966B (en) | 2023-03-28 |
Family
ID=67818283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910453337.0A Active CN110222966B (en) | 2019-05-28 | 2019-05-28 | Distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110222966B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111163376B (en) * | 2019-12-16 | 2021-05-14 | 天津大学 | Distribution network synchronous measurement and communication equipment configuration method for distributed state estimation |
CN111983386B (en) * | 2020-08-19 | 2023-05-09 | 剑科云智(深圳)科技有限公司 | Distributed state estimation system of power distribution network |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106921156A (en) * | 2015-12-25 | 2017-07-04 | 中国电力科学研究院 | A kind of active distribution network method for estimating state based on many sampling period hybrid measurements |
CN107563550A (en) * | 2017-08-23 | 2018-01-09 | 武汉大学 | A kind of Optimal Configuration Method of the real-time distributed state estimation of power distribution network based on PMU and PMU |
CN107843810A (en) * | 2017-11-01 | 2018-03-27 | 东南大学 | A kind of active power distribution network fault section tuning on-line method based on state estimation |
CN108879708A (en) * | 2018-08-28 | 2018-11-23 | 东北大学 | A kind of the reactive voltage partition method and system of active distribution network |
CN109193652A (en) * | 2018-10-30 | 2019-01-11 | 贵州电网有限责任公司 | A kind of distribution network failure self-healing system containing distributed generation resource based on Situation Awareness |
CN109412145A (en) * | 2018-10-16 | 2019-03-01 | 国网上海市电力公司 | A kind of active distribution network dynamic characteristic appraisal procedure based on synchro measure data |
-
2019
- 2019-05-28 CN CN201910453337.0A patent/CN110222966B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106921156A (en) * | 2015-12-25 | 2017-07-04 | 中国电力科学研究院 | A kind of active distribution network method for estimating state based on many sampling period hybrid measurements |
CN107563550A (en) * | 2017-08-23 | 2018-01-09 | 武汉大学 | A kind of Optimal Configuration Method of the real-time distributed state estimation of power distribution network based on PMU and PMU |
CN107843810A (en) * | 2017-11-01 | 2018-03-27 | 东南大学 | A kind of active power distribution network fault section tuning on-line method based on state estimation |
CN108879708A (en) * | 2018-08-28 | 2018-11-23 | 东北大学 | A kind of the reactive voltage partition method and system of active distribution network |
CN109412145A (en) * | 2018-10-16 | 2019-03-01 | 国网上海市电力公司 | A kind of active distribution network dynamic characteristic appraisal procedure based on synchro measure data |
CN109193652A (en) * | 2018-10-30 | 2019-01-11 | 贵州电网有限责任公司 | A kind of distribution network failure self-healing system containing distributed generation resource based on Situation Awareness |
Non-Patent Citations (3)
Title |
---|
Optimal PMU placement considering topology constraints;Yuanyuan Zhao等;《ES》;20150520;第240-248页 * |
基于定制遗传算法考虑配电网多种拓扑可观性的PMU优化配置;孔祥玉 等;《电力自动化设备》;20200130;第40卷(第1期);第66-72页 * |
基于新型PMU配置的局部配电网状态估计;吕前程 等;《南方电网技术》;20190430;第13卷(第4期);第54-59页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110222966A (en) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hendrickx et al. | Open multi-agent systems: Gossiping with random arrivals and departures | |
CN110222966B (en) | Distribution network distributed state estimation-oriented synchronous phasor measurement configuration partitioning method | |
Marzano et al. | A kalman filter for quasi-dynamic od flow estimation/updating | |
CN110138620B (en) | Distribution network synchronous measurement and communication link configuration method facing distributed state estimation | |
CN110784883B (en) | Base station construction evaluation method, device, equipment and storage medium | |
CN106780064A (en) | A kind of region partitioning method, device and the network equipment | |
CN111291145A (en) | Mapping method, device and storage medium of wireless hotspot and interest point | |
CN102779142A (en) | Quick community discovery method based on community closeness | |
CN110794263B (en) | Method for positioning fault section of power distribution network line with distributed power supply | |
CN111105005B (en) | Wind power prediction method | |
CN116796904A (en) | Method, system, electronic equipment and medium for predicting new line passenger flow of rail transit | |
CN111163376B (en) | Distribution network synchronous measurement and communication equipment configuration method for distributed state estimation | |
CN110991562B (en) | Animal group geographic division method based on species composition characteristics | |
da Silva et al. | Combining k-means method and complex network analysis to evaluate city mobility | |
CN104507095B (en) | Indoor AP dispositions methods based on user location distribution | |
CN103957530B (en) | The cellular network frequency spectrum distributing method of the meta-heuristic based on figure | |
CN111008730A (en) | Crowd concentration degree prediction model construction method and device based on urban space structure | |
CN115296295A (en) | Low-voltage distribution area topology identification method and system based on cluster analysis and graph theory | |
CN115329814A (en) | Low-voltage user link identification method and device based on image signal processing | |
CN115588086A (en) | Map dividing method, map dividing device, computer readable storage medium and processor | |
CN110504689B (en) | Power system load flow calculation method based on cloud calculation technology | |
CN110312260B (en) | Method and system for optimizing physical identification of overlapped cells of LTE (Long term evolution) power wireless private network | |
CN111210088B (en) | Traffic state index prediction method based on space-time factors | |
Ma et al. | Modelling mobile traffic patterns using a generative adversarial neural networks | |
CN107704540B (en) | Method for rapidly extracting road network in designated administrative district from OSM data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |