Reliability Evaluation of Smart Microgrids Considering Cyber Failures and Disturbances under Various Cyber Network Topologies and Distributed Generation’s Scenarios

: Smart microgrids (SMGs), as cyber–physical systems, are essential parts of smart grids. The SMGs’ cyber networks facilitate efﬁcient system operation. However, cyber failures and interferences might adversely affect the SMGs. The available studies about SMGs have paid less attention to SMGs’ cyber–physical features compared to other subjects. Although a few current research works have studied the cyber impacts on SMGs’ reliability, there is a research gap about reliability evaluation simultaneously concerning all cyber failures and interferences under various cyber network topologies and renewable distributions scenarios. This article aims to ﬁll such a gap by developing a new Monte Carlo simulation-based reliability assessment method considering cyber elements’ failures, data/information transmission errors, and routing errors under various cyber network topologies. Considering the microgrid control center (MGCC) faults in comparion to other failures and interferences is one of the major contributions of this study. The reliability evaluation of SMGs under various cyber network topologies, particularly based on an MGCC’s redundancy, highlights this research’s advantages. Moreover, studying the interactions of uncertainties for cyber systems and distributed generations (DGs) under various DG scenarios is another contribution. The proposed method is applied to a test system using actual historical data. The comparative test results illustrate the advantages of the proposed method.


Introduction
The microgrid (MG) concept has received a great deal of attention because of its contribution to emission mitigation and sustainable development [1][2][3]. The cyber system, including communication networks, information technologies, and smart monitoring/control systems, boosts the MGs' reliability, efficiency, and technoeconomic aspects [4,5]. However, the cyber system impacts, such as other uncertainties, might influence smart grids and MGs' reliability [6][7][8]. Although the impacts of the uncertainties in the output power of renewable distributed generations (DGs) on MGs' reliability have been studied in the literature, such as in [9][10][11], less attention has been paid to the eventual cyber-power uncertainties and interdependencies of cyber-power MGs (CPMGs).
The Monte Carlo simulation (MCS) methods could help evaluate smart grids/MGs' reliability considering cyber-power interdependencies and other stochastic behaviors [12]. The sequential MCS has been adopted to quantify the MGs' reliability due to frequency control issues, while cyber failures have been of concern in [13]. In [14,15], the MCSbased smart grid's reliability evaluation methods considering the impacts of direct and indirect cyber-power interdependencies have been reported, respectively. Wang et al. [16] evaluated the effects of data transmission errors, latency, cyber element defects, and routing errors on MG reliability using the MCS. Although the impacts of most essential parameters have been addressed in [16], the defects of the MG control center (MGCC) have not been considered. In addition, the CPMGs' reliability and interaction of cyber and renewable DG uncertainties under various DG scenarios has not been studied in [16].
The computation time of MCS-based methods affects their effectiveness, particularly for reliability-oriented short-term operational decisions. Therefore, some research works have been devoted to developing non-MCS methods to speed up the smart grid/CPMG's reliability assessment. Falahati et al. [17] reported a new probability table (P- Table)-based reliability assessment approach that considers direct cyber-power interdependencies, while other system uncertainties, such as stochastic behaviors of renewable DGs' output power, have not been considered. The indirect cyber-power interdependencies due to malfunctions of monitoring and protection systems have been examined in [18]. Zhu et al. [19] presented an analytical method that uses complex network theory to investigate the effect of communication network interruptions and failures on smart grid reliability. The reliability of each communication device and packet transmission delay has been researched in [19]. In most non-MCS reliability evaluation methods for CPMGs, the renewable DGs uncertainties have not been researched due to their limitations compared to MCS-based ones. In a few available studies, such as [20], both the DGs and cyber-power uncertainties have been considered using non-MCS methods. Barani et al. [20] have studied the uncertainty of renewable DGs along with communication technologies and cyber networks in MGs' reliability evaluation. Although cyber component failure impacts have been concerned in [20], various parts of networks and communication technologies, such as data transmission delay and data transfer errors, have not been considered.
In Table 1, a summary of the literature review in the area of smart grid and MGs' reliability evaluation considering the cyber-physical challenges is presented. In light of the literature review, there is a research gap about developing a reliability evaluation method for CPMGs, which simultaneously concerns all cyber challenges, particularly the failures of MGCC. This study aims to respond to such a research gap by proposing a new MCS-based reliability evaluation method. Proposing a new method considering both uncertainties of the cyber system and renewable DGs is the major contribution of this research. Developing a new MCS-based reliability assessment method for CPMGs using the theory graph concepts based on cyber-physical interdependencies is another contribution of this study. Moreover, the CPMGs and the interactions of different cyber uncertainties with DGs' uncertainties under various DG scenarios would be studied.
Comparing available MCS-based and non-MCS reliability assessment methods infers that the MCS-based ones could be more effective for precisely studying all uncertainties. Hence, the MCS has been selected in this research. However, the non-MCS methods that considers all cyber and physical uncertainties could be developed in future works.
The major contributions of this research are listed as follows: • Proposing a new MCS-based reliability evaluation of CPMGs that considers all cyber faults and interruptions, i.e., information transmission errors, routing errors, and delay of information channels; • Considering the MGCC faults' impacts on CPMG's reliability; • Simulating the uncertainties of the physical system of CPMG, e.g., stochastic behaviors of wind turbine (WT) and photovoltaic (PV)outputs; • Studying the CPMG's reliability under different cyber network topologies, including cyber network topologies with the redundancies of the MGCC; • An analysis of the CPMG's reliability under various DG scenarios; • Sensitivity analyses to gain insight into how the cyber failures and interruptions might affect the CPMG's reliability.
The proposed method is compared to available ones such as in [16], which has not addressed the MGCC faults. The test results illustrate the advantages of the proposed method. In addition, the proposed method is studied under different DG scenarios, including PV and WT DGs. The sensitivity analyses are performed to obtain insight into how the changes in physical and cyber characteristics might affect the CPMGs' reliability.
The organization of the rest of this paper is as follows: Section 2 presents the CPMG's control/energy management structure and different cyber-power interdependencies on system reliability, Section 3 introduces the proposed reliability evaluation method and reliability modeling of CPMGs' subsystems, Section 4 presents the simulation results and discussions are given, and, finally, the conclusion is reported in Section 5.

CPMG's Architecture
In Figure 1, the architecture of the understudy CPMGs, including cyber and physical subsystems, is shown [16,20]. The subsystems of CPMGs could be categorized into two separate general networks: physical and cyber networks, as shown in Figure 1b. Moreover, the cyber network could be divided into different layers, e.g., interface, communication, and decision layers. The different layers of the CPMGs have been depicted in Figure 1c.
In this paper, the central control mode of the CPMG has been considered. MGCC makes dissension for the MG's stability in the centralized control mode by collecting data from control devices. Accordingly, the MGCC is the most important subsystem of the CPMG, and its eventual failures might affect the CPMG dramatically. The centralized control system of the CPMG includes the MGCC and control units of sources and loads (microunit control (MC) and load control (LC)). (c) different layers and interfaces [16,20]. (c) different layers and interfaces [16,20]. The MGCC must control the power balance between the supply side and demand side at each moment. In addition, the MGCC is in charge of frequency and voltage stability. The MCs and LCs are responsible for collecting DGs and loads information such as output power, voltage, and frequency, and sending these data to the MGCC. The complete MG control process could be described as follows [16]: 1.
The gathering of DGs and loads information by the MCs and LCs, respectively; 2.
Transmitting the information through the information channel connecting the MGCC, MCs, and LCs; 3.
Checking of the power balance condition and other operation aspects by the MGCC; 4.
Adopting orders based on information received, such as increasing power generation or decreasing power in the supply side, and load management strategies on the demand side; 5.
Transmitting the commands generated by the MGCC using the information channel between the MGCC, MCs, and LCs; and 6.
Executing the MGCCs' generated commands by local controllers.
In CPMGs, in addition to physical failures and interruptions, the cyber system might affect the system's reliability [25]. In this paper, the following cyber-power interdependencies are considered: • Failure of the cyber system's elements: The appropriate and stable operation of CPMGs requires available power and information channels. The failures in cyber layers, e.g., in the MGCC, MCs, and LCs result in energy not being supplied, similar to failures in power elements [26,27]. • Information transmission faults: It is essential to receive the correct information from the sensors and control units for the MGCC to make correct decisions [28]. Hence, the information transmission faults might lead to incorrect decisions. Consequently, the appropriate operation of the CPMG and its stability might be threatened.

•
Routing faults: The CPMG equipment's real-time operation needs to receive the correct information package. Therefore, routing faults in the sending or receiving of the information package might affect the instantaneous function of the control and cyber systems [29]. Hence, the CPMG is influenced by routing faults.

•
Delay in data transfer: The fast performance of cyber elements and information channels is crucial to meet power balance conditions in CPMGs. Therefore, a delay in data transfer could cause a lack of timely communication between the control units, resulting in a loss of power balance [30].
In Table 2, the uncertainties in both cyber and power networks of CPMGs have been listed. Comprehensive consideration of all cyber and power uncertainties is one of the major contributions of this research.
As seen in Table 2, physical/power and cyber subsystem failures are major sources of a CPMG's uncertainties. The time to failures (TTF) and time to repairs (TTR) of a CPMG's subsystems are stochastic parameters, and their probabilistic behaviors should be concerned in reliability evaluations [31].
Environmental parameters, e.g., wind speed and solar irradiance, influence renewable DGs' output power. Hence, the uncertainties of environmental parameters should be considered in the physical/power network.
Moreover, the uncertainties and faults of information transmission delay and information transmission faults (payload errors and hear errors) are concerned in this article. The probability distributions of discussed stochastic parameters have also been presented in Table 2.

Cyber-Physical Interdependencies Due to Cyber Failures
The state of the CPMG should be determined to evaluate the system's reliability, such as in [17,37]. The CPMG's state, including both cyber and power elements, could be modeled using the following equation: MS i = ME 1 (t), . . . , ME j (t), . . . , ME N P +N C (t) = ME 1 (t), . . . , ME k (t), . . . , ME N P (t), . . . , ME N P +l (t), . . . , ME N P +N C (t) i = 1 : N MS , j = 1 : where the CPMG's state at each time step (MS i (t)) is distinguished based on the state of its physical elements ME k (t) and cyber elements MS l (t). The two-state reliability model is one of the most common models for modeling the CPMGs, as depicted in Figure 2 [38]. The binary variables represent the CPMG elements' states. The zero value of each element state represents that the corresponding element is unavailable. Otherwise, the discussed element is in service [39]. The availability state of any subsystem could be represented using the time to failure (TTF) and time to repair (TTR), as shown in Figure 3 [40,41]. In this paper, the two-state reliability model has been used for both the cyber and physical networks' elements. Hence, the failures in cyber and physical networks could be determined using the two-state reliability model. In Figure 2, λ and µ denote the failure rate and repair rate, respectively. It means that the system is transited from the Up sate to the Down state based on the failure rate (λ) and similarly the system transition from the Down state to the Up state occurs according to the repair rate (µ). [40,41]. In this paper, the two-state reliability model has been used for both the cyber and physical networks' elements. Hence, the failures in cyber and physical networks could be determined using the two-state reliability model. In Figure 2, λ and μ denote the failure rate and repair rate, respectively. It means that the system is transited from the Up sate to the Down state based on the failure rate ( λ ) and similarly the system transition from the Down state to the Up state occurs according to the repair rate ( μ ).  In this paper, the exponential distribution has been used to simulate the Up and Down states [43]. The availability and unavailability probabilities have been mathematically expressed in (2) and (3).
( ) 1 exp 1: To simulate the stochastic time-series of CPMG elements' availability , the TTFs and TTRs should be generated using the MCS, as shown in (4) and (5). As seen, the mean time to failure (MTTF), mean time to repair (MTTR), and random variables uniformly distributed in (0 1) are utilized to generate the TTF and TTR.
( ) Figure 2. Two-state reliability model for CPMG's elements [16]. [40,41]. In this paper, the two-state reliability model has been used for both the cyber and physical networks' elements. Hence, the failures in cyber and physical networks could be determined using the two-state reliability model. In Figure 2, λ and μ denote the failure rate and repair rate, respectively. It means that the system is transited from the Up sate to the Down state based on the failure rate ( λ ) and similarly the system transition from the Down state to the Up state occurs according to the repair rate ( μ ).

Figure 2.
Two-state reliability model for CPMG's elements [16]. In this paper, the exponential distribution has been used to simulate the Up and Down states [43]. The availability and unavailability probabilities have been mathematically expressed in (2) and (3).
To simulate the stochastic time-series of CPMG elements' availability , the TTFs and TTRs should be generated using the MCS, as shown in (4) and (5). As seen, the mean time to failure (MTTF), mean time to repair (MTTR), and random variables uniformly distributed in (0 1) are utilized to generate the TTF and TTR.
In this paper, the exponential distribution has been used to simulate the Up and Down states [43]. The availability and unavailability probabilities have been mathematically expressed in (2) and (3).
To simulate the stochastic time-series of CPMG elements' availability MS i (t) t = 1 : 8760 × T , the TTFs and TTRs should be generated using the MCS, as shown in (4) and (5). As seen, the mean time to failure (MTTF), mean time to repair (MTTR), and random variables uniformly distributed in (0 1) are utilized to generate the TTF and TTR.
The cyber-physical interdependencies matrix (CPIM) at each time step could be defined using (6). As shown, the CPIM consists of power interconnections, cyber interconnections, and cyber-power interdependencies. In this paper, this matrix is used to create the graph of the CPMG.
Mapping the cyber failures to the physical/power network is essential for the reliability assessment of CPMGs. Hence, in the proposed method, the cyber-physical interdependencies are first concerned with mapping the unavailability of interconnected cyber elements to the physical layer by (7). If there are some cyber-power links between cyber elements and a physical element, all interconnected cyber elements' appropriate performance and availability are necessary. The physical element state should be changed due to any failure in interconnected cyber elements.
Afterward, the CPIM should be updated due to failures in power and cyber elements because the system graph and its interconnections are affected by these failures. In (8), updating the interconnections of the CPMGs is shown.
Finally, the power channels and information channels should be determined using (9)- (11). The shortest path function (SPF) based on theory graph concepts has been utilized to distinguish power and information channels. If there is a path between the power sending terminal (PST) and the power receiving terminal (PRT), the summation of the shortest path would not be zero, and it means that there is a power channel (Pch), as shown in (10). Otherwise, the power channel index would be assigned to zero, and no power channel exists between discussed terminals. Similarly, the state of the information channel (Ich) between the cyber sending terminal (CST) and the cyber receiving terminal (CRT) is determined using (11).

Information Transmission Errors and Delays
To assess the reliability of CPMGs, it is essential to concern the routing errors, data transmission delays, and dynamic transmission errors, as shown in Figure 4.

Information Transmission Errors and Delays
To assess the reliability of CPMGs, it is essential to concern the routing errors, data transmission delays, and dynamic transmission errors, as shown in Figure 4. The ideal information channel could be presented using (12).
, 1 : The transmission with delay in an information channel should be simulated according to (13). The delay might affect the information transmission if it exceeds the maximum allowed delay margin.
, 1 : The time delay in information transmission is probabilistic. The exponential probability distribution is one of the most well-known distributions, and has been used for simulating the statistical behaviors of time delays. The time delay at each time step could be simulated using the MCS based on (14). As demonstrated in (14), a random variable uniformly distributed in (0 1) is generated, and the mean value of time delays based on measured actual historical data is used to simulate the stochastic time delays by exponential probability distribution. The ideal information channel could be presented using (12).
The transmission with delay in an information channel should be simulated according to (13). The delay might affect the information transmission if it exceeds the maximum allowed delay margin.
The time delay in information transmission is probabilistic. The exponential probability distribution is one of the most well-known distributions, and has been used for simulating the statistical behaviors of time delays. The time delay at each time step could be simulated using the MCS based on (14). As demonstrated in (14), a random variable uniformly distributed in (0 1) is generated, and the mean value of time delays based on measured actual historical data is used to simulate the stochastic time delays by exponential probability distribution.
The data packet from the LC/MC to the MGCC, such as other data packets in other applications, consists of payload and header [44,45]. The first part of the data packet (payload) contains commands and information that must be transmitted, such as current, voltage, frequency, and so on. The second part (header) contains the destination address.
For instance, in a typical data packet (PAC sent PAC f send (t), DES f sent ), the transmission destination and the data contents are represented by PAC f send (t) and DES f sent , respectively. The information transmission without any errors in payloads and headers of information channels is shown in (13). However, the errors in payloads E f Payload (t) and headers E f Header (t) should be addressed to precisely calculate the reliability of CPMGs. The Gaussian distribution of white noise has been reported to study the statistical behaviors of payload and header errors in transmission information. The information transmission considering time delay, payload errors, and header errors should be simulated by the MCS. The white noise and a random variable uniformly distributed in (0 1) should be generated to simulate the payload and header errors, as mathematically expressed in (15) and (16) [16].

Stochastic Simulation of Renewable-Based DGs
The wind speed and technical characteristics influence the output energy of WTs. Since the wind speed is an uncertain environmental parameter, its stochastic behavior should be simulated using the MCS. The Weibull distribution has been selected in this paper to model the statistical behaviors of wind speed, which is one of the most common distributions for simulating the wind speed's stochastic nature [46]. As expressed in (17), the wind speed is simulated by the Weibull cumulative density function (CDF) inverse and a random variable in the proposed MCS-based method. Moreover, the WT output power is computed by (18).
The output power of PV units depends on the amount of solar radiation and the solar clearness index (SCI). Since the distribution of SCI mainly follows the Beta probability distribution [47,48], the MCS simulates the SCI by the Beta CDF inverse based on historical data. The Beta PDF and simulating the SCI using the MCS have been presented in (19) and (20) [42].
Finally, the output power of PV units should be distinguished using (21)-(24) based on environmental parameters and technical characteristics [48,49].

Stochastic Simulation of Loads
In this paper, the uncertainties of loads have been concerned. The normal probability distribution is one of the most common probability distributions to model the uncertainties of load values. The probabilistic load values are simulated using the normal CDF inverse by the MCS. Modeling the uncertainties of the CPMG's load is one of this study's contributions, which has not received a great deal of attention in current studies.

Reliability Calculations
The power balance condition should be met as depicted in (25). The adequacy assessment should be performed to check the values of the supply side, demand side, and their difference according to (26). Consequently, the load curtailment and at-risk state at each time step could be distinguished by (27) and (28) [50].
The reliability indices are determined using the load curtailment of all time steps, according to (29)- (31).
The flowchart of the proposed reliability assessment of CPMGs considering cyberphysical interdependencies, cyber failures, routing errors, time delays in information transmission, and other physical uncertainties, such as the uncertainties of renewable power output, is shown in Figure 5. In this paper, satisfying the number of MCS iterations and the accuracy level have been defined for the stopping criteria of the proposed MCSbased reliability evaluation method.
The flowchart of the proposed reliability assessment of CPMGs considering cyberphysical interdependencies, cyber failures, routing errors, time delays in information transmission, and other physical uncertainties, such as the uncertainties of renewable power output, is shown in Figure 5. In this paper, satisfying the number of MCS iterations and the accuracy level have been defined for the stopping criteria of the proposed MCSbased reliability evaluation method.

Test Results and Discussions
In this paper, the physical structure of the understudy CPMG is shown in Figure 6 [16].

Test Results and Discussions
In this paper, the physical structure of the understudy CPMG is shown in Figure 6   The selected CPMG has been equipped with a 500 kWh ESS. The maximum charging and discharging rate of the ESS have been assumed to be 10 kW. The MGT's maximum and minimum allowed output power and operating range are 100 kW and 10 kW, respectively.
The stochastic behaviors of load value at each time step have been simulated based on 15 min load values for three years, measured by the Kashan Electrical Distribution Company. The Sanjesh Afzar Asia (a manufacturer in Tehran, Iran) measurement instruments have been used for collecting the historical data of load values. The historical wind speed and SCI data used in this study are based on actual measurements in a city located in the center of Iran [51]. In Figure 7, typical one-year historical data of load, solar irradiance, and wind speed have been shown.
The reliability specifications of the understudy CPMG's elements are given in Table  3 [36,[52][53][54] and Table 4 [13,16,55,56]. Also, in Table 5, the information transmission errors based on data of [36,[52][53][54] has been demonstrated.  The selected CPMG has been equipped with a 500 kWh ESS. The maximum charging and discharging rate of the ESS have been assumed to be 10 kW. The MGT's maximum and minimum allowed output power and operating range are 100 kW and 10 kW, respectively.
The stochastic behaviors of load value at each time step have been simulated based on 15 min load values for three years, measured by the Kashan Electrical Distribution Company. The Sanjesh Afzar Asia (a manufacturer in Tehran, Iran) measurement instruments have been used for collecting the historical data of load values. The historical wind speed and SCI data used in this study are based on actual measurements in a city located in the center of Iran [51]. In Figure 7, typical one-year historical data of load, solar irradiance, and wind speed have been shown.  The reliability specifications of the understudy CPMG's elements are given in Table 3 [36,[52][53][54] and Table 4 [13,16,55,56]. Also, in Table 5, the information transmission errors based on data of [36,[52][53][54] has been demonstrated. In addition, to evaluate the different cyber network topologies' impacts on CPMG reliability, the four types of cyber network topologies are considered, as depicted in Figure 8. Topologies 1, 2, and 3 have been extracted from [16], and the fourth topology has been suggested to study the impacts of MGCC and its redundancy. kW and a WT with a maximum capacity of 25 kW are considered as the hybrid renewable DGs of the CPMG.
In addition, to evaluate the different cyber network topologies' impacts on CPMG reliability, the four types of cyber network topologies are considered, as depicted in Figure  8. Topologies 1, 2, and 3 have been extracted from [16], and the fourth topology has been suggested to study the impacts of MGCC and its redundancy.  The first cyber network topology shown in Figure 8 is the base network without any redundancy. In the second topology, the redundancy for the switches and optic fibers has been considered. In the third cyber network topology, in addition to redundant switches and optic fibers, the redundant LCs and MCs have been considered. In current research works such as [16], the redundancy of CPMGs has not been considered, and is studied in the fourth suggested topology.
In the selected test system, the battery is responsible for the initial control, and the MGT is responsible for generating the power shortage in addition to the renewable DGs. The control and operation cycle has been simulated continuously for 1 min time steps in the proposed MCS-based reliability evaluation method. The iteration numbers of 1000 MCS have been simulated to guarantee the convergence of the simulation-based calculations.
In Table 6, the proposed method's reliability evaluation results under various scenarios and cyber network topologies are presented. The reliability indices without considering the cyber failures and other cyber interruptions are demonstrated.  The first cyber network topology shown in Figure 8 is the base network without any redundancy. In the second topology, the redundancy for the switches and optic fibers has been considered. In the third cyber network topology, in addition to redundant switches and optic fibers, the redundant LCs and MCs have been considered. In current research works such as [16], the redundancy of CPMGs has not been considered, and is studied in the fourth suggested topology.
In the selected test system, the battery is responsible for the initial control, and the MGT is responsible for generating the power shortage in addition to the renewable DGs. The control and operation cycle has been simulated continuously for 1 min time steps in the proposed MCS-based reliability evaluation method. The iteration numbers of 1000 MCS have been simulated to guarantee the convergence of the simulation-based calculations.
In Table 6, the proposed method's reliability evaluation results under various scenarios and cyber network topologies are presented. The reliability indices without considering the cyber failures and other cyber interruptions are demonstrated.
Moreover, the convergence diagrams of the MCS-based methods under Scenarios 1, 2, and 3 (PV and WT-based and hybrid renewable-based DGs scenarios) are shown in Figures 9-11. As depicted, the expected energy not-supplied (EENS) of the proposed MCS-based method under different cyber network topologies and DG scenarios has been converged completely.

MCS-based method under different cyber network topologies and DG scenarios has been converged completely.
Test results infer that the CPMG under the first scenario (WT-100) has better reliability indices. The WT generation at all time steps and suitable wind speeds results in better adequacy and reliability indices under Scenarios 1 and 3 compared to Scenario 2 (PV-100). Selecting an appropriate CPMG structure based on studies under different scenarios leads to an improvement in system reliability.
Regardless of the renewable DG scenarios, the most impacts of the cyber network have occurred in Topology 1 without any redundancy. Test results highlight the importance of cyber impacts on CPMG's reliability. The 55.21%, 51.24%, and 50.3% EENS increments have occurred due to the cyber impacts for Topology 1 under scenarios 1, 2, and 3, respectively. It could be concluded that the cyber impacts for Topology 1 are highlighted, while the CPMG reliability and adequacy levels are satisfying (such as Scenario 1) because of the effectiveness of renewable DGs.
Moreover, test results shown in Table 6 imply that the appropriate cyber network (such as cyber network with Topology 4) mitigates the eventual cyber impacts. As revealed by the obtained test results, the cyber impacts under Scenario 1 and Topology 4 have been reduced to 1.5% compared to 55.21% (corresponding to Topology 1).   Test results infer that the CPMG under the first scenario (WT-100) has better reliability indices. The WT generation at all time steps and suitable wind speeds results in better adequacy and reliability indices under Scenarios 1 and 3 compared to Scenario 2 (PV-100). Selecting an appropriate CPMG structure based on studies under different scenarios leads to an improvement in system reliability.
Regardless of the renewable DG scenarios, the most impacts of the cyber network have occurred in Topology 1 without any redundancy. Test results highlight the importance of cyber impacts on CPMG's reliability. The 55.21%, 51.24%, and 50.3% EENS increments have occurred due to the cyber impacts for Topology 1 under scenarios 1, 2, and 3, respectively. It could be concluded that the cyber impacts for Topology 1 are highlighted, while the CPMG reliability and adequacy levels are satisfying (such as Scenario 1) because of the effectiveness of renewable DGs.
Moreover, test results shown in Table 6 imply that the appropriate cyber network (such as cyber network with Topology 4) mitigates the eventual cyber impacts. As revealed by the obtained test results, the cyber impacts under Scenario 1 and Topology 4 have been reduced to 1.5% compared to 55.21% (corresponding to Topology 1).  The EENS increments due to different types of cyber failures and interruptions under various scenarios and topologies are depicted in Figure 12. The EENS value of the CPMG without considering the cyber impacts has been considered the base condition, and the EENS increases due to all cyber impacts under various scenarios have been determined. Finally, the percentage of each cyber impact has been evaluated. The MGCC impacts under all scenarios, while Topologies 1, 2, and 3 have been selected, are significantly affected by other cyber failures and interruptions. In contrast, under Topology 4 of the cyber network with the redundancy of the MGCC, the MGCC failure did not dramatically affect the system reliability. It is therefore concluded that the MGCC failures are important. However, it is possible to mitigate or decrease the negative impacts of MGCC failures on the CPMG reliability by selecting an appropriate cyber network topology.
Another notable result is the EENS increments due to routing errors and information terminal failures. Therefore, it is essential to consider these types of errors and failures in any reliability evaluation of CPMGs that considers different cyber failures and interruptions. The EENS increments due to different types of cyber failures and interruptions under various scenarios and topologies are depicted in Figure 12. The EENS value of the CPMG without considering the cyber impacts has been considered the base condition, and the EENS increases due to all cyber impacts under various scenarios have been determined. Finally, the percentage of each cyber impact has been evaluated. The MGCC impacts under all scenarios, while Topologies 1, 2, and 3 have been selected, are significantly affected by other cyber failures and interruptions. In contrast, under Topology 4 of the cyber network with the redundancy of the MGCC, the MGCC failure did not dramatically affect the system reliability. It is therefore concluded that the MGCC failures are important. However, it is possible to mitigate or decrease the negative impacts of MGCC failures on the CPMG reliability by selecting an appropriate cyber network topology.
Another notable result is the EENS increments due to routing errors and information terminal failures. Therefore, it is essential to consider these types of errors and failures in any reliability evaluation of CPMGs that considers different cyber failures and interruptions. Sustainability 2021, 13, x FOR PEER REVIEW 20 of 32 (a) Scenario 1.
(c) Scenario 3. In Table 7, to illustrate the EENS decrements that modify the cyber network topologies, the EENS and LOLP increments due to cyber impacts for Topology 1 have been considered as the base condition, and the EENS and LOLP decrements for other network topologies compared to Topology 1 have been identified. Comparative test results infer that In Table 7, to illustrate the EENS decrements that modify the cyber network topologies, the EENS and LOLP increments due to cyber impacts for Topology 1 have been considered as the base condition, and the EENS and LOLP decrements for other network topologies compared to Topology 1 have been identified. Comparative test results infer that the cyber network topology affects the CPMG reliability effectively. Regardless of the DG scenarios, approximately 97% of negative cyber impacts could be omitted using an appropriate cyber network topology, such as Topology 4. The comparison of test results under Topologies 2 and 3 with Topology 4 highlights the importance of the MGCC failures. As observed, the redundancy of the MGCC to mitigate its eventual failures under Topology 4 results in 33%, 28%, and 35% EENS improvement. To illustrate the contributions of the proposed research, the reliability results of the understudy CPMG that do not consider the MGCC faults, but rather consider other cyber impacts, such as in [16], are presented in Table 8. The differences in test results based on the proposed method and other available ones (without considering the MGCC faults) have been demonstrated in Table 9. The test results infer that under cyber network topologies 1, 2, and 3, considerable inaccuracy occurs in the reliability evaluation of the CPMG without considering the MGCC faults, and the advantages of the proposed method are highlighted. Moreover, the results under Topology 4 show that the impacts of MGCC faults could be completely omitted by selecting an appropriate cyber network.  To further evaluate the impacts of each cyber element's failure on the CPMG reliability, a sensitivity analysis has been performed to obtain insight into how the changes in failures of the cyber elements affect the CPMG reliability, as shown in Figure 13. In this analysis, the failure rates of cyber elements have gradually increased up to ten times under various scenarios.
As revealed by test results shown in Figure 13, regardless of the DG scenarios, the EENS under Topology 1 is very sensitive against the changes in MGCC failure rate. In addition, the switch and information terminal failures could influence the CPMG reliability if they would change. The failure increment of optical fibers under Topology 1 cannot negatively influence the CPMG reliability because its failure rate is very low, and its impacts are negligible compared to other element faults.
As seen, the sensitivity of the CPMG's EENS due to the faults of switches has been solved under Topology 2. However, the changes in the MGCC and information terminal failure rates lead to a significant increment of the CPMG's EENS. Accordingly, it is essential to modify the cyber network topology. Although the sensitivity of the CPMG's EENS due to failures of switches and information terminals has been mitigated under Topology 3, the MGCC's failure increment might affect the CPMG's EENS. This means that it is not possible to calculate the CPMG's reliability accurately without considering the MGCC faults such as in [16]. Moreover, the test results under Topology 4 emphasize the effectiveness of using an appropriate cyber network topology by the MGCC redundancy.
In order to evaluate the impact of data transmission quality on reliability indices, sensitivity analyses based on changes in the probability of routing errors have been performed, as depicted in Figure 14. As seen, the probability of routing errors has been gradually increased from 0.0001 to 0.001, and the CPMG's EENS under different cyber network topologies for Scenario 3 (PV and WT-based and hybrid renewable-based DGs) has been assessed.
The test results show that the impacts of routing errors are higher than payload errors under all cyber network topologies. In addition, the EENS increase due to changes in routing and payload errors does not dramatically depend on the cyber network topology.
Furthermore, to evaluate the effect of tolerable latency, a sensitivity analysis has been performed on the CPMG information channel's tolerable latency. In Figure 15, the maximum tolerable latency rate has been gradually increased from 1 to 10 s under different cyber network topologies and DG scenarios. Test results show that a maximum tolerable latency rate of the information channels higher than 2 s under Topologies 2, 3, and 4 is satisfactiory.
However, the impacts of the maximum tolerable latency rate of the information channels on CPMG's reliability are considerable. As seen, if it is desired to mitigate the impacts of delays of information channels on the CPMG's reliability, the maximum tolerable latency rate of the information channels should be greater than 6 s for Topology 1. Simulation results imply that the impacts of the maximum tolerable latency rate of the information channels depend on cyber network topology.

Conclusions
SMGs have received a great deal of attention in recent years because the advantages of DGs, particularly renewable DGs, are more achievable using SMGs. Although the cyber As seen, if it is desired to mitigate the impacts of delays of information channels on the CPMG's reliability, the maximum tolerable latency rate of the information channels should be greater than 6 s for Topology 1. Simulation results imply that the impacts of the maximum tolerable latency rate of the information channels depend on cyber network topology.

Conclusions
SMGs have received a great deal of attention in recent years because the advantages of DGs, particularly renewable DGs, are more achievable using SMGs. Although the cyber systems facilitate the energy management schemes in CPMGs, cyber faults, cyber interferences, and cyber-physical interdependencies might affect the CPMG's reliability. A research gap exists in the field of CPMG's reliability evaluation, simultaneously considering all cyber failures and interferences under various cyber network topologies and renewable distributions scenarios. This research tried to fill such a gap by developing a new MCS-based reliability assessment method considering cyber elements' failures, data/information transmission errors, and routing errors under various cyber network topologies. Developing the proposed method by applying the theory graph concepts is one of the contributions of this paper. It should be noted that MGCC faults, aside from other failures and interferences, have been concerned in the proposed method, which is one of the essential contributions of this study. The test results highlighted the advantages of the proposed method considering the MGCC faults' impacts on the CPMG's reliability. It has been concluded based on obtained simulation results that the negative impacts of cyber faults and interruptions could be decreased effectively by selecting an appropriate cyber network topology. The simulation results inferred that the negative impacts of cyberphysical interdependencies could be limited to less than 3% using a suitable cyber network topology with the redundancy of MGCC.

PP
Power-power interconnection CP Cyber-power interdependency PC Power-cyber interconnection CC Cyber-cyber interconnection PP Updated power-power interconnection after mapping the cyber failures to physical system CP Updated cyber-power interdependency after mapping the cyber failures to physical system PC Updated power-cyber interconnection after mapping the cyber failures to physical system CC Updated cyber-cyber interconnection after mapping the cyber failures to physical system PP Updated power-power interconnection after considering the destroyed interconnections due to cyber/physical failures CP Updated cyber-power interdependency after considering the destroyed interconnections due to cyber/physical failures PC Updated power-cyber interconnection after considering the destroyed interconnections due to cyber/physical failures CC Updated cyber-cyber interconnection α, β Beta probability density function (PDF) parameters