CN112487658A - Method, device and system for identifying key nodes of power grid - Google Patents

Method, device and system for identifying key nodes of power grid Download PDF

Info

Publication number
CN112487658A
CN112487658A CN202011468470.2A CN202011468470A CN112487658A CN 112487658 A CN112487658 A CN 112487658A CN 202011468470 A CN202011468470 A CN 202011468470A CN 112487658 A CN112487658 A CN 112487658A
Authority
CN
China
Prior art keywords
power
node
nodes
power grid
key
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.)
Granted
Application number
CN202011468470.2A
Other languages
Chinese (zh)
Other versions
CN112487658B (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.)
Chongqing University of Post and Telecommunications
Original Assignee
Chongqing University of Post and Telecommunications
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 Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN202011468470.2A priority Critical patent/CN112487658B/en
Publication of CN112487658A publication Critical patent/CN112487658A/en
Application granted granted Critical
Publication of CN112487658B publication Critical patent/CN112487658B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid distribution networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention belongs to the field of electric power data processing, and relates to a method, a device and a system for identifying key nodes of a power grid; the method comprises the steps of constructing a power-efficient power grid topological structure according to validity rules of power node injection power and line transmission power in a power grid; calculating the average electrical distance and the electrical betweenness center of the power nodes in the power grid topological structure, respectively distributing weight vectors, and constructing a criticality evaluation model of the nodes; and calculating a power node core value according to the criticality evaluation model, carrying out recursive network decomposition according to the interval to which the power node core value belongs, and taking the power node contained in the last layer of sub-network as an identified key node set. According to the method, the rule characteristics of the large-scale power system are excavated according to a series of rules, the standard power grid topological structure is constructed, the network characteristics and the electrical characteristics of the power system are comprehensively considered, a multi-factor-based power grid key node identification model is established, and the efficiency and the accuracy of key node identification in the complex power grid system are improved.

Description

Method, device and system for identifying key nodes of power grid
Technical Field
The invention belongs to the field of electric power data processing, and relates to a method, a device and a system for identifying a key node of a power grid.
Background
In recent years, rapid development of various emerging technologies makes people have a qualitative leap in life level. The power system is an important infrastructure for the development of all scientific and technical technologies, and makes remarkable contribution in the development process of other social networks such as the internet, a communication network, a traffic network and the like. With the continuous expansion of the scale of the power grid and the improvement of the integration level of the components, the structural characteristics of the power grid become more and more complex, and the improvement of the complexity plays a significant role in maintaining the efficient power transmission. Although the modern power grid forms a cross-regional large power grid pattern with high voltage level and long-distance transmission, and interconnection and intercommunication of network architectures are basically realized, safe and stable operation of the power grid is still a big problem which is painful and difficult to solve. The safe and stable operation of the power system depends on a plurality of factors, such as hardware devices of the power grid, the climate around the power grid, and the like, wherein the safety of some key parts in the power grid is more important to be ignored. Therefore, accurately identifying the key nodes or lines in the power grid system and taking certain protective measures on the key nodes or lines are important means for maintaining the robustness of the power system.
In the past, most of researches on large-area power failure accidents are based on a theory of restoration, the essence of which is that firstly, a rule is discovered by simplifying a research object, then, a system rule is restored, and the theoretical research is usually carried out by limiting the physical properties of elements. The reduction theory method generally describes a power grid into a group of multidimensional differential algebraic equations, and then solves the equations through computer simulation, so that the evolution characteristics in the whole system are ignored to a certain extent, and therefore a new thought is needed to research the characteristics of the power grid. The complex network theory is that a bus is abstracted into a network node, and a transmission line is abstracted into a network link, so that a power grid is regarded as a network with interaction between units or individuals, and a new thought is provided for researching the structure of the power grid and the cascade reaction process.
Based on the prior art research, it can be found that many large-scale grid fault accidents are caused by the fault of some special nodes or lines in the system. Therefore, the method accurately identifies the key nodes or lines in the power grid system and takes targeted protective measures, and is an effective means for preventing the large-scale cascading faults of the power grid. At present, most of the key nodes of the power grid system are identified mainly from the perspective of a power grid topological structure by combining a general network model (such as a small-world network model, a scale-free network model, a regular network model and the like) and common indexes (such as a clustering coefficient, degree and the like) in a complex network theory to identify the key nodes in the power system. However, in a complex power system, the criticality of a node in a power grid cannot be analyzed only from the viewpoint of topology, and more physical properties of the power system, such as directional transmission of power of a power grid line, influence of a change of a state of a certain node on the whole power grid flow, and the like, need to be considered.
In addition, because the conventional sorting algorithm (e.g., k _ shell algorithm, toposis algorithm, pageank algorithm, etc.) has certain limitations in many practical application scenarios, on one hand, the prior art only considers the basic characteristics of the network topology, and is not combined with the practical power grid scenario, for example, power transmission is not considered, so that certain improvement needs to be made on the prior art to identify the key nodes in the power grid system more accurately.
Disclosure of Invention
In view of the above, the invention aims to provide a power grid key node identification method based on an improved k _ shell algorithm, which combines the network topology characteristics and the electrical characteristics of a power system to establish a new model for multi-factor evaluation of node criticality, and provides a power grid key node identification method, device and system by improving the k _ shell algorithm.
The scheme adopted by the invention for solving the technical problems comprises the following steps:
in a first aspect of the present invention, the present invention provides a method for identifying a key node of a power grid, where the method includes the following steps:
step 1) constructing a power-efficient power grid topological structure according to validity rules of power node injection power and line transmission power in a power grid;
step 2) calculating the average electrical distance and the electrical permittivity center of the power node in the power grid topological structure, respectively distributing weight vectors, and constructing a criticality evaluation model of the node;
and 3) calculating a power node core value according to the criticality evaluation model, performing recursive network decomposition according to the interval to which the power node core value belongs, and taking the power node contained in the last layer of sub-network as an identified key node set.
In a second aspect of the present invention, the present invention further provides an identification apparatus for a key node of a power grid, the identification apparatus comprising:
the topological structure constructing module is used for constructing a power-efficient power grid topological structure according to validity rules of power node injection power and line transmission power in a power grid;
the criticality evaluation module is used for calculating the average electrical distance and the electrical betweenness center of the power nodes, respectively distributing weight vectors and constructing a criticality evaluation model;
the core value calculation module is used for calculating the core value of each power node according to the criticality evaluation model;
the interval dividing module is used for dividing a plurality of intervals according to the core values of the power nodes and determining the minimum interval;
the topology updating module is used for deleting the power nodes in the minimum interval, updating the power grid topology structure and returning to the core value calculating module;
and the node output module is used for outputting the power node of the last interval, namely the power node is the key point set.
In a third aspect of the present invention, the present invention further provides an identification system for a key node of a power grid, where the identification system includes:
the data acquisition module is used for acquiring the power data of the power nodes in the power grid;
and the identification device of the key nodes of the power grid is used for identifying the collected power data of the power nodes and outputting a key point set.
The invention has the beneficial effects that:
1. according to the power grid topological structure, the power grid topological structure with power effectiveness is constructed according to the effectiveness rules of power node injection power and line transmission power in a power grid, the electrical characteristics of a complex power system can be reflected by the power grid topological structure constructed by the method, other interference factors are eliminated, and compared with the existing power grid topological structure, the method improves the operation efficiency of the power grid structure to a certain extent, so that key nodes in the power grid can be extracted quickly and accurately.
2. The method is characterized in that a new node criticality evaluation model is provided according to the average path length and the betweenness center of the power nodes and in combination with electrical characteristics, and a weight vector of the criticality evaluation model is obtained by utilizing an analytic hierarchy process for training; the evaluation model considers the network characteristics and the electrical characteristics of the power system at the same time, and can reflect the key degree of the power nodes more objectively;
3. according to the method, the core values of the nodes are redefined according to the criticality evaluation model, the core values of all the nodes are subjected to interval division, recursive network decomposition is performed according to the space to which the core values of the nodes belong, and the nodes contained in the last layer of sub-network are used as a final key point set. The method improves the limitation of the k _ shell algorithm in a power grid application scene, and can improve the accuracy of identifying the key nodes of the power grid.
Drawings
For the purposes of promoting a better understanding of the objects, aspects and advantages of the invention, reference will now be made to the following detailed description taken in conjunction with the accompanying drawings in which:
FIG. 1 is a general structure diagram of a power grid key node identification model constructed by the invention;
fig. 2 is a flowchart of a method for identifying a key node of a power grid according to an embodiment of the present invention;
FIG. 3 is a flow chart of an analytic hierarchy process for assigning weights used in an embodiment of the present invention;
FIG. 4 is a flow chart of a key node set screening process using an improved k _ shell algorithm according to an embodiment of the present invention;
fig. 5 is a structural diagram of an identification apparatus of a key node of a power grid according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The method is used for researching the characteristics of the large-scale power system based on the complex network theory, comprehensively considering the network characteristics and the electrical characteristics of the power system and combining with the ranking algorithm based on multi-factor improvement, so that the key nodes in the power system can be accurately identified.
Fig. 1 is a node identification model architecture diagram adopted in the present invention, as shown in fig. 1, after power data acquired by a power system is abstracted and normalized, a power grid topology structure is formed according to validity rules of power node injection power and line transmission power, wherein nodes in the power grid topology structure include power generation nodes, intermediate nodes (transmission nodes) and load nodes, in the power grid topology structure, an average electrical aggregation and electrical permittivity center of the power nodes are calculated, and weights are assigned to them according to an analytic hierarchy process, so as to construct a new criticality evaluation model; and decomposing the network according to a k-shell algorithm improved by a section division method, thereby extracting a key node set.
Fig. 2 is a flowchart of an identification method for a key node of a power grid according to the present invention, and as shown in fig. 2, the identification method includes the following steps:
step 1) constructing a power-efficient power grid topological structure according to validity rules of power node injection power and line transmission power in a power grid;
in this embodiment, a power station, a substation, or an intermediate connection device in a power system is abstracted as a power node of a power network, and a transmission line is used as a power edge of the network; classifying the power nodes in the power grid based on validity rules of power node injection power and line transmission power in the power grid, and setting that at most one edge is connected between every two power nodes, thereby constructing a power validity power grid topological structure.
The node injection power validity rule is that for the whole power grid system, the nodes are subdivided into power generation nodes, transmission nodes and load nodes according to the net injection power values of the nodes; specifically, for the whole power grid system, if the net injection power of the node is a positive value, the node is a power generation node for injecting active power into the power grid; if the net injection power of the node is a negative value, the node is a load node for receiving power to the power grid; if the net injected power of the node is zero, the node is a transmitting node.
The power validity rule of circuit transmission means that a strategy of forward superposition and reverse cancellation is adopted, and only one edge with non-negative transmission power is connected between the specified node pairs; specifically, if a plurality of lines are connected between pairs of power nodes, the transmission power of the plurality of lines is processed by adopting a forward superposition and reverse cancellation strategy because the power transmitted by each line may have differences in magnitude and direction, and finally, it is specified that at most one edge with non-negative transmission power is connected between every two power nodes.
Step 2) calculating the average electrical distance and the electrical permittivity center of the power node in the power grid topological structure, respectively distributing weight vectors, and constructing a criticality evaluation model of the node;
firstly, for the average electrical distance of the power nodes, on the basis of obtaining the power grid topological structure, the average path length of the power nodes is considered, and the calculation method is referred to as follows by combining the electrical characteristics of the power grid:
solving the shortest path i, k between any power node pair ij based on the shortest path algorithm1,k2,…,kmJ, and the length of the shortest path, the sum of the transmission power P is solved based on the load of the power node itself and the line transmission powerijExpressed as:
Figure BDA0002835387690000061
wherein, wijRefers to the actual load between power node i and power node j, defined as wij=min(Si,Sj) And w isii=0,SiRefers to the rated generating capacity, S, of the power node ijIs the rated generating capacity of power node j.
The shortest path algorithm can adopt a floyd-warshall algorithm, and the Flouard algorithm is an algorithm for finding the shortest path between multiple source points in a given weighted graph by using a dynamic programming idea, and can accurately find out the shortest path between any power node pair ij.
Electrical distance d between any pair of nodesijThe ratio of the sum of power to the length of the shortest path between two nodes is defined as follows:
Figure BDA0002835387690000062
wherein M isijThe length of the shortest or most effective path between the power node i and the power node j is defined, and if the shortest or most effective path between the power node i and the power node j cannot be reached, M is definedijInfinity, and has Mii=1。
On the basis of the shortest path between the power nodes, the invention provides the definition of the average electrical distance of the power nodes, which is expressed as:
Figure BDA0002835387690000063
wherein D isiRepresents the average electrical distance of power node i; v represents a power node set in the power grid; and N refers to the total number of power nodes in the power grid.
Secondly, after the average electrical distance of the power nodes is solved, the electrical node betweenness centers of the three power nodes are respectively defined by combining the electrical characteristics of the power grid for the electrical node betweenness centers on the basis of the power grid topological structure;
the betweenness center of a node in a general complex network is generally used to characterize the criticality of the node in the whole network, depending on the number of times that the shortest path or the most effective path in the network system passes through the node, and therefore, the betweenness center in a conventional power grid is defined as follows:
Figure BDA0002835387690000071
after normalization processing is carried out on the data, the data are expressed as follows:
Figure BDA0002835387690000072
wherein σijIs the number of shortest or most efficient paths, σ, between power node i and power node jij(k) G, L refers to the set of power generation nodes and the set of load nodes, respectively, as the number of passing power nodes k in the shortest or most efficient path between power node i and power node j.
On the basis of normalizing the betweenness central index, the invention also combines the electrical and directional transmission characteristics of the power grid, respectively defines the electrical betweenness central index for three conditions of a power generation node, a load node, an intermediate node and the like, and is represented as follows:
Figure BDA0002835387690000073
wherein, PijRefers to the sum of the power transmitted through the shortest or most efficient path between power node i and power node j, and G, L refers to the power generation node set and the load node set respectively; wherein in the upper formula
Figure BDA0002835387690000075
Then it indicates that power node k belongs to the middleA node set; considering the directionality of the power nodes in the grid, PkiRepresents the sum of the power transmitted over the shortest or most efficient path between the generating node k and the other power nodes i; pikRepresenting the sum of the power transmitted over the shortest or most efficient path between the other power node i and the load node k.
After the average electrical distance and the electrical betweenness center of the power node are solved, the invention also distributes weight vectors to the average electrical distance and the electrical betweenness center respectively based on an analytic hierarchy process
Figure BDA0002835387690000074
And constructing a new node criticality evaluation model according to an objective principle of node criticality setting.
Fig. 3 is a flow chart of assigning weights by an analytic hierarchy process according to an embodiment of the present invention, where as shown in fig. 3, the assigning process includes the following steps:
establishing a hierarchical structure model, namely a decision criterion of an average electrical distance and an electrical betweenness center of an input power node;
constructing a judgment matrix, and comparing the decision criteria in pairs by adopting relative scales to obtain a judgment matrix A:
Figure BDA0002835387690000081
wherein a isijWhich refers to the scale value of pairwise comparisons between decision criteria. In order to verify the reasonability of the setting of the scale value, whether the constructed judgment matrix has consistency is checked by defining a consistency ratio, so that a weight vector of the index is determined; the matrix a needs to be checked for consistency next.
If A satisfies
Figure BDA0002835387690000082
A is a consistency matrix, the rank is 1, the only nonzero eigenvalue is lambda, and the eigenvector corresponding to the eigenvalue is normalized and then is used as a weight vector psi; if A is not oneThe consistency matrix is checked to see if the inconsistency of A is within a reasonable range. Taking the maximum eigenvalue lambda of A, defining a consistency index CI:
Figure BDA0002835387690000083
in order to measure the CI size more objectively, m consistency matrixes are constructed randomly, and a random consistency index RI is introduced:
Figure BDA0002835387690000084
then a consistency ratio can be defined:
Figure BDA0002835387690000085
only when CR <0.1 is satisfied, a is considered to have satisfactory consistency, and the eigenvector corresponding to the maximum eigenvalue λ of a is taken to perform normalization operation, thereby obtaining the weight vector ψ. If not, resetting the scale value of the judgment matrix A until obtaining a judgment matrix meeting the condition.
And 3) calculating a power node core value according to the criticality evaluation model, performing recursive network decomposition according to the interval to which the power node core value belongs, and taking the power node contained in the last layer of sub-network as an identified key node set.
FIG. 4 is a flow chart of a key node set screening process employed by the embodiments of the present invention; as shown in fig. 4, the method performs network decomposition based on an improved k _ shell algorithm, and finally identifies a key node set of the power system, and mainly includes the following steps:
1) on the basis of using an analytic hierarchy process, the weight vector of the average electrical distance of the power node and the electrical dielectric constant center is obtained
Figure BDA0002835387690000091
According to nodeThe smaller the average electrical distance is and the larger the electrical medium center is, the more critical the node is to the power grid system is shown, and a new node criticality evaluation model is constructed:
Figure BDA0002835387690000092
wherein k (i) represents criticality of power node i;
Figure BDA0002835387690000093
representing a first weight vector, DiRepresents the average electrical distance of the power node i,
Figure BDA0002835387690000094
representing a second weight vector; b ise(i) Represents the electrical betweenness center of the power node i.
According to the criticality evaluation model, defining a new node core value calculation formula as follows:
Figure BDA0002835387690000095
wherein ks (i) represents a core value of power node i; kminRepresenting the minimum value of the criticality of all nodes in the power grid; kmaxRepresenting the maximum value of criticality of all power nodes in the power grid; k (i) represents the criticality of power node i.
In the embodiment, a maximum and minimum normalization mode is adopted, the key values are scaled in an equal proportion to serve as new node core values, and after the core values are fixed in a [0,1] interval, the centralized distribution degree of the core values can be enhanced, so that a space can be divided for recursive network decomposition in the follow-up process.
2) Calculating the core value of each power node, and dividing all the core values into l according to reasonable difference values1Individual interval (pi)12…,Πl1) Removing those kernels belonging to the minimum interval
Figure BDA0002835387690000096
After the nodes are removed, a subgraph is left in the power grid topological structure, the core value of each power node in the subgraph is recalculated, and if the core value of the power node still belongs to the subgraph
Figure BDA0002835387690000097
Within the interval, the power nodes are continuously deleted until a subgraph G is left1All the power nodes in the system have different core values
Figure BDA0002835387690000098
Within the interval, the deleted power nodes belong to an S (1) set;
in some embodiments, the key value calculated by the criticality evaluation model is directly used as an interval division index of the power node core value; the division standard of the space of the kernel value comprises the steps of calculating key values of all nodes, and carrying out grade division according to the distribution condition of the key values, wherein the range of the interval depends on the distribution density, if the distribution is dense, the space range is divided into small spaces, and otherwise, a large space can be divided.
In other embodiments, in addition to the above division according to the key value, referring to the above embodiments, the present embodiment may also divide the kernel space according to the new kernel calculation formula provided in the present invention.
3) Similar to the above steps, the erasure of a core value belongs to
Figure BDA0002835387690000101
A power node within the interval, wherein
Figure BDA0002835387690000102
Is full of
Figure BDA0002835387690000103
Finally, subgraph G is obtained2,G2The core values of all the power nodes in the network are larger than the interval
Figure BDA0002835387690000104
A value of;
4) and in the same way, until all the last power node core values are located in an interval, the last power node core values are the key power node set in the power grid topological structure.
It can be understood that the recursive network and the last layer of sub-network of the invention both refer to that the power grid topology structure comprises a power grid topology structure before updating and a sub-network after updating, and in addition, the improvement of the k-shell algorithm of the invention mainly lies in dividing a space for recursive network decomposition and improving the traditional kernel value calculation.
Fig. 5 is a structural diagram of an identification device of a key node of a power grid according to the present invention, and as shown in fig. 5, the identification device includes:
the topological structure constructing module is used for constructing a power-efficient power grid topological structure according to validity rules of power node injection power and line transmission power in a power grid;
the criticality evaluation module is used for calculating the average electrical distance and the electrical betweenness center of the power nodes, respectively distributing weight vectors and constructing a criticality evaluation model;
the core value calculation module is used for calculating the core value of each power node according to the criticality evaluation model;
the interval dividing module is used for dividing a plurality of intervals according to the core values of the power nodes and determining the minimum interval;
the topology updating module is used for deleting the power nodes in the minimum interval, updating the power grid topology structure and returning to the core value calculating module;
and the node output module is used for outputting the power node of the last interval, namely the power node is the key point set.
In this embodiment, the core value calculation module and the topology update module are directly connected except for using the interval division module as intermediate transition, after the core value of each power node in the power grid topology structure is calculated, different intervals are divided for the power nodes according to the interval division module, in the topology update module, the power node in the minimum interval is removed according to a k-shell mode, the power grid topology structure is updated, the updated power grid topology structure is used as a sub-graph, the core value of each power node in the sub-graph is calculated again through the core value calculation module, and at first, the power node in the sub-graph which still exists in the minimum interval is removed; and then, continuously removing the power nodes in the secondary cell according to a k-shell mode, repeating the iteration process until all the last power node core values are located in the last interval, namely the last subgraph, and taking all the power node sets in the subgraph as key power node sets in the power grid topological structure.
In some embodiments, the present invention further provides an identification system for a critical node of a power grid, the identification system comprising:
the data acquisition module is used for acquiring the power data of the power nodes in the power grid;
and the identification device of the key nodes of the power grid is used for identifying the collected power data of the power nodes and outputting a key point set.
It can be understood that the method, the device and the system for identifying key nodes in the invention belong to the same inventive concept, and corresponding features thereof can be mutually cited, which is not described in detail herein.
In the description of the present invention, it is to be understood that the terms "coaxial", "bottom", "one end", "top", "middle", "other end", "upper", "one side", "top", "inner", "outer", "front", "center", "both ends", and the like, indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplicity of description, and do not indicate or imply that the devices or elements referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, are not to be construed as limiting the present invention.
In the present invention, unless otherwise expressly stated or limited, the terms "mounted," "disposed," "connected," "fixed," "rotated," and the like are to be construed broadly, e.g., as meaning fixedly connected, detachably connected, or integrally formed; can be mechanically or electrically connected; the terms may be directly connected or indirectly connected through an intermediate, and may be communication between two elements or interaction relationship between two elements, unless otherwise specifically limited, and the specific meaning of the terms in the present invention will be understood by those skilled in the art according to specific situations.
Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that changes, modifications, substitutions and alterations can be made in these embodiments without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.

Claims (9)

1. A method for identifying a key node of a power grid is characterized by comprising the following steps:
step 1) constructing a power-efficient power grid topological structure according to validity rules of power node injection power and line transmission power in a power grid;
step 2) calculating the average electrical distance and the electrical permittivity center of the power node in the power grid topological structure, respectively distributing weight vectors, and constructing a criticality evaluation model of the node;
and 3) calculating a power node core value according to the criticality evaluation model, performing recursive network decomposition according to the interval to which the power node core value belongs, and taking the power node contained in the last layer of sub-network as an identified key node set.
2. The method for identifying key nodes of the power grid according to claim 1, wherein the step 1) comprises dividing the power nodes into power generation nodes, transmission nodes and load nodes according to the net injection power values of the power nodes, and adopting a forward superposition and reverse cancellation strategy to connect only one edge with non-negative transmission power between the power nodes, so as to construct the power grid topology.
3. The method for identifying the key node of the power grid according to claim 1, wherein the key degree evaluation model of the node is represented as:
Figure FDA0002835387680000011
wherein k (i) represents criticality of power node i;
Figure FDA0002835387680000012
representing a first weight vector, DiRepresents the average electrical distance of the power node i,
Figure FDA0002835387680000013
representing a second weight vector; b ise(i) Represents the electrical betweenness center of the power node i.
4. The method for identifying the key nodes of the power grid according to claim 1 or 3, wherein the calculation mode of the weight vector comprises the steps of constructing a hierarchical structure model of decision criteria of average electrical distance and electrical medium center of the nodes by adopting an analytic hierarchy process; comparing the two decision criteria pairwise by adopting relative scale to obtain a judgment matrix; and detecting whether the judgment matrixes have consistency according to the consistency ratio, and taking indexes of the judgment matrixes with consistency as corresponding weight vectors.
5. The method for identifying key nodes of a power grid according to claim 1 or 3, wherein the calculation mode of the electrical betweenness center comprises:
Figure FDA0002835387680000021
wherein, PkiRepresenting the shortest or the longest passage between power node k and power node iThe sum of the powers of the active path transmissions; pikRepresents the sum of the power transmitted over the shortest or most efficient path between power node i and power node k, G represents the set of power generation nodes; l represents a load node set; cB(k) Representing the traditional betweenness center of the normalized power node k.
6. The method for identifying key nodes of a power grid according to claim 1, wherein the power node is checked in a manner that includes:
Figure FDA0002835387680000022
wherein ks (i) represents a core value of power node i; kminRepresenting the minimum value of the criticality of all nodes in the power grid; kmaxRepresenting the maximum value of criticality of all power nodes in the power grid; k (i) represents the criticality of power node i.
7. The method for identifying key nodes of a power grid according to claim 1, wherein the core values of all the power nodes are divided into a plurality of intervals, a minimum interval is determined, the power nodes in the minimum interval are removed, the power grid topology is updated, the core values of the remaining power nodes in the updated power grid topology are calculated, the steps of removing the power nodes and updating the calculation are repeated until the core values of all the power nodes are in the same interval, and all the power nodes in the interval are used as the identified key node set.
8. An identification device for a key node of a power grid, the identification device comprising:
the topological structure constructing module is used for constructing a power-efficient power grid topological structure according to validity rules of power node injection power and line transmission power in a power grid;
the criticality evaluation module is used for calculating the average electrical distance and the electrical betweenness center of the power nodes, respectively distributing weight vectors and constructing a criticality evaluation model;
the core value calculation module is used for calculating the core value of each power node according to the criticality evaluation model;
the interval dividing module is used for dividing a plurality of intervals according to the core values of the power nodes and determining the minimum interval;
the topology updating module is used for deleting the power nodes in the minimum interval, updating the power grid topology structure and returning to the core value calculating module;
and the node output module is used for outputting the power node of the last interval, namely the power node is the key point set.
9. An identification system for a critical node of a power grid, the identification system comprising:
the data acquisition module is used for acquiring the power data of the power nodes in the power grid;
a power grid key node identification apparatus as claimed in claim 8, configured to identify the collected power data of the power node and output a key point set.
CN202011468470.2A 2020-12-14 2020-12-14 Method, device and system for identifying key nodes of power grid Active CN112487658B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011468470.2A CN112487658B (en) 2020-12-14 2020-12-14 Method, device and system for identifying key nodes of power grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011468470.2A CN112487658B (en) 2020-12-14 2020-12-14 Method, device and system for identifying key nodes of power grid

Publications (2)

Publication Number Publication Date
CN112487658A true CN112487658A (en) 2021-03-12
CN112487658B CN112487658B (en) 2022-09-16

Family

ID=74916837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011468470.2A Active CN112487658B (en) 2020-12-14 2020-12-14 Method, device and system for identifying key nodes of power grid

Country Status (1)

Country Link
CN (1) CN112487658B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113507389A (en) * 2021-07-08 2021-10-15 华北电力大学 Power grid key node identification method and system
CN113536508A (en) * 2021-07-30 2021-10-22 齐鲁工业大学 Method and system for classifying manufacturing network nodes
CN113591937A (en) * 2021-07-09 2021-11-02 国家电网有限公司 Power system key node identification method based on local distance coding
CN114943430A (en) * 2022-05-13 2022-08-26 平顶山学院 Smart grid key node identification method based on deep reinforcement learning

Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110106321A1 (en) * 2009-11-03 2011-05-05 Spirae, Inc. Dynamic distributed power grid control system
US20140219133A1 (en) * 2013-02-05 2014-08-07 Cisco Technology, Inc. Proactive and selective time-stamping of packet headers based on quality of service experience and node location
CN104362638A (en) * 2014-12-04 2015-02-18 国家电网公司 Method for regulating voltage of key node of grid polymerization based on measurement of phasor measurement unit
GB201500943D0 (en) * 2015-01-20 2015-03-04 Tumbling Dice Ltd Improvements in and relating to network analysis
CN105303460A (en) * 2015-10-30 2016-02-03 国家电网公司 Identification method of key nodes and key branches in power grid
CN105426685A (en) * 2015-11-27 2016-03-23 云南电网有限责任公司电力科学研究院 Lightning flashover risk evaluation method for power system
CN105893740A (en) * 2016-03-29 2016-08-24 中国人民解放军国防科学技术大学 Method used for mining vulnerable node and circuit in power network
CN106779309A (en) * 2016-11-24 2017-05-31 华北电力大学(保定) The multi-level discrimination method of critical circuits multi-angle
CN106850254A (en) * 2016-12-20 2017-06-13 国网新疆电力公司信息通信公司 Key node recognition methods in a kind of power telecom network
CN107453357A (en) * 2017-08-24 2017-12-08 天津大学 A kind of State Estimation for Distribution Network based on hierarchical solving
CN107482626A (en) * 2017-08-17 2017-12-15 广东电网有限责任公司惠州供电局 A kind of regional power grid key node recognition methods
CN107622360A (en) * 2017-10-20 2018-01-23 广东电网有限责任公司电力调度控制中心 A kind of critical circuits recognition methods for considering subjective and objective factor
CN107947151A (en) * 2017-11-02 2018-04-20 广西电网有限责任公司电力科学研究院 A kind of recognition methods of electric system key node
US20180260561A1 (en) * 2017-03-08 2018-09-13 General Electric Company Generic framework to detect cyber threats in electric power grid
US20180262525A1 (en) * 2017-03-09 2018-09-13 General Electric Company Multi-modal, multi-disciplinary feature discovery to detect cyber threats in electric power grid
CN108879708A (en) * 2018-08-28 2018-11-23 东北大学 A kind of the reactive voltage partition method and system of active distribution network
CN109687438A (en) * 2018-12-19 2019-04-26 西南交通大学 It is a kind of meter and high-speed rail impact load effect under power grid vulnerable line discrimination method
CN109830963A (en) * 2019-01-08 2019-05-31 国网山东省电力公司滨州供电公司 A kind of power grid key node recognition methods counted and voltage class influences
CN110034581A (en) * 2019-04-23 2019-07-19 湘潭大学 The electrical betweenness vulnerability assessment method in the section of electric system under wind-electricity integration
CN110570094A (en) * 2019-08-13 2019-12-13 国网四川省电力公司成都供电公司 Method and system for evaluating strength of grid structure of medium-voltage distribution network
US20200233956A1 (en) * 2019-01-23 2020-07-23 General Electric Company Framework for cyber-physical system protection of electric vehicle charging stations and power grid
CN111697590A (en) * 2020-06-19 2020-09-22 上海交通大学 Entropy weight method-based power system key node identification method and system

Patent Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110106321A1 (en) * 2009-11-03 2011-05-05 Spirae, Inc. Dynamic distributed power grid control system
US20140219133A1 (en) * 2013-02-05 2014-08-07 Cisco Technology, Inc. Proactive and selective time-stamping of packet headers based on quality of service experience and node location
CN104362638A (en) * 2014-12-04 2015-02-18 国家电网公司 Method for regulating voltage of key node of grid polymerization based on measurement of phasor measurement unit
GB201500943D0 (en) * 2015-01-20 2015-03-04 Tumbling Dice Ltd Improvements in and relating to network analysis
CN105303460A (en) * 2015-10-30 2016-02-03 国家电网公司 Identification method of key nodes and key branches in power grid
CN105426685A (en) * 2015-11-27 2016-03-23 云南电网有限责任公司电力科学研究院 Lightning flashover risk evaluation method for power system
CN105893740A (en) * 2016-03-29 2016-08-24 中国人民解放军国防科学技术大学 Method used for mining vulnerable node and circuit in power network
CN106779309A (en) * 2016-11-24 2017-05-31 华北电力大学(保定) The multi-level discrimination method of critical circuits multi-angle
CN106850254A (en) * 2016-12-20 2017-06-13 国网新疆电力公司信息通信公司 Key node recognition methods in a kind of power telecom network
US20180260561A1 (en) * 2017-03-08 2018-09-13 General Electric Company Generic framework to detect cyber threats in electric power grid
US20180262525A1 (en) * 2017-03-09 2018-09-13 General Electric Company Multi-modal, multi-disciplinary feature discovery to detect cyber threats in electric power grid
CN107482626A (en) * 2017-08-17 2017-12-15 广东电网有限责任公司惠州供电局 A kind of regional power grid key node recognition methods
CN107453357A (en) * 2017-08-24 2017-12-08 天津大学 A kind of State Estimation for Distribution Network based on hierarchical solving
CN107622360A (en) * 2017-10-20 2018-01-23 广东电网有限责任公司电力调度控制中心 A kind of critical circuits recognition methods for considering subjective and objective factor
CN107947151A (en) * 2017-11-02 2018-04-20 广西电网有限责任公司电力科学研究院 A kind of recognition methods of electric system key node
CN108879708A (en) * 2018-08-28 2018-11-23 东北大学 A kind of the reactive voltage partition method and system of active distribution network
CN109687438A (en) * 2018-12-19 2019-04-26 西南交通大学 It is a kind of meter and high-speed rail impact load effect under power grid vulnerable line discrimination method
CN109830963A (en) * 2019-01-08 2019-05-31 国网山东省电力公司滨州供电公司 A kind of power grid key node recognition methods counted and voltage class influences
US20200233956A1 (en) * 2019-01-23 2020-07-23 General Electric Company Framework for cyber-physical system protection of electric vehicle charging stations and power grid
CN110034581A (en) * 2019-04-23 2019-07-19 湘潭大学 The electrical betweenness vulnerability assessment method in the section of electric system under wind-electricity integration
CN110570094A (en) * 2019-08-13 2019-12-13 国网四川省电力公司成都供电公司 Method and system for evaluating strength of grid structure of medium-voltage distribution network
CN111697590A (en) * 2020-06-19 2020-09-22 上海交通大学 Entropy weight method-based power system key node identification method and system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ZIO,E等: "Identifying groups of critical edges in a realistic electrical network by multi-objective genetic algorithms", 《RELIABILITY ENGINEERING & SYSTEM SAFETY》 *
徐毅非: "基于复杂网络理论的无功电压分区与关键线路辨识", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113507389A (en) * 2021-07-08 2021-10-15 华北电力大学 Power grid key node identification method and system
CN113507389B (en) * 2021-07-08 2022-05-03 华北电力大学 Power grid key node identification method and system
CN113591937A (en) * 2021-07-09 2021-11-02 国家电网有限公司 Power system key node identification method based on local distance coding
CN113591937B (en) * 2021-07-09 2023-09-26 国家电网有限公司 Power system key node identification method based on local distance coding
CN113536508A (en) * 2021-07-30 2021-10-22 齐鲁工业大学 Method and system for classifying manufacturing network nodes
CN113536508B (en) * 2021-07-30 2023-11-21 齐鲁工业大学 Manufacturing network node classification method and system
CN114943430A (en) * 2022-05-13 2022-08-26 平顶山学院 Smart grid key node identification method based on deep reinforcement learning

Also Published As

Publication number Publication date
CN112487658B (en) 2022-09-16

Similar Documents

Publication Publication Date Title
CN112487658B (en) Method, device and system for identifying key nodes of power grid
CN106850254B (en) Method for identifying key nodes in power communication network
CN102594909B (en) Multi-objective community detection method based on spectrum information of common neighbour matrix
CN107482626A (en) A kind of regional power grid key node recognition methods
CN104992266A (en) Method of determining power grid node importance degree and system thereof
CN110969347A (en) Power transmission network structure form evaluation method
CN110111024A (en) Scientific and technological achievement market valuation method based on AHP model of fuzzy synthetic evaluation
CN110417011A (en) A kind of online dynamic secure estimation method based on mutual information Yu iteration random forest
CN107947151B (en) Identification method for key nodes of power system
CN107133684A (en) A kind of random matrix construction method towards GA for reactive power optimization
Xiao et al. Network security situation prediction method based on MEA-BP
Li et al. An effective reliability evaluation method for power communication network based on community structure
CN105228185A (en) A kind of method for Fuzzy Redundancy node identities in identification communication network
CN112103950B (en) Power grid partitioning method based on improved GN splitting algorithm
CN111464327A (en) Spatial information network survivability evaluation method based on graph convolution network
Chabert‐Liddell et al. Impact of the mesoscale structure of a bipartite ecological interaction network on its robustness through a probabilistic modeling
CN115577292A (en) Vulnerable line identification method based on abstract dual network and cascading failure super network
Shen et al. A maintainability evaluation method for wireless sensor networks based on AHP and fuzzy comprehensive evaluation
Lu et al. EMH: Extended Mixing H-index centrality for identification important users in social networks based on neighborhood diversity
Sachan et al. Multiple correspondence analysis to study failures in a diverse population of a cable
CN108268982B (en) Large-scale active power distribution network decomposition strategy evaluation method and device
Zhang et al. Tensor-based spectral analysis of cascading failures over multilayer complex systems
Scardapane et al. Distributed spectral clustering based on Euclidean distance matrix completion
US20240054259A1 (en) Method for identifying fragile lines in power grids based on electrical betweenness
CN109784723A (en) Power transmission and transformation Project Risk Evaluation and terminal device

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