A method for modeling and simulation analysis of power CPS interlocking failure and related equipment

Abstract

The application relates to the technical field of power information physical system, in particular to a power CPS interlocking fault modeling simulation analysis method and related equipment, which comprises the following steps: collecting voltage data, active power, reactive power of each power node in the power CPS at each collection time and data exchange amount between each information node; obtaining related power nodes of each power node; obtaining active power difference between each power node and other power nodes at each collection time; calculating the trend similarity of each power node at each collection time; then obtaining each power community and each information community; calculating the overlap degree of each information node at each collection time; and then re-dividing the power community and analyzing the interlocking fault of the power CPS. The application improves the analysis and processing efficiency of the cascading fault.

Description

Technical Field

[0001] This application relates to the field of power cyber-physical systems technology, specifically to a modeling and simulation analysis method and related equipment for power CPS interlocking faults. Background Technology

[0002] Power Cyber-Physical Systems (CPS) are systems that deeply integrate computing, communication, and control capabilities. Through the coordination of virtual networks in the information space and physical networks in the physical space, they enable the physical system to achieve higher reliability and security. Interlocking failures in power CPS arise from the coupling between power nodes and information nodes within the system. A failure in one power node can trigger failures in other power nodes, leading to large-scale power outages.

[0003] Because power nodes in a power CPS network possess community characteristics, interlocking faults can be analyzed and addressed based on these community characteristics to prevent large-scale power outages. However, traditional community partitioning methods are based on the degree of power nodes and the weights of edges in the network. Due to the varying distribution of power flows in a power CPS network, directly applying these methods can lead to deviations in community partitioning, resulting in inaccuracies in the analysis and handling of interlocking faults. Summary of the Invention

[0004] To address the aforementioned technical problems, the purpose of this application is to provide a modeling and simulation analysis method and related equipment for power CPS interlocking faults. The specific technical solution adopted is as follows:

[0005] In a first aspect, embodiments of this application provide a modeling and simulation analysis method for power CPS interlocking faults, the method comprising the following steps:

[0006] The system collects voltage data, active power, reactive power, and data exchange volume between information nodes at each data collection time for each power node in the power CPS. Each power node corresponds one-to-one with each information node, and each power node contains all devices within a region. Information nodes are responsible for transmitting information to each other.

[0007] Obtain the relevant power nodes for each power node;

[0008] Based on the difference between current flow direction and active power, the active power difference between each power node and other power nodes is obtained at each acquisition time.

[0009] Based on the active power difference between a power node and its related power nodes, and the differences in voltage data between power nodes, the similarity of the movement of each power node at each data acquisition time is obtained.

[0010] The community discovery algorithm is combined with the heading similarity to obtain each power community and each information community.

[0011] Based on the heading similarity, the distance and the exchange amount between information nodes, the overlap degree of each information node at each collection time is obtained.

[0012] Based on the overlap degree of each information node in each information community, the power community is re-divided, and the interlocking fault of the power CPS is analyzed based on the re-divided power community.

[0013] Further, the method for obtaining the data exchange amount between the information nodes is that the data exchange amount between any two information nodes within one second at each collection time is counted as the data exchange amount at each collection time.

[0014] Further, the method for obtaining the related power nodes is that the other power nodes whose number of power nodes passed on the position connection line between each power node and the other power nodes is less than or equal to a preset number are taken as the related power nodes of each power node.

[0015] Further, the method for obtaining the active difference value is that:

[0016] For each power node at each collection time, if the current direction flows from each power node to other power nodes, the difference value between the active power of each power node and the other power nodes is taken as the active difference value between each power node and the other power nodes; if the current direction flows from the other power nodes to each power node, the difference value between the active power of each power node and the other power nodes is taken as the active difference value between each power node and the other power nodes.

[0017] Further, the method for obtaining the heading similarity is that:

[0018] For each power node at each collection time, the ratio of the active difference value between the power node and its related power nodes to the sum of the active power of all power nodes in the power grid is taken as the first ratio between the power node and its related power nodes, and the sum of all first ratios is taken as the centrality coefficient of each power node at each collection time.

[0019] The voltage similarity of each power node at each collection time is obtained based on the difference between the voltage data of the power node.

[0020] The normalized value of the product of the centrality coefficient and the voltage similarity is taken as the heading similarity of each power node at each collection time.

[0021] Further, the calculation formula of the voltage similarity is: ; in the formula, a voltage similarity of the i-th power node at each collection time point; N represents a number of relevant power nodes of the i-th power node in the power grid, 、 respectively represent voltage data of the i-th and j-th power nodes at each collection time point, 、 respectively represent a maximum value of voltage data and a minimum value of voltage data among all relevant power nodes of the i-th power node at each collection time point.

[0022] Further, the acquisition method of each power community and each information community is as follows:

[0023] The degree of a node in a community discovery algorithm is replaced by the trend similarity of a power node, and the community discovery algorithm is used to divide the power nodes into communities, to obtain each community of the power nodes at each collection time point as each power community, and to divide the information nodes corresponding to the power nodes into communities according to the division of the power node communities, to obtain each community of the information nodes at each collection time point as each information community.

[0024] Further, the acquisition method of the overlap degree is as follows:

[0025] For each power community at each collection time point, the power node with the maximum trend similarity in the power community is obtained as the first power node of each power community at each collection time point, and the corresponding information node of the first power node is obtained as the first information node of the information community in which the corresponding information node is located at each collection time point.

[0026] For each collection time point, the two first information nodes closest to each information node are obtained, the distance between each information node and one of the two first information nodes closest to the information node is obtained as the first nearest distance of each information node at each collection time point, and the distance between each information node and the other first information node is obtained as the second nearest distance of each information node at each collection time point.

[0027] ; in the formula, represents the overlap degree of the i-th information node at each collection time point; 、 respectively represent the first nearest distance and the second nearest distance of the i-th information node at each collection time point, and M represents a number of all information nodes directly connected to the i-th information node, where the direct connection means that the information nodes are connected to the i-th information node and there is no other information node between the i-th information node and the information nodes, represents the data exchange amount between the i-th information node and the j-th information node at each collection time point.

[0028] Further, the power community is re-divided, and interlocking failures of the power CPS are analyzed based on the re-divided power community, including:

[0029] For each information node at each collection time, a distance between the two first information nodes closest to the information node is calculated as a first distance, when the first distance is less than a first preset threshold, an information node with the maximum overlap in the two information communities where the two first information nodes closest to the information node are located is taken as an overlapping information node, and the overlapping information node is taken as an information node in the two information communities at the same time, the information communities at each collection time are re-divided, and the re-divided information communities are mapped to the power nodes to obtain the re-divided power communities at each collection time.

[0030] For each power node, if the voltage of the power node drops to a preset value or the voltage is greater than a preset multiple of the rated voltage, the connection between the power community where the fault power node is located and other power nodes is disconnected, and if the overlapping power node fails, the connection between the overlapping power node and all power nodes in the power community where the overlapping power node is located is disconnected.

[0031] In a second aspect, an embodiment of the present application provides a related device for modeling and simulation analysis of interlocking failures of a power CPS, and the related device comprises:

[0032] A data collection module is configured to collect voltage data, active power, reactive power of each power node in the power CPS at each collection time, and data exchange amount between each information node, each power node corresponds to each information node one by one, the power node contains all devices in a region, and the information nodes are responsible for transmitting information;

[0033] An overlap degree acquisition module is configured to acquire related power nodes of each power node, acquire active difference values between each power node and other power nodes at each collection time based on the difference between the current flow direction and the active power, acquire trend similarity of each power node at each collection time based on the active difference values between the power node and the related power nodes of the power node and the difference between the voltage data of the power node, acquire each power community and each information community by using a community discovery algorithm combined with the trend similarity, and acquire overlap degrees of each information node at each collection time based on the trend similarity, the distance between the information nodes, and the exchange amount.

[0034] A power CPS interlocking failure analysis module is configured to re-divide the power community based on the overlap degrees of each information node in each information community, and analyze interlocking failures of the power CPS based on the re-divided power community.

[0035] The present application has at least the following beneficial effects:

[0036] Firstly, the application obtains the relevant power nodes of each power node according to the number of other power nodes between the power nodes in the power CPS model, obtains the active power difference value according to the difference of the active power between the power nodes, obtains the centrality coefficient of each power node according to the active power difference value between each power node and the relevant power nodes, reflects the key degree of the position of each power node, obtains the voltage similarity of each power node according to the difference between the voltage data of each power node, reflects the similarity degree of each power node and the nearby power nodes, and then obtains the trend similarity of each power node according to the centrality coefficient and the voltage similarity, improves the Louvain algorithm according to the trend similarity, and considers the relationship between power and voltage when community division; further, the overlap degree of each information node is obtained according to the trend similarity, the community is re-divided considering the community overlap, the interlocking fault of the power CPS is analyzed according to the re-divided community, the calculation of the modularity in the traditional Louvain algorithm is improved, the distribution of the power CPS network is considered at the same time based on the degree of the power node in the network, so that the community divided by the improved method is more in line with the actual situation; the community division is further optimized considering the overlap in the power CPS network in the actual situation, and the analysis and processing efficiency of the cascading failure is improved. BRIEF DESCRIPTION OF DRAWINGS

[0037] In order to more clearly illustrate the technical solutions and advantages of the embodiments of the present application or the prior art, the drawings needed to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are only some embodiments of the present application, and those skilled in the art can obtain other drawings according to these drawings without creating any creative labor.

[0038] Figure 1 A step flow chart of a power CPS interlocking fault modeling simulation analysis method provided by an embodiment of the present application is shown in the figure.

[0039] Figure 2 A power CPS model schematic diagram provided by an embodiment of the present application is shown in the figure. DETAILED DESCRIPTION

[0040] In order to further illustrate the technical means and effects adopted by the present application to achieve the predetermined invention purpose, the following describes the specific implementation, structure, features and effects of a power CPS interlocking fault modeling simulation analysis method and related equipment according to the present application in combination with the drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. In addition, the specific features, structures or characteristics in one or more embodiments can be combined in any suitable form.

[0041] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.

[0042] The application provides a modeling simulation analysis method for power CPS interlocking failure and related equipment.

[0043] Please refer to Figure 1 , which shows a step flow chart of a modeling simulation analysis method for power CPS interlocking failure provided by an embodiment of the application, which includes the following steps:

[0044] Step S001, collect the voltage data, active power, reactive power of each power node in the power CPS at each collection time, and the data exchange amount between each information node.

[0045] The power CPS includes a power network and an information network. The power network mainly includes power stations, distribution stations, substations, users, power transmission lines, sensors, and communication equipment. All devices in a region are regarded as a power node, and the connection between each power node represents the power transmission line between each region. The information network is a virtual network mainly composed of a dispatching center and each information node. The connection between each information node represents an information transmission line, is responsible for transmitting the collected information of the power node and the control signal issued by the control center, and the connection can be a wireless network or a wired network. Each power node corresponds to an information node.

[0046] The power CPS model schematic diagram is shown in Figure 2 , wherein the solid circle represents the power node, the connection between the power nodes represents the power transmission line, and the set of the power nodes and the connection represents the power network; the hollow circle represents the information node, the connection between the information nodes represents the information transmission line, and the set of the information nodes and the connection represents the information network. The connection between the power node and the information node represents the coupling relationship therebetween, and the power node and the information node are one-to-one coupled, but the power transmission line and the information transmission line are not necessarily corresponding.

[0047] The voltage sensor is used to measure the voltage data of each power node at each collection time, and the power sensor is used to measure the active power and the reactive power of each power node at each collection time. The measurement interval of the collection time is 1s.

[0048] Further, the data exchange amount between any two information nodes within one second at each collection time is counted as the data exchange amount at each collection time. In this embodiment, the incremental data is compared by time stamp, and the data exchange amount is recorded.

[0049] In step S002, the related power nodes of each power node are obtained; based on the difference between the flow direction of the current and the active power, the active difference between each power node and other power nodes at each collection time is obtained; based on the active difference between the power node and its related power nodes and the difference between the voltage data of the power node, the trend similarity of each power node at each collection time is obtained; the community discovery algorithm is used in combination with the trend similarity to obtain each power community and each information community; based on the trend similarity, the distance and exchange amount between information nodes, the overlap degree of each information node at each collection time is obtained; and the power community is re-divided based on the overlap degree of each information node in each information community.

[0050] The interlocking fault of the power CPS refers to a series of chain reactions triggered by a single or multiple initial power node faults, causing the failure of multiple power nodes, propagating along the power transmission line, and eventually possibly leading to a large-area power outage accident, also known as a cascading failure. Some power nodes are closely related to other power nodes, and such power nodes are more likely to cause interlocking faults. Some power nodes are sparsely connected to other power nodes and are not easily affected by other power nodes. Therefore, the community structure between power nodes has an important influence on the propagation of interlocking faults, that is, interlocking faults are more likely to occur between community power nodes, and interlocking faults across communities need a long time to propagate.

[0051] Since the power grid is constructed based on the construction principle of hierarchical zoning, the community structure presents significant community structure characteristics from the regional perspective. The community structure refers to an organization composed of power nodes with similar functions or structures, and the connection between the power nodes in the community is close, while the connection between the power nodes in different communities is loose.

[0052] In order to prevent interlocking faults from causing large-area power outage accidents, the community where the power node is located is identified at an early stage of the failure of the power node, and the power nodes in the community are isolated and disconnected from other communities, so that the interlocking fault can be prevented from propagating to other communities through key power nodes, causing more power node failures.

[0053] Traditional community division algorithms, such as the Louvain method, have the advantages of high efficiency suitable for large-scale networks and can identify hierarchical community structures. This method identifies communities through modularity, and the calculation of modularity depends on the degree of the power nodes in the network. However, the trend of each branch in the power grid is not the same, so directly using the Louvain method to divide the community will differ from the actual situation. The present application divides the community based on the degree of the power node, which is closer to the actual situation.

[0054] The connectivity of each power node in the power grid is different. Some power nodes are connected with multiple power nodes, are in a more critical position, and are more suitable as the central power node of the community. Some power nodes are connected with only a few power nodes, and are more likely to be the edge power node in the community. In the power grid, the degree of the power node connected with more other power nodes is larger, but due to the different power distribution of each power node and the line between the power nodes, the power node with the largest degree is not necessarily the central power node. In order to make the division of the community structure of the power grid power node more in line with the actual situation, in addition to considering the degree of the power node, the distribution of the direction is considered to calculate the centrality coefficient, and the voltage similarity is calculated by considering the relationship between the power nodes.

[0055] The direction central power node in the power grid means that the importance of the power node in the power transmission process is larger, and is more likely to be the central power node of a community.

[0056] For each power node in the power grid, the more active power between the other power nodes connected with the power node, the more likely the power node is a direction center; further, in order to ensure the community nature and avoid too large calculation amount, the other power nodes with the number of power nodes between each power node less than or equal to 5 are regarded as the related power nodes of each power node.

[0057] For each power node at each collection time, if the current direction flows from each power node to other power nodes, the difference between the active power of other power nodes and each power node is regarded as the active difference between each power node and other power nodes; if the current direction flows from other power nodes to each power node, the difference between the active power of each power node and other power nodes is regarded as the active difference between each power node and other power nodes.

[0058] The calculation formula of the centrality coefficient of each power node is: ; in the formula, Ci represents the centrality coefficient of the i th power node at each collection time; j represents the j th related power node of the i th power node, and N represents the number of related power nodes of the i th power node in the power grid, represents the active difference between the i th power node and the j th related power node at each collection time, represents the sum of the active power of all power nodes in the power grid at each collection time.

[0059] It should be noted that for each power node, the larger the centrality coefficient, the more critical the position of the power node, and the more likely to be located in the center of the community. Conversely, the power node is more likely to be located at the edge of the community. The calculation of the modularity according to the centrality coefficient is more in line with the division of the community in the actual situation.

[0060] At least two power nodes are contained in a community structure in a power grid, and the power nodes in the community generally have similar properties, the voltage of two power nodes closer in distance is more similar, and they are more likely to be in the same community. The voltage similarity between power nodes represents the similarity between the two power nodes, and the calculation formula of the voltage similarity between each power node and other power nodes is: ; wherein, represents the voltage similarity of the i th power node at each collection time; N represents the number of relevant power nodes of the i th power node in the power grid, , respectively represent the voltage data of the i th and j th power nodes at each collection time, , respectively represent the maximum value and the minimum value of the voltage data of all relevant power nodes of the i th power node at each collection time.

[0061] It should be noted that the greater the voltage similarity, the higher the similarity between the power node and the surrounding power nodes, and the more likely they belong to a community. On the contrary, the lower the similarity between the power node and the surrounding power nodes, the less likely they belong to a community.

[0062] Further, according to the centrality coefficient and the voltage similarity of each power node, the trend similarity of each power node is obtained, and the calculation formula is: ; wherein, represents the trend similarity of the i th power node at each collection time; norm() represents a normalization function, , respectively represent the centrality coefficient and the voltage similarity of the i th power node at each collection time.

[0063] It should be noted that when the centrality coefficient and the voltage similarity of the power node are greater, the trend similarity of the power node is greater; otherwise, the trend similarity of the power node is smaller.

[0064] Further, the calculation formula of the trend similarity is brought into the calculation of the modularity, and the degree of the node in the Louvain algorithm is replaced by the trend similarity of the power node, so that the power and voltage relationship between the power nodes is considered when the community is divided.

[0065] The improved modularity calculation method is used to divide the power nodes into communities by using the Louvain algorithm to obtain each community of the power nodes at each collection time as each power community. According to the division of the power node communities, the information nodes corresponding to the power nodes are divided into communities to obtain each community of the information nodes at each collection time as each information community. It should be noted that the process of dividing the communities by using the Louvain algorithm is a known technology, and the specific steps are not repeated.

[0066] The communities divided by using the Louvain method are isolated, that is, there is no power node belonging to two communities at the same time, but in the actual situation, there is a phenomenon of community overlap. For example, the substation power node located between two communities can communicate with two communities at the same time, but it is not closely connected with other power nodes in the community. At this time, it can be considered that the two communities overlap, and the substation power node belongs to two communities at the same time, that is, the substation power node is an overlapping power node. When isolating the community where the fault occurs, the connection of the overlapping power node should not be disconnected, because the overlapping power node is a key power node connecting the community and other power nodes, and disconnecting the overlapping power node may cause the community where the fault does not occur to lose connection.

[0067] For each power community at each collection time, the power node with the maximum similarity in the community is obtained as the first power node of each power community at each collection time, and the corresponding information node of the first power node is obtained as the first information node of the information community where the corresponding information node is located at each collection time.

[0068] For each collection time, the two first information nodes closest to each information node are obtained, the distance between each information node and the closest first information node is obtained as the first nearest distance of each information node at each collection time, and the distance between each information node and the other first information node is obtained as the second nearest distance of each information node at each collection time.

[0069] In order to determine the overlapping power node, the community divided by the Louvain method is further optimized, and the information of the information node in the information network coupled with the power network is used for confirmation. Since the overlapping power node is usually located at the edge of the community far from the community center, and because of the communication between multiple communities, the data exchange amount is large. According to the above analysis, for each information node, the overlap degree is calculated to represent the probability that the information node is an overlapping power node. The calculation formula of the overlap degree of each information node is: ; in the formula, represents the overlap degree of the i-th information node at each collection time; , respectively represent the first nearest distance and the second nearest distance of the i-th information node at each collection time, M represents the number of all information nodes directly connected to the i-th information node, and directly connected means that the information node is connected to the i-th information node and there is no other information node between the i-th information node, represents the data exchange amount of the i-th information node and the j-th information node at each collection time.

[0070] It should be noted that in the calculation formula of the overlap degree Therefore The closer to 1, the more likely the i-th information node is at the edge of the community, i.e. the greater the possibility of belonging to the overlapping power node; and the overlapping power node communicates multiple communities at the same time, so the greater the average data exchange amount of the i-th information node, the more likely it is an overlapping power node.

[0071] For each information node, the greater the overlap degree of the information node, the more likely an information node is at the edge of the community and the greater the data exchange amount, the more likely the information node is an overlapping power node and should belong to multiple communities at the same time; on the contrary, the less likely the information node is an overlapping power node.

[0072] For each information node at each collection time, the distance between the two first information nodes closest to the information node is calculated as the first distance, and when the first distance is less than a first preset threshold, the information node with the maximum overlap degree in the two information communities where the two first information nodes closest to the information node are located is taken as an overlapping information node, and the overlapping information node is taken as an information node in the two information communities at the same time. In this way, the redivision of the information community at each collection time is completed; further, the redivided information community is mapped to the power node to obtain the redivided power community at each collection time.

[0073] At this point, the redivision of the power community at each collection time is completed.

[0074] Step S003, based on the redivided power community, the interlocking fault of the power CPS is analyzed.

[0075] When the power grid is in normal operation, the trend between each power node and the transmission line should be balanced, i.e. the generated power is equal to the sum of the power consumption and the transmission loss, and when a certain power node fails, the trend balance will be broken, causing other power nodes or transmission lines to operate overload, causing interlocking failure.

[0076] Further, when the sensor monitors that the voltage of the power node drops to 0 or the voltage is greater than 1.2 times of the rated voltage, the connection between the power community where the fault power node is located and other power nodes is disconnected, the generation power is adjusted through the information network to recalculate the trend to maintain balance and avoid the failure of the power nodes in other communities.

[0077] According to the power community division in S002, when the power node fails, the power connection between the power node and the power nodes in the power community where the power node is located except the overlapping power nodes is disconnected, because the overlapping power nodes communicate multiple communities; if the overlapping power node fails, the connection between it and the power nodes in all the power communities where it is located is disconnected, to avoid interlocking failure and cause large-area power outage. Thus, a modeling simulation analysis method for power CPS interlocking failure is completed.

[0078] Based on the same inventive concept as the above method, the embodiments of the present application provide a related device for modeling simulation analysis of power CPS interlocking failure, the related device comprises:

[0079] a data acquisition module, configured to acquire voltage data, active power, reactive power of each power node in the power CPS at each acquisition time, and data exchange amount between each information node; each power node corresponds to each information node one by one, the power node contains all devices in a region, and the information nodes are responsible for transmitting information;

[0080] an overlapping degree acquisition module, configured to acquire related power nodes of each power node; based on the difference between the current flow direction and the active power, an active difference value between each power node and other power nodes at each acquisition time is acquired; based on the active difference value between the power node and its related power nodes and the difference between the voltage data of the power node, a trend similarity of each power node at each acquisition time is acquired; a community discovery algorithm is combined with the trend similarity to acquire each power community and each information community; based on the trend similarity, the distance and the exchange amount between the information nodes, an overlapping degree of each information node at each acquisition time is acquired;

[0081] a power CPS interlocking failure analysis module, configured to redivide the power communities based on the overlapping degree of each information node in each information community, and analyze the interlocking failure of the power CPS based on the redivided power communities.

[0082] It should be noted that the above-mentioned embodiments of the present application are in the order of description only, and do not represent the advantages and disadvantages of the embodiments. Moreover, the above describes specific embodiments of the present application. In addition, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multi-task processing and parallel processing are possible or can be advantageous.

[0083] The various embodiments in the present application are described in a progressive manner, and the same or similar parts between the various embodiments can be referred to each other, and each embodiment focuses on the difference from other embodiments.

[0084] The above description is merely the preferred embodiments of the present application, and is not intended to limit the present application. Any modification, equivalent replacement, improvement, etc. made within the principles of the present application shall be included in the protection scope of the present application.

Claims

1. A method for modeling and simulation analysis of power CPS interlocking failure, characterized in that, The method comprises the following steps: Collecting voltage data, active power, reactive power of each power node in the power CPS at each collection time, and data exchange amount between each information node; each power node corresponds to each information node, and the power node comprises all devices in a region, and the information nodes are responsible for transmitting information; Obtaining relevant power nodes of each power node; Based on the difference between the flow direction of the current and the active power, obtaining the active difference value between each power node and other power nodes at each collection time; Based on the active difference value between the power node and its relevant power node, and the difference between the voltage data of the power node, obtaining the trend similarity of each power node at each collection time; Using a community discovery algorithm combined with the trend similarity to obtain each power community and each information community; Based on the trend similarity, the distance and the exchange amount between the information nodes, obtaining the overlap degree of each information node at each collection time; Based on the overlap degree of each information node in each information community, re-dividing the power community, and analyzing the interlocking failure of the power CPS based on the re-divided power community; The method for obtaining the trend similarity is as follows: For each power node at each collection time, the ratio of the active difference value between the power node and its relevant power node to the sum of the active power of all power nodes in the power grid is taken as the first ratio between the power node and its relevant power node, and the sum of all first ratios is taken as the centrality coefficient of each power node at each collection time; Based on the difference between the voltage data of the power node, obtaining the voltage similarity of each power node at each collection time; The normalized value of the product of the centrality coefficient and the voltage similarity is taken as the trend similarity of each power node at each collection time; The method for obtaining the overlap degree is as follows: For each power community at each collection time, the power node with the maximum trend similarity in each power community is taken as the first power node of each power community at each collection time, and the corresponding information node of the first power node is taken as the first information node of the information community where the corresponding information node is located at each collection time; For each collection time, the two first information nodes closest to each information node are obtained, the distance between each information node and the closest first information node is taken as the first nearest distance of each information node at each collection time, and the distance between each information node and the other first information node is taken as the second nearest distance of each information node at each collection time; ; wherein represents the overlap degree of the i-th information node at each collection time point; , respectively represent the first nearest distance and the second nearest distance of the i-th information node at each collection time point, and M represents the number of all information nodes directly connected to the i-th information node, wherein directly connected means connected to the i-th information node and without other information nodes between the i-th information node and the directly connected information node, represents the data exchange amount between the i-th information node and the j-th information node at each collection time point.

2. The method of claim 1, wherein the power CPS interlock fault is modeled as a series of voltage and current sources. The method for obtaining the data exchange amount between the information nodes is as follows: the data exchange amount between any two information nodes within one second before each collection time is taken as the data exchange amount at each collection time.

3. The method of claim 1, wherein the power CPS interlock fault is modeled as a series of voltage and current sources. The method for obtaining the relevant power node is as follows: other power nodes with the number of power nodes passing through the position connection line between each power node being less than or equal to a preset number are taken as the relevant power nodes of each power node.

4. The method of claim 1, wherein the power CPS interlock fault is modeled as a series of voltage and current sources. The method for obtaining the active difference value is as follows: If the current direction is from each power node to other power nodes, the difference between the active power of each power node and other power nodes is taken as the active difference between each power node and other power nodes; If the current direction is from other power nodes to each power node, the difference between the active power of each power node and other power nodes is taken as the active difference between each power node and other power nodes.

5. The method of claim 1, wherein, The calculation formula of the voltage similarity is: ; in the formula, represents the voltage similarity of the i-th power node at each collection moment. N represents the number of power nodes related to the i-th power node in the power grid. , Let represent the voltage data of the i-th and j-th power nodes at each acquisition time, respectively. , These represent the maximum and minimum voltage values ​​of the i-th power node at each acquisition time, respectively.

6. The method of claim 1, wherein, The acquisition method of the power community and the information community is: The degree of the node in the community discovery algorithm is replaced by the trend similarity of the power node, the community division of the power node is performed by using the community discovery algorithm, each community of the power node at each collection time is acquired as each power community, and the information node corresponding to the power node is divided into communities according to the division of the community of the power node, and each community of the information node at each collection time is acquired as each information community.

7. The method of claim 1, wherein the power CPS interlock fault is modeled as a series of voltage and current sources. The power community is re-divided, and the interlocking fault of the power CPS is analyzed based on the re-divided power community, including: For each information node at each collection time, the distance between the two first information nodes closest to the information node is taken as the first distance, when the first distance is less than a first preset threshold, the information node with the maximum overlap degree in the two information communities where the two first information nodes closest to the information node are located is taken as an overlapping information node, the overlapping information node is taken as an information node in the two information communities at the same time, the information community at each collection time is re-divided, and the re-divided information community is mapped to the power node to acquire the re-divided power community at each collection time; For each power node, if the voltage of the power node drops to a preset value or the voltage is greater than a preset multiple of the rated voltage, the connection between the power community where the fault power node is located and other power nodes is disconnected, and if the overlapping power node fails, the connection between the overlapping power node and the power nodes in all power communities where the overlapping power node is located is disconnected.

8. A device for modeling and simulation analysis of power CPS interlocking faults, realizing the method for modeling and simulation analysis of power CPS interlocking faults according to any one of claims 1-7, characterized in that, The related device includes: A data collection module is configured to collect voltage data, active power, reactive power of each power node in the power CPS at each collection time, and data exchange amount between each information node; each power node corresponds to each information node one by one, the power node includes all devices in a region, and the information nodes are responsible for transmitting information; An overlap degree acquisition module is configured to acquire related power nodes of each power node, acquire active difference between each power node and other power nodes at each collection time based on the difference between the current direction and the active power, acquire trend similarity of each power node at each collection time based on the active difference between the power node and the related power node of the power node and the difference between the voltage data of the power node, acquire each power community and each information community by using a community discovery algorithm combined with the trend similarity, and acquire overlap degree of each information node at each collection time based on the trend similarity, the distance between the information nodes and the exchange amount. The power CPS interlocking fault analysis module is used for re-dividing the power community based on the overlapping degree of each information node in each information community, and analyzing the interlocking fault of the power CPS based on the re-divided power community.

Images