CN111697579B - Method, system and medium for determining key power transmission section of power system - Google Patents
Method, system and medium for determining key power transmission section of power system Download PDFInfo
- Publication number
- CN111697579B CN111697579B CN202010611996.5A CN202010611996A CN111697579B CN 111697579 B CN111697579 B CN 111697579B CN 202010611996 A CN202010611996 A CN 202010611996A CN 111697579 B CN111697579 B CN 111697579B
- Authority
- CN
- China
- Prior art keywords
- power transmission
- node
- power
- nodes
- level
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
Abstract
The invention discloses a method, a system and a medium for determining a key power transmission section of a power system, wherein the method for determining the key power transmission section of the power system comprises the following steps: inputting a topology structure diagram G of a power grid; partitioning a power grid based on the topological structure diagram G to obtain a power transmission section between partitions; and screening out key power transmission sections from all the power transmission sections, wherein an input topological structure diagram G of the power grid is formed by a node set V and an edge set E formed by edges between the nodes, each node in the node set V represents one power transmission node of the power grid, and each edge in the edge set E represents a power transmission line between two power transmission nodes of the power grid. The method can adapt to the frequent change condition of the power system, can quickly, accurately and automatically generate the power transmission section and the key power transmission section thereof, can meet the quick and accurate safety analysis requirements of the power grid, and has the advantages of high detection accuracy and high detection speed.
Description
Technical Field
The invention relates to a power grid management scheduling technology of a power system, in particular to a method, a system and a medium for determining a key power transmission section of the power system.
Background
With the continuous enlargement of the scale of the power system, the scale of the power grid is more and more huge, and the structure is more and more complex. Therefore, a new challenge is also presented for grid management scheduling of power systems, because the grid management scheduling system of a single power system has been difficult to implement to monitor and analyze all nodes and elements in the grid. In order to solve the problems, an effective method is to perform partition management on the power grid, so that the large-scale power grid is divided into a plurality of independent partitions, and independent management, scheduling and analysis are performed by taking the partitions as units, thereby solving the problem of expansibility of the power grid.
The power grid is divided into a plurality of subareas, a plurality of power transmission sections are formed, and each subarea can be divided into subsystems which can be independently managed, dispatched and analyzed through the power transmission sections. Therefore, how to scientifically determine the transmission section has great significance for the management, the dispatching and the analysis of the power grid. The current method for determining the power transmission section generally comprises the steps of dividing according to regions and then screening out key power transmission sections according to expert experience. However, the method does not consider the actual condition of the power grid across regions, so that the division of the power transmission section is unreasonable, the determination of the key power transmission section is further influenced, the method cannot adapt to the frequent change condition of the power system, the missing selection of the power transmission section and the key power transmission section is easily caused, the requirement of fast and accurate power grid safety analysis cannot be met, and the development of power grid scheduling and operation work is restricted.
Disclosure of Invention
The technical problems to be solved by the invention are as follows: in view of the above problems in the prior art, a method, system and medium for determining a critical power transmission section of a power system are provided. The method can adapt to the frequent change condition of the power system, can quickly, accurately and automatically generate the power transmission section and the key power transmission section thereof, can meet the requirements of quick and accurate power grid safety analysis, and has the advantages of high detection accuracy and high detection speed.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a method for determining a key power transmission section of a power system comprises the following implementation steps:
1) Inputting a topology structure diagram G of the power grid;
2) Partitioning a power grid based on the topological structure diagram G to obtain a power transmission section between partitions;
3) And screening out key power transmission sections from all power transmission sections.
Optionally, the topology structure diagram G of the power grid input in step 1) is formed by an edge set E formed by a node set V and edges between nodes, where each node in the node set V represents one power transmission node of the power grid, and each edge in the edge set E represents a power transmission line between two power transmission nodes of the power grid.
Optionally, the power grid is a power grid including multiple voltage classes, and a side of the edge set E has power transmission line voltage information, and the detailed step of step 2) includes:
2.1A) acquiring voltage information of all the transmission lines in the edge set E;
2.2A) selecting a power transmission node connected with the edge with the highest voltage grade as a first-grade power transmission node to obtain a first-grade power transmission node set;
2.3A) aiming at each primary power transmission node in the primary power transmission node set, respectively establishing a corresponding process or thread and suspending the process or thread so that each primary power transmission node corresponds to an independent process or thread, and initializing the tree structure of each primary power transmission node to only contain the primary power transmission node as a root node;
2.4A) initializing the level variable i to 1;
2.5A) controlling the synchronous start of the process or the thread of each first-level power transmission node, gradually searching the ith-level power transmission node which has a direct or indirect connection relation with the current first-level power transmission node through the process or the thread traversing edge set E, if the ith-level power transmission node is found and is not additionally provided with a use mark by other first-level power transmission nodes in the node set V, attaching the ith-level power transmission node to the ith layer of the tree structure of the current first-level power transmission node according to the hierarchy of the ith-level power transmission node, and enabling a father node of the ith-level power transmission node to be the current first-level power transmission node or a previous-level power transmission node which is connected with the current first-level power transmission node and additionally providing the use mark for the ith-level power transmission node in the node set V; suspending the process or the thread;
2.6A) judging whether power transmission nodes which are not additionally marked by using still exist in the node set V, if so, adding 1 to the hierarchy variable i, and skipping to execute the step 2.5A); otherwise, taking the obtained tree structure of each primary power transmission node as a partition;
2.7A) searching leaf nodes of the tree structure of each first-level power transmission node, and forming a leaf node pair by one leaf node of the tree structure of any first-level power transmission node and one leaf node of the tree structure of different first-level power transmission nodes to obtain a leaf node pair set;
2.8A) selecting a set of current leaf node pairs from the set of leaf node pairs;
2.9A) judging whether the matching of one edge of two power transmission nodes corresponding to the current leaf node in the edge set E is successful, and if so, judging that the current leaf node pair is a power transmission section between the two corresponding subareas;
2.10A) judging whether the leaf node pair set is selected in a traversing way, if not, skipping to execute the step 2.8A); otherwise, judging that all the power transmission sections in the power grid are determined to be finished, and skipping to execute the step 3).
Optionally, the power grid is a single-voltage-class power grid, and the nodes in the node set V have geographical location information of the power transmission nodes, and the detailed step of step 2) includes:
2.1B) determining the number of edges connected with each power transmission node in the node set V to obtain an edge number binary array consisting of various edge numbers and the corresponding node numbers;
2.2B) sorting the edge number binary arrays in a descending order according to the edge numbers;
2.3B) taking out an element from the head end of the edge number binary array, judging whether the ratio of the number of the nodes of all the taken-out elements to the number of all the nodes in the node set V exceeds a preset threshold value or not, and skipping to execute the next step if the ratio exceeds the preset threshold value; otherwise, skipping to execute the step 2.3B again);
2.4B) taking the power transmission nodes corresponding to all the taken-out elements as first-level power transmission nodes to obtain a first-level power transmission node set;
2.5B) generating a plane layout according to the geographical positions of all power transmission nodes in the node set V, and recording the reference position of each first-stage power transmission node in the first-stage power transmission node set in the plane layout;
2.6B) initialization step size k, radius r of each reference position i Wherein i represents the number of the reference position;
2.7B) taking each reference position as a circle center in a plane layout drawing, and taking the radius r of each reference position as a circle center i Increasing the step size k as a new radius r i Generating circles, judging whether the circles at two reference positions form an intersection, and if the circles at any two reference positions form an intersection, determining the radius r of the two reference positions i The original radius r before the step length k is increased is fixed i ；
2.8B) determining whether there is still a radius r i If the reference position with unfixed value still has radius r i Skipping to execute the step 2.7B) when the reference position with the unfixed value is taken; otherwise, skipping to execute the next step;
2.9B) judging whether the radius r which does not fall into each reference position still exists in the node set V i If the scattered power transmission nodes in the covered area still exist, the geographical position distance between each scattered power transmission node and each reference position is calculated for each scattered power transmission node, and the scattered power transmission nodes are added into the radius r of the reference position with the nearest geographical position distance i Within the area of coverage; finally, the radius r of each reference position is obtained i The covered area and the scattered power transmission nodes are used as initial subareas;
2.10B) forming a power transmission node pair aiming at the power transmission nodes which are in connection relation with other partitions in each partition, wherein the power transmission node pair is a power transmission section between the two corresponding partitions; jump execution step 3).
Optionally, the detailed steps of step 3) include: and calculating the power flow betweenness for all power transmission sections, and screening the power transmission sections with the power flow betweenness larger than a preset threshold value from all power transmission sections as finally obtained key power transmission sections.
Optionally, the function expression of the trend argument is:
in the above formula, F ij Representing the tidal current betweenness of power transmission sections i-j between the partitions i and j, G representing all power generation node sets providing power for the power transmission sections i-j, L representing a load node set to which the power transmission sections i-j flow, and min (S) m ，S n ) Is a weight coefficient of S m 、S n The smaller of the two, where S m Representing the actual output, S, of the generating node m n Representing the actual load of the load node n, P ij,m For the part of the power transmission line ij load flow between transmission sections i-j originating from the power generation node m, P ij,n Part of the flow of the power flow to the load n, P, of the transmission line ij n Representing the power flow of the load node n, P ij The i side active power of the transmission line ij; p Ln To track the power of the available load n by reverse order, A unm For the reverse-order distribution of the matrix, P, in n rows and m columns Gm The active output of the power generation node m.
In addition, the invention also provides a system for determining the key power transmission section of the power system, which is characterized by comprising the following steps:
the data input program unit is used for inputting a topology structure chart G of the power grid;
the power grid partitioning program unit is used for partitioning a power grid based on the topological structure chart G to obtain a power transmission section between partitions;
and the section screening program unit is used for screening key power transmission sections from all power transmission sections.
The invention also provides a system for determining a critical power transmission section of an electric power system, which comprises computer equipment and is characterized in that the computer equipment is programmed or configured to execute the steps of the method for determining the critical power transmission section of the electric power system.
The invention also provides a system for determining the critical power transmission section of the power system, which comprises a computer device, and is characterized in that a computer program which is programmed or configured to execute the method for determining the critical power transmission section of the power system is stored on a memory of the computer device.
The present invention also provides a computer-readable storage medium, characterized in that the computer-readable storage medium has stored thereon a computer program programmed or configured to perform the method for determining a critical power transmission section of an electrical power system.
Compared with the prior art, the invention has the following advantages: the method for determining the key power transmission section of the power system comprises the following steps: inputting a topology structure diagram G of a power grid; partitioning a power grid based on the topological structure diagram G to obtain a power transmission section between partitions; and screening out key power transmission sections from all power transmission sections, wherein an input topology structure diagram G of the power grid is composed of a node set V and an edge set E composed of edges between the nodes, each node in the node set V represents one power transmission node of the power grid, and each edge in the edge set E represents a power transmission line between two power transmission nodes of the power grid. The invention can adapt to the frequent change of the power system by adopting the topological structure diagram G of the power grid as input, and can quickly, accurately and automatically generate a power transmission section and a key power transmission section thereof based on the power grid partitioning of the topological structure diagram G of the power grid, can meet the requirements of quick and accurate power grid safety analysis, and has the advantages of high detection accuracy and high detection speed.
Drawings
FIG. 1 is a schematic diagram of a basic flow of a method according to an embodiment of the present invention.
Fig. 2 is a flowchart of an alternative implementation of step 2) in the embodiment of the present invention.
FIG. 3 is a flowchart of an alternative implementation of step 2) of the embodiment of the present invention.
Detailed Description
As shown in fig. 1, the implementation steps of the method for determining the key power transmission section of the power system in this embodiment include:
1) Inputting a topology structure diagram G of the power grid;
2) Partitioning a power grid based on the topological structure diagram G to obtain a power transmission section between partitions;
3) And screening out key power transmission sections from all power transmission sections.
The method for determining the key power transmission section of the power system in the embodiment adopts the topological structure diagram G of the power grid as input, so that the method for determining the key power transmission section of the power system in the embodiment can adapt to the frequent change condition of the power system.
In this embodiment, the topology structure diagram G of the power grid input in step 1) is formed by a node set V and an edge set E formed by edges between nodes, where each node in the node set V represents one power transmission node (substation) of the power grid, and each edge in the edge set E represents a power transmission line between two power transmission nodes of the power grid. Wherein the node set V is composed of nodes V 1 ～v m Forming m nodes in total; edge set E is defined by edge E 1 ～e n Make up, a total of n edges.
As an alternative embodiment, as shown in fig. 2, the power grid is a power grid including multiple voltage levels, and the sideband in the edge set E has transmission line voltage information, and the detailed steps of step 2) include:
2.1A) acquiring voltage information of all the transmission lines in the edge set E; for example, in this embodiment, voltages of all power transmission lines in a certain power grid include 750kV, 330kV, 220kV, 10kV, and the like;
2.2A) selecting a power transmission node connected with the edge with the highest voltage grade as a primary power transmission node to obtain a primary power transmission node set; in the embodiment, a power transmission node connected with a 750kV power transmission line is selected as a primary power transmission node, and an obtained primary power transmission node set is recorded as F _ nodes; by selecting the power transmission node connected with the edge with the highest voltage grade as the first-grade power transmission node, the partition based on the power transmission node connected with the edge with the highest voltage grade as the center node can be constructed, so that the partition is more accurate;
2.3A) aiming at each primary power transmission node in the primary power transmission node set F _ nodes, respectively establishing a corresponding process or thread and suspending the process or thread, so that each primary power transmission node corresponds to an independent process or thread, and initializing the tree structure of each primary power transmission node to only contain the primary power transmission node as a root node;
2.4A) initializing the level variable i to 1;
2.5A) controlling the processes or threads of all the first-level power transmission nodes to be started synchronously, searching the i-level power transmission nodes which are in direct or indirect connection with the current first-level power transmission node step by step through the process or thread traversal edge set E, if one i-level power transmission node is found and is not additionally provided with a use mark by other first-level power transmission nodes in the node set V, adding the i-level power transmission node to the i-level of the tree structure of the current first-level power transmission node according to the level of the i-level power transmission node, enabling a father node of the i-level power transmission node to be the current first-level power transmission node or the previous-level power transmission node which is connected with the current first-level power transmission node and additionally providing a use mark for the i-level power transmission node in the node set V; suspending the process or the thread;
2.6A) judging whether a power transmission node which is not additionally marked by use still exists in the node set V, if so, adding 1 to the level variable i, and skipping to execute the step 2.5A); otherwise, taking the obtained tree structure of each primary power transmission node as a partition; through parallel execution of processes or threads, a synchronous expansion spanning tree structure based on the first-level power transmission nodes can be realized for the node set V, the processing efficiency is high, and the tree structures of different first-level power transmission nodes can be prevented from being overlapped through synchronous expansion of a plurality of tree structures and additional use of marks;
2.7A) searching leaf nodes of the tree structure of each primary power transmission node, and forming a leaf node pair by one leaf node of the tree structure of any primary power transmission node and one leaf node of the tree structure of different primary power transmission nodes to obtain a leaf node pair set; for example, the partition 1 includes leaf nodes a, b, and c, and the partition 2 includes leaf nodes d, e, and f, so that leaf node pairs formed between the partition 1 and the partition 2 include ad, ae, af, bd, be, bf, cd, ce, and cf; leaf nodes between the same partitions belong to the same partition, and a power transmission section is not formed, so that leaf node pairs do not need to be formed;
2.8A) selecting a set of current leaf node pairs from the set of leaf node pairs;
2.9A) judging whether the matching of one edge of two power transmission nodes corresponding to the current leaf node in the edge set E is successful, and if so, judging that the current leaf node pair is a power transmission section between the two corresponding subareas;
2.10A) judging whether the traversing and the selection of the leaf node pair set are finished, if not, skipping to execute the step 2.8A); otherwise, judging that all the power transmission sections in the power grid are determined to be finished, and skipping to execute the step 3).
As another alternative, as shown in fig. 3, the power grid is a single voltage class power grid (e.g. 220kV power grid, etc.), and the nodes in the node set V have geographical location information of the transmission nodes, and the detailed step of step 2) includes:
2.1B) determining the number of edges connected with each power transmission node in the node set V to obtain an edge number binary array consisting of various edge numbers and the corresponding node numbers; for example, the number of edges connected to each power transmission node in a node set V in a certain power grid includes four types, i.e., 5, 1, 2, 4, and 3, and the number of each type is 8, 56, 34, 22, and 24, respectively, so that the obtained edge number binary array is:
{(5,8),(1,56),(2,34),(4,22),(3,24)}
2.2B) sorting the edge number binary arrays in a descending order according to the edge number; still taking the edge number binary array as an example, the result after sorting is: { (5, 8), (4, 22), (3, 24), (2, 34), (1, 56) };
2.3B) taking out an element from the head end of the edge number binary array, judging whether the ratio of the number of the nodes of all the taken-out elements to the number of all the nodes in the node set V exceeds a preset threshold value or not, and if the ratio exceeds the preset threshold value, skipping to execute the next step; otherwise, skipping to execute the step 2.3B again); the purpose of comparing with the preset threshold value is to control the number of the first-level power transmission nodes so as to avoid generating excessive partitions, and the value of the partition can be manually set according to the requirement;
2.4B) taking the power transmission nodes corresponding to all the taken-out elements as first-level power transmission nodes to obtain a first-level power transmission node set;
2.5B) generating a plane layout according to the geographical positions of all the power transmission nodes in the node set V, and recording the reference position of each first-stage power transmission node in the first-stage power transmission node set in the plane layout;
2.6B) initialization step size k, radius r of each reference position i Wherein i represents the number of the reference position;
2.7B) taking each reference position as a circle center in a plane layout drawing, and taking the radius r of each reference position as a circle center i Increasing the step size k as a new radius r i Generating circles, judging whether the circles at two reference positions form an intersection, and if the circles at any two reference positions form an intersection, determining the radius r of the two reference positions i The original radius r before the step length k is increased is fixed i ；
2.8B) determining whether there is still a radius r i If the reference position with unfixed value still has radius r i Skipping to execute the step 2.7B) when the reference position with the value not fixedly taken; otherwise, skipping to execute the next step; in the embodiment, by adopting a mode of expanding a plane layout diagram, the partition generation of each reference position can be simply and quickly realized, and the partitions of each reference position can be prevented from being overlapped;
2.9B) judging whether the radius r which does not fall into each reference position still exists in the node set V i If the scattered power transmission nodes of the covered area still exist, the geographical position distances of the scattered power transmission nodes and each reference position are respectively calculated for each scattered power transmission node, and the scattered power transmission nodes are added into the radius r of one reference position with the nearest geographical position distance i Within the area of coverage; finally, the radius r of each reference position is obtained i The covered area and the scattered power transmission nodes are used as initial subareas; by processing scattered power transmission nodes, full coverage of the power transmission nodes in a partitioned mode can be achieved;
2.10B) forming a power transmission node pair aiming at the power transmission nodes which are in connection relation with other partitions in each partition, wherein the power transmission node pair is a power transmission section between the two corresponding partitions; jump execution step 3).
The partitions connected with the key power transmission section are weak in physical connection, the transmission flow is large, large-scale transfer of the flow is easily caused when the key power transmission section breaks down, and the operation of the whole power grid is threatened, so that the key power transmission section is an object needing important monitoring. The detailed steps of step 3) in this embodiment include: and calculating the power flow betweenness for all power transmission sections, and screening the power transmission sections with the power flow betweenness larger than a preset threshold value from all power transmission sections as finally obtained key power transmission sections.
In this embodiment, the function expression of the power flow betweenness is:
in the above formula, F ij Representing the current betweenness of a power transmission section i-j between the subareas i and j, G representing all power generation node sets for providing power for the power transmission section i-j, L representing a load node set to which the power transmission section i-j flows, and min (S) m ，S n ) Is a weight coefficient of S m 、S n The smaller of the two, where S m Representing the actual output, S, of the generating node m n Representing the actual load of the load node n, P ij,m For the part of the power flow of the transmission line ij between transmission sections i-j originating from the power generation node m, P ij,n Part n, P of the flow direction load in the power flow of the transmission line ij n Representing the power flow of the load node n, P ij The active power of the i side of the transmission line ij; p Ln To track the power of the available load n by reverse order, A unm For the reverse-order distribution of the matrix, P, in n rows and m columns Gm The active output of the power generation node m.
And (3) adopting a power flow tracking method, respectively carrying out forward flow and reverse flow tracking calculation on the basis of the sequence/reverse sequence distribution matrix, and obtaining the power flow power composition in the power transmission line ij. The calculation formula of the part flowing to the load n in the power flow of the power transmission line ij is as follows:
in the above formula, P ij Is the i-side active power of line ij; p is i A node flow (a value equal to the sum of the injected or outgoing flows) for node i; p Ln Is the active load of the load node n; a. The d For sequential assignment of matrices, in denotes i rows and n columns.
The calculation formula of the part m from the generator in the power flow of the transmission line ij is as follows:
in the above formula, P ij Is the i-side active power of line ij; p is i A node flow (a value equal to the sum of the injected or outgoing flows) for node i; p Ln Is the active load of the load node n; p Gm The active output of the generator node m; a. The u The matrix is assigned in reverse order, im denotes i rows and m columns.
Sequential assignment matrix A d The ith row and j column elements of (1) are:
in the above formula, P j-1 Node flow, P, for node j-1 j For the node flow of node j,
is the outgoing set of node i.
Reverse order distribution matrix A u The ith row and j column elements of (1) are:
in the above formula, P j-1 Node flow, P, for node j-1 j For the node flow of node j,
is the incoming set of node i.
The power of the load n available by the reverse order tracking can be expressed as:
the contribution share of the power generation node m, namely the transmission power from the power generation node m to the load node n, is as follows:
the part of the transmission line ij flowing from the generator m to the load n is P ij (m, n) is:
in the formula, P ij,m Is the part of the line ij flow originating from the generator m; p ij,n Part of the flow in the power flow of the line ij towards the load n; g ij Set of all generators, L, providing power to line ij ij Is the set of all loads to which line ij flows. Therefore, the function expression of the tidal current betweenness in the embodiment can be derived as follows:
in the above formula, F ij Representing the tidal current betweenness of power transmission sections i-j between the partitions i and j, G representing all power generation node sets providing power for the power transmission sections i-j, L representing a load node set to which the power transmission sections i-j flow, and min (S) m ，S n ) Is a weight coefficient of S m 、S n The smaller of the two, where S m Represents the actual output, S, of the power generation node m n Representing the actual load of the load node n, P ij,m For the part of the power flow of the transmission line ij between transmission sections i-j originating from the power generation node m, P ij,n Part of the flow of the power flow to the load n, P, of the transmission line ij n Representing the power flow of the load node n, P ij The i side active power of the transmission line ij; p is Ln To track the power of the available load n by reverse order, A unm Allocating a matrix, P, for n rows and m columns in reverse order Gm The active output of the power generation node m.
In addition, this embodiment also provides a system for determining a critical transmission section of an electric power system, including:
the data input program unit is used for inputting a topology structure chart G of the power grid;
the power grid partitioning program unit is used for partitioning a power grid based on the topological structure diagram G to obtain a power transmission section between partitions;
and the section screening program unit is used for screening key power transmission sections from all power transmission sections.
In addition, the power grid partition program unit is programmed or configured to perform the foregoing steps 2.1A) -2.10A), and/or the foregoing steps 2.1B) -2.10B), which will not be described herein again; the section filter unit is programmed or configured to perform the aforementioned steps 3.1) -3.2).
In addition, the present embodiment also provides a system for determining a critical power transmission section of an electric power system, which includes a computer device programmed or configured to execute the steps of the method for determining a critical power transmission section of an electric power system.
In addition, the present embodiment also provides a system for determining a critical power transmission section of an electric power system, which includes a computer device, where a memory of the computer device stores a computer program programmed or configured to execute the method for determining a critical power transmission section of an electric power system.
Furthermore, the present embodiment also provides a computer-readable storage medium having stored thereon a computer program programmed or configured to execute the aforementioned method for determining a critical power transmission section of a power system.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein. The present application is directed to methods, apparatus (systems), and computer program products according to embodiments of the present application, wherein the instructions that execute on the flowcharts and/or processors of the computer program product create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks. These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only a preferred embodiment of the present invention, and the protection scope of the present invention is not limited to the above embodiments, and all technical solutions belonging to the idea of the present invention belong to the protection scope of the present invention. It should be noted that modifications and embellishments within the scope of the invention may occur to those skilled in the art without departing from the principle of the invention, and are considered to be within the scope of the invention.
Claims (7)
1. A method for determining a critical power transmission section of a power system is characterized by comprising the following implementation steps:
1) Inputting a topology structure diagram G of a power grid;
2) Partitioning a power grid based on the topological structure diagram G to obtain a power transmission section between partitions;
3) Screening out key power transmission sections from all power transmission sections;
the topology structure diagram G of the power grid input in the step 1) is composed of a node set V and an edge set E formed by edges between the nodes, wherein each node in the node set V represents one power transmission node of the power grid, and each edge in the edge set E represents a power transmission line between two power transmission nodes of the power grid;
the power grid is a power grid comprising a plurality of voltage levels, the sideband in the edge set E has power transmission line voltage information, and the detailed steps of the step 2) comprise:
2.1A) acquiring voltage information of all the transmission lines in the edge set E;
2.2A) selecting a power transmission node connected with the edge with the highest voltage grade as a first-grade power transmission node to obtain a first-grade power transmission node set;
2.3A) respectively establishing corresponding processes or threads for each primary power transmission node in the primary power transmission node set and suspending the processes or threads so that each primary power transmission node corresponds to an independent process or thread, and initializing the tree structure of each primary power transmission node to only contain the primary power transmission node as a root node;
2.4A) initializing the level variable i to 1;
2.5A) controlling the synchronous start of the process or the thread of each first-level power transmission node, gradually searching the ith-level power transmission node which has a direct or indirect connection relation with the current first-level power transmission node through the process or the thread traversing edge set E, if the ith-level power transmission node is found and is not additionally provided with a use mark by other first-level power transmission nodes in the node set V, attaching the ith-level power transmission node to the ith layer of the tree structure of the current first-level power transmission node according to the hierarchy of the ith-level power transmission node, and enabling a father node of the ith-level power transmission node to be the current first-level power transmission node or a previous-level power transmission node which is connected with the current first-level power transmission node and additionally providing the use mark for the ith-level power transmission node in the node set V; suspending the process or the thread;
2.6A) judging whether a power transmission node which is not additionally marked by use still exists in the node set V, if so, adding 1 to the level variable i, and skipping to execute the step 2.5A); otherwise, taking the obtained tree structure of each primary power transmission node as a partition;
2.7A) searching leaf nodes of the tree structure of each first-level power transmission node, and forming a leaf node pair by one leaf node of the tree structure of any first-level power transmission node and one leaf node of the tree structure of different first-level power transmission nodes to obtain a leaf node pair set;
2.8A) traversing and selecting a group of current leaf node pairs from the leaf node pair set;
2.9A) judging whether the matching of one edge of two power transmission nodes corresponding to the current leaf node in the edge set E is successful, and if so, judging that the current leaf node pair is a power transmission section between the two corresponding subareas;
2.10A) judging whether the traversing and the selection of the leaf node pair set are finished, if not, skipping to execute the step 2.8A); otherwise, judging that all the power transmission sections in the power grid are finished, and skipping to execute the step 3);
or the power grid is a single-voltage-class power grid, and nodes in the node set V have geographical position information of power transmission nodes, and the detailed step of the step 2) comprises the following steps:
2.1B) determining the number of edges connected with each power transmission node in the node set V to obtain an edge number binary array consisting of various edge numbers and the corresponding node numbers;
2.2B) sorting the edge number binary arrays in a descending order according to the edge number;
2.3B) taking out an element from the head end of the edge number binary array, judging whether the ratio of the number of the nodes of all the taken-out elements to the number of all the nodes in the node set V exceeds a preset threshold value or not, and if the ratio exceeds the preset threshold value, skipping to execute the next step; otherwise, skipping to execute the step 2.3B again);
2.4B) taking the power transmission nodes corresponding to all the taken-out elements as first-level power transmission nodes to obtain a first-level power transmission node set;
2.5B) generating a plane layout according to the geographical positions of all the power transmission nodes in the node set V, and recording the reference position of each first-stage power transmission node in the first-stage power transmission node set in the plane layout;
2.6B) initialization step size k, radius r of each reference position i Wherein i represents the number of the reference position;
2.7B) taking each reference position as a circle center in the plan layout drawing, and taking the radius r of each reference position as a circle center i Increasing the step k as the new radius r i Generating circles, judging whether the circles at two reference positions form an intersection, and if the circles at any two reference positions form an intersection, determining the radius r of the two reference positions i The original radius r before the step length k is increased is fixed i ；
2.8B) determining whether there is still a radius r i If the reference position with unfixed value still has radius r i Skipping to execute the step 2.7B) when the reference position with the unfixed value is taken; otherwise, skipping to execute the next step;
2.9B) judging whether the radius r which does not fall into each reference position still exists in the node set V i If the scattered power transmission nodes in the covered area still exist, the geographical position distance between each scattered power transmission node and each reference position is calculated for each scattered power transmission node, and the scattered power transmission nodes are added into the radius r of the reference position with the nearest geographical position distance i Within the area of coverage; finally, the radius r of each reference position is obtained i The covered area and the scattered power transmission nodes are used as initial subareas;
2.10B) forming a power transmission node pair aiming at the power transmission nodes which are in connection relation with other partitions in each partition, wherein the power transmission node pair is a power transmission section between the two corresponding partitions; jump execution step 3).
2. The method for determining the key transmission section of the power system according to claim 1, wherein the detailed step of step 3) comprises: and calculating the power flow betweenness for all power transmission sections, and screening the power transmission sections with the power flow betweenness larger than a preset threshold value from all power transmission sections as finally obtained key power transmission sections.
3. The method for determining a critical transmission profile of an electric power system according to claim 2, wherein the function expression of the power flow betweenness is as follows:
in the above formula, F ij Representing the tidal current betweenness of power transmission sections i-j between the partitions i and j, G representing all power generation node sets providing power for the power transmission sections i-j, L representing a load node set to which the power transmission sections i-j flow, and min (S) m ,S n ) Is a weight coefficient of S m 、S n The smaller of the two, where S m Representing the actual output, S, of the generating node m n Representing the actual load of the load node n, P ij,m For the part of the power transmission line ij load flow between transmission sections i-j originating from the power generation node m, P ij,n Part n, P of the flow direction load in the power flow of the transmission line ij n Representing the power flow of the load node n, P ij The i side active power of the transmission line ij; p Ln To track the power of the available load n by reverse order, A unm For the reverse-order distribution of the matrix, P, in n rows and m columns Gm The active output of the power generation node m.
4. A system for determining a critical transmission profile of an electrical power system, comprising:
the data input program unit is used for inputting a topology structure chart G of the power grid;
the power grid partitioning program unit is used for partitioning a power grid based on the topological structure chart G to obtain a power transmission section between partitions;
the section screening program unit is used for screening key power transmission sections from all power transmission sections;
the topology structure diagram G of the power grid input in the data input program unit is composed of a node set V and an edge set E formed by edges between the nodes, each node in the node set V represents one power transmission node of the power grid, and each edge in the edge set E represents a power transmission line between two power transmission nodes of the power grid;
the power grid is a power grid comprising a plurality of voltage levels, the voltage information of the power transmission line is arranged on one side of the side set E, and the detailed steps of partitioning the power grid by the power grid partitioning program unit based on the topology structure diagram G comprise:
2.1A) acquiring voltage information of all the power transmission lines in the edge set E;
2.2A) selecting a power transmission node connected with the edge with the highest voltage grade as a first-grade power transmission node to obtain a first-grade power transmission node set;
2.3A) aiming at each primary power transmission node in the primary power transmission node set, respectively establishing a corresponding process or thread and suspending the process or thread so that each primary power transmission node corresponds to an independent process or thread, and initializing the tree structure of each primary power transmission node to only contain the primary power transmission node as a root node;
2.4A) initializing the level variable i to 1;
2.5A) controlling the synchronous start of the process or the thread of each first-level power transmission node, gradually searching the ith-level power transmission node which has a direct or indirect connection relation with the current first-level power transmission node through the process or the thread traversing edge set E, if the ith-level power transmission node is found and is not additionally provided with a use mark by other first-level power transmission nodes in the node set V, attaching the ith-level power transmission node to the ith layer of the tree structure of the current first-level power transmission node according to the hierarchy of the ith-level power transmission node, and enabling a father node of the ith-level power transmission node to be the current first-level power transmission node or a previous-level power transmission node which is connected with the current first-level power transmission node and additionally providing the use mark for the ith-level power transmission node in the node set V; suspending the process or the thread;
2.6A) judging whether a power transmission node which is not additionally marked by use still exists in the node set V, if so, adding 1 to the level variable i, and skipping to execute the step 2.5A); otherwise, taking the obtained tree structure of each primary power transmission node as a partition;
2.7A) searching leaf nodes of the tree structure of each first-level power transmission node, and forming a leaf node pair by one leaf node of the tree structure of any first-level power transmission node and one leaf node of the tree structure of different first-level power transmission nodes to obtain a leaf node pair set;
2.8A) traversing and selecting a group of current leaf node pairs from the leaf node pair set;
2.9A) judging whether the matching of one edge of two power transmission nodes corresponding to the current leaf node in the edge set E is successful, and if so, judging that the current leaf node pair is a power transmission section between the two corresponding subareas;
2.10A) judging whether the traversing and the selection of the leaf node pair set are finished, if not, skipping to execute the step 2.8A); otherwise, judging that all the power transmission sections in the power grid are finished;
or the power grid is a single voltage-class power grid, the nodes in the node set V have geographical position information of power transmission nodes, and the detailed steps of the power grid partitioning program unit for partitioning the power grid based on the topology structure diagram G comprise:
2.1B) determining the number of edges connected with each power transmission node in the node set V to obtain an edge number binary array consisting of various edge numbers and the corresponding node numbers;
2.2B) sorting the edge number binary arrays in a descending order according to the edge numbers;
2.3B) taking out an element from the head end of the edge number binary array, judging whether the ratio of the number of the nodes of all the taken-out elements to the number of all the nodes in the node set V exceeds a preset threshold value or not, and if the ratio exceeds the preset threshold value, skipping to execute the next step; otherwise, skipping to execute the step 2.3B again);
2.4B) taking the power transmission nodes corresponding to all the taken-out elements as first-level power transmission nodes to obtain a first-level power transmission node set;
2.5B) generating a plane layout according to the geographical positions of all the power transmission nodes in the node set V, and recording the reference position of each first-stage power transmission node in the first-stage power transmission node set in the plane layout;
2.6B) initialization step size k, radius r of each reference position i Wherein i represents the number of the reference position;
2.7B) taking each reference position as a circle center in the plan layout drawing, and taking the radius r of each reference position as a circle center i Increasing the step k as the new radius r i Generating circles, judging whether circles at two reference positions form an intersection, and if circles at any two reference positions form an intersection, determining the radius r of the two reference positions i The original radius r before the step length k is increased is fixed i ；
2.8B) determining whether there is still a radius r i If the reference position with unfixed value still has radius r i Skipping to execute the step 2.7B) when the reference position with the value not fixedly taken; otherwise, skipping to execute the next step;
2.9B) judging whether the radius r which does not fall into each reference position still exists in the node set V i If the scattered power transmission nodes in the covered area still exist, the geographical position distance between each scattered power transmission node and each reference position is calculated for each scattered power transmission node, and the scattered power transmission nodes are added into the radius r of the reference position with the nearest geographical position distance i Within the area of coverage; finally, the radius r of each reference position is obtained i The covered area and the scattered power transmission nodes are used as initial subareas;
2.10B) forming a power transmission node pair aiming at the power transmission nodes which have connection relations with other partitions in each partition, wherein the power transmission node pair is a power transmission section between two corresponding partitions.
5. A system for determining a critical power transmission section of an electrical power system, comprising a computer device, characterized in that the computer device is programmed or configured to perform the steps of the method for determining a critical power transmission section of an electrical power system according to any one of claims 1 to 3.
6. A system for determining a critical power transmission section of an electrical power system, comprising a computer device, characterized in that a computer program programmed or configured to perform a method for determining a critical power transmission section of an electrical power system according to any one of claims 1 to 3 is stored on a memory of the computer device.
7. A computer-readable storage medium having stored thereon a computer program programmed or configured to perform a method of determining a critical power transmission profile of a power system according to any one of claims 1 to 3.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010611996.5A CN111697579B (en) | 2020-06-30 | 2020-06-30 | Method, system and medium for determining key power transmission section of power system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010611996.5A CN111697579B (en) | 2020-06-30 | 2020-06-30 | Method, system and medium for determining key power transmission section of power system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111697579A CN111697579A (en) | 2020-09-22 |
| CN111697579B true CN111697579B (en) | 2022-11-29 |
Family
ID=72484610
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010611996.5A Active CN111697579B (en) | 2020-06-30 | 2020-06-30 | Method, system and medium for determining key power transmission section of power system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111697579B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112989129A (en) * | 2021-03-18 | 2021-06-18 | 国家电网有限公司 | Power transmission network key perception section identification method and system based on graph theory |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018049737A1 (en) * | 2016-09-18 | 2018-03-22 | 国电南瑞科技股份有限公司 | Safe correction calculation method based on partition load control |
| CN109638838A (en) * | 2019-01-21 | 2019-04-16 | 广东电网有限责任公司 | The recognition methods of power grid key sections, device and electronic equipment |
-
2020
- 2020-06-30 CN CN202010611996.5A patent/CN111697579B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2018049737A1 (en) * | 2016-09-18 | 2018-03-22 | 国电南瑞科技股份有限公司 | Safe correction calculation method based on partition load control |
| CN109638838A (en) * | 2019-01-21 | 2019-04-16 | 广东电网有限责任公司 | The recognition methods of power grid key sections, device and electronic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111697579A (en) | 2020-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Musharavati et al. | Enhanced simulated-annealing-based algorithms and their applications to process planning in reconfigurable manufacturing systems | |
| CN110609759B (en) | Fault root cause analysis method and device | |
| CN108009016A (en) | A kind of balancing resource load control method and colony dispatching device | |
| CN108345544B (en) | Software defect distribution influence factor analysis method based on complex network | |
| CN110135716B (en) | Power grid infrastructure project dynamic early warning identification method and system | |
| CN104915897B (en) | A kind of computer implemented method of Electric Power Network Planning evaluation assignment | |
| CN109034562A (en) | A kind of social networks node importance appraisal procedure and system | |
| CN110147632A (en) | A kind of topology matching route bus method considering non-uniform track and barrier | |
| JP4327783B2 (en) | Node group sorting apparatus and method for power system | |
| CN111697579B (en) | Method, system and medium for determining key power transmission section of power system | |
| CN107343164A (en) | The distribution method and device of vision control task | |
| Sergeeva et al. | Dynamic airspace configurations generated by evolutionary algorithms | |
| Meincheim et al. | Combining process mining with trace clustering: manufacturing shop floor process-an applied case | |
| Guo et al. | Identification of power system dynamic signature using hierarchical clustering | |
| CN106780747A (en) | A kind of method that Fast Segmentation CFD calculates grid | |
| CN106844932A (en) | Micro-capacitance sensor minimal cut set based on BFS quickly seeks method | |
| CN111695583A (en) | Feature selection method based on causal network | |
| Rodrigues et al. | Evaluation of nosql dbms in private cloud environment: an approach based on stochastic modeling | |
| Li et al. | ChatTwin: Toward Automated Digital Twin Generation for Data Center via Large Language Models | |
| KR20140006491A (en) | Effective graph clustering apparatus and method for probabilistic graph | |
| Wagner et al. | Automated generation of material flow network models with cluster structures using random walks | |
| CN111667149A (en) | System efficiency evaluation method for large-scale scientific and technological engineering, simulation and expert evaluation | |
| Maghsoodi et al. | Predicting cascading failure with machine learning methods in the interdependent networks | |
| CN113111523B (en) | Meta-heuristic method for solving intelligent layout of multi-constraint double-layer aisle | |
| CN110955521A (en) | Power distribution network wide area distributed type sinking calculation system and method based on big data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |