CN107196290B - Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization - Google Patents

Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization Download PDF

Info

Publication number
CN107196290B
CN107196290B CN201710243292.5A CN201710243292A CN107196290B CN 107196290 B CN107196290 B CN 107196290B CN 201710243292 A CN201710243292 A CN 201710243292A CN 107196290 B CN107196290 B CN 107196290B
Authority
CN
China
Prior art keywords
branch
bus
measurement
relationship
buses
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
Application number
CN201710243292.5A
Other languages
Chinese (zh)
Other versions
CN107196290A (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.)
Quanzhou Power Supply Co of State Grid Fujian Electric Power Co Ltd
Quanzhou Yixing Electric Power Co Ltd
Original Assignee
Quanzhou Power Supply Co of State Grid Fujian Electric Power Co Ltd
Quanzhou Yixing Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Quanzhou Power Supply Co of State Grid Fujian Electric Power Co Ltd, Quanzhou Yixing Electric Power Co Ltd filed Critical Quanzhou Power Supply Co of State Grid Fujian Electric Power Co Ltd
Priority to CN201710243292.5A priority Critical patent/CN107196290B/en
Publication of CN107196290A publication Critical patent/CN107196290A/en
Application granted granted Critical
Publication of CN107196290B publication Critical patent/CN107196290B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/39Circuit design at the physical level
    • G06F30/394Routing
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

The invention discloses a branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization, which is established on the basis of error measurement detection, identification and correction, defines a branch and bus association relation island by analyzing the association relation of branches and buses, enables the measurement balance index of each node in the island to be qualified by searching the branch and bus connection relation with the actually acquired equipment state with the highest similarity as possible based on the characteristic that the measurement of each branch in the island meets the balance relation, and gives an identification result. Simulation analysis and practical application prove the effectiveness of the method.

Description

Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization
Technical Field
The invention relates to a branch and bus connection relation identification method, in particular to a branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization.
Background
For the dispatching automation system, the correctness of the network topology is very important because the network topology directly influences the state estimation calculation, and the calculation result of the state estimation calculation is a data source of a series of network analysis software including power flow calculation, reactive voltage sensitivity calculation, static safety analysis and the like. The correctness of the network topology is established on the basis of correct remote signaling data, so that the detection, identification and correction of errors of the network structure are important contents for real-time information processing of the dispatching automation system.
The common network topology detection and identification method is to utilize the remote signaling of switch and the remote measurement of line load flow to check if they are correspondent to find out suspicious remote signaling, then use search method to check and detect according to the possible network structure form by means of state estimation
Figure BDA0001269793190000011
Whether the value is smaller than the threshold value or not is judged to judge whether the network structural form is correct or not. The method has the advantages of simplicity and low requirement on measurement configuration. But failures may occur when multiple structural errors are simultaneously erroneous with bad data. In the automatic system of most dispatching mechanisms in China, widely adopted network structure identification is only completed independently by comparing whether telemetering signals and telesignaling signals are consistent, and only in the coarse detection stage of basic two teledata, prompt alarm information is given, but no further identification and correction measures are taken during state estimation and calculation, but network topology analysis is carried out based on actually acquired equipment states, so that the calculation result has larger deviation with the actual operation condition of a power grid。
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and provides a branch and bus connection relation identification method which is used for identifying the result and optimizing based on the measurement balance index and the equipment state similarity by analyzing the correlation relation between branches and buses, defining a branch and bus correlation relation island, and searching the branch and bus connection relation which enables the measurement balance index of each node in the island to be qualified and has the highest similarity with the actually acquired equipment state based on the characteristic that the measurement of each branch in the island meets the balance relation.
A branch and bus connection relationship identification method based on measurement balance index and equipment state similarity optimization comprises the following steps:
(1) Error measurement identification and correction
Searching for a branch associated with a bus that is not qualified for measuring balance indexes, if the branch satisfies the following conditions: a. the number of the ports is not less than 2; b. if the ports of the branch are configured with measurement, the measurement balance index of the branch can be detected and still used as P' L Denotes a measurement balance index of the line, denoted by P L Indicating absence of home port T k The line measurement balance index;
if P' L Satisfies a balance criterion, and P ″) L Also satisfies the balance criterion, then considers the port T k Is not connected to a node; if P' L Satisfies the balance criterion, and P L If the balance criterion is not satisfied, the port T is considered to be k Connecting with the node;
if P' L Do not meet the balance criterion, need to be on T k And measuring other branch circuit ports, and detecting when the following conditions are met: 1. for a branch port connected with a node, the measurement balance index of the connected node meets a balance criterion; 2. for the branch port not connected with the node, the measurement is 0; available P ″) L Alternative port T k And considering port T k Connecting with the node;
(2) Bus tie relationship identification
If communication paths exist among the buses, the communication relation among the buses is identified. And identifying the bus connection relation according to the bus connection relation reflected by the dispatching automation system and the measurement configuration on the bus communication path.
a. The bus connection relation reflected by the dispatching automation system is connection
If the bus communication path is not configured with measurement, the bus communication relation is doubtful; if the bus communication path is configured with measurement, checking whether the measurement is 0, if the measurement is not 0, the bus communication relationship is identified as connection, otherwise, the bus communication relationship is doubted.
b. The bus bar connection relation reflected by the dispatching automation system is unconnected
If the bus communication path is not configured with measurement, the bus communication relation is doubtful; if the bus communication path is configured with measurement, then checking whether the measurement is 0, if the measurement is 0, the identification result is that the bus communication relationship is not connected, otherwise, the identification result is that the bus communication relationship is doubted.
(3) Incidence relation of branch and bus
The association between the branch and the bus bar may be a one-to-one relationship, that is, the branch has only a communication path to one bus bar; a one-to-many relationship is also possible, i.e. a branch having a communication path to a plurality of busbars. The relationship between the branch and the bus can be expressed by a graph, wherein the vertex of the graph corresponds to the branch and the bus, and the arc between the vertex and the bus corresponds to a communication path from the branch to the bus. Therefore, a plurality of subgraphs without connection relation are formed, and each subgraph corresponds to the topological connection relation between a group of branches and the bus, and is called a topological island. Assuming no communication between topological islands, the power flow of all branches within an island should satisfy a balance relationship for each topological island.
If all the faulty branch measurements are corrected, all the relevant branch measurements should satisfy the measurement balance relationship at a certain voltage level, the measurement balance index of each island can be recorded, if the measurement balance index of a certain island is qualified, the island is not in contact with other islands, otherwise, the island is necessarily in contact with other islands. Based on the point, the inter-island contact relationship can be further checked, and meanwhile, for the island with unqualified measurement balance indexes, an inter-island contact path can be tried to be found, so that a new island is formed to meet the measurement balance relationship.
(4) Branch and bus connection relation identification
When the branch measurement of an island is balanced, the reasons for the failure of the node measurement balance index in the island can be classified into the following three categories: (1) the contact relation between the buses is wrong (the contact relation is questioned); (2) measuring branch circuits which are not 0 and are not connected with the connecting bus; (3) the connection relation between the branch and the bus is wrong.
1. Single bus connection relationship identification
If the branch boundary connection point only has a communication path to 1 bus and the communication path to the bus is not only 1, then since the path communicated with the same bus is not different in topology, in order to judge which path meets the actual situation, the following judging method is introduced:
searching whether measurement is configured on a port through which the path passes, if so, comparing the measurement with the measurement of the branch, and if the difference value between the two is small enough, determining that the path is established;
if all paths are not provided with ports for measurement, selecting the paths by introducing a similarity discrimination method based on the small probability event that the equipment state error belongs to, wherein the method comprises the following steps: assuming that paths i and j are both paths leading from the branch boundary node to a certain bus, assuming that the paths are true, the number of elements in the path i with correct state is n it The number of elements with a state error is n if (ii) a The number of elements in path j that are in the correct state is n jt The number of elements with a state error is n jf . Then the similarity of the device states of the two paths is:
similar i =n it /(n it +n if ),similar j =n jt /(n jt +n jf )
by comparing similar i And similar j A path with greater device state similarity is selected.
2. Multi-bus connection relationship identification
a. Branch and bus connection path coding
If the branch boundary connection point has communication paths to different buses, the connection relationship between one branch and at most two buses is realized according to the actual power grid connection and operation mode. According to the characteristic, the connection relation between the branch circuits and the buses and the communication relation between the buses are coded: by using
Figure BDA0001269793190000051
To represent the connection relationship between the branch and the bus, which is "1" or "0"; by a i To express the connection relationship between the buses, for one island, the connection relationship between each branch and the bus can be coded, where the number of branches is m, the number of bus connections in question is n, and the coding is:
Figure BDA0001269793190000052
b. method for optimizing similarity based on measurement balance index and equipment state
The reasonable identification result should satisfy the following 2 principles: (1) the measurement balance index of each node in the island must be qualified; (2) it must have as high a similarity as possible to the actual collected branch and bus connection. Therefore, a method for optimizing based on measurement balance indexes and equipment state similarity is introduced, and the following equipment state similarity function is defined:
if the identification result is consistent with the actual situation, the identification result is used for checking the state of the equipment obtained by actual acquisition, and m is used t N is used for representing the number of connecting paths between the correct branch and the bus t Representing the correct number of bus connections, the similarity is represented as follows:
similar=(m t +n t )/(m+n)
and when the identification result meets the following conditions and the measurement balance indexes of all the nodes in the island are qualified, the identification result is considered to be reasonable.
similar≥(m+n-2)/(m+n)
Step 1:
to substring { l 1 ,l 2 …,l n And (5) carrying out small probability mutation operation, and randomly changing 1 bus connection relation. If the substring is mutated, the connection relation between the branch and the bus is not changed, the measurement balance index of each node in the island is calculated, and if the jth suspected bus connection relation is mutated, the code of the connection relation between the branch and the bus is as follows:
Figure BDA0001269793190000061
note: />
Figure BDA0001269793190000062
Is to be paired with j And (4) taking the opposite.
If the substring is not mutated, the following operations are carried out:
determining two buses related to the branch circuit according to the branch circuit bus association relationship diagram, and when the following two conditions are simultaneously met: (1) the two buses belong to different nodes; (2) if the branch measurement value is in the same direction as the measurement balance index of the connected node (the branch measurement and the measurement balance index of the connected node are both positive or negative), the branch is taken as the branch of the alternative 1 adjustable path. Searching all adjustable path branches to form an adjustable path branch set B 1 If branch set B 1 If not, the slave branch set B 1 Randomly selecting 1 branch circuit, changing the connection path between the branch circuit and the bus, if the branch circuit is collected B 1 If the bus is empty, the two buses are searched for the branch belonging to the same node, and the branch is classified into an adjustable path branch set B 2 Set of slave branches B 2 Randomly selecting 1 branch circuit, and changing the connection path of the branch circuit and the bus;
setting an adjusting branch b i And if the connection path is connected with the bus, the connection relation between the branch and the bus is coded as follows:
Figure BDA0001269793190000063
note: />
Figure BDA0001269793190000064
Is paired with>
Figure BDA0001269793190000065
The result is obtained by taking the inverse of the formula,
step 2:
according to the method, a population with the scale of m can be obtained after repeating for m times, each individual in the population corresponds to the connection relation between each branch and a bus in an island, the measurement balance index and the equipment state similarity of each node in the island corresponding to the individual can be respectively calculated, and the individual which can enable the measurement balance indexes of all the nodes to be qualified and the similarity to meet the requirement is recorded without further operation; otherwise, repeating the operation of the step 1. After the step is finished, the population of the 1 st generation can be obtained;
and step 3:
and (3) repeating the steps 1 and 2 on the basis of the population of the 1 st generation to obtain a new population of the 1 st generation. After several generations of operations, all reasonable branch and bus connection relations can be obtained, and finally, an identification result is given.
Drawings
Fig. 1 is a diagram of an electrical main wiring having a plurality of islands of branch and bus bar associations.
Fig. 2 is a diagram of an electrical main wiring having only 1 branch associated with a bus bar island.
Fig. 3 is a diagram showing the relationship between branch busbars corresponding to the main wiring shown in fig. 1.
Fig. 4 is a diagram showing the relationship between branch busbars corresponding to the main wiring shown in fig. 2.
Detailed Description
The present invention will be described in more detail with reference to examples.
Example 1
1. Error side data identification and correction
The abnormal measurement and state data on the power grid dispatching control system can cause the abnormal measurement and state data, so the abnormal measurement and state data can be detected and identified by detecting and analyzing the abnormal measurement and state data of the measurement and balance index.
1.1 detection basis of bus measurement balance index
Setting the CT transformation ratio of the interval measurement on the bus as I CT PT ratio of/5 to U PT 100, the reference value of the interval power is therefore:
Figure BDA0001269793190000071
the practical standard of the dispatching automation system stipulates that the allowable error of the active/reactive power quantity side is 1% of the reference value, that is, if the power reference value of a certain interval is S B Then, according to the above standard, if N intervals are connected to the bus in total (according to the result of the network topology analysis), the allowable bus active/reactive power measurement balance index is:
Figure BDA0001269793190000081
if the absolute value of the bus power measurement balance index exceeds the detected value
Figure BDA0001269793190000082
It can be determined that the bus bar measurement balance index is abnormal.
1.2 detection basis of line measurement balance index
If the line totally associates N stations (according to the analysis result of the network topology), the active measurement of each station is P' 1 、P′ 2 、…P′ N (considering that the positive power reference direction is that the outgoing bus is positive and the incoming bus is negative), the active power measurement balance index of the line is as follows:
Figure BDA0001269793190000083
of formula (II) to P' L And the line active power measurement balance index is represented.
Setting the power reference value of each station end as S B1 、S B2 、…S BN Then, the allowed line active power measurement balance index is:
Figure BDA0001269793190000084
in the formula-P L Representing the line active power loss.
Assuming that the statistical line loss rate does not exceed Kl, one can obtain:
Figure BDA0001269793190000085
in the formula-P + Representing the line input active sum, P - Representing the line output active sum.
If the number is 1, …, N + Is the input end of the line, numbered N + The interval of +1, …, N is the output end of the line, then the equations (1-4), (1-5) can be obtained:
Figure BDA0001269793190000091
of formula (II) to P' + Representing the total line input power measurement sum.
Figure BDA0001269793190000092
Of formula (II) to P' - And the total sum of the line output power measurement is represented.
Get
Figure BDA0001269793190000093
And &>
Figure BDA0001269793190000094
The smaller of which is identified as Min P And substituting the following formula (1-2) to obtain the following line active power measurement balance index judgment data:
Figure BDA0001269793190000095
and when the active power measurement balance index of the line exceeds the numerical range shown in the formula (1-6), judging that the active power measurement balance index is abnormal.
The detection basis of the active power measurement balance index of the transformer, the derivation method and the detection basis of the line active power measurement balance index are similar, and are not repeated.
1.3 identifying and correcting error measurements
Searching for a branch associated with a bus that is not qualified for measuring balance indexes, if the branch satisfies the following conditions: 1. the number of the ports is not less than 2; 2. measurements are allocated to each port of the branch. The measurement balance index of the branch can be detected. Still from P' L Denotes a measurement balance index of the line, denoted by P L Indicating absence of home port T k The line measurement balance index.
If P' L Satisfies a balance criterion, and P ″) L Also satisfies the balance criterion, then considers the port T k Are not connected to a node; if P' L Satisfies the balance criterion, and P L If the balance criterion is not satisfied, the port T is considered to be k Connecting with the node;
if P' L Do not meet the balance criterion, need to be on T k And measuring other branch ports. When the following conditions are satisfied: 1. for a branch port connected with a node, the measurement balance index of the connected node meets a balance criterion; 2. for a branch port not connected to a node, its measurement is 0; available P ″) L Alternative port T k And considering port T k Connected to the node.
Because the number of detected balanced nodes and balanced branches in an actual dispatching automation system is always large, in most cases, the error measurement data in the system can be identified and corrected by using an alternating iterative calculation method of 'correcting unbalanced branches by balanced nodes and correcting unbalanced branches by balanced branches'.
2. The invention relates to a method for identifying the connection relation between a power grid branch and a bus based on measurement balance indexes and equipment state similarity optimization
2.1 bus tie relationship identification
TABLE 2-1 bus bar connection relation table
Figure RE-GDA0001342733000000101
The connection relation between the buses is judged according to the attached table 2-1, and the bus connection relation with doubtful state is kept in the current state.
2.2 Association of Branch and bus
The association relationship between the branch circuits and the bus bars may be a one-to-one relationship, that is, the branch circuits only have communication paths to one bus bar; a one-to-many relationship is also possible, i.e. the branches have communication paths to a plurality of busbars. The relationship between the branch and the bus can be expressed by a graph, wherein the vertex of the graph corresponds to the branch and the bus, and the arc between the vertex and the bus corresponds to a communication path from the branch to the bus. Therefore, a plurality of subgraphs without connection relation are formed, and each subgraph corresponds to the topological connection relation between a group of branches and the bus, and is called a topological island. Assuming no communication between topological islands, the power flow of all branches within an island should satisfy a balance relationship for each topological island.
If all the faulty branch measurements are corrected, and all the relevant branch measurements should satisfy the measurement balance relationship at a certain voltage level, the measurement balance index of each island can be recorded, if the measurement balance index of a certain island is qualified, the island is not in a connection relationship with other islands, otherwise, the island is necessarily in communication with other islands. Based on the point, the inter-island contact relationship can be further checked, and meanwhile, for the island with unqualified measurement balance indexes, an inter-island contact path can be tried to be found, so that a new island is formed to meet the measurement balance relationship.
2.3 Branch and bus connection relationship identification
2.3.1 Single bus connection identification
If the branch boundary connection point only has communication paths to 1 bus and the communication paths to the bus are not only 1, then since the paths communicated with the same bus are not different in topology, in order to judge which path meets the actual situation, the following judgment method is introduced:
searching whether measurement is configured on a port through which the path passes, if so, comparing the measurement with the measurement of the branch, and if the difference value between the two is small enough, determining that the path is established;
if all paths are not provided with ports for measurement, selecting the paths by introducing a similarity discrimination method based on the small probability event that the equipment state error belongs to, wherein the method comprises the following steps: assuming that paths i and j are both paths leading from the branch boundary node to a certain bus, assuming that the paths are true, the number of elements in the path i with correct state is n it The number of elements with a state error is n if (ii) a The number of elements in path j that are in the correct state is n jt The number of elements with a state error is n jf . Then the similarity of the device states of the two paths is:
similar i =n it /(n it +n if ),similar j =n jt /(n jt +n jf )
by comparing similar i And similar j A path with greater device state similarity is selected.
2.3.2 Multi-bus connection relationship identification
a. Branch and bus connection path coding
If the branch boundary connection point has communication paths to different buses, the connection relationship between one branch and at most two buses is realized according to the actual power grid connection and operation mode. According to the characteristic, the connection relation between the branch and the bus and the communication relation between the buses are coded: by using
Figure BDA0001269793190000121
To represent the connection relation between the branch and the bus, which is '1' or '0'; by a i To express the connection relationship between the buses, for one island, the connection relationship between each branch and the bus can be coded, where the number of branches is m, the number of bus connections in question is n, and the coding is:
Figure BDA0001269793190000122
b. method for optimizing similarity based on measurement balance index and equipment state
Reasonable identification results should satisfy the following 2 principles: (1) the measurement balance index of each node in the island must be qualified; (2) it must have as high a similarity as possible to the actual collected branch and bus connection. Therefore, a method for optimizing based on measurement balance indexes and equipment state similarity is introduced, and the following equipment state similarity function is defined:
if the identification result is consistent with the actual situation, the identification result is used to check the state of the equipment actually acquired, and m is used t N is used for representing the number of correct branch and bus connecting paths t Representing the correct number of bus connections, the similarity is represented as follows:
similar=(m t +n t )/(m+n)
and when the identification result meets the following conditions and the measurement balance indexes of all the nodes in the island are qualified, the identification result is considered to be reasonable.
similar≥(m+n-2)/(m+n)
Step 1:
to substring { l 1 ,l 2 …,l n And (5) carrying out small probability mutation operation, and randomly changing 1 bus connection relation. If the substring is mutated, the connection relation between the branch and the bus is not changed, the measurement balance index of each node in the island is calculated, and if the j suspected bus connection relation is mutated, the branch and the bus areThe coding of the bus connection relation is as follows:
Figure BDA0001269793190000131
note: />
Figure BDA0001269793190000132
Is a pair of j And taking the opposite.
If the substring is not mutated, the following operations are carried out:
determining two buses related to the branch according to the branch bus association relationship diagram, and when the following two conditions are simultaneously met: (1) the two buses belong to different nodes; (2) if the branch measurement value is in the same direction as the measurement balance index of the connected node (the branch measurement and the measurement balance index of the connected node are both positive or negative), the branch is taken as the branch of the alternative 1 adjustable path. Searching all adjustable path branches to form an adjustable path branch set B 1 If branch set B 1 If not, the slave branch set B 1 Randomly selecting 1 branch circuit, changing the connection path between the branch circuit and the bus, if the branch circuit set B 1 If the bus is empty, the two buses are searched for the branch belonging to the same node, and the branch is classified into an adjustable path branch set B 2 Set of slave branches B 2 Randomly selecting 1 branch circuit, and changing the connection path of the branch circuit and the bus;
setting an adjusting branch b i And if the connection path is connected with the bus, the connection relation between the branch and the bus is coded as follows:
Figure BDA0001269793190000133
note: />
Figure BDA0001269793190000134
Is paired with>
Figure BDA0001269793190000135
The result is obtained by taking the inverse of the formula,
step 2:
according to the method, a population with the scale of m can be obtained after repeating for m times, each individual in the population corresponds to the connection relation between each branch and the bus in the island, the measurement balance index and the equipment state similarity of each node in the island corresponding to the individual can be respectively calculated, and the individual which can enable all the node measurement balance indexes to be qualified and the similarity to meet the requirement is recorded without further operation; otherwise, repeating the operation of the step 1. After the step is finished, the population of the 1 st generation can be obtained;
and step 3:
and (3) repeating the steps 1 and 2 on the basis of the population of the 1 st generation to obtain a new population of the 1 st generation. After several generations of operations, all reasonable branch and bus connection relations can be obtained, and finally, an identification result is given.
3. Specific applications of the invention
Taking 220 kv clean Mongolian to 110 kv system in the quanza power grid as an example, the following is an initial data section reflected by the scheduling control system:
TABLE 3-1 Primary model and real-time section of regional grid
Figure RE-GDA0001342733000000141
Figure RE-GDA0001342733000000151
In this example, the 110 kv system is a typical double-bus connection, and the measurement of all the branches satisfies the balance relationship, so the 110 kv system forms a single island, and the equipment state identification method based on the measurement balance index and similarity optimization described herein can be adopted.
(1) Changing the buses connected with the No. 1 main transformer 18A and the No. 2 main transformer 18B, simulating the connection relation between 1 wrong branch and the buses, recording as case 1, and calculating the following results:
TABLE 3-1 case 1 optimization procedure calculation results for each generation
Figure RE-GDA0001342733000000152
Figure RE-GDA0001342733000000161
Table 3-2 case 1 optimization procedure calculation time (seconds) for each generation
Figure RE-GDA0001342733000000162
The optimization process finally provides a unique feasible solution individual, the feasible solution is obtained after the 1 st generation of calculation, and no new feasible solution appears after the 2 nd generation.
(2) Changing the bus connected with the clear flag line 185 and the clear shop II 193, simulating the connection relation between 1 wrong branch and the bus, and recording as case 2, wherein the calculation result is as follows:
table 3-3 case 2 optimization procedure calculation results of each generation
Figure RE-GDA0001342733000000163
Figure RE-GDA0001342733000000171
Table 3-4 case 2 optimization procedure calculation time (seconds) for each generation
Figure RE-GDA0001342733000000172
As can be seen from tables 3-3, the optimization process finally provides two feasible solution individuals, which are obtained after the 1 st generation calculation, and no new feasible solution appears after the 2 nd generation. However, a single individual needs 2 generations at least to obtain all feasible solutions, and part of individuals still obtain only 1 feasible solution after 5 generations of calculation. As can be seen from tables 3-4, the algorithm takes only about 1 minute to obtain all feasible solutions.
The parts not described in the present embodiment are the same as those in the prior art.

Claims (1)

1. A branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization is characterized by comprising the following steps:
(1) Error measurement identification and correction:
searching for a branch associated with a bus with unqualified measurement balance index, if the branch meets the following conditions: a. the number of the ports is not less than 2; b. if the branch is configured with measurement at each port, the measurement balance index of the branch is detected and still P is used L ' indicating the measurement balance index of the line, P L "means not including the port T of the home side k The line measurement balance index;
if P is L ' satisfy the balance criterion, and P L "the balance criterion is also met, then Port T is considered k Is not connected with the node; if P is L ' satisfy the balance criterion, and P L "not meeting the balance criterion, then consider the port T k Connecting with the node;
if P is L ' not meeting the balance criterion, requiring a pair of T k The measurement of other branch ports is detected when the following conditions are met: 1. for a branch port connected with a node, the measurement balance index of the connected node meets a balance criterion; 2. for a branch port not connected to a node, its measurement is 0; by P L "replacement Port T k And considering port T k Connecting with the node;
(2) And (3) identifying the bus connection relation:
if there is a communication path between buses, the communication relationship between buses should be identified, the communication relationship between buses is identified according to the bus communication relationship reflected by the dispatching automation system and the measurement configuration on the bus communication path,
a. the bus connection relation reflected by the dispatching automation system is connection:
if the bus communication path is not configured with measurement, the bus communication relation is doubtful; if the bus communication path is configured with measurement, checking whether the measurement is 0, if the measurement is not 0, identifying that the bus communication relationship is connection, otherwise, identifying that the bus communication relationship is doubtful;
b. the bus connection relation reflected by the dispatching automation system is unconnected:
if the bus communication path is not configured with measurement, the bus connection relation is questioned; if the bus communication path is configured with measurement, then checking whether the measurement is 0, if the measurement is 0, the bus communication relationship is identified as unconnected, otherwise, the bus communication relationship is doubted,
(3) The incidence relation of the branch and the bus is as follows:
when the correlation between the branch and the bus is one-to-one, the branch only has a communication path to one bus; when the association relationship between the branch circuits and the buses is a one-to-many relationship, namely the branch circuits are provided with communication paths to a plurality of buses, the association relationship between the branch circuits and the buses is expressed by a graph mode, the top points of the graph correspond to the branch circuits and the buses, and arcs between the top points and the buses correspond to one communication path from the branch circuits to the buses, so that a plurality of sub-graphs without connection relationship with each other are formed, each sub-graph corresponds to a group of topological connection relationship between the branch circuits and the buses and is called a topological island, and when the topological islands are not connected, the tide of all the branch circuits in the island should satisfy a balance relationship;
if all the wrong branch measurements are corrected, all the related branch measurements should meet the measurement balance relationship under a certain voltage level, recording the measurement balance index of each island, if the measurement balance index of a certain island is qualified, not having the contact relationship with other islands, otherwise having the contact with other islands, further checking the contact relationship between the islands based on the point, and trying to find a contact path between the islands for the island with the unqualified measurement balance index, thereby forming a new island to meet the measurement balance relationship;
(4) Identifying the connection relation between the branch and the bus:
when the branch measurement of an island is balanced, the reasons for the disqualification of the node measurement balance index in the island are divided into the following three categories: (1) the contact relation between the buses is wrong, namely the contact relation is questioned; (2) measuring the branch not 0 and the connecting bus; (3) the connection relationship between the branch and the bus is wrong,
A. single bus connection relationship identification
If the branch boundary connection point only has communication paths to 1 bus and the communication paths to the bus are not only 1, then since the paths communicated with the same bus are not different in topology, in order to judge which path meets the actual situation, the following judgment method is introduced:
searching whether measurement is configured on a port through which the path passes, if so, comparing the measurement with the measurement of the branch, and if the difference value between the two is small enough, determining that the path is established;
if all paths are not provided with ports for measurement, selecting the paths by introducing a similarity discrimination method based on the small probability event that the equipment state error belongs to, wherein the method comprises the following steps: assuming that paths i and j are both paths leading from the branch boundary node to a certain bus, assuming that the paths are true, the number of elements in the path i with correct state is n it The number of elements with a state error is n if (ii) a The number of correctly-conditioned elements on path j is n jt The number of elements with a state error is n jf Then, the similarity of the device states of the two paths is:
similar i =n it /(n it +n if ),similar j =n jt /(n jt +n jf )
by comparing similar i And similar j Selecting a path with greater equipment state similarity;
B. multi-bus connection relationship identification
a. Branch and bus connection path coding
If the branch boundary connection point has communication paths to different buses, at most two buses have a communication relation for one branch according to the actual power grid connection and operation mode, and according to the characteristic, the branch is connected with the busesAnd coding the connection relation of the buses and the communication relation among the buses: by using
Figure FDA0004058446000000031
To represent the connection relationship between the branch and the bus, which is "1" or "0"; by using i To express the connection relationship between the buses, for one island, the connection relationship between each branch and the bus is coded, and the number of branches is set as m, the number of bus connections in question is set as n, and the coding is:
Figure FDA0004058446000000041
b. method for optimizing similarity based on measurement balance index and equipment state
Reasonable identification results should satisfy the following 2 principles: (1) the measurement balance index of each node in the island must be qualified; (2) the method has the following advantages that the method has the similarity with the connection relation between actually collected branches and buses as high as possible, and therefore, a method for optimizing the similarity based on measurement balance indexes and equipment states is introduced, and the following equipment state similarity function is defined:
if the identification result is consistent with the actual situation, the identification result is used to check the equipment state obtained by actual collection, and m is used t N is used for representing the number of correct branch and bus connecting paths t Representing the correct bus bar tie coefficient, the similarity is represented as follows:
similar=(m t +n t )/(m+n)
when the identification result satisfies the following conditions and the measurement balance indexes of each node in the island are all qualified, the identification result is considered to be reasonable,
similar≥(m+n-2)/(m+n)
step 1:
for substring { l } 1 ,l 2 …,l n Adopting small probability variation operation to randomly change 1 bus connection relation, if substrings are varied, not changing connection relation of branch and bus, calculating measurement balance index of each node in islandIn the third embodiment, if the jth suspected bus tie relationship is mutated, the codes of the branch and bus connection relationships are:
Figure FDA0004058446000000042
annotating>
Figure FDA0004058446000000043
Is a pair of j The result is obtained by taking the inverse of the formula,
if the substring is not mutated, the following operations are carried out:
determining two buses related to the branch circuit according to the branch circuit bus association relationship diagram, and when the following two conditions are simultaneously met: (1) the two buses belong to different nodes; (2) the branch measurement value is in the same direction with the measurement balance index of the connected node, that is, the branch measurement and the measurement balance index of the connected node are both positive or negative, then the branch is used as the branch of the alternative 1 adjustable paths, and the branches of all adjustable paths are searched to form an adjustable path branch set B 1 If branch set B 1 If not, the slave branch set B 1 Randomly selecting 1 branch circuit, changing the connection path between the branch circuit and the bus, if the branch circuit set B 1 If the bus is empty, the two buses are searched for the branch belonging to the same node, and the branch is classified into an adjustable path branch set B 2 Set of slave branches B 2 Randomly selecting 1 branch circuit, and changing the connection path of the branch circuit and the bus;
setting an adjusting branch b i And the connection path with the bus, the connection relation between the branch and the bus is coded as follows:
Figure FDA0004058446000000051
annotating>
Figure FDA0004058446000000052
Is paired with>
Figure FDA0004058446000000053
The process is carried out in a reverse manner,
step 2:
according to the method, a population with the scale of m is obtained after repeating for m times, each individual in the population corresponds to the connection relation between each branch and a bus in the island, the measurement balance index and the equipment state similarity of each node in the island corresponding to the individual are respectively calculated, and the individual which can enable the measurement balance indexes of all the nodes to be qualified and the similarity to meet the requirement is recorded without further operation; otherwise, repeating the operation of the step 1, and obtaining the population of the 1 st generation after the step is finished;
and step 3:
and (3) repeating the steps 1 and 2 on the basis of the population of the generation 1 to obtain a new generation 1 population, obtaining all reasonable branch and bus connection relations after a plurality of generations of operations, and finally giving an identification result.
CN201710243292.5A 2017-04-14 2017-04-14 Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization Active CN107196290B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710243292.5A CN107196290B (en) 2017-04-14 2017-04-14 Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710243292.5A CN107196290B (en) 2017-04-14 2017-04-14 Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization

Publications (2)

Publication Number Publication Date
CN107196290A CN107196290A (en) 2017-09-22
CN107196290B true CN107196290B (en) 2023-04-07

Family

ID=59870992

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710243292.5A Active CN107196290B (en) 2017-04-14 2017-04-14 Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization

Country Status (1)

Country Link
CN (1) CN107196290B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110867906B (en) * 2019-12-02 2021-11-02 南京工程学院 Power grid topological relation identification method based on power balance index

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102005757A (en) * 2010-10-29 2011-04-06 江西电力调度中心 Method for recognizing grid topology based on measurement spanning tree
KR101396131B1 (en) * 2013-11-18 2014-05-19 한국과학기술정보연구원 Apparatus and method for measuring relation similarity based pattern
CN104466959A (en) * 2014-12-31 2015-03-25 广东电网有限责任公司电力调度控制中心 Power system key line identification method and system
CN104901308A (en) * 2015-06-26 2015-09-09 太原理工大学 Power system critical circuit identification method
CN105140919A (en) * 2015-09-18 2015-12-09 国网山东省电力公司菏泽供电公司 Method for rapid topology analysis and establishment of topological island based on graph theory
CN106505557A (en) * 2016-11-14 2017-03-15 中国电力科学研究院 A kind of remote measurement misidentification method and device

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102005757A (en) * 2010-10-29 2011-04-06 江西电力调度中心 Method for recognizing grid topology based on measurement spanning tree
KR101396131B1 (en) * 2013-11-18 2014-05-19 한국과학기술정보연구원 Apparatus and method for measuring relation similarity based pattern
CN104466959A (en) * 2014-12-31 2015-03-25 广东电网有限责任公司电力调度控制中心 Power system key line identification method and system
CN104901308A (en) * 2015-06-26 2015-09-09 太原理工大学 Power system critical circuit identification method
CN105140919A (en) * 2015-09-18 2015-12-09 国网山东省电力公司菏泽供电公司 Method for rapid topology analysis and establishment of topological island based on graph theory
CN106505557A (en) * 2016-11-14 2017-03-15 中国电力科学研究院 A kind of remote measurement misidentification method and device

Also Published As

Publication number Publication date
CN107196290A (en) 2017-09-22

Similar Documents

Publication Publication Date Title
Peppanen et al. Distribution system low-voltage circuit topology estimation using smart metering data
CN108898311A (en) A kind of data quality checking method towards intelligent distribution network repairing dispatching platform
CN111999692B (en) Multi-table error calibration method and device
CN110867906B (en) Power grid topological relation identification method based on power balance index
CN108199891B (en) Cps network attack identification method based on artificial neural network multi-angle comprehensive decision
CN104618064B (en) Static information monitoring method based on CRC check and device
CN108919044B (en) Active identification method for unit distribution power grid faults based on mutual verification mechanism
CN103869192A (en) Smart power grid line loss detection method and system
CN105404269B (en) A kind of remote signalling Acceptance Test System of intelligent substation three
WO2020177706A1 (en) Method for determining fault section and fault time of power distribution network
CN104834305B (en) Distribution automation terminal remote measurement exception analysis system and method based on DMS systems
CN107196290B (en) Branch and bus connection relation identification method based on measurement balance index and equipment state similarity optimization
CN111999691B (en) Error calibration method and error calibration device for metering sensor device
CN109714788A (en) The sensor node fault diagnosis algorithm of online fault-tolerant and neighbours cooperation
CN108594075B (en) Power distribution network power failure fault positioning method based on improved ant colony algorithm
CN110443481A (en) Distribution power automation terminal state evaluation system and method based on mixing K- nearest neighbor algorithm
CN110909774B (en) Power transformation equipment heating defect reason distinguishing method based on Bayesian classification
CN111208464A (en) Online evaluation system and method for measurement accuracy of primary and secondary power distribution complete equipment
CN115951123A (en) Electric energy metering method and system based on wireless communication
CN110866615A (en) Transformer substation, fault-free time calculation method thereof and loop state maintenance method
CN113283881B (en) Automatic auditing method and system for telecontrol information source
CN113191062B (en) Power distribution network fault section positioning method and system based on multisource incomplete information
CN115954865A (en) Power distribution network low-voltage topology analysis method and system based on artificial intelligence
CN113391256B (en) Electric energy meter metering fault analysis method and system of field operation terminal
CN112684300B (en) Active power distribution network fault diagnosis method and device using bidirectional monitoring information

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