CN111141992B - Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data - Google Patents

Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data Download PDF

Info

Publication number
CN111141992B
CN111141992B CN201910966007.1A CN201910966007A CN111141992B CN 111141992 B CN111141992 B CN 111141992B CN 201910966007 A CN201910966007 A CN 201910966007A CN 111141992 B CN111141992 B CN 111141992B
Authority
CN
China
Prior art keywords
node
port
value
data
section
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
CN201910966007.1A
Other languages
Chinese (zh)
Other versions
CN111141992A (en
Inventor
肖博文
肖祥
刘强
李小可
朱绍杰
刘随阳
胡小青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Jingmen Power Supply Co of State Grid Hubei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Jingmen Power Supply Co of State Grid Hubei 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 State Grid Corp of China SGCC, Jingmen Power Supply Co of State Grid Hubei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201910966007.1A priority Critical patent/CN111141992B/en
Publication of CN111141992A publication Critical patent/CN111141992A/en
Application granted granted Critical
Publication of CN111141992B publication Critical patent/CN111141992B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08—Locating faults in cables, transmission lines, or networks
    • G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08—Locating faults in cables, transmission lines, or networks
    • G01R31/088—Aspects of digital computing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

An adaptive layered sector location method using FTU telemetry and telemetry data, comprising the steps of: step S1: separating a plurality of two ports of the whole power distribution network by taking the multi-branch node as a boundary, and obtaining related data of each two port; step S2: obtaining incidence matrix data about the two ports according to the topological structure of the power distribution network; step S3: reading remote signaling data of a port node acquired by an FTU (fiber to the Unit), and establishing a switching function and a fitness function of a first-layer positioning model; step S4: reading remote signaling data of a common node inside a fault port acquired by an FTU; step S5: reading telemetry data of a common node inside a fault port acquired by an FTU (fiber to the Unit), and positioning a fault section in a second layer positioning model; step S6: the invention sends out alarm to the distorted position of the remote signaling data, and has the advantages that: the method has higher positioning accuracy, and the positioning fault-tolerant capability is greatly improved compared with a section positioning method only using remote signaling data.

Description

Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data
Technical Field
The invention relates to the technical field of distribution network fault positioning, isolation and power supply recovery methods, in particular to the technical field of FTU-based distribution network section positioning methods.
Background
With the continuous improvement of the degree of distribution automation, the distribution network section positioning method based on the FTU becomes a research hotspot. The existing FTU-based section positioning method is roughly divided into two types, one is a matrix algorithm, and the other is an artificial intelligence algorithm. The section positioning method based on the intelligent optimization algorithm is a research focus due to the good fault-tolerant capability.
According to different improvement targets, the power distribution network section positioning method based on the intelligent optimization algorithm at the present stage can be roughly divided into three types. The first type is that a single-layer positioning model is constructed on a positioning model only according to remote signaling data; on the positioning algorithm, a single intelligent algorithm or an improved algorithm thereof or an algorithm of a plurality of populations evolving in parallel is adopted. The method can improve the fault-tolerant capability and the positioning accuracy of the section positioning on a certain level, but the fault-tolerant capability of the remote signaling data is limited when multipoint distortion occurs, and a positioning model needs to be reconstructed when the topological structure of the power distribution network changes. The second method is to introduce the idea of partitioning or layering into fault location, and reduce the dimension of the intelligent algorithm through region partitioning, thereby improving the search efficiency of fault location to a certain extent. But the positioning process only considers remote signaling data, and the fault tolerance capability and the positioning model have weak self-adaption capability to the topological structure. The third method is to use an algebraic method to replace a logic method to construct a switching function, use a complementary constraint condition to relax a discrete variable into a continuous variable, and use the method to provide a fault location model smoothing solving algorithm based on disturbance factors. The method has good numerical stability, particularly greatly improves the positioning speed, but is only suitable for single faults, only takes remote signaling data into consideration in positioning, and the self-adaption capability of a positioning model to a topological structure is weak.
In summary, these methods, although some achievements are achieved on different improvement targets, still have the following disadvantages:
1) the method only utilizes the remote signaling data of the FTU, has less utilization of the remote signaling data, and can generally accurately position the result when the remote signaling data has a small amount of distortion in the communication process, but the probability of remote signaling error codes is higher than that of remote signaling error codes in the communication error code range, and the positioning accuracy is greatly reduced when information distortion occurs at multiple points; 2) the fitness function and the switching function constructed in the positioning model are closely related to the topological structure of the power distribution network, and when the topological structure of the power distribution network changes or the section positioning is carried out on a plurality of different power distribution networks, the fitness function and the switching function cannot be updated or transplanted in a self-adaptive mode.
Disclosure of Invention
In order to solve the problems, the invention provides a self-adaptive layered section positioning method which can improve the positioning fault tolerance by using the auxiliary criterion of section positioning formed by telemetering data when the section positioning is judged wrongly under the high telesignaling distortion rate, and can adaptively update or transplant the fitness function and the switching function when the topological structure of the power distribution network changes or the section positioning is carried out aiming at a plurality of different power distribution networks, thereby avoiding the complex and fussy modeling process by using FTU telesignaling and telemetering data.
The technical solution adopted by the invention to solve the technical problems is as follows:
an adaptive layered sector location method using FTU telemetry and telemetry data, comprising the steps of:
step S1: before a fault occurs, separating a plurality of two ports of the whole power distribution network by taking a multi-branch node as a boundary, and obtaining relevant data of each two port, wherein the two ports comprise: a port node and a section that a port line contains; all the two ports form a first layer positioning model, and each two port forms a second layer positioning model;
step S2: before a fault occurs, obtaining incidence matrix data about two ports according to a topological structure of a power distribution network, wherein the incidence matrix data comprises the following steps: port node-port node incidence matrix AFPort node-port segment association matrix BFPort road-port segment matrix of main power supply
Figure GDA0003346061120000031
Port road-port section matrix of distributed power supply
Figure GDA0003346061120000032
Meanwhile, according to the internal topological structure of the two ports, the method obtainsAnd the incidence matrix data inside the two ports comprises: port internal node-port internal node incidence matrix ASPort internal node-port internal segment association matrix BSRoad-port internal segment matrix of main power supply
Figure GDA0003346061120000033
Road-port segment matrix for distributed power
Figure GDA0003346061120000034
Step S3: after the fault occurs, reading remote signaling data of the port node acquired by the FTU, and then according to the incidence matrix data A of the two portsF、BF、
Figure GDA0003346061120000035
Adaptively establishing a switch function and a fitness function of a first-layer positioning model, then solving the positioning model by using a Binary Particle Swarm Optimization (BPSO) to position a failed two-port;
step S4: reading the remote signaling data of the common node in the fault port acquired by the FTU, and then according to the incidence matrix data A in the two portsS、BS、
Figure GDA0003346061120000036
Adaptively establishing a switching function and a fitness function of a second layer positioning model; then, solving the second layer positioning model by using an exhaustion method in 0-1 planning, and positioning a fault section;
step S5: reading telemetering data of common nodes in a fault port acquired by FTU (fiber to the Unit), and utilizing the provided auxiliary criterion and port incidence matrix
Figure GDA0003346061120000037
Adaptively locating a fault section in a second layer location model;
step S6: if the two positioning results in the second layer positioning model are consistent, outputting the positioning result, otherwise, utilizing the zone positioning result obtained by the auxiliary criterion based on the telemetering data, simultaneously checking the zone positioning result based on the telecommand data in the second layer positioning model and the port positioning result based on the telecommand data in the first layer positioning model, and sending an alarm to the position where the telecommand data is distorted.
In the step S2, the node-node correlation matrix (A)F、AS) Node-segment association matrix (B)F、BS) Road-segment association matrix with power supply
Figure GDA0003346061120000038
The following definition is adopted for the elements in (1):
Figure GDA0003346061120000041
Figure GDA0003346061120000042
Figure GDA0003346061120000043
in step S3, the first layer positioning model adopts the following switching function:
Figure GDA0003346061120000044
in the formula Ii(s) represents a switching function of the switching device,
Figure GDA0003346061120000045
respectively representing the power supply from node i to upstream suNode i to downstream power supply sdThe state of the middle section, M 'and N' are the number of the upstream power supplies and the number of the downstream power supplies respectively; si,d、si,uRespectively representing the states of all sections from node i to downstream and from node i to upstream, M, N being the number of all sections upstream and downstreamThe number of segments; Π denotes the logical OR, ". denotes the numerical multiplication, Ku、KdRespectively representing the power supply coefficients of an upstream power supply and a downstream power supply, wherein the power supply access is 1, and the power supply exit is 0;
in the function of the switch-on/off,
Figure GDA0003346061120000046
the self-adaptive construction process comprises the following steps:
1) at M' number
Figure GDA0003346061120000047
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming M' sets [ i, s ]u](ii) a 2) For M's [ i, s ]u]Inner zone status
Figure GDA0003346061120000048
Respectively taking logical OR and then inverting to obtain M' s
Figure GDA0003346061120000049
3) According to the turn-on condition K of M' power suppliesuM' pieces of
Figure GDA00033460611200000410
Taking a logical OR;
in the function of the switch-on/off,
Figure GDA00033460611200000411
the self-adaptive construction process comprises the following steps:
1) at the number of N
Figure GDA00033460611200000412
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming N' sets [ i, s ]d](ii) a 2) For N' pieces of [ i, s ]d]Inner zone status
Figure GDA00033460611200000413
Respectively taking logical OR and then inverting to obtain N' ones
Figure GDA0003346061120000051
3) According to the turn-on condition K of N' power suppliesdN' pieces of
Figure GDA0003346061120000052
Taking a logical OR;
in the function of the switch-on/off,
Figure GDA0003346061120000053
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in AFSearching for an element having a value of 1 in the ith row of (1); 2) finding upstream section s of node ikDownstream section s of node jlAnd combining the segments sk、slLoad set [ i, u ]]: in BFSearching for an element having a value of-1 in the ith row of B, and searching for an element having a value of 1 in the jth row of B; 3) find section skUpstream node m, segment slA downstream node n: in BFSearch for an element with a value of 1 in the k-th column of (1), and search for an element with a value of 1 in BFSearching for an element having a value of-1 in the l column; 4) respectively making i ═ m and j ═ n, returning to 1), repeating the above three steps until the set [ i, u ═ n]No longer changed; 5) all [ i, u ]]All sectors within a sector are logically ORed;
in the function of the switch-on/off,
Figure GDA0003346061120000054
the self-adaptive construction process comprises the following steps:
1) finding the downstream segment s of node iiAnd put in the set [ i, d ]]: in BFSearch for the element with the value of 1 in the ith row of (1) and store the position of the column in which the element is located in [ i, d ]](ii) a 2) Find section siA downstream node q: in BFSearching for an element with a value of-1 in the ith column and recording the position q of the row where the element is located; 3) let i ═ q, return to 1), repeat the above two steps until the set [ i, d ═ q-]No longer changed; 4) all of [ i, d]All sectors within a sector are logically ORed;
the first layer positioning model adopts the following fitness function:
Figure GDA0003346061120000055
wherein f (X) represents the fitness value of particle X; in the first layer positioning model, D represents the number of equivalent two ports; the number of nodes of the whole power distribution network is T; i isjIs the real value of the node state;
Figure GDA0003346061120000056
is the expected value of the node state calculated by the switching function; siRepresenting the state of the port segment in a first layer localization model; eta is a weight coefficient.
In step S4, the second layer positioning model adopts the following switching function:
Figure GDA0003346061120000061
in the function of the switch-on/off,
Figure GDA0003346061120000062
the self-adaptive construction process comprises the following steps:
1) in that
Figure GDA0003346061120000063
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming the set i, su](ii) a 2) To [ i, s ]u]Inner zone status
Figure GDA0003346061120000064
Respectively taking logical OR and inverting
Figure GDA0003346061120000065
3) According to the power-on condition KuTo obtain
Figure GDA0003346061120000066
In the function of the switch-on/off,
Figure GDA0003346061120000067
the self-adaptive construction process comprises the following steps:
1) in that
Figure GDA0003346061120000068
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming the set i, sd](ii) a 2) To [ i, s ]d]Inner zone status
Figure GDA0003346061120000069
Respectively taking logical OR and inverting
Figure GDA00033460611200000610
3) According to the power-on condition KdTo obtain
Figure GDA00033460611200000611
In the function of the switch-on/off,
Figure GDA00033460611200000612
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in ASSearching for an element having a value of 1 in the ith row of (1); 2) find the upstream segment s of node jl: in BSSearch for an element with a value of 1 in line j and combine slLoad set [ i, u ]](ii) a 3) Repeating the above two steps until l ═ M; all [ i, u ]]All sectors within a sector are logically ORed;
in the function of the switch-on/off,
Figure GDA00033460611200000613
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in ASSearching for an element having a value of 1 in the ith row of (1); 2) find the downstream segment s of node jq: in BSSearch for an element with a value of-1 in line j and combine sqLoad set [ i, d ]](ii) a 3) Repeat the aboveTwo steps, until d ═ N; all of [ i, d]All sectors within a sector are logically ORed;
the fitness function adopted by the second layer positioning model is the same as that adopted by the first layer positioning model, and the difference is that: in the second layer of the localization model, D represents the number of nodes inside the port, siIndicating the sector state.
In step S5, the positioning assistance criterion constructed based on the telemetry data is:
Figure GDA0003346061120000071
wherein, IiIndicating the detected current at the node upstream of the faulty section, Ii+1Indicating the current detected at the node downstream of the faulty section, KsetIs a set threshold value;
the self-adaptive fault discrimination process comprises the following steps: according to
Figure GDA0003346061120000072
The current ratios of adjacent nodes are calculated from upstream to downstream in sequence, and the first section meeting the threshold is determined as a fault section.
Compared with the prior art, the invention has the following 2 outstanding advantages:
1. when a large amount of node information is distorted in a power distribution network fault acquisition system, the method still has high positioning accuracy, and compared with a section positioning method only using remote signaling data, the positioning fault tolerance of the method is greatly improved;
2. by analyzing the construction characteristics of the fitness function and the switch function, the self-adaptive modeling of section positioning is carried out based on the node-node incidence matrix A, the node-section incidence matrix B and the road-section incidence matrix T of the power supply, when the topological structure of the power distribution network changes or the section positioning is carried out aiming at a plurality of different power distribution networks, the fitness function and the switch function can be updated or transplanted in a self-adaptive manner, and the complex and tedious modeling process is avoided.
Drawings
FIG. 1 is a diagram of a 14-node distribution network architecture of the present invention;
FIG. 2 is a flow chart of the entire sector location of the present invention;
FIG. 3 is a two-port partition of a 14-node distribution network of the present invention;
FIG. 4 is an iteration convergence curve of the BPSO iteration of the present invention;
FIG. 5 is a curve of the fitness value change in port four of the present invention;
fig. 6 is a curve of the variation of the fitness value in port seven of the present invention.
Detailed Description
The invention is further explained below with reference to the drawings and the embodiments.
Taking a 14-node power distribution network model as an embodiment of the invention, as shown in fig. 1, the model comprises 16 nodes and 14 sections, the nodes are represented by 1-16, the sections are represented by (1) - (14), the middle node 1 is a circuit breaker, the nodes 3 and 10 are tie switches, the rest nodes are section switches, distributed power supplies DG1 and DG2 are connected to the ends of the sections 6 and 14, the transformer capacity is 20MVA, the capacities of DG1 and DG2 are set to be 0.3MVA, and the flow of positioning the whole section is shown in fig. 2.
An adaptive layered sector location method using FTU telemetry and telemetry data, comprising the steps of:
step S1: before a fault occurs, separating a plurality of two ports of the whole power distribution network by taking a multi-branch node as a boundary, and obtaining relevant data of each two port, wherein the two ports comprise: a port node and a section that a port line contains; all the two ports form a first layer positioning model, and each two port forms a second layer positioning model;
in this embodiment, the multi-branch node is taken as a boundary method, and the whole power distribution network can be divided into seven two-port R1-R7As shown in fig. 3, the section and port node included in two ports are shown in table 1.
TABLE 1 nodes and segments encompassed by two ports
Figure GDA0003346061120000081
Step S2: before a fault occurs, obtaining incidence matrix data about two ports according to a topological structure of a power distribution network, wherein the incidence matrix data comprises the following steps: port node-port node incidence matrix AFPort node-port segment association matrix BFPort road-port segment matrix of main power supply
Figure GDA0003346061120000091
Port road-port section matrix of distributed power supply
Figure GDA0003346061120000092
Meanwhile, obtaining incidence matrix data inside the two ports according to the topology structure inside the two ports, comprising the following steps: port internal node-port internal node incidence matrix ASPort internal node-port internal segment association matrix BSRoad-port internal segment matrix of main power supply
Figure GDA0003346061120000093
Road-port segment matrix for distributed power
Figure GDA0003346061120000094
Node-node incidence matrix (A)F、AS) Node-segment association matrix (B)F、BS) Road-segment association matrix with power supply
Figure GDA0003346061120000095
The following definition is adopted for the elements in (1):
Figure GDA0003346061120000096
Figure GDA0003346061120000097
Figure GDA0003346061120000098
in the embodiment, before the fault occurs, the incidence matrix data A about the two ports is obtained according to the topological structure of the power distribution networkF、BF、
Figure GDA0003346061120000099
Comprises the following steps:
Figure GDA00033460611200000910
Figure GDA0003346061120000101
Figure GDA0003346061120000102
incidence matrix data A inside portS、BS、
Figure GDA0003346061120000103
(only data for port four and port seven are given) are:
Figure GDA0003346061120000104
Figure GDA0003346061120000105
Figure GDA0003346061120000106
step S3: after the fault occurs, reading remote signaling data of the port node acquired by the FTU, and then according to the incidence matrix data A of the two portsF、BF、
Figure GDA0003346061120000107
Adaptively establishing a switch function and a fitness function of a first-layer positioning model, then solving the positioning model by using a Binary Particle Swarm Optimization (BPSO) to position a failed two-port;
the first layer positioning model adopts the following switching functions:
Figure GDA0003346061120000111
in the formula Ii(s) represents a switching function of the switching device,
Figure GDA0003346061120000112
respectively representing the power supply from node i to upstream suNode i to downstream power supply sdThe state of the middle section, M 'and N' are the number of the upstream power supplies and the number of the downstream power supplies respectively; si,d、si,uRespectively representing the states of all sections from node i to downstream and from node i to upstream, wherein M, N is the number of all sections at upstream and the number of all sections at downstream; Π denotes the logical OR, ". denotes the numerical multiplication, Ku、KdRespectively representing the power supply coefficients of an upstream power supply and a downstream power supply, wherein the power supply access is 1, and the power supply exit is 0;
in the function of the switch-on/off,
Figure GDA0003346061120000113
the self-adaptive construction process comprises the following steps:
1) at M' number
Figure GDA0003346061120000114
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming M' sets [ i, s ]u](ii) a 2) For M's [ i, s ]u]Inner zone status
Figure GDA0003346061120000115
Respectively taking logical OR and then inverting to obtain M' s
Figure GDA0003346061120000116
3) According to the turn-on condition K of M' power suppliesuM' pieces of
Figure GDA0003346061120000117
Taking a logical OR;
in the function of the switch-on/off,
Figure GDA0003346061120000118
the self-adaptive construction process comprises the following steps:
1) at the number of N
Figure GDA0003346061120000119
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming N' sets [ i, s ]d](ii) a 2) For N' pieces of [ i, s ]d]Inner zone status
Figure GDA00033460611200001110
Respectively taking logical OR and then inverting to obtain N' ones
Figure GDA00033460611200001111
3) According to the turn-on condition K of N' power suppliesdN' pieces of
Figure GDA00033460611200001112
Taking a logical OR;
in the function of the switch-on/off,
Figure GDA00033460611200001113
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in AFSearching for an element having a value of 1 in the ith row of (1); 2) finding upstream section s of node ikDownstream section s of node jlAnd combining the segments sk、slLoad set [ i, u ]]: in BFSearching for an element having a value of-1 in the ith row of B, and searching for an element having a value of 1 in the jth row of B; 3) find section skUpstream node m, segment slDownstreamAnd (3) node n: in BFSearch for an element with a value of 1 in the k-th column of (1), and search for an element with a value of 1 in BFSearching for an element having a value of-1 in the l column; 4) respectively making i ═ m and j ═ n, returning to 1), repeating the above three steps until the set [ i, u ═ n]No longer changed; 5) all [ i, u ]]All sectors within a sector are logically ORed;
in the function of the switch-on/off,
Figure GDA0003346061120000121
the self-adaptive construction process comprises the following steps:
1) finding the downstream segment s of node iiAnd put in the set [ i, d ]]: in BFSearch for the element with the value of 1 in the ith row of (1) and store the position of the column in which the element is located in [ i, d ]](ii) a 2) Find section siA downstream node q: in BFSearching for an element with a value of-1 in the ith column and recording the position q of the row where the element is located; 3) let i ═ q, return to 1), repeat the above two steps until the set [ i, d ═ q-]No longer changed; 4) all of [ i, d]All sectors within a sector are logically ORed;
the first layer positioning model adopts the following fitness function:
Figure GDA0003346061120000122
wherein f (X) represents the fitness value of particle X; in the first layer positioning model, D represents the number of equivalent two ports; the number of nodes of the whole power distribution network is T; i isjIs the real value of the node state;
Figure GDA0003346061120000123
is the expected value of the node state calculated by the switching function; siRepresenting the state of the port segment in a first layer localization model; η is a weight coefficient, and is usually set to 0.5.
In the positioning of the section of the power distribution network, the fitness function reflects the accumulation of the difference between the actual acquisition value and the expected value of each node, and only has a large relation with the number of the nodes but has a small relation with the topological structure of the power distribution network, so that the self-adaptive construction process of the first-layer positioning model fitness function is the self-adaptive statistical nodeThe number of the processes is as follows: computing a port node-port node incidence matrix AFAnd assigning the value of the matrix to D.
In this embodiment, when all DG are connected, a three-phase metallic short circuit occurs in the segment 12 on the dual-source branch, and the remote signaling data of the port nodes 1, 2, 4, 5, 7, 9, and 11 collected by the FTU are: [ 10111-11 ], it is clear that false positives occurred for port nodes 5, 9; adaptively constructing a first-layer positioning model by using incidence matrix data about two ports, and performing optimization iteration by using the BPSO, wherein an iteration convergence curve of the iteration of the BPSO is shown in FIG. 4, and it can be found from the figure that immature convergence occurs at the moment, and the corresponding result is as follows: [ 0001001 ]; thus, it is determined that the ports four and seven have failed, and it is apparent that the port four is misjudged.
Step S4: reading the remote signaling data of the common node in the fault port acquired by the FTU, and then according to the incidence matrix data A in the two portsS、BS、
Figure GDA0003346061120000131
Adaptively establishing a switching function and a fitness function of a second layer positioning model; then, solving the second layer positioning model by using an exhaustion method in 0-1 planning, and positioning a fault section;
the second layer positioning model adopts the following switching functions:
Figure GDA0003346061120000132
in the function of the switch-on/off,
Figure GDA0003346061120000133
the self-adaptive construction process comprises the following steps:
1) in that
Figure GDA0003346061120000134
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming the set i, su](ii) a 2) To [ i, s ]u]Inner zone status
Figure GDA0003346061120000135
Respectively taking logical OR and inverting
Figure GDA0003346061120000136
3) According to the power-on condition KuTo obtain
Figure GDA0003346061120000137
In the function of the switch-on/off,
Figure GDA0003346061120000138
the self-adaptive construction process comprises the following steps:
1) in that
Figure GDA0003346061120000141
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming the set i, sd](ii) a 2) To [ i, s ]d]Inner zone status
Figure GDA0003346061120000142
Respectively taking logical OR and inverting
Figure GDA0003346061120000143
3) According to the power-on condition KdTo obtain
Figure GDA0003346061120000144
In the function of the switch-on/off,
Figure GDA0003346061120000145
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in ASSearching for an element having a value of 1 in the ith row of (1); 2) find the upstream segment s of node jl: in BSSearch for an element with a value of 1 in line j and combine slLoad collections[i,u](ii) a 3) Repeating the above two steps until l ═ M; all [ i, u ]]All sectors within a sector are logically ORed;
in the function of the switch-on/off,
Figure GDA0003346061120000146
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in ASSearching for an element having a value of 1 in the ith row of (1); 2) find the downstream segment s of node jq: in BSSearch for an element with a value of-1 in line j and combine sqLoad set [ i, d ]](ii) a 3) Repeating the two steps until d is equal to N; all of [ i, d]All sectors within a sector take a logical or.
The fitness function adopted by the second layer positioning model is the same as that adopted by the first layer positioning model, and the difference is that: in the second layer of the localization model, D represents the number of nodes inside the port, siIndicating the sector state.
Because the fitness function adopted by the second layer positioning model is the same as the fitness function adopted by the first layer positioning model, and the only difference is that the number of nodes is different, the self-adaptive construction of the fitness function of the second layer positioning model is realized by the following steps: computing port internal node-port internal node incidence matrix ASAnd assigning the value of the matrix to D.
In this embodiment, according to the positioning result of the first-layer positioning model, the remote signaling data [1, -1, -1] of the four internal nodes 5, 6, 15 of the port are read, and the second-layer positioning model is adaptively constructed, and is solved by using an exhaustion method, the change curve of the fitness value is shown in fig. 5, the corresponding solving result is [1,0], and then the fault of the section (5) is determined, and it is obvious that the erroneous judgment occurs here.
Similarly, the remote signaling data [1,1, -1, -1, -1] of the seven internal nodes 11, 12, 13, 14, 16 of the port are read, the solution is carried out by using an exhaustion method, a change curve of the fitness value is shown in fig. 6, the solution result is [0100], and then the fault of the section (12) is determined.
Step S5: reading telemetering data of common nodes in fault ports acquired by FTU (fiber to the Unit), and utilizingProposed auxiliary criteria and port association matrix
Figure GDA0003346061120000151
Adaptively locating a fault section in a second layer location model;
the positioning auxiliary criterion constructed based on the telemetering data is as follows:
Figure GDA0003346061120000152
wherein, IiIndicating the detected current at the node upstream of the faulty section, Ii+1Indicating the current detected at the node downstream of the faulty section, KsetFor the set threshold, the invention is set to 5;
the self-adaptive fault discrimination process comprises the following steps: according to
Figure GDA0003346061120000153
Sequentially calculating the current ratio of adjacent nodes from upstream to downstream, and judging a first section meeting a threshold value as a fault section; it is noted that for end sectors where there are no FTU nodes downstream, the default is to set to a fractional number towards zero, the invention takes 10-5。
In this embodiment, the port four is started with the auxiliary criterion, and according to the positioning result of the first layer positioning model, the telemetry data of the port four internal nodes 5, 6, and 15 are read and are respectively 14.42A, 15.74A, and 17.31A, and the current ratio is sequentially calculated as: 0.91, 0.90, and no section fault is determined.
For port seven, according to the positioning result of the first layer positioning model, reading the remote signaling data of internal nodes 11, 12, 13, 14 and 16 of port seven, which are 183.45a, 171.15a, 23.01A, 30.34A and 38.21A respectively, and sequentially calculating the current ratio according to an auxiliary criterion as follows: 1.071, 7.441, 0.758, 0.794, and thus it is determined that section 12 is malfunctioning.
Step S6: if the two positioning results in the second layer positioning model are consistent, outputting the positioning result, otherwise, utilizing the zone positioning result obtained by the auxiliary criterion based on the telemetering data, simultaneously checking the zone positioning result based on the telecommand data in the second layer positioning model and the port positioning result based on the telecommand data in the first layer positioning model, and sending an alarm to the position where the telecommand data is distorted.
In this embodiment, for port seven, the two positioning results of the second layer positioning model are consistent and output. For port four, the two positioning results of the second layer positioning model are inconsistent, so that result verification is performed, and the verification result finds that: and the port nodes 5 and 9 generate false alarm, and alarm is given to the positions 5 and 9 of the two points.

Claims (5)

1. An adaptive layered sector location method using FTU telemetry and telemetry data, comprising the steps of:
step S1: before a fault occurs, separating a plurality of two ports of the whole power distribution network by taking a multi-branch node as a boundary, and obtaining relevant data of each two port, wherein the two ports comprise: a port node and a section that a port line contains; all the two ports form a first layer positioning model, and each two port forms a second layer positioning model;
step S2: before a fault occurs, obtaining incidence matrix data about two ports according to a topological structure of a power distribution network, wherein the incidence matrix data comprises the following steps: port node-port node incidence matrix AFPort node-port segment association matrix BFPort road-port segment matrix of main power supply
Figure FDA0003346061110000011
Port road-port section matrix of distributed power supply
Figure FDA0003346061110000012
Meanwhile, obtaining incidence matrix data inside the two ports according to the topology structure inside the two ports, comprising the following steps: port internal node-port internal node incidence matrix ASPort internal node-port internal segment association matrix BSRoad-port internal segment matrix of main power supply
Figure FDA0003346061110000013
Road-port segment matrix for distributed power
Figure FDA0003346061110000014
Step S3: after the fault occurs, reading remote signaling data of the port node acquired by the FTU, and then according to the incidence matrix data A of the two portsF、BF、
Figure FDA0003346061110000015
Adaptively establishing a switch function and a fitness function of a first-layer positioning model, then solving the positioning model by using a Binary Particle Swarm Optimization (BPSO) to position a failed two-port;
step S4: reading the remote signaling data of the common node in the fault port acquired by the FTU, and then according to the incidence matrix data A in the two portsS、BS、
Figure FDA0003346061110000016
Adaptively establishing a switching function and a fitness function of a second layer positioning model; then, solving the second layer positioning model by using an exhaustion method in 0-1 planning, and positioning a fault section;
step S5: reading telemetering data of common nodes in a fault port acquired by FTU (fiber to the Unit), and utilizing the provided auxiliary criterion and port incidence matrix
Figure FDA0003346061110000017
Adaptively locating a fault section in a second layer location model;
step S6: if the two positioning results in the second layer positioning model are consistent, outputting the positioning result, otherwise, utilizing the zone positioning result obtained by the auxiliary criterion based on the telemetering data, simultaneously checking the zone positioning result based on the telecommand data in the second layer positioning model and the port positioning result based on the telecommand data in the first layer positioning model, and sending an alarm to the position where the telecommand data is distorted.
2. The method of adaptive layered sector location using FTU telemetry and signaling data of claim 1, wherein: in the step S2, the node-node correlation matrix (A)F、AS) Node-segment association matrix (B)F、BS) Road-segment association matrix with power supply
Figure FDA0003346061110000021
The following definition is adopted for the elements in (1):
Figure FDA0003346061110000022
Figure FDA0003346061110000023
Figure FDA0003346061110000024
3. the method of adaptive layered sector location using FTU telemetry and signaling data of claim 2, wherein: in step S3, the first layer positioning model adopts the following switching function:
Figure FDA0003346061110000025
in the formula Ii(s) represents a switching function of the switching device,
Figure FDA0003346061110000026
respectively representing the power supply from node i to upstream suNode i to downstream power supply sdIn the middle section, M ', N' being upstream power supply respectivelyThe number and the number of downstream power supplies; si,d、si,uRespectively representing the states of all sections from node i to downstream and from node i to upstream, wherein M, N is the number of all sections at upstream and the number of all sections at downstream; Π denotes the logical OR, ". denotes the numerical multiplication, Ku、KdRespectively representing the power supply coefficients of an upstream power supply and a downstream power supply, wherein the power supply access is 1, and the power supply exit is 0;
in the function of the switch-on/off,
Figure FDA0003346061110000031
the self-adaptive construction process comprises the following steps:
1) at M' number
Figure FDA0003346061110000032
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming M' sets [ i, s ]u](ii) a 2) For M's [ i, s ]u]Inner zone status
Figure FDA0003346061110000033
Respectively taking logical OR and then inverting to obtain M' s
Figure FDA0003346061110000034
3) According to the turn-on condition K of M' power suppliesuM' pieces of
Figure FDA0003346061110000035
Taking a logical OR;
in the function of the switch-on/off,
Figure FDA0003346061110000036
the self-adaptive construction process comprises the following steps:
1) at the number of N
Figure FDA0003346061110000037
Searching for an element with a value of 1 in the ith row of (2) and recording the position of the column in which the element is located, wherein all the positions form NSet [ i, s ]d](ii) a 2) For N' pieces of [ i, s ]d]Inner zone status
Figure FDA0003346061110000038
Respectively taking logical OR and then inverting to obtain N' ones
Figure FDA0003346061110000039
3) According to the turn-on condition K of N' power suppliesdN' pieces of
Figure FDA00033460611100000310
Taking a logical OR;
in the function of the switch-on/off,
Figure FDA00033460611100000311
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in AFSearching for an element having a value of 1 in the ith row of (1); 2) finding upstream section s of node ikDownstream section s of node jlAnd combining the segments sk、slLoad set [ i, u ]]: in BFSearching for an element having a value of-1 in the ith row of B, and searching for an element having a value of 1 in the jth row of B; 3) find section skUpstream node m, segment slA downstream node n: in BFSearch for an element with a value of 1 in the k-th column of (1), and search for an element with a value of 1 in BFSearching for an element having a value of-1 in the l column; 4) respectively making i ═ m and j ═ n, returning to 1), repeating the above three steps until the set [ i, u ═ n]No longer changed; 5) all [ i, u ]]All sectors within a sector are logically ORed;
in the function of the switch-on/off,
Figure FDA00033460611100000312
the self-adaptive construction process comprises the following steps:
1) finding the downstream segment s of node iiAnd put in the set [ i, d ]]: in BFSearch for the element with the value of 1 in the ith row of (1) and store the position of the column in which the element is located in [ i, d ]](ii) a 2) Find section siA downstream node q: in BFSearching for an element with a value of-1 in the ith column and recording the position q of the row where the element is located; 3) let i ═ q, return to 1), repeat the above two steps until the set [ i, d ═ q-]No longer changed; 4) all of [ i, d]All sectors within a sector are logically ORed;
the first layer positioning model adopts the following fitness function:
Figure FDA0003346061110000041
wherein f (X) represents the fitness value of particle X; in the first layer positioning model, D represents the number of equivalent two ports; the number of nodes of the whole power distribution network is T; i isjIs the real value of the node state;
Figure FDA0003346061110000042
is the expected value of the node state calculated by the switching function; siRepresenting the state of the port segment in a first layer localization model; eta is a weight coefficient.
4. The method of claim 3 for adaptive layered sector location using FTU telemetry and telemetry data, wherein: in step S4, the second layer positioning model adopts the following switching function:
Figure FDA0003346061110000043
in the function of the switch-on/off,
Figure FDA0003346061110000044
the self-adaptive construction process comprises the following steps:
1) in that
Figure FDA0003346061110000045
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming the set i, su](ii) a 2) To [ i, s ]u]Inner zoneSegment state
Figure FDA0003346061110000046
Respectively taking logical OR and inverting
Figure FDA0003346061110000047
3) According to the power-on condition KuTo obtain
Figure FDA0003346061110000048
In the function of the switch-on/off,
Figure FDA0003346061110000049
the self-adaptive construction process comprises the following steps:
1) in that
Figure FDA0003346061110000051
Search for the element with value 1 in the ith row of (2) and record the position of the column in which it is located, all positions forming the set i, sd](ii) a 2) To [ i, s ]d]Inner zone status
Figure FDA0003346061110000052
Respectively taking logical OR and inverting
Figure FDA0003346061110000053
3) According to the power-on condition KdTo obtain
Figure FDA0003346061110000054
In the function of the switch-on/off,
Figure FDA0003346061110000055
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in ASSearching for an element having a value of 1 in the ith row of (1); 2) find the upstream segment s of node jl: in BSSearch for an element with a value of 1 in line j and combine slLoad set [ i, u ]](ii) a 3) Repeating the above two steps until l ═ M; all [ i, u ]]All sectors within a sector are logically ORed;
in the function of the switch-on/off,
Figure FDA0003346061110000056
the self-adaptive construction process comprises the following steps:
1) find the neighbor node j of node i: in ASSearching for an element having a value of 1 in the ith row of (1); 2) find the downstream segment s of node jq: in BSSearch for an element with a value of-1 in line j and combine sqLoad set [ i, d ]](ii) a 3) Repeating the two steps until d is equal to N; all of [ i, d]All sectors within a sector are logically ORed;
the fitness function adopted by the second layer positioning model is the same as that adopted by the first layer positioning model, and the difference is that: in the second layer of the localization model, D represents the number of nodes inside the port, siIndicating the sector state.
5. The method of claim 4 for adaptive layered sector location using FTU telemetry and telemetry data, wherein: in step S5, the positioning assistance criterion constructed based on the telemetry data is:
Figure FDA0003346061110000057
wherein, IiIndicating the detected current at the node upstream of the faulty section, Ii+1Indicating the current detected at the node downstream of the faulty section, KsetIs a set threshold value;
the self-adaptive fault discrimination process comprises the following steps: according to
Figure FDA0003346061110000061
The current ratios of adjacent nodes are calculated from upstream to downstream in sequence, and the first section meeting the threshold is determined as a fault section.
CN201910966007.1A 2019-10-12 2019-10-12 Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data Active CN111141992B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910966007.1A CN111141992B (en) 2019-10-12 2019-10-12 Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910966007.1A CN111141992B (en) 2019-10-12 2019-10-12 Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data

Publications (2)

Publication Number Publication Date
CN111141992A CN111141992A (en) 2020-05-12
CN111141992B true CN111141992B (en) 2022-01-18

Family

ID=70516842

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910966007.1A Active CN111141992B (en) 2019-10-12 2019-10-12 Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data

Country Status (1)

Country Link
CN (1) CN111141992B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112034301A (en) * 2020-05-13 2020-12-04 江苏方天电力技术有限公司 Power distribution network fault section positioning method based on circuit breaker
CN111551825B (en) * 2020-05-28 2023-04-18 中国矿业大学(北京) Self-adaptive power distribution network fault positioning method based on fault current path
CN112485587B (en) * 2020-11-11 2024-04-19 国网福建省电力有限公司宁德供电公司 Layered positioning method for fault section of distribution network containing distributed photovoltaic
CN113191062B (en) * 2021-04-13 2024-01-09 云南电网有限责任公司昆明供电局 Power distribution network fault section positioning method and system based on multisource incomplete information

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108037414A (en) * 2017-12-11 2018-05-15 福州大学 A kind of electrical power distribution network fault location method based on hierarchical mode and intelligent checking algorithm
CN109635411A (en) * 2018-12-06 2019-04-16 湖北鄂电德力电气有限公司 A kind of distribution network failure Hierarchical Location method counted and FTU is failed to report and reported by mistake

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7342878B2 (en) * 2002-12-23 2008-03-11 International Business Machines Corporation Input port routing circuit that performs output port filtering
US10534027B2 (en) * 2017-01-17 2020-01-14 Fluke Corporation Phase coherent main and remote units of a network analyzer

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108037414A (en) * 2017-12-11 2018-05-15 福州大学 A kind of electrical power distribution network fault location method based on hierarchical mode and intelligent checking algorithm
CN109635411A (en) * 2018-12-06 2019-04-16 湖北鄂电德力电气有限公司 A kind of distribution network failure Hierarchical Location method counted and FTU is failed to report and reported by mistake

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"An_accurate_fault_location_method_of_smart_distribution_network";Tan Zhihai等;《2014 China International Conference on Electricity Distribution》;20140926;第916-920页 *
"基于0_1整数规划的配电网区段定位方法";刘严等;《电网与清洁能源》;20181031;第06-11页 *
"基于分层模型和智能校验算法的配电网故障定位技术";王秋杰等;《电工技术学报》;20181130;第5327-5337页 *

Also Published As

Publication number Publication date
CN111141992A (en) 2020-05-12

Similar Documents

Publication Publication Date Title
CN111141992B (en) Self-adaptive layered section positioning method utilizing FTU remote signaling and remote measuring data
CN110994612A (en) Power distribution network fault rapid recovery method based on network topology partition layering
CN106229964A (en) A kind of based on the electrical power distribution network fault location method improving binary particle swarm algorithm
CN106532689B (en) Power distribution network topological structure optimization method and system
CN111551825B (en) Self-adaptive power distribution network fault positioning method based on fault current path
CN103902775A (en) Multilayer obstacle-avoiding Steiner minimal tree construction method for very large scale integration
CN104578427A (en) Fault self-healing method for power distribution network containing microgrid power source
CN109470997A (en) A kind of distribution network segment positioning method using multifactor dimensionality reduction
CN116223973A (en) Distributed power distribution network fault positioning method based on improved gray wolf optimization algorithm
CN113468745B (en) Method and system for rapidly evaluating reliability of power distribution network based on historical faults
CN112103950B (en) Power grid partitioning method based on improved GN splitting algorithm
CN110687397A (en) Active power distribution network fault positioning method based on improved artificial fish swarm algorithm
CN112014687A (en) Layered positioning method for fault section of active power distribution network containing distributed power supply
CN109345155B (en) Model-based hierarchical diagnosis method for power distribution network faults
CN109033603B (en) Intelligent substation secondary system simulation method based on source flow path chain
Srivastava et al. Parallel self-organising hierarchical neural network-based fast voltage estimation
CN113507116B (en) Power distribution network load transfer method, device, equipment and storage medium
CN110703725B (en) Path optimization method suitable for aerospace attitude orbit control system
CN110197305B (en) Relay protection data model searching and optimizing method and system based on shortest path algorithm
CN111369052B (en) Simplified road network KSP optimization algorithm
CN113009274B (en) Power distribution network fault section positioning method and system based on IELM algorithm
CN106227696B (en) Method for rapidly reconstructing high-performance target array
CN112462187A (en) Power distribution network fault layered positioning method and device considering FTU (fiber to the Unit) missing report and false report
CN109713665A (en) A kind of minimal hitting set algorithm suitable for the multiple multiphase failure of power distribution network
Zhao et al. Fault section location for distribution network containing DG based on IBQPSO

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