Method, system, medium and device for evaluating invulnerability of power system dynamic damage

Abstract

The application discloses a kind of power system dynamic damage invulnerability evaluation method, system, medium and equipment, including according to power grid real-time operation data, constructing power grid directed graph and adjacency matrix;According to several kinds of power grid damaged mode and the centrality index of each power grid node, the power grid directed graph and adjacency matrix are sequentially updated, until the invulnerability index corresponding to current power grid damaged mode does not comply with preset index condition, output current invulnerability index;According to the invulnerability index corresponding to each power grid damaged mode, output power system dynamic damage invulnerability evaluation result.The application introduces a kind of invulnerability index, and the power grid directed graph and the invulnerability index are dynamically updated according to different power grid damaged mode, realizes the invulnerability of current power grid in the case of power grid dynamic damage, dynamically evaluates.

Description

Technical Field

[0001] This invention belongs to the field of power system engineering and relates to a method, system, medium and equipment for evaluating the resilience of power system dynamic damage. Background Technology

[0002] With the development of large power grids, power grids are characterized by high voltage, large capacity, AC / DC hybrid operation, and long-distance power supply, resulting in increasingly complex grid structures. Under extreme events, power grids face the risk of short circuits and out-of-circuit failures. Faults at local nodes can lead to rapid load loss, triggering a surge in power at surrounding nodes, further expanding the fault range, and causing large-scale cascading power grid failures.

[0003] In existing technologies, complex network theory is used to graphically represent the power system, equating power generation, substations, and users as nodes, and transmission and distribution as edges. This allows for a macroscopic assessment of the power grid's damage status, or the construction of typical large-scale dynamic power flow models, such as AC optimal power flow models and hybrid search models combining transient stability and AC power flow, to achieve dynamic equivalence of power flow across large power grids. However, currently, there is still a lack of evaluation methods that can dynamically reflect the current resilience of the power grid under dynamic damage conditions. Summary of the Invention

[0004] To address the shortcomings of existing technologies, this application provides a method, system, medium, and equipment for evaluating the resilience of a power system under dynamic damage conditions, solving the problem of dynamically evaluating the current resilience of the power grid under dynamic damage conditions.

[0005] To achieve the above objectives, in a first aspect, the present invention provides a method for evaluating the resilience of a power system to dynamic damage, comprising:

[0006] Based on real-time power grid operation data, construct a directed graph and adjacency matrix for the power grid;

[0007] Based on several power grid damage modes and the centrality indices of each power grid node, the directed graph and adjacency matrix of the power grid are updated sequentially until the resilience index corresponding to the current power grid damage mode does not meet the preset index conditions, at which point the current resilience index is output. During each update, key nodes corresponding to the current power grid damage mode are removed from the current directed graph. The centrality indices of each power grid node are calculated based on the current directed graph and the adjacency matrix, and are used to determine the key nodes in the current directed graph. The resilience index is calculated based on the electrical parameters of each node in the updated directed graph and the maximum connected subgraph ratio.

[0008] Based on the resilience indicators corresponding to the damage modes of each power grid, the resilience evaluation results of dynamic damage to the power system are output.

[0009] Compared to existing technologies, the embodiments of this application have the following beneficial effects: By constructing a directed graph and adjacency matrix of the power grid, the topology and connection relationships of the power grid can be intuitively represented, providing basic data support for subsequent calculations and analyses; calculating centrality indices can help identify key nodes in the power grid, and screening out key nodes helps distinguish the importance of different nodes. After removing key nodes, the resilience indices are calculated, ensuring the relevance and effectiveness of the assessment, while also reducing the amount of computation and improving assessment efficiency; by simulating different damage modes and removing key nodes before calculating resilience indices, the resilience performance of the power grid under different impact conditions can be comprehensively assessed. At the same time, the directed graph and adjacency matrix of the power grid are dynamically updated to reflect changes in the power grid structure in a timely manner, ensuring the accuracy and timeliness of the assessment results; based on each resilience index, the resilience evaluation results of dynamic damage to the power system are output, which can comprehensively assess the resilience performance of the power grid and provide a scientific basis for the maintenance and optimization of the power grid.

[0010] In some embodiments of the first aspect of this application, constructing a directed graph and adjacency matrix of the power grid based on real-time power grid operation data includes:

[0011] Define a directed power graph as G = (V, E), where G represents the set of nodes V and line edges E, and satisfies the condition that node V... i ∈V, edge e i =(v i ,v j )∈E;

[0012] Construct the adjacency matrix:

[0013] Compared with existing technologies, the above embodiments have the following advantages: by constructing a directed graph and adjacency matrix of the power grid, the topology and connection relationship of the power grid can be represented intuitively, providing basic data support for subsequent calculations and analysis.

[0014] In some embodiments of the first aspect of this application, the centrality index of each power grid node is calculated based on the current directed graph of the power grid and the adjacency matrix, and is used to determine the key nodes in the current directed graph of the power grid, including:

[0015] The centrality index of each node is calculated based on the degree centrality index, proximity centrality index, center centrality index, and eigenvector centrality index of each node.

[0016] The centrality indices of each node are sorted from high to low, and the nodes corresponding to the top centrality indices are taken as the key nodes of the current directed graph of the power grid.

[0017] Compared with existing technologies, the above embodiments have the following beneficial effects: calculating centrality indices can help identify key nodes in the power grid, and screening out key nodes helps to distinguish the importance of different nodes.

[0018] In some embodiments of the first aspect of this application, the step of removing the corresponding critical node in the current directed graph of the power grid according to the current power grid damage mode during each update includes:

[0019] Based on the way the power grid is damaged, the corresponding critical nodes are removed from the current directed power grid graph. The removal method is as follows:

[0020] Represents node v i The way the damage occurs, when At that time, remove node v i ,when At that time, remove node v i and the corresponding edge e i .

[0021] Compared with existing technologies, the above embodiments have the following beneficial effects: based on the characteristics of different damage modes, the directed graph and adjacency matrix of the power grid are dynamically updated, reflecting changes in the power grid structure in a timely manner, and ensuring the accuracy and timeliness of the assessment results.

[0022] In some embodiments of the first aspect of this application, the output of the current survivability index includes:

[0023] The survivability index is calculated as follows:

[0024]

[0025] Where ψ represents the real-time transmission capacity ratio. Represents node v i The load loss ratio, Represents node v i Real-time transmission capacity, Represents node v i Historical transmission capacity Represents node v i The rated maximum load, T represents the number of nodes in the current directed graph of the power grid, P{k} represents the resilience index calculated in the Kth time, MCP represents the maximum connected subgraph proportion, and ω1, ω2, ω3 are the weights of the resilience index, determined by the damage mode. Sure.

[0026] Compared with existing technologies, the above embodiments have the following beneficial effects: outputting the current resilience index, gradually assessing the resilience of the power grid under different node damage conditions, and helping to understand the vulnerability and resilience of the power grid at different stages.

[0027] In some embodiments of the first aspect of this application, the resilience index corresponding to the current power grid damage mode does not meet the preset index conditions, including:

[0028] The conditions for the indicators are as follows:

[0029]

[0030] Where p{1→k} is the resilience index sequence obtained by merging the resilience indexes from the first to the Kth records, Ex(P{1→k}) represents the expected value of P{1→k}, and Std(P{1→k}) represents the standard deviation of P{1→k}.

[0031] If the current resilience index does not meet the specified conditions, it is determined that the current power grid's resilience does not meet the requirements after being hit, and the calculation is stopped.

[0032] Compared with the prior art, the above embodiments have the following beneficial effects: by setting preset index conditions, the resilience of the power grid under different node damage conditions can be gradually evaluated, the critical state of the power grid can be determined, and a reference basis can be provided for the formulation of protection measures for the power grid.

[0033] In some embodiments of the first aspect of this application, the centrality index of each power grid node is calculated based on the current directed graph of the power grid and the adjacency matrix, including:

[0034] The centrality index of each node is calculated based on its degree centrality, proximity centrality, center centrality, and eigenvector centrality; the calculation method is as follows:

[0035]

[0036] Where C d (v i ), J d (v i Z d (v i ), T d (v i ), They are nodes v i The degree centrality index, proximity centrality index, center centrality index, eigenvector centrality index, and centrality index are given, where N represents the number of nodes, and g(v j ,v k ) represents node vj and v k The number of shortest paths, n(v j ,v k ) represents node v j and v k The number of shortest paths passing through node v i The number of shortest paths, (n-1)(n-2) / 2 represents the number of paths that pass through node v. i The maximum possible value of the number of shortest paths, d(v i ,v j ) represents node v i and v j Euclidean distance, λ j θ represents the characteristic value, and θ represents the proportionality coefficient.

[0037] Compared with existing technologies, the above embodiments have the following advantages: calculating the centrality index can help identify key nodes in the power grid, screening out key nodes helps to distinguish the importance of different nodes, and then calculating the resilience index after removing key nodes ensures the relevance and effectiveness of the assessment, while also reducing the amount of calculation and improving the assessment efficiency.

[0038] Secondly, the present invention also provides a resilience evaluation system for dynamic damage of power systems, comprising: a graph construction module, an index output module, and a result output module.

[0039] The graph construction module is used to construct a directed graph and an adjacency matrix of the power grid based on real-time power grid operation data.

[0040] The indicator output module is used to sequentially update the directed graph of the power grid and the adjacency matrix based on several power grid damage modes and the centrality indicators of each power grid node, until the resilience indicator corresponding to the current power grid damage mode does not meet the preset indicator conditions, and then outputs the current resilience indicator. During each update, the corresponding key node in the current directed graph of the power grid is removed based on the current power grid damage mode. The centrality indicators of each power grid node are calculated based on the current directed graph of the power grid and the adjacency matrix, and are used to determine the key nodes in the current directed graph of the power grid. The resilience indicator is calculated based on the electrical parameters of each node in the updated directed graph of the power grid and the maximum connected subgraph ratio.

[0041] The result output module is used to output the resilience evaluation results of dynamic damage to the power system based on the resilience index corresponding to each power grid damage mode.

[0042] Compared with existing technologies, the above embodiments of this application have the following beneficial effects: By constructing a directed graph and adjacency matrix of the power grid, the topology and connection relationships of the power grid can be intuitively represented, providing basic data support for subsequent calculations and analyses; calculating centrality indices can help identify key nodes in the power grid, and screening out key nodes helps distinguish the importance of different nodes. After removing key nodes, the resilience indices are calculated, ensuring the relevance and effectiveness of the assessment, while also reducing the amount of computation and improving assessment efficiency; by simulating different damage modes and removing key nodes before calculating resilience indices, the resilience performance of the power grid under different impact conditions can be comprehensively assessed. At the same time, the directed graph and adjacency matrix of the power grid are dynamically updated to reflect changes in the power grid structure in a timely manner, ensuring the accuracy and timeliness of the assessment results; based on each resilience index, the resilience evaluation results of dynamic damage to the power system are output, which can comprehensively assess the resilience performance of the power grid and provide a scientific basis for the maintenance and optimization of the power grid.

[0043] Thirdly, the present invention also provides a device for evaluating the resilience of dynamic damage to a power system, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the computer program, when loaded onto the processor, implements the steps of the aforementioned method for evaluating the resilience of dynamic damage to a power system.

[0044] Fourthly, embodiments of this application also provide a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the aforementioned method for evaluating the resilience of dynamic damage to a power system. Attached Figure Description

[0045] Figure 1 This is a flowchart illustrating a method for evaluating the resilience of a power system to dynamic damage, provided in some embodiments of the present invention.

[0046] Figure 2 This is a schematic diagram of the structure of a power system dynamic damage resilience evaluation system provided in some embodiments of the present invention.

[0047] Figure 3 : This is a structural diagram of a power system dynamic damage resilience evaluation device provided in some embodiments of the present invention. Detailed Implementation

[0048] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0049] Example 1:

[0050] Please refer to Figure 1 The present invention provides a method for evaluating the resilience of a power system to dynamic damage, comprising steps S1 to S3:

[0051] Step S1: Construct a directed graph and adjacency matrix of the power grid based on real-time power grid operation data.

[0052] The real-time power grid operation data may include parameter information and actual operation information of each node and line, such as the real-time and historical transmission capacity of the nodes, node load loss, rated maximum load, etc.

[0053] In some embodiments of this application, step S1 includes steps S11 to S12, the specific steps of which are as follows:

[0054] S11: Define a directed power graph as G = (V, E), where G represents the set of nodes V and line edges E, and satisfies the condition that node V... i ∈V, edge e i =(v i ,v j )∈E;

[0055] S12: Construct the adjacency matrix:

[0056] Furthermore, when constructing a directed graph of the power grid, we can focus only on the transmission network above 35kV voltage level and directly connected nodes, ignoring the distribution network nodes and ground-connected branches below 35kV voltage level. This approach can pre-divide the backbone network, i.e. the transmission network above 35kV voltage level, before the subsequent key node screening, reducing the amount of calculation during key node screening and improving efficiency.

[0057] In this embodiment, by constructing a directed graph of the power grid and an adjacency matrix, the topology and connection relationships of the power grid can be intuitively represented, providing basic data support for subsequent calculations and analysis.

[0058] Step S2: Based on several power grid damage modes and the centrality index of each power grid node, sequentially update the directed graph of the power grid and the adjacency matrix until the resilience index corresponding to the current power grid damage mode does not meet the preset index conditions, and output the current resilience index; wherein, during each update, according to the current power grid damage mode, remove the corresponding key node in the current directed graph of the power grid; the centrality index of each power grid node is calculated based on the current directed graph of the power grid and the adjacency matrix, and is used to determine the key node in the current directed graph of the power grid; the resilience index is calculated based on the electrical parameters of each node in the updated directed graph of the power grid and the maximum connected subgraph ratio.

[0059] Preferably, in some embodiments of this application, the centrality index of each power grid node is calculated based on the current directed graph of the power grid and the adjacency matrix, and is used to determine the key nodes in the current directed graph of the power grid, specifically including steps S21 to S22:

[0060] S21: Calculate the centrality index of each node based on its degree centrality index, proximity centrality index, center centrality index, and eigenvector centrality index.

[0061] In step S21, after calculating the degree centrality index, proximity centrality index, center centrality index, and eigenvector centrality index of each node, the average value of these four indices can be taken to obtain the centrality index of the corresponding node.

[0062] S22: Sort the centrality indices of each node from high to low, and select the nodes corresponding to the top centrality indices as the key nodes of the current directed graph of the power grid.

[0063] In specific implementation, when selecting key nodes in step S22, the nodes ranked in the top 5% can be used as key nodes in the current directed graph of the power grid.

[0064] Furthermore, the centrality index of each node can be calculated according to the following method:

[0065]

[0066] Where C d (v i ), J d (v i Z d (v i ), T d (v i ), They are nodes v iThe degree centrality index, proximity centrality index, center centrality index, eigenvector centrality index, and centrality index are given, where N represents the number of nodes, and g(v j ,v k ) represents node v j and v k The number of shortest paths, n(v j ,v k ) represents node v j and v k The number of shortest paths passing through node v i The number of shortest paths, (n-1)(n-2) / 2 represents the number of paths that pass through node v. i The maximum possible value of the number of shortest paths, d(v i ,v j ) represents node v i and v j Euclidean distance, λ j θ represents the characteristic value, and θ represents the proportionality coefficient.

[0067] In this preferred embodiment, calculating the centrality index can help identify critical nodes in the power grid. Screening out critical nodes helps distinguish the importance of different nodes. After removing critical nodes, the resilience index is calculated to ensure the relevance and effectiveness of the assessment, while also reducing the amount of computation and improving the assessment efficiency.

[0068] Preferably, the removal of the corresponding critical node in the current directed graph of the power grid based on the current power grid damage mode during each update can be implemented through the following methods:

[0069] Based on the way the power grid is damaged, the corresponding critical nodes are removed from the current directed power grid graph. The removal method is as follows:

[0070] Represents node v i The way the damage occurs, when At that time, remove node v i ,when At that time, remove node v i and the corresponding edge e i .

[0071] In this preferred embodiment, the method of removing nodes is accurately characterized according to the characteristics of different damage modes, and the directed graph and adjacency matrix of the power grid are dynamically updated to reflect changes in the power grid structure in a timely manner, ensuring the accuracy and timeliness of the evaluation results.

[0072] Preferably, the current survivability index can be calculated in the following way:

[0073]

[0074] Where ψ represents the real-time transmission capacity ratio. Represents node v i The load loss ratio, Represents node v i Real-time transmission capacity, Represents node v i Historical transmission capacity Represents node v i The rated maximum load, T represents the number of nodes in the current directed graph of the power grid, P{k} represents the resilience index calculated in the Kth time, MCP represents the maximum connected subgraph proportion, and ω1, ω2, ω3 are the weights of the resilience index, determined by the damage mode. Sure.

[0075] For example, the weights of the resilience index can be allocated as follows:

[0076] The values ​​of ω1, ω2, and ω3 are all within the range of 0 to 1, and the sum of the three weights is 1.

[0077] In this preferred embodiment, after each power grid is attacked, critical nodes are removed, the directed graph and adjacency matrix of the power grid are updated, and the current resilience index is calculated and output. The resilience of the power grid under different node damage conditions is evaluated step by step, which helps to understand the vulnerability and resilience of the power grid at different stages.

[0078] Furthermore, the resilience index corresponding to the current power grid damage mode does not meet the preset index conditions, including:

[0079] The conditions for the indicators are as follows:

[0080]

[0081] Where p{1→k} is the resilience index sequence obtained by merging the resilience indexes from the first to the Kth records, Ex(P{1→k}) represents the expected value of P{1→k}, and Std(P{1→k}) represents the standard deviation of P{1→k}.

[0082] If the current resilience index does not meet the specified conditions, it is determined that the current power grid's resilience does not meet the requirements after being hit, and the calculation is stopped.

[0083] In this preferred embodiment, by setting preset index conditions, the resilience of the power grid under different node damage conditions is gradually evaluated, which can determine the critical state of the power grid and provide a reference for formulating protective measures for the power grid.

[0084] Step S3: Based on the resilience index corresponding to each power grid damage mode, output the resilience evaluation result of dynamic damage to the power system.

[0085] In this embodiment, based on various resilience indicators, the resilience evaluation results of dynamic damage to the power system are output, which can comprehensively assess the resilience performance of the power grid and provide a scientific basis for the maintenance and optimization of the power grid.

[0086] In summary, compared with the prior art, the above embodiments of this application have the following beneficial effects: By constructing a directed graph and adjacency matrix of the power grid, the topology and connection relationship of the power grid can be intuitively represented, providing basic data support for subsequent calculations and analyses; calculating the centrality index can help identify key nodes in the power grid, and screening out key nodes helps to distinguish the importance of different nodes. After removing key nodes, the resilience index is calculated to ensure the pertinence and effectiveness of the assessment, while also reducing the amount of computation and improving the assessment efficiency; by simulating different damage modes and removing key nodes before calculating the resilience index, the resilience performance of the power grid under different impact conditions can be comprehensively assessed. At the same time, the directed graph and adjacency matrix of the power grid are dynamically updated to reflect changes in the power grid structure in a timely manner, ensuring the accuracy and timeliness of the assessment results; based on each resilience index, the resilience evaluation results of dynamic damage to the power system are output, which can comprehensively assess the resilience performance of the power grid and provide a scientific basis for the maintenance and optimization of the power grid.

[0087] Example 2:

[0088] Please refer to Figure 2 Based on the same inventive concept, the present invention discloses a power system dynamic damage resilience evaluation system, comprising: a graph construction module M1, an index output module M2, and a result output module M3.

[0089] The graph construction module M1 is used to construct a directed graph and adjacency matrix of the power grid based on real-time power grid operation data.

[0090] The graph construction module M1 includes a directed graph construction unit and an adjacency matrix construction unit; the directed graph construction unit is used to construct a directed graph of the power grid, and the construction method is as follows:

[0091] Define a directed power graph as G = (V, E), where G represents the set of nodes V and line edges E, and satisfies the condition that node V... i ∈V, edge e i =(v i ,v j )∈E;

[0092] The adjacency matrix construction unit is used to construct the adjacency matrix, and the construction method is as follows:

[0093]

[0094] In this embodiment, the graph construction module M1 can intuitively represent the topology and connection relationships of the power grid by constructing a directed graph and adjacency matrix, providing basic data support for subsequent calculations and analysis.

[0095] The index output module M2 is used to sequentially update the directed graph of the power grid and the adjacency matrix based on several power grid damage modes and the centrality index of each power grid node, until the resilience index corresponding to the current power grid damage mode does not meet the preset index conditions, and then output the current resilience index. During each update, the corresponding key node in the current directed graph of the power grid is removed based on the current power grid damage mode. The centrality index of each power grid node is calculated based on the current directed graph of the power grid and the adjacency matrix, and is used to determine the key nodes in the current directed graph of the power grid. The resilience index is calculated based on the electrical parameters of each node in the updated directed graph of the power grid and the maximum connected subgraph ratio.

[0096] The indicator output module M2 includes: a centrality indicator update unit and a filtering unit.

[0097] The centrality index update unit is used to calculate the centrality index of each node based on the degree centrality index, proximity centrality index, intermediate centrality index and eigenvector centrality index of each node.

[0098] The filtering unit is used to sort the centrality index of each node from high to low, and the nodes corresponding to the top centrality indexes are used as the key nodes of the current directed graph of the power grid.

[0099] Furthermore, the centrality index of each node can be calculated according to the following method:

[0100]

[0101]

[0102] Where C d (v i ), J d (v i Z d (v i ), T d (v i ), They are nodes v i The degree centrality index, proximity centrality index, center centrality index, eigenvector centrality index, and centrality index are given, where N represents the number of nodes, and g(v j ,v k ) represents node v j and vk The number of shortest paths, n(v j ,v k ) represents node v j and v k The number of shortest paths passing through node v i The number of shortest paths, (n-1)(n-2) / 2 represents the number of paths that pass through node v. i The maximum possible value of the number of shortest paths, d(v i ,v j ) represents node v i and v j Euclidean distance, λ j θ represents the characteristic value, and θ represents the proportionality coefficient.

[0103] In this embodiment, the index output module M2 can help identify key nodes in the power grid by calculating the centrality index. Screening out key nodes helps to distinguish the importance of different nodes. After removing key nodes, the resilience index is calculated again to ensure the relevance and effectiveness of the assessment, while also reducing the amount of calculation and improving the assessment efficiency.

[0104] The indicator output module M2 also includes a graph update unit.

[0105] The graph update unit is used to remove corresponding key nodes from the current directed power grid graph according to the power grid damage method. The removal method is as follows:

[0106] Represents node v i The way the damage occurs, when At that time, remove node v i ,when At that time, remove node v i and the corresponding edge e i .

[0107] In this embodiment, the index output module M2 dynamically updates the directed graph and adjacency matrix of the power grid according to the characteristics of different damage modes, so as to reflect the changes in the power grid structure in a timely manner and ensure the accuracy and timeliness of the evaluation results.

[0108] The index output module M2 also includes: a survivability index calculation unit.

[0109] The survivability index calculation unit is used to calculate the survivability index, and the calculation method is as follows:

[0110]

[0111] Where ψ represents the real-time transmission capacity ratio. Represents node v i The load loss ratio, Represents node v i Real-time transmission capacity, Represents node v i Historical transmission capacity Represents node v i The rated maximum load, T represents the number of nodes in the current directed graph of the power grid, P{k} represents the resilience index calculated in the Kth time, MCP represents the maximum connected subgraph proportion, and ω1, ω2, ω3 are the weights of the resilience index, determined by the damage mode. Sure.

[0112] In this embodiment, the index output module M2 outputs the current resilience index to progressively assess the resilience of the power grid under different node damage conditions, which helps to understand the vulnerability and resilience of the power grid at different stages.

[0113] The indicator output module M2 also includes an indicator judgment unit.

[0114] The indicator judgment unit is used to determine whether the current survivability indicator meets the preset indicator conditions, wherein the indicator conditions are:

[0115]

[0116] Where p{1→k} is the resilience index sequence obtained by merging the resilience indexes from the first to the Kth records, Ex(P{1→k}) represents the expected value of P{1→k}, and Std(P{1→k}) represents the standard deviation of P{1→k}.

[0117] If the current resilience index does not meet the specified conditions, it is determined that the current power grid's resilience does not meet the requirements after being hit, and the calculation is stopped.

[0118] In this embodiment, the index output module M2 evaluates the resilience of the power grid under different node damage conditions by setting preset index conditions. It can determine the critical state of the power grid and provide a reference for formulating protective measures for the power grid.

[0119] The result output module M3 is used to output the resilience evaluation results of dynamic damage to the power system based on the resilience index corresponding to each power grid damage mode.

[0120] The result output module M3 in this embodiment outputs the resilience evaluation results of dynamic damage to the power system based on various resilience indicators, which can comprehensively evaluate the resilience performance of the power grid and provide a scientific basis for the maintenance and optimization of the power grid.

[0121] In summary, compared with existing technologies, the embodiments of this application have the following beneficial effects: By constructing a directed graph and adjacency matrix of the power grid, the topology and connection relationships of the power grid can be intuitively represented, providing basic data support for subsequent calculations and analyses; calculating centrality indices can help identify key nodes in the power grid, and screening out key nodes helps distinguish the importance of different nodes. After removing key nodes, the resilience indices are calculated again, ensuring the relevance and effectiveness of the assessment, while also reducing computational load and improving assessment efficiency; by simulating different damage modes and removing key nodes before calculating resilience indices, the resilience performance of the power grid under different impact conditions can be comprehensively assessed. Simultaneously, dynamically updating the directed graph and adjacency matrix of the power grid reflects changes in the power grid structure in a timely manner, ensuring the accuracy and timeliness of the assessment results; based on each resilience index, the resilience evaluation results of dynamic damage to the power system are output, which can comprehensively assess the resilience performance of the power grid and provide a scientific basis for power grid maintenance and optimization.

[0122] The above-described division of modules is only a logical functional division. In actual implementation, there may be other division methods, such as multiple modules being combined or integrated into another system.

[0123] Example 3:

[0124] Figure 3 A structural diagram of a power system dynamic damage resilience evaluation device according to this application is presented. (See diagram below.) Figure 3 As shown, the survivability assessment device may include: a processor N1, a memory N2, a data interface N3, and a communication bus N4.

[0125] Wherein: processor N1, memory N2, and data interface N3 communicate with each other through communication bus N4; data interface N3 is used for data communication with other devices such as input devices or output devices; processor N1 is used to execute program N5, specifically to execute the relevant steps in the above embodiment of the power system dynamic damage resilience evaluation method.

[0126] Specifically, program N5 may include program code, which includes computer-executable instructions.

[0127] The processor N1 may be a central processing unit (CPU), an application-specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of this application. The power system dynamic damage resilience assessment device includes one or more processors, which may be processors of the same type, such as one or more CPUs, or processors of different types, such as one or more CPUs and one or more ASICs.

[0128] Memory N2 is used to store program N5. Memory N2 may include high-speed RAM memory and may also include non-volatile memory, such as at least one disk storage device.

[0129] The algorithms or displays provided herein are not inherently related to any particular computer, virtual system, or other device. Furthermore, the embodiments in this application are not directed to any particular programming language.

[0130] Example 4:

[0131] This invention also provides a computer-readable storage medium storing at least one executable instruction that, when executed on a power system dynamic damage resilience assessment device / system, causes the power system dynamic damage resilience assessment device / system to perform the power system dynamic damage resilience assessment method in any of the above method embodiments.

[0132] Numerous specific details are set forth in the specification provided herein. However, it will be understood that embodiments of this application may be practiced without these specific details. Similarly, for the purpose of simplification and aiding understanding of one or more aspects of the invention, in the above description of exemplary embodiments of this application, various features of the embodiments are sometimes grouped together in a single embodiment, figure, or description thereof. The claims, which follow the detailed description, are hereby expressly incorporated into that detailed description, wherein each claim itself is a separate embodiment of this application.

[0133] Those skilled in the art will understand that the modules in the device of the embodiment can be adaptively changed and placed in one or more devices different from that embodiment. Modules, units, or components in the embodiment can be combined into a single module, unit, or component, and further, they can be divided into multiple sub-modules, sub-units, or sub-components, except that at least some of such features and / or processes or units are mutually exclusive.

Claims

1. A method for evaluating the resilience of a power system to dynamic damage, characterized in that, include: Based on real-time power grid operation data, construct a directed graph and adjacency matrix for the power grid; Based on several power grid damage modes and the centrality indices of each power grid node, the directed graph and adjacency matrix of the power grid are updated sequentially until the resilience index corresponding to the current power grid damage mode does not meet the preset index conditions, at which point the current resilience index is output. During each update, key nodes corresponding to the current power grid damage mode are removed from the current directed graph. The centrality indices of each power grid node are calculated based on the current directed graph and the adjacency matrix, and are used to determine the key nodes in the current directed graph. The resilience index is calculated based on the electrical parameters of each node in the updated directed graph and the maximum connected subgraph ratio. Based on the resilience indicators corresponding to the damage modes of each power grid, the resilience evaluation results of dynamic damage to the power system are output. The output of the current survivability index includes: The survivability index is calculated as follows: Where ψ represents the real-time transmission capacity ratio. Represents node v i The load loss ratio, Represents node v i Real-time transmission capacity, Represents node v i Historical transmission capacity Represents node v i The rated maximum load, T represents the number of nodes in the current directed graph of the power grid, P{k} represents the resilience index calculated in the Kth time, MCP represents the maximum connected subgraph proportion, and ω1, ω2, ω3 are the weights of the resilience index, which are determined by the damage mode.

2. The method for evaluating the resilience of a power system to dynamic damage as described in claim 1, characterized in that, The process of constructing a directed graph and adjacency matrix of the power grid based on real-time power grid operation data includes: Define a directed power graph as G = (V, E), where G represents the set of nodes V and line edges E, and satisfies the condition that node V... i ∈V, edge e i =(v i ,v j )∈E; Construct the adjacency matrix:

3. The method for evaluating the resilience of a power system to dynamic damage as described in claim 2, characterized in that, The centrality index of each power grid node is calculated based on the current directed graph of the power grid and the adjacency matrix, and is used to determine the key nodes in the current directed graph of the power grid, including: The centrality index of each node is calculated based on the degree centrality index, proximity centrality index, center centrality index, and eigenvector centrality index of each node. The centrality indices of each node are sorted from high to low, and the nodes corresponding to the top centrality indices are taken as the key nodes of the current directed graph of the power grid.

4. The method for evaluating the resilience of a power system to dynamic damage as described in claim 3, characterized in that, During each update, based on the current power grid damage method, the corresponding critical node in the current directed graph of the power grid is removed, including: Based on the way the power grid is damaged, the corresponding critical nodes are removed from the current directed power grid graph. The removal method is as follows: Represents node v i The way the damage occurs, when At that time, remove node v i ,when At that time, remove node v i and the corresponding edge e i .

5. The method for evaluating the resilience of a power system to dynamic damage as described in claim 4, characterized in that, The failure of the resilience index corresponding to the current power grid damage mode to meet the preset index conditions includes: The conditions for the indicators are as follows: Where p{1→k} is the resilience index sequence obtained by merging the resilience indexes from the first to the Kth records, Ex(P{1→k}) represents the expected value of P{1→k}, and Std(P{1→k}) represents the standard deviation of P{1→k}. If the current resilience index does not meet the specified conditions, it is determined that the current power grid's resilience does not meet the requirements after being hit, and the calculation is stopped.

6. A method for evaluating the resilience of a power system to dynamic damage as described in any one of claims 1 to 5, characterized in that, The centrality index of each power grid node is calculated based on the current directed graph of the power grid and the adjacency matrix, including: The centrality index of each node is calculated based on its degree centrality index, proximity centrality index, center centrality index, and eigenvector centrality index.

7. A resilience evaluation system for dynamic damage in a power system, characterized in that, include: The module includes a graph construction module, an indicator output module, and a result output module. The graph construction module is used to construct a directed graph and an adjacency matrix of the power grid based on real-time power grid operation data. The indicator output module is used to sequentially update the directed graph of the power grid and the adjacency matrix based on several power grid damage modes and the centrality indicators of each power grid node, until the resilience indicator corresponding to the current power grid damage mode does not meet the preset indicator conditions, and then outputs the current resilience indicator. During each update, the corresponding key node in the current directed graph of the power grid is removed based on the current power grid damage mode. The centrality indicators of each power grid node are calculated based on the current directed graph of the power grid and the adjacency matrix, and are used to determine the key nodes in the current directed graph of the power grid. The resilience indicator is calculated based on the electrical parameters of each node in the updated directed graph of the power grid and the maximum connected subgraph ratio. The result output module is used to output the resilience evaluation results of dynamic damage to the power system based on the resilience index corresponding to each power grid damage mode; The indicator output module includes: a survivability indicator calculation unit; The survivability index calculation unit is used to calculate the survivability index, and the calculation method is as follows: Where ψ represents the real-time transmission capacity ratio. Represents node v i The load loss ratio, Represents node v i Real-time transmission capacity, Represents node v i Historical transmission capacity Represents node v i The rated maximum load, T represents the number of nodes in the current directed graph of the power grid, P{k} represents the resilience index calculated in the Kth time, MCP represents the maximum connected subgraph proportion, and ω1, ω2, ω3 are the weights of the resilience index, which are determined by the damage mode.

8. A device for evaluating the resilience to dynamic damage of a power system, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the computer program is loaded into the processor, it implements the steps of the method for evaluating the resilience of dynamic damage to a power system according to any one of claims 1-6.

9. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by the processor, it implements the steps of the method for evaluating the resilience of dynamic damage to a power system according to any one of claims 1-6.

Images