Research Article  Open Access
XinRui Liu, Yuan Meng, Peng Chang, "Node Importance Evaluation of CyberPhysical System under CyberAttacks Spreading", Complexity, vol. 2021, Article ID 6641030, 15 pages, 2021. https://doi.org/10.1155/2021/6641030
Node Importance Evaluation of CyberPhysical System under CyberAttacks Spreading
Abstract
The study of cyberattacks, and in particular the spread of attack on the power cyberphysical system, has recently attracted considerable attention. Identifying and evaluating the important nodes under the cyberattack propagation scenario are of great significance for improving the reliability and survivability of the power system. In this paper, we improve the closeness centrality algorithm and propose a compound centrality algorithm based on adaptive coefficient to evaluate the importance of singlelayer network nodes. Moreover, we quantitatively calculated the decouple degree of cascading failures caused by exposed nodes formed by attack propagation. At last, experiments based on the IEEE 57 test system show that the proposed compound centrality algorithm can match the cyberattack propagation scenario well, and we give the importance values of the nodes in a specific attack scenario.
1. Introduction
In recent years, electric power systems are facing unprecedented threat from cyberattacks due to the rapid development of the information network and the higher integration of critical infrastructure and IED (intelligent electronic device) equipment in power CPS (cyberphysical system) [1–3]. The widely deployed ICT (information and communications technology) system makes the interaction more complex among multiple systems, especially more vulnerable under cyberattacks. The failure of a singlelayer system after being attacked by a cyberattack will spread through the interdependent network, causing fragmentation and cascading failures [4, 5] and induce a largescale power flow abnormal transferring. The abovementioned cascading failure process and cyberattacks will eventually lead to a blackout that led to the collapse of the power system such as large blackouts which occurred in North America in 2003 [6], Rome in 2004 [7], and Ukraine in 2015 [8].
A large number of studies have shown that the power grid has smallworld effect and scalefree property [9]. It shows strong robustness under random attacks but very fragile to deliberate attacks. Wang et al. in [10, 11] proposed methods to improve the stability of the power system. Albert et al. in [12] found that the power grid can maintain stable under most disturbances, but when the key power nodes are attacked, the synchronization ability of the grid will be greatly reduced. Therefore, identifying and assessing important nodes in the power grid and performing prevention and control are of great significance for improving the reliability and survivability of the power system. The current node importance evaluation methods include evaluation methods based on local information such as Kshell decomposition method [13], based on node path such as closeness centrality [14], based on feature vectors such as PageRank algorithm [15], and based on node removal [16] and contraction [17]. The authors in [17] applied the node contraction method to the power grid and verified the feasibility of identifying important nodes based on the topological structure of the power grid. A comprehensive evaluation index that takes into account both electrical characteristics and topological structure characteristics is proposed to identify important nodes in the power grid [18]. Based on the node link strength defined by power flow tracking, the author in [19] identifies important nodes in the power system from the perspective of global energy transmission. In [20], the authors analyze the benefits, losses, costs, and other factors of network attacks and use dynamic Bayesian networks to comprehensively evaluate the attack effects of network nodes. In [21], the improved threat propagation tree is used to evaluate the situation and a CPPS security situation assessment model that considers threat propagation is proposed. The above node importance evaluation method only focuses on the characteristics of a singlelayer network and does not consider the method of identifying key nodes of the power system under the interdependent network.
About interdependent networks, Buldyrev first analyzed the cascading failure in 2010 and proposed a cascading failure model based on network topology [22]. The author in [23] proposed electrical distance and node electrical coupling connectivity metric to identify key nodes in complex power grids. The average load balance of adjacent nodes and the parameter of network load rate are combined to measure the impact of the disabled node on the network and judge the importance of nodes in the interdependent network [24]. The above node importance evaluation considers the impact of disabled nodes on the interdependent network from different aspects but does not effectively analyze the composite value of nodes in the power CPS under the cyberattack propagation scenario.
Whether the above node importance evaluation is based on a singlelayer network or an interdependent network, the importance value is only a fixed value calculated based on a certain characteristic of the system. Without considering the potential dangers brought by the spread of cyberattacks, such fixed indexes cannot meet the system’s need to distinguish important nodes under the attack spread scenarios. In this paper, the analysis of potential dangers caused by the spread of cyberattacks is shown in Figure 1. Suppose that at a certain time during the interval between two detections, a certain device in the information network is attacked by a cyberattack (such as worms, and Trojan Horse).(1)First, before the attacked device is successfully detected by the power CPS and the countermeasures are executed, the spread of the cyberattack is in the first stage. In this stage, the attack may have spread from the initially attacked device to the remaining devices that have topological connections or information interactions with it. Take the information layer network as an example: the data detection frequency cannot match the information exchange rate, which will cause the cyberattack to be spread to the new information equipment due to the information exchange before the next detection. If the attack makes the transmitted device in exposed state (the device has been successfully attacked, but the attacker did not perform any attack operations), then the power system cannot successfully detect such exposed device in the next detection.(2)Subsequently, the power system detects the initial attacked device and executes countermeasures. Although the current detection failure rate can meet the expectations of the dispatch center, it cannot be reduced to zero. Therefore, the response strategy made by the dispatch center can only reduce and block the harm caused by cyberattacks to a certain extent. And, this strategy does not consider the cyberattack propagation described in (1). In fact, after the first propagation stage, there are most likely exposed nodes in the power system.(3)The period from the execution of the strategy to the next detection, the spread of cyberattacks has reached the second stage. Although the previous strategies dealt with the initially attacked device, the exposed nodes formed in the first propagation stage will continue deepen the scope of cyberattacks over time.(4)After the above two propagation stages, most devices have the potential to be attacked. This potential risk may expose the power system to great danger during the interval between two detections.
In order to improve the proposed node importance evaluation method of power CPS and effectively reduce the potential risks caused by the cyberattack propagation between the detection gap, in this paper, from the information network and cyberphysical interdependent network two aspects, we propose a node importance evaluation method suitable for cyberattack scenarios. The main contributions of this paper are as follows:(1)At the information network level, considering the propagation of cyberattacks among information devices, a compound centrality index and its calculation method of nodes based on adaptive coefficients are proposed. This index improves the onesidedness of the existing centrality indexes, and the timevarying calculation result better matches the importance of the power CPS node under the cyberattack propagation scenario.(2)Establish a cyberphysical interdependent network model to analyze the cascading failure behavior characteristics of the system and give a quantitative analysis of the degree of power grid decoupling under different attack scenarios. Calculate the degree of decoupling caused by cascading failures which caused by exposed nodes under the spread of cyberattacks and comprehensively reflect the importance of different nodes in power CPS.
In the second section, we proposed the node risk and established a network attack propagation model. The third section improves closeness centrality and proposes a compound centrality algorithm based on adaptive coefficients. In the fourth section, the power system cyberphysical interdependent network is established to quantitatively calculate the degree of cascading failure that the exposed node formed by the spread of cyberattacks may cause. In the fifth section, a case study is given to prove the advantages of our proposed algorithm and report the conclusion.
2. CyberAttack Propagation Model
Cyberattacks occur more commonly in power communication networks due to the deep integration of communication and measurement devices such as phasor measurement units (PMUs). Attackers use new viruses and worms with autonomous propagation capabilities to falsify or control the resources of the information layer to endanger the security and reliability of the power grid. Therefore, in this section, establish a cyberattack propagation model for power CPS information layer devices which can apply to all equipped information devices with communication capabilities, such as information acquisition system SCADA (Supervisory Control and Data Acquisition), energy management system EMS (Energy Management System), WAMS (Wide Area Measurement System), and so on.
Since SCADA, EMS, and WAMS have a certain time interval in collecting, monitoring, processing information, and transmitting instructions, while considering the characteristics of information systems based on discrete events and the propagation mechanism of new viruses such as worms, the unit propagation time of the attack is defined as the data interaction interval of power CPS. According to the propagation stage of the cyberattack, the information equipment is divided into 5 states. A single device can only be in one state within a unit of propagation time . Under certain conditions, the current state can be transformed into other states. The situation transformation process is shown in Figure 2.(1)Infected. The device has been successfully infected by a cyberattack and has been manipulated by the attacker to perform some actions (for example, tampering with operating data and modifying switch status). It has the ability to infect other devices. The detection system can detect the device.(2)Exposed. The device successfully infected by a cyberattack, but the attacker has not performed any attack operations. It has the ability to infect other devices. The detection system cannot find the device.(3)Susceptible. The device directly connected to the infected or exposed device in the communication topology. There are security vulnerabilities that can be infected by malicious code such as worms but have not been infected yet and do not have the ability to infect other devices.(4)Isolated. The device that has been attacked is automatically isolated by the detection software or manually by the system operator. No longer have hardware connections or data interactions with other devices. No longer has the ability to infect other devices.(5)Normal. Devices are in a normal state of untouched cyberattacks.
2.1. Probability of Propagation
In this section, we calculate the cyberattack propagation probability between pairs of individuals who have a connection in the communication topology. Consider a pair of devices which are connected, one of which is infective or exposed and the other susceptible. Attack spreads through them by contact from to . Then, the rate that not being attacked in the continuous time system is as follows [25]:where is the probability that the cyberattack successfully spreads from infective or exposed device to susceptible device under the selfprotection mechanism of the device and information system. The infective device remains infective for a time that has not been detected and disabled by the detection system.
Due to the cyber systems based on discrete events, use discrete timesteps rather than continuous time, in which case instead of taking the limit in equation (1) we simply set , givingwhere is the average detection time of the online detection system that is measured in data transmission timesteps.
2.2. Propagation Model
For modeling the attack propagation clearly, we abstract different types of communication devices in the information layer as nodes. According to the communication topology of the information network, the connection matrix is defined. The matrix elements are binary variables which equals 1 if is kept connected to in communication and 0 otherwise. We define the probability that a node may be attacked and turn into exposed state as the node risk . Let us assume at the initial time between two detection gaps, one or several nodes in the CPS information layer are compromised by cyberattacks and turn into infected while the remaining nodes are uncompromised. In the following period of time, the remaining nodes will turn into exposed with a certain probability, accompanied by a timevarying node risk.
We define as the set of information nodes; is the set of infected nodes; and is the time when the attack spreads through n units of propagation time . The attack spreads in stages, and the risk of each node changes with time as follows [26]:(1)At the initial time , the risk of infected and other nodes are as follows: where is a binary variable which equals 1 () if the node has the possibility of being attacked at time , and ; otherwise, is the node risk of at time .(2)By time , cyberattacks may spread from infected nodes to susceptible nodes via routers and communication connections, turning susceptible nodes into exposed nodes, and the detection software has not detected them yet. The risk of each node is given by the following equation: where is the probability that the attack propagates from node to node during the time , and it is given in equation (2).(3)By time , the range and volume of cyberattacks further expand over time, and the risk of exposed nodes deepens. The risk of each node is an iterative process over time. According to equations (4) and (5), the node risk by time is given by the following equation: where is the probability that the node is not attacked at time but is attacked at time . In equation (7), only is unknown during the iteration process. This variable is discussed as follows: Substituting equation (8) into equation (7) to obtain the risk of some nodes at time :(4)By time (n > 2), the node risk is given as follows:(5)During the above propagation process, if the initially infected node is detected and processed by the power system at time , then in the time after that, the risk of the node is
Notice. One of the ways the power system responds to cyberattacks is as follows. After successfully detecting the cyberattack, the power system will disable infected nodes automatically by the detection software or manually by the system operator. These nodes are no longer connected to the information network (no more data interaction with other nodes), while the system no longer trusts its uploaded data. Namely, this type of node turns into the isolated state. Then, they may be recovered to normal by the system and reconnect to the network, but in the propagation discussed in this paper, such nodes no longer participate in the propagation of cyberattacks, that is, they cannot be reinfected by cyberattacks.
3. Node Centrality Algorithm Based on Adaptive Coefficient
The connection form of the node and its position in the topological structure have a vital influence on the promotion or interruption of the spread of cyberattacks [4]. At present, a variety of centrality algorithms have been applied to node importance evaluation. In this section, we propose a node importance evaluation algorithm that considers cyberattacks propagation and the potential threats it may cause in the detection gap, which is more flexible and more suitable for attack scenarios.
3.1. Centrality Analysis
According to the connection form of the node and its position in the topology, common node centrality algorithms include degree centrality [27], closeness centrality [28], betweenness centrality [28], and so on, all describe the importance of nodes in the network from different respects, and they are given as follows:where N is the number of nodes in the network, is a binary decision variable that judges the topological structure or information interaction between nodes and , is the shortest path between nodes and ; in this paper, we use the minimum number of nodes in the path from to as the value of the shortest path. is the number of shortest paths from to ; is the number of shortest paths from to which passing by node .
Degree centrality indicates the sum of the number of nodes directly connected to the designated node. Closeness centrality is the reciprocal of the average shortest path from the designated node to all other reachable nodes. In general, the closer a node is to other nodes, the greater its closeness centrality. Conversely, the smaller the node closeness centrality, the more the node is at the edge of the network. Betweenness centrality is the number of times that the designated node is located on the shortest path between any two other nodes.
The three algorithms describe the importance of nodes from the perspectives of local characteristics, global characteristics, and propagation characteristics, while they also have certain limitations. Degree centrality can only onesidedly reflect the closeness between the designated node and its surrounding nodes. For example, the connection node between two partitioned networks has a smaller degree of centrality but a higher degree of importance. Or for some nodes with high degree centrality, the clustering network it embeds may be at the edge of the system. At this time, the degree centrality is relatively high but the degree of importance is average. However, as the nodes in a complex network, especially in largenode systems such as power systems, the position of the node in the system and some of its functions are far more important than the number of nodes around it. The descending order curve of the normalized values of DC, CC, and BC of each node in the IEEE14 system is shown in Figure 3. The comparative analysis is as follows:(1)The relative DC values of all nodes are arranged in descending order as shown in. It descends in steps, and there are 3 steps and some steps are larger in width. It descends in steps, and there are 3 steps and some steps are larger in width. For example, the DC of 6 nodes in the graph are all 0.4. The importance of these nodes in the same step cannot be distinguished.(2)The descending arrangement curve of relative CC is relatively smooth, but there is little difference among the values. The variable interval of the value is 0.585–1, which is too narrow.(3)The relative BC descending order curve is also a smoother curve with values distributed in the interval 01. However, there is a platform segment with a value of 0, which accounts for 28.57% of the total system nodes, at its end. The node in this segment cannot effectively distinguish the importance.
3.2. Centrality Algorithm Based on Adaptive Coefficient
According to the analysis results of the three centrality algorithms in 3.1, we can get the following. The three algorithms have their respective advantages and disadvantages. If only one centrality algorithm is selected as the basis for judging the importance of nodes, the sorting result is not appropriate and accurate. These centrality algorithms make it impossible to effectively distinguish the importance of each node in the system under a specific attack scenario based on the fixed value calculated by the topology structure. Moreover, if the importance of nodes is only distinguished based on the node risk P that proposed in 2.2, the position of the node in the system is ignored, and the mutual influence between nodes in the process of attack propagation is discarded. Therefore, after comprehensively considering the global attributes of CC and the propagation attributes of BC, this section improves the CC algorithm according to the propagation characteristics of cyberattacks and proposes a compound centrality algorithm that considers both improved CC and BC to weaken the differences and shortcomings of existing centrality algorithms.
Improved closeness centrality algorithm:where is the shortest path from to , is the node initially infected by the cyberattack, and is the remaining nodes in the system.
The compound centrality algorithm proposed in this paper performs a weighted summation of the two indexes of improved CC and BC, but the two are not completely independent indexes. For example, nodes with a high improved CC will have a higher BC under a certain probability. Summing the two indexes with a weight of 0.5, reasons such as duplication of information and redundant factors will cause the calculated composite index to be inaccurate. To solve this problem, this section proposes an adaptive coefficient as the weight of the two centralities.
In the propagation of the power CPS cyberattack, the propagation probability between nodes , as a key factor, has always been defined as a fixed value in previous studies. Considering the differences in the importance of each node in the network, a correction equation is proposed to modify to a certain extent:where is the correction coefficient and is the centrality index of node . The of node changes with the centrality index . The higher the , the lower the . That is, when increases, the probability of a cyberattack spreading from node to node decreases.
The iterative solution process of the adaptive coefficient is given by equations (17)–(23). Equation (22) is the termination condition of the iteration, and equation (23) is the objective function. The appropriate initial values of are selected and substituted into formula (17) to start the iteration. In equations (18) and (19), the composite centrality index of the node and the propagation probability are updated with the update of . According to the risk model proposed in 2.2, use the iteratively updated to calculate the risk of each node and the average of the overall node risk, which is shown in equation (21). Substituting the average of the overall node risk in the current and previous iterations into equation (17) again will start a new round of iterative correction process. Until satisfies the iteration termination condition , the iteration process stops and jumps out of the iteration. Take the minimum the average of the overall node risk in all rounds of iterations as the objective function and find the optimal solution that satisfies the objective function in the iteration process. That is, at the target time ( is an integer multiple of ), after correcting with the optimal , the average of the overall node risk in the system is the smallest. This optimal solution is the adaptive coefficient of this paper. It can be seen from the above algorithm flow that if we want to start the iterative process of , we need to give the first two initial values of . The rationality of the initial value selection and its influence on the number of iterations and the judgment of the optimal solution are discussed in the case study in Section 5.where is the adaptive coefficient of the kth iteration, is the compound centrality of node when , is the probability of successful propagation of the attack from to modified according to , is when , the risk of node at time T, is when , the average risk of all nodes in the system at time T, and is the decision function of . A certain in the iterative process minimizes the average risk of nodes in the system at the target time T.
4. Cascading Failure of Power Interdependent Network
4.1. Interdependent Network Model
Electric power CPS is a multidimensional heterogeneous system that deeply embeds perception, information processing, and control platforms into the power system to meet realtime monitoring and achieve commanddriven of the power system. The power network provides power for the information network, and the control and analysis module in information network reversely drives the power network. The interdependence between the two networks enables the power system to be modeled as a power cyberphysical interdependent network. Based on the “undirected” and “disordered” characteristics of general types of cyberattacks (nondirected attacks) spreading in information systems, both the power network and the information network are equivalent to nonweighted undirected networks. Combining the actual situation of China’s power system [5], only considering the characteristics of interconnection between nodes, using complex network theory to simplify the power and information network is as follows:(1)Ignore the functional differences among the plants and stations, and regard the power generation nodes and substation nodes, dispatching nodes and routing nodes in the information network as equivalent nodes, regardless the difference in the types, quantities, and deployment modes of devices in various sites(2)Ignoring the differences in information protocols and hierarchical structures among information nodes at all levels, it is considered that the lines between nodes can transmit bidirectional information, and multiple information lines in the same direction are combined to eliminate multiple edges and selfloops.
Use and to represent a power network with n nodes and k branches and an information network with m nodes and branches, where is the set of edges and is the set of nodes in the network, respectively. Each node in the power network is connected to an information node in . The power node provides electrical energy support to the information node, and the information node receives the status information sent by the power node and feeds back control instructions. Naturally, each node in connects and depends on the corresponding node in , and vice versa. In addition, according to the important status of the scheduling nodes, it is set as an autonomous node independent of the power grid. It deploys a complete backup power supply and power generation equipment, which is not affected by power grid energy fluctuations.
China’s power line information network is a dedicated resource for the power system, and most of the communication lines are laid along with highvoltage transmission lines. The geographical similarity of the two layouts makes the topological structure between the information network and the power network highly similar. On the other hand, in order to meet the needs of control and dispatch, the information network also has dispatching nodes. So, the information network has more stations than the power network. At the same time, the optical information network needs to be formed into a ring to protect its selfhealing ability, and the structure of the information network is more complicated [29]. Therefore, the dependent network model selects the “parttoone correspondence” coupling mode, and the established power cyberphysical dependent network model is shown in Figure 4.
Different from the singlelayer network with only connectivity link, there are two types of edges in the interdependent network: the connectivity link and the dependency link. The nodes in the singlelayer network rely on the internal connectivity link (the black solid lines in Figure 4) to achieve corresponding functions. For example, the power generation nodes, substation nodes, and load nodes in the power network realize the generation, transmission, and consumption of electric energy through the transmission line. The topological structure of the information network is more complicated than that of the power network. In Figure 5, there are connecting links that do not exist in the power network between the information node 2 and the information node 8. The dependency link between the two networks (the red dotted line in Figure 4) is used as a medium for energy or information exchange to realize the mutual influence between the twolayer networks. The dependency matrix is defined according to the connection relationship of the power cyberphysical dependent network, and the matrix element uses logic elements of “1” and “0” to indicate whether there is a dependency edge between the power and information nodes.
(a)
(b)
(c)
(d)
In summary, establish a power cyberphysical interdependence network that includes power grids, information networks, and the interdependency links , where , is the set of nodes in a singlelayer network, is the set of connectivity links in a singlelayer network, and is the set of the dependency links between the power cyberphysical dependent network.
4.2. Cascading Failure Behavior
According to the operating principle of the dependent network, the failure of the power node or information node will lead to the failure of the dependent node in the other network. The power network or information network will be broken into several fragmented networks, and the scope of the attack will expand with the expansion of the network fragmentation and eventually cause cascading failures in the dependent network [29]. Describe the cascading failure behavior caused by the attack as follows:(1)When one or some information nodes (power nodes) in the information network (or power network ) are attacked and fail, the connectivity links and the dependency links on these nodes fail too;(2)Corresponding nodes in the power network fail due to the interdependence with the failure node in . Corresponding, the connectivity links and the dependency links on the nodes also fail. After removing all the abovementioned faulty nodes, the connectivity links, and the dependency links, the power and information network is decomposed into several fragmented networks.(3)(On the basis of 2), identify the connected subgraph in and , and judge the nodes that do not belong to the connected subgraph as failed nodes. Remove the newly determined failed node and its connectivity links and dependency links. Based on this process, the system will reach stability after a certain number of propagations of the failure in the dependent network.(4)Identify the set of nodes in the maximal connected subgraph in the stable network, and finally determine the stable state of the decoupled network.
Notice. According to the main research content of this paper, the information layer network attack and its propagation mechanism demand for network connectivity use “the set of nodes in the maximal connected subgraph” [29] to judge the stable state of the dependent network. The set element is the powercyber node group , and the dependency links between , where belongs to the maximal connected subgraph in the power network and also belongs to the maximal connected subgraph in the information network. In order to ensure the effectiveness of connectivity, after the fragmentation of the network, only the nodes in the set of nodes in the maximal connected subgraph are available, and the remaining nodes fail.
According to the cascading failure behavior, taking the IEEE14 node as an example, the failure decoupling process of the dependent network is shown in Figure 6 when the initial attack node is 4 and 11. The upper network in the figure is the power system , and the lower network is the information system . Figure 6(a) shows that the connectivity of is similar to that of , but compared with , has one more independent green “scheduling node” and green connectivity links between the scheduling nodes and each information node, and three more green connectivity links between 4–14, 3–11, and 5–12. In Figure 6(b), the nodes 4 and 11 failed due to cyberattacks and turned to red. Correspondingly, the connectivity link connected to it fails and turns to a dashed line, and the dependency link fails and turns to a red line, and the corresponding dependent nodes and connectivity links in are treated in the same way. In Figure 6(c), we delete the failed node and the failed edge in Figure 6(b). The remaining nodes that do not belong to the connected subgraph are turned into blue, and the connecting edges on these nodes turn to blue dashed lines. At this time, it can be seen that both networks are fragmented since 5–12 connecting edge in , and the number of nodes in the connected subgraph is greater than that in , and the degree of fragmentation is not as serious as . Figure 7(d) filters out the set of nodes in the maximal connected subgraph from Figure 7(c). At this time, the power interdependent network finally reaches a stable state.
4.3. Decoupling Degree Based on Cascading Failure
According to the behavior of cascading failures, it can be seen that the cascading failure propagation process in this paper is suitable for random cyberattacks and deliberate cyberattacks of power CPS. The attacker can choose one or some nodes as the attack target to form different combinations of cyberattacks. Different combinations of attacks will cause different decoupling processes. The exposed nodes formed in the process of cyberattack propagation will not only cause the attack propagation between the information layer devices during the detection gap but also cause the cascading failure of the power cyberphysical dependent network. The node survival rate is established to describe the degree of decoupling of power CPS when the system reaches a steady state after the cascading failure. The larger the value of , the more the number of nodes remaining in the stable state of the system and the lower the degree of system decoupling.where and are the number of effective nodes in the power network before and after the cascade failure and and are the number of effective nodes in the information network before and after the cascade failure.
5. Experimental Results
We use the connection relationship of the IEEE57bus test system to simulate the information layer network connection relationship, ignoring the weight of each link; it is shown in Figure 5.
5.1. Compound Centrality Algorithm
Notice. In the information layer cyberattack propagation model described in Section 1, the initial time should be the time when the cyberattack occurred. However, at present, the power system cannot accurately analyze which time the cyberattack occurred between the two detection gaps. Therefore, in the experiment, the moment when the cyberattack is detected is regarded as the initial moment in the risk calculation. And, set the power system to make a response strategy to disable the attacked device at time .
Establish the attack scenario with attack the information device equipped on node 18. The BC and improved CC indexes of the remaining nodes calculated according to formulas (13) and (14) are shown in Table 1, and the two indexes are ranked according to importance. It can be seen from the table:(1)The improved CC index mainly serves the attack scenario. The closer the node to the initial infected node, the lower the node’s improved CC value. At the same time, the indexes have a certain degree of repeatability; for example, the value of the improved CC index for 7 nodes is equal to 4.(2)The importance ranking results of BC and improved CC are obviously different. Some nodes have higher improved CC and lower BC, such as node 9. Some nodes have lower improved CC and higher BC, such as node 37.
 
Notice. According to the improved CC, the lower the CC, the closer the node is to the attacked node and the higher the node importance. 
Perform log normalization processing () on the BC and improved CC indexes in Table 1, and use them as the input data to execute the compound centrality algorithm proposed in this paper. The iterative results and iteration times N of the adaptive coefficients under different correction coefficients and different target times are shown in Table 2.

According to the iterative results in Table 2, it can be seen that in the initial stage of attack propagation, there is a certain proportion of improved CC in the compound centrality. That is, in the initial stage of propagation, mainly modifying of node which is near the infected node has a better effect on reducing the overall risk. However, as the attack spreads, the proportion of BC in the compound centrality gradually increases. That is, the propagation characteristics possessed by nodes in the topological structure in the later stage of propagation have a more obvious effect on reducing risk which caused by the spread of the attack. In addition, it can be seen that the iteration times of this algorithm are less than 10, and the convergence speed of the algorithm is fast, which can meet the demand for timeliness of power system scheduling. And, we verified that the selection of the initial value of has no effect on the convergence value of the iteration and the optimal solution. It has a slight impact on the iteration times, but the number of iterations can also be guaranteed within 10 times.
In the case of , and , , using compound centrality, improved CC, BC, and DC to modify , the risk of the remaining nodes of the system is shown in Figures 7 and 8.
Figure 7 shows the following:(1)In the initial stage of attack propagation, the risk of several nodes that are topologically close to the initial infection node 18 rises fastest (for example, node 4, 19, 5, 6, 20)(2)In the initial stage of attack propagation, in view of the strong correlation between the improved CC and the cyberattack scenario, the minimum risk after the use of the improved CC indicator to modify is slightly better than the effect of using the composite centrality indicator to modify.
Figure 8 shows the following:(1)Node 4 is located in the direct connection position of the infected node 18, and node 6 is in a special position in the connection relationship (with a high BC and DC index). When the conventional centrality is used to modify , because the index does not consider the characteristics of each attack scenario, the risk of node 4 and node 6 rises faster, which cannot effectively reduce the risk in the system. Note that the correction effect of the BC is better than that of the DC; in other words, during the attack propagation process, the propagation attribute of the node is more important than the local attribute.(2)When using the improved CC and compound centrality to modify, since the characteristics of the attack scenario is considered in the index, the attack spread in the system is effectively curbed and the risk of node 4 and node 6 is reduced. Note that the improved CC only considers the actual situation of the attack, which is really effective in reducing the risk of node 4 and node 19 which are close to infected node 18, but does not consider the role of nodes in the subsequent attack propagation. Therefore, the effect of risk reduction for other nodes is worse than that of the compound centrality index.
Figures 7 and 8 show the following:(1)As the attack spreads, the compound centrality is used to modify , which effectively reduces the impact of the attack on the 19th node with a higher initial risk.(2)Although in the initial stage of attack propagation, the effect of using the improved CC to modify is slightly better than that of the compound centrality modification; considering the overall attack propagation, the effect of compound centrality correction is better.
In summary, the compound centrality considers both the characteristics of attack scenarios and the role of nodes in attack propagation, and using it to modify can effectively reduce the risk of the remaining nodes in the system. The ranking result of compound centrality is more suitable for identification requirements of key nodes in the power system under cyberattack scenarios. In the case of and , the importance ranking of the remaining nodes in the system is shown in Figure 9 when the compound centrality, improved CC, BC, and DC indexes are used to modify .
5.2. Cascading Failure Analysis
The node importance evaluation in this paper focuses on the distinction of the importance of a single node. Therefore, according to the cascading failure model established in Section 3, each time a single node is selected for attack and the maximal connected subgraph and the node survival rate under the stable state are obtained. Analyze the difference of decoupling results when attacking different single nodes. The experiment uses the IEEE57 system “onetoone correspondence” coupling method. Sequentially attacking a single node of the information network, the power CPS decoupling steady state results are shown in Table 3.

It can be seen that the decoupling caused by the attack on the 24^{th} node is the greatest. After the decoupling, there are only 45 nodes left in the stable system, accounting for 78.9% of the total number of nodes. In the case of cascading failures caused by a cyberattack on the power CPS, the importance of nodes based on the node survival rate is ranked as follows: {24,29,,, , , , 33}. There are some nodes with the same in the sorting result, so the importance of these nodes cannot be accurately distinguished. The centrality of each node calculated by the compound centrality algorithm of the onesided network can just make up for this shortcoming. Considering two indexes at the same time, the importance of nodes can be distinguished smoothly.
In summary, the experiment gives indexes that can describe the importance of IEEE57 nodes from two aspects: the compound centrality of the singlelayer network and the value that reflects the degree of cascading failure decoupling, as shown in Table 3. Both indexes can be used as a reference basis for operators and operation centers in the operation of the power system. The dispatch center can select a certain index according to the actual situation of the cyberattack to distinguish the importance of the node and guide the followup strategy.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Key R&D Program of China (2018YFA0702200) and the Fundamental Research Funds for the Central Universities (N2004013).
References
 Q. Guo, S. Xin, H. Sun et al., “Power system cyberphysical modelling and security assessment:motivation and ideas,” Proceedings of the CSEE, vol. 36, no. 6, pp. 1481–1489, 2016. View at: Google Scholar
 Y. Liu and P. Ning, False Data Injection Attacks against State Estimation in Electric Power Grids, Association for Computing Machinery, New York, NY, USA, 2009.
 S. Cui, Z. Han, S. Kar et al., “Coordinated datainjection attack and detection in the smart grid: a detailed look at enriching detection solutions,” IEEE Signal Processing Magazine, vol. 29, no. 5, pp. 106–115, 2012. View at: Google Scholar
 Y. Wang, K. Gao, T. Zhao et al., “Assessing the harmfulness of cascading failures across space in electric cyberphysical system based on improved attack graph,” Proceedings of the CSEE, vol. 36, no. 6, pp. 1490–1499, 2016. View at: Google Scholar
 X. Yi, B. Wang, D. Chen et al., “Review on interdependent networks theory and its applications in the structural vulnerability analysis of electrical cyberphysical system,” Proceedings of the CSEE, vol. 36, no. 17, pp. 4521–4532, 2016. View at: Google Scholar
 A. Vespignani, “The fragility of interdependency,” Nature, vol. 464, no. 7291, pp. 984985, 2010. View at: Publisher Site  Google Scholar
 X. Yi, Research on Structural Vulnerability of Power InformationPhysical System Based on Interdependence Network Theory, Wuhan University, Wuhan, China, 2016.
 Q. Guo, S. Xin, J. Wang et al., “Comprehensive security assessment for a cyber physical energy system a lesson from Ukraine’s blackout,” Automation of Electric Power Systems, vol. 40, no. 5, pp. 1–3, 2016. View at: Google Scholar
 Z. Meng, Z. Lu, and J. Song, “Comparison analysis OF the smallworld topological model of chinese and american power grids,” Automation of Electric Power Systems, vol. 28, no. 15, pp. 21–24, 2004. View at: Google Scholar
 R. Wang, Q. Sun, P. Zhang et al., “Reducedorder transfer function model of the droopcontrolled inverter via Jordan continuedfraction expansion,” IEEE Transactions on Energy Conversion, vol. 35, no. 3, pp. 1585–1595, 2020. View at: Google Scholar
 W. Hu, C. Ruan, H. Nian, and D. Sun, “Zerosequence current suppression strategy with common mode voltage control for openend winding PMSM drives with common DC bus,” IEEE Transactions on Industrial Electronics, vol. 99, 2020. View at: Publisher Site  Google Scholar
 R. Albert, I. Albert, and G. L. Nakarado, “Structural vulnerability of the North American power grid,” Physical Review E, vol. 69, no. 2, Article ID 25103, 2004. View at: Publisher Site  Google Scholar
 M. Kitsak, L. K. Gallos, S. Havlin et al., “Identification of influential spreaders in complex networks,” Nature Physics, vol. 6, no. 11, pp. 888–893, 2010. View at: Publisher Site  Google Scholar
 S. Dolev, Y. Elovici, and R. Puzis, “Routing betweenness centrality,” Journal of the ACM, vol. 57, no. 4, pp. 1–27, 2010. View at: Publisher Site  Google Scholar
 A. N. Langville and C. D. Meyer, “Deeper inside page rank,” Internet Mathematics, vol. 1, no. 3, 2003. View at: Google Scholar
 Y. Chen, A. Hu, and X. Hu, “Evaluation method for node importance in communication networks,” Journal of China Institute of Communications, vol. 25, no. 8, pp. 129–134, 2004. View at: Google Scholar
 Y. Tan, J. Wu, and H. Deng, “Evaluation method for node importance based on node contraction in complex networks,” Systems Engineering Theory& Practice, vol. 26, no. 11, pp. 79–83, 2006. View at: Google Scholar
 Q. Xie, C. Deng, H. Zhao et al., “Evaluation method for node importance of power grid based on the weighted network model,” Automation of Electric Power Systems, vol. 33, no. 4, pp. 21–24, 2009. View at: Google Scholar
 J. Wang, X. Gu, T. Wang et al., “Power system critical node identification based on power tracing and link analysis method,” Power System Protection and Control, vol. 45, no. 6, pp. 22–29, 2017. View at: Google Scholar
 Z. P. Network, “Security situation analysis based on a dynamic bayesian network and phase space reconstruction,” Journal of Supercomputing, vol. 76, no. 2, pp. 1342–1357, 2020. View at: Google Scholar
 G. Li, P. Huang, Y. Chen et al., “Security situation assessment method for cyber physical power system considering threat propagation characteristics,” Electric Power Construction, vol. 40, no. 5, pp. 29–37, 2019. View at: Google Scholar
 V. Buldyrev Sergey, R. Parshani, G. Paul, H. Stanley, and S. Havlin, “Catastrophic cascade of failures in interdependent networks,” Nature, vol. 464, 2010. View at: Google Scholar
 Y. Tan, X. Li, Y. Cai et al., “Critical node identification for complex power grid based on electrical distance,” Proceedings of the CSEE, vol. 34, no. 1, pp. 146–152, 2014. View at: Google Scholar
 R. Wu, B. Zhang, and L. Tang, “A cascading failure based nodal importance evaluation method applied in dual network coupling model,” Power System Technology, vol. 4, pp. 1053–1058, 2015. View at: Google Scholar
 M. E. J. Newman, “The spread of epidemic disease on networks,” Physical Review E Statal Nonlinear & Soft Matter Physics, vol. 66, no. 1, Article ID 16128, 2002. View at: Publisher Site  Google Scholar
 S. Mousavian, J. Valenzuela, and J. Wang, “A probabilistic risk mitigation model for cyberattacks to PMU networks,” IEEE Transactions on Power Systems, vol. 30, no. 1, pp. 156–165, 2015. View at: Publisher Site  Google Scholar
 D. Krackhardt, “Assessing the political landscape: structure, cognition, and power in organizations,” Administrative Science Quarterly, vol. 35, no. 2, pp. 342–369, 1990. View at: Publisher Site  Google Scholar
 L. C. Freeman, “Centrality in social networks: conceptual clarification,” Social Networks, vol. 1, no. 3, 1979. View at: Google Scholar
 X. Ji, B. Wang, D. Liu et al., “Improving interdependent networks robustness by adding connectivity links,” Physical A: Statistical Mechanics and its Applications, vol. 444, pp. 9–19, 2016. View at: Google Scholar
Copyright
Copyright © 2021 XinRui Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.