Towards Avoiding Cascading Failures in Transmission Expansion Planning of Modern Active Power Systems Using Hybrid Snake-Sine Cosine Optimization Algorithm

: In this paper, a transmission expansion planning (TEP) model is proposed to guarantee the resilience of power systems and mitigate cascading failures’ impacts. The energy storage systems and fault current limiters’ planning models are integrated into the TEP problem to minimize cascading outages and comply with short-circuit current reliability constraints. Most studies in the literature adopt a single strategy to simulate power systems’ cascading failures that may not be enough to guarantee networks’ resilience. This work elaborates on two scenarios for initiating cascading failures to study the impact of various initiating events on the planned system’s strength and the projects required. The TEP problem is formulated as a non-linear, non-convex large-scale problem. To avoid linearization issues and enhance meta-heuristics performance, a hybridization of two meta-heuristic techniques, namely snake optimizer and sine cosine algorithm (SO-SCA), is proposed to solve the problem. Two hybridization strategies are suggested to improve the exploration and exploitation stages. Deﬁning future loads growth is essential for TEP. Hence, a load forecasting technique based on SO-SCA is investigated and compared with some methods reported in the literature. The results obtained proved the efﬁciency of the proposed approach in predicting load growth. TEP’s calculations were carried out on the Garver and the IEEE 24-bus system. The results demonstrated the superiority of the hybrid SO-SCA in solving the TEP problem. Moreover, the projects required to expand networks differed according to the type of cascading failures’ initiating scenario.


Introduction
The transmission expansion planning (TEP) problem corresponds to allocating different power facilities, such as transmission lines, energy storage systems (ESSs), thyristorcontrolled series compensators (TCSC), and fault current limiters (FCLs), to satisfy the expected load growth and fulfill operational and reliability constraints [1]. TEP is a wellstudied problem in the literature; however, there are still many challenges related to solution methods and offering a mathematical model that can accommodate the power system's uncertainties and ensure reliability constraints [2].
Several mathematical, heuristic and meta-heuristic methods were reported in the literature to solve TEP models [2]. Some of these works are introduced below as follows. In Hamidpour et al. [3], the AC power flow-based model (ACTEP) was adopted. The non-linear and non-convex constraints were linearized through Taylor's expansion and polygon approaches. The CPLEX solver solved the linearized model within GAMS. The suggested strategy was conducted on the IEEE 6-bus and 30-bus test systems. The numerical results verified the efficiency of the proposed solution. Esmaili et al. [4] presented a space. The results showed the effectiveness of the suggested approach in enhancing the exploration and exploitation stages and getting high-quality solutions.
A heuristic approach based on fuzzy logic was developed by Sousa et al. [13] for TEP. The simulation results carried out on different test systems confirmed the proposed technique's effectiveness in finding high-quality solutions. In Armaghani et al. [14], an analytical multi-stage TEP approach combined with a forward pseudo-dynamic planning method was proposed for finding critical transmission lines that may initiate and propagate the cascading failure.
Recently, cascading failures' impacts during TEP have been taken into account in several studies [15][16][17]. Mehrtash et al. [15] contributed to a risk-based planning model that minimizes cascading failures' effects. The suggested problem was modeled as a mixedinteger second-order cone programming problem. Risk indices were considered in the developed model to measure the non-identical probability and individual contingencies' consequences. The results obtained from the IEEE 24-bus system demonstrated that the system planned was more robust against high-risk contingencies that may have resulted in cascading failures. Qorbani et al. [16] developed a resilient TEP model to mitigate cascading outages' effects and ensure N-1 security. The outage of two-connected lines to each node was the triggering event for cascading failures. Gjorgiev et al. [17] proposed a risk-informed TEP method. An AC-based cascading failure model was developed to simulate post contingencies power systems' responses. A multi-objective optimization was investigated and solved using a meta-heuristic algorithm. The suggested procedure succeeded in mitigating the effects of cascading failures.
Calculation of future peak loads is essential for TEP; hence the accuracy of the load forecasting tool implemented is necessary to guarantee the accuracy of the planning process. The hyper-parameters of the forecasting models play an essential role in improving the model's performance [18][19][20][21]. A long-term load prediction-based particle swarm optimization was implemented by AlRashidi et al. [21]. The results demonstrated the superiority of the applied optimizer in estimating forecasting models' parameters compared to the least error square method. An adaptive neural fuzzy system was recommended by Fathy et al. [22] and Refaat et al. [23] as an efficient tool for load forecasting. However, it had a problem selecting the type and the number of membership functions. Hong et al. [24] applied the support vector machines method to forecast future loads growth; however, it has parameters selection problems. In this work, a load forecasting approach based on a hybrid of two-meta heuristic techniques is proposed and compared with some methods presented in the literature.
Based on all the above discussion, the existing research gaps are illustrated as follows: • The mathematical methods efficiently solve linear and simple TEP problems. However, they need a high computational burden for solving large-scale non-linear TEP problems with a vast search space. These methods have difficulty finding the optimal solution, mainly if non-convexities exist with the search space. Most of the existing works apply a linearization scheme and a decomposing strategy to simplify the problem, but it may result in linearization errors that affect the accuracy of obtained solutions.

•
As known, the meta-heuristics are not transparent algorithms. They are time-consuming and challenging to reach the optimal solution in a single run [2]. Much effort is still needed for improving the performance of meta-heuristic techniques in solving TEP problems.

•
Although heuristic methods need a lower computational burden compared to mathematical and meta-heuristic algorithms, they cannot guarantee the optimal solution or high-quality solutions.

•
The hybridization of several optimization algorithms is recommended for solving the TEP problems. However, most current studies proposed it for solving the DCTEP problem that is simple compared to the AC models. Testing the performance of the hybrid algorithms in solving the ACTEP problems is rarely discussed in the literature. This work corresponds to introducing a robust TEP model and investigating a longterm load forecasting method to address the mentioned shortcomings. The ACTEP is adopted in this work. It handles short-circuit current constraints and cascading failures impacts. The proposed model is non-linear and non-convex. To avoid linearization, a hybrid of two meta-heuristic algorithms is developed. The main contributions of this study can be outlined as follows: • The meta-heuristic algorithm, Snake optimizer (SO) [25], is combined with the sine cosine algorithm (SCA) [26] to solve the proposed TEP problem and calculate the long-term load forecasting. Two combination approaches are adopted to enhance both algorithms' exploration and exploitation phases.

•
The use of renewable energy sources (RESs) threatens the modern power systems' resilience and may be a trigger to cascading failures and rolling blackouts [27,28]. Therefore, the adopted model deals with RESs' stochastic behavior.

•
Two scenarios for initiating cascading failures are suggested to ensure the planned system's strength and study the impact of various initiating events on the projects required.

•
The impact of ESSs and FCLs in minimizing cascading outages and fulfilling shortcircuit current constraints is also investigated. In this regard, Sodium sulfur batteries (NaS) and FCLs' sizing and placing problems are included in the TEP problem.
The remainder of this paper is systematized as follows. The problem's formulation is presented in Section 2, while the hybrid of SO-SCA is explained in Section 3. Section 4 shows and discusses the results obtained. Finally, in Section 5, the main remarks derived from this work are presented.

Problem Formulation
The TEP problem is modeled as an optimization problem to expand power networks, increase networks' resilience and achieve reliability constraints in the presence of RESs. The proposed objective function targets to minimize the cost of installed circuits, FCLs, generation units and ESSs.
The AC model is adopted in this study to formulate the proposed problem. It is superior to the DC model in terms of the accuracy of the results obtained and the ability to integrate many problems into a single model, such as combining the problem of reactive power units' allocation and the TEP into a single problem [27]. A clustered data of wind speed and load consumed over the year simulate power systems' stochastic behavior.
As mentioned before, TEP is always linked to future loads; therefore, implementing an accurate load forecasting technique is essential to improve the planning process. In this regard, the load forecasting problem is formulated as an estimation process, and the hybrid of SO-SCA is applied to minimize the error associated with the estimated model.

Mathematical Model of TEP
The proposed TEP model has a single objective function and is subjected to normal and security constraints, as illustrated in (1)- (29). Equation (1) provides the proposed cost function. The first term concerns the cost of installed circuits, while the second and third terms cover the cost of FCLs and NaS batteries required to achieve security and short-circuit current constraints. The fourth term covers the cost of new generation units installed. The reminder terms present the operating cost of NaS batteries and thermal units and RESs, respectively.
where the recovery factor λ(1+λ) y (1+λ) y −1 is inserted in the cost function to consider the money's time value [29], where λ and y are the discount rate and the lifetime of the project in years, respectively. The problem constraints are provided in (2)- (29). The steady-state constraints are shown in (2)-(24), whereas reliability constraints are given in (25)- (29). Constraints (2) and (3) are related to active and reactive powers injected at each node. P h ESS,i is negative if an ESS is in the discharging mode, and it is positive if an ESS is in the charging mode.
So that P h i and Q h i i are calculated as follows [30]: Constraints (6) and (7) ensure that active and reactive powers generated from thermal stations do not violate their limits. For RESs, it is assumed that they generate active power only, and the output power must be within the permissible bounds, as shown in (8). Constraints (9) and (10) guarantee that the buses' voltage and angle do not surpass their limits.
Equation (11) explains that any circuit built in any previous scenario will be present in the following scenarios. It also ensures that the number of circuits between any two nodes should be less than the maximum. Equations (12) and (13) show that the power flows in new and existing circuits do not exceed the thermal limits.
where S s,h ij and S r,h ij are calculated as follows: The other parameters are calculated as follows [30]: Constraints related to ESSs' planning are given in (20)- (24). Equation (20) restricts ESSs' number at each bus, whilst (21)-(23) maintain ESSs' power and capacity during charging and discharging modes to be lower than nominal values. ESSs' storage capacity at each hour is calculated using (24).
Two security scanning matrices, L1 and L2, are suggested as triggering events for the initiation of cascading failures. L1 simulates a power system under a single circuit's outage between two nodes. In comparison, L2 regards the outage of all circuits between any two nodes as the triggering event for cascading failures. N n represents total circuits installed Mathematics 2022, 10, 1323 7 of 25 at route number n. Systems planned must be able to fulfill the constraints (2)-(27) after contingencies L1 and L2.

Long-Term Load Forecasting Model
The proposed load forecasting problem is modeled as an estimation process, and the hybrid of SO-SCA is applied to find the optimal coefficients to reduce estimation errors. Linear, quadratic models are commonly used in the literature for forecasting load growth in the long term, and they are adopted in this work.
The historical data used in the forecasting process is a matrix (D) that comprises m rows and two columns. The first column represents the input data, whilst the second column represents the output data. This work uses the input data of historical years (y), and peak loads (P) are the output data. For m rows, the discrete state-space form represents the forecasting system is formulated as follows [21]: where P(y) is a vector (m × 1) that represents calculated outputs, while ϕ is a vector (n × 1) that represents estimated parameters. H(y) is a row vector representing the relation between P(y) and ϕ. ε(y) is the error vector results in the estimation process. In this study, the linear and the quadratic models are adopted. For the linear model: while for the quadratic models:

The Hybrid Snake-Sine & Cosine-Based Optimization Algorithm
The proposed problem is a large-scale non-linear optimization problem, and to solve it, a hybrid of two meta-heuristic techniques, namely SO-SCA, is developed. SO is recommended for solving several complex and challenging engineering problems. It is capable of converging fast to the global optimum. While, SCA is superior in solving unimodal, multi-modal, and composite test functions compared to other algorithms. SO and SCA have some issues finding optimal or high-quality solutions in one run. Moreover, it is challenging to decide on suitable design parameters for all problems because well-defined algorithm parameters for one situation may not be ideal for solving another problem. The hybrid SO-SCA is employed in this work to overcome these issues. Two combination approaches are proposed to enhance the exploration and exploitation stages and get high-quality solutions in as few runs as possible.

Snake Optimizer
SO was developed by Hashim et al. [25]. It is based on the idea of imitating snakes' mating behavior. Snakes commonly fight to get the best partners if they have enough food and a low ambient temperature. Like any meta-heuristic technique, SO initializes with a randomly generated population, and then the updating process is carried out in two phases: exploration and exploitation. During the updating process, it is assumed that the number of males equals the number of females. This process goes on until the maximum number of iterations (T) is reached.
In the exploration phase, it is assumed that there is not enough food, and snakes randomly move toward searching for food. The food's quantity (Q) is calculated by The exploration's behavior is modeled as follows: where X t+1 i,m and X t+1 i, f indicate to the ith male and female positions. X t rand,m and X t rand, f are the positions of random males and females selected from the population. t is the current iteration. rand is a random value and C 2 is a constant that set to 0.05. A m and A f measure respectively the male and female abilities to get food, and they are calculated ; respectively. f rand,m and f rand, f are the fitness values of X t rand,m and X t rand, f ; respectively. While f i,m and f i, f are the fitness values of ith male and female.
In the exploitation phase, it is assumed that the food is available, and the next positions are defined according to the ambient temperature (TEMP). TEMP is defined by exp −t T . If the temperature is high and exceeds the threshold value, the new position for males and females is located using (35). Where X t f ood is the best individuals position and C 3 equals 2.
Otherwise, the snakes are in the fighting or mating modes and are randomly switched between both modes. The new positions are calculated using (36) and (37) for the fighting mode and using (38) and (39) for the mating mode.
where FM and FF are males' and females' fighting abilities and calculated by exp

Sine Cosine Algorithm
SCA was proposed by Mirjalili et al. [26]. It is based on the use of sine and cosine functions to vary new positions outward and toward the best positions in the hope of enhancing the exploitation and explorations phases. SCA starts with a randomly generated population, and then the updating scheme given in (40) is employed to define the new position of each individual in the population. Where r 2 , r 3 and r 4 are randomly defined, and r 1 is calculated by C − C t T . C is constant.

Proposed Hybrid Scheme
Two combination approaches for SO-SCA are proposed. The pseudo-code of the first approach is presented in Algorithm 1. SO and SCA operate in parallel to calculate new variables in this approach. According to Algorithm 1, candidate solutions are randomly generated, and their corresponding finesses are calculated. Individuals who violate the steady-state and security constraints are penalized by adding additional costs. After that, the best solution and best fitness value are defined. In the beginning, the population is divided into two groups that have the same number of individuals. New solutions for the first group are calculated using SO's updating scheme provided in (33)-(39). While new positions of individuals in the second group are obtained using SCA's updating scheme given in (40). The best solution for both groups is decided. If the SO gives the best solution, more individuals will be assigned for SO in the next iteration, relying on the switching rate (N rate ). Otherwise, more individuals will be assigned for SCA. Set SO and SCA parameters.

5:
Calculate the cost function and check the problem's constraints. 6: Add a penalty cost to solutions that violate the constraints. 7: Define the best and the worst solutions.  In the second approach, a series combination is applied, as given in Algorithm 2. The population is first updated using SCA's updating scheme. Then, the position of the best solution and its fitness is updated regarding the output of SCA. After that, the SO updating scheme is applied to calculate the position of new variables, and the best solution's position is defined based on the output of SO. Set SO and SCA parameters.

4:
Calculate the cost function and check the problem's constraints.

5:
Add a penalty cost to solutions that violate the constraints. 6: Define the best and the worst solutions. 7: For t = 2:T 8: For i = 1 : N p 9: Define the new position of each individual (X t+1 i,SCA ) in X SCA using (40). 10: Calculate the fitness value of each new solution ( f SCA i ).

11:
Check the problem constraints and add a penalty cost if the new solution violates the constraints. 12: end for 13: Find the best and the worst solutions obtained using SCA (X SCA Best , f SCA Best ). 14: Update the positions of the candidate solution (X t i = X t+1 i,SCA ); 15: Update the best solution and the best fitness. 16:

Test Systems
The proposed model and solution method were examined on the Garver network and the IEEE 24-bus system. A data set introduced in [21] was used to investigate the adopted load forecasting algorithm. The wind speed and load consumed depicted in Figure 1 were used to validate the proposed TEP model. Investment and operating factors of different facilities used are given in Table 1. The results were obtained using the MATLAB r2017a platform.

Test Systems
The proposed model and solution method were examined on the Garver network and the IEEE 24-bus system. A data set introduced in [21] was used to investigate the adopted load forecasting algorithm. The wind speed and load consumed depicted in Figure 1 were used to validate the proposed TEP model. Investment and operating factors of different facilities used are given in Table 1. The results were obtained using the MATLAB r2017a platform.    Wind unit [31] 0.139 × 10 6 $/MW 0.232 0 Thermal units [31] 0.536 × 10 6 $/MW 0.123 0.144 FCL 0.5 × 10 6 $/p.u NA NA NaS [29] 0.446 × 10 6 $/MW NA 0.2981 The Garver network is a benchmark system and is widely used for investigating the efficiency of TEP's models. It has six buses, three-generation units, and five load buses [22], [32]. Bus 5 was suggested as a candidate position for new wind units, and buses 2, 5, and 6 were suggested for new ESSs. The network's initial configuration is shown in Figure 2. New candidate routes are presented in dash lines. The IEEE-24 bus system has 24 buses, eleven generation units, and seventeen load buses. Some modifications to the system were made in [32] to make congestions. The candidate locations for new wind units were buses 3, 10, and 19. While the candidate locations for ESSs were buses 2, 6, 8, 14, 20, 21, 22, and 24. The initial configuration of the system is depicted in Figure 3. The IEEE-24 bus system has 24 buses, eleven generation units, and seventeen load buses. Some modifications to the system were made in [32] to make congestions. The candidate locations for new wind units were buses 3, 10, and 19. While the candidate locations for ESSs were buses 2, 6, 8, 14, 20, 21, 22, and 24. The initial configuration of the system is depicted in Figure 3. Mathematics 2022, 10, x FOR PEER REVIEW 13 of 26

Testing the Performance of SO-SCA
Figures 4-7 show the convergence curves of the hybrid algorithm using the two approaches, SO and SCA, over ten runs. Noteworthy, all the simulations in this section were carried out without embedding L1 and L2 securities. The population size was 30 for SO-SCA #2 and 40 for SO-SCA #1, SO, and SCA. It is shown that the series approach effectively solved the Garver network in terms of standard deviation and quality of solutions compared to the parallel approach, as shown in Table 2. For the IEEE 24-bus system, SO-SCA #1 gave a lower standard deviation than SO-SCA #2, as illustrated in Table 3. Regarding SO and SCA, SCA achieved better performance than SO for the 24-bus system, but SO was superior to SCA in solving the Garver network.  Figures 4-7 show the convergence curves of the hybrid algorithm using the two approaches, SO and SCA, over ten runs. Noteworthy, all the simulations in this section were carried out without embedding L1 and L2 securities. The population size was 30 for SO-SCA #2 and 40 for SO-SCA #1, SO, and SCA. It is shown that the series approach effectively solved the Garver network in terms of standard deviation and quality of solutions compared to the parallel approach, as shown in Table 2. For the IEEE 24-bus system, SO-SCA #1 gave a lower standard deviation than SO-SCA #2, as illustrated in Table 3. Regarding SO and SCA, SCA achieved better performance than SO for the 24-bus system, but SO was superior to SCA in solving the Garver network.         Based on all results, it can be concluded that the hybrid algorithm is efficient at solving TEP problems compared to the individual algorithms like SO and SCA. Both hybridization frameworks moved SO and SCA in each iteration to a better search area exploiting the best solution obtained by one of them. The results revealed that the series approach is recommended for solving small-scale systems, while the parallel process is preferred in solving larger systems.

Testing the Performance of SO-SCA-Based Load Forecasting Technique
The performance of SO-SCA-the based load forecasting technique was compared with PSO-, SCA-, SO, and least error squares (LES) estimation-based methods. The estimated parameters are given in Tables 4 and 5. The predicted load using all algorithms is depicted in Figure 8.  Tables 4 and 5 show that SO-SCA#1 and SO-SCA#2 gave the least average error of 3.803074% for the linear model and 3.323598% for the quadratic model. The results also illustrated that SCA was superior to SO, PSO and LES in terms of average error. The linear model was 3.832615% for SCA, 3.8336% for SO, and 3.8336% and 4.2103 for PSO-and LES-based techniques. While for the quadratic model, it was 3.3375%, 3.3488%, 3.3488%, and 3.3873% for SCA-, SO-, PSO-, LES-based techniques.
The results demonstrated that the quadratic model was more suitable for fitting the testing data than the linear model because the minimum average error was 3.803074% for the linear model and 3.323598% for the quadratic model. zation frameworks moved SO and SCA in each iteration to a better search area exploiting the best solution obtained by one of them. The results revealed that the series approach is recommended for solving small-scale systems, while the parallel process is preferred in solving larger systems.

Testing the Performance of SO-SCA-Based Load Forecasting Technique
The performance of SO-SCA-the based load forecasting technique was compared with PSO-, SCA-, SO, and least error squares (LES) estimation-based methods. The estimated parameters are given in Tables 4 and 5. The predicted load using all algorithms is depicted in Figure 8.   5 show that SO-SCA#1 and SO-SCA#2 gave the least average error of 3.803074% for the linear model and 3.323598% for the quadratic model. The results also illustrated that SCA was superior to SO, PSO and LES in terms of average error. The linear model was 3.832615% for SCA, 3.8336% for SO, and 3.8336% and 4.2103 for PSO-and LES-  Table 6 details the results obtained using the proposed planning scheme. The results showed a significant change for projects required to achieve the L1 and L2 security constraints. The maximum number of circuits built in each route did not exceed 2 when the triggering event L2 was applied. In contrast, it did not surpass 3 for the triggering event L1. The total number of circuits installed was 7 for the event L1 and 10 for the event L2. The results emphasized that to achieve L1 security, there was a need for a new circuit at routes 3-5 and 5-6, two circuits at route 4-6, and three circuits at route 2-6. While meeting L2 security required new circuits at routes 2-3, 1-6, 3-4, 3-6, 4-6, and 5-6.  The results showed that the use of ESSs was essential to maintain power continuity and avoid cascading failures. The results indicated that NaS's capacities required for the events L2 and L1 were 198.35 and 192.79, respectively. The best locations for the ESSs, to achieve L1 security were at nodes 2 and 5. Whilst the locations of ESSs needed to achieve L2 security were at buses 2, 5, and 6. Figure 9 shows the short circuit current at scenario number24 as an example to verify that it did not exceed the maximum limit (6.5 p.u) and that the short-circuit constraints were achieved in each scenario. It is worth mentioning that the total size of FCLs required ranged between 0.59 + 5.9i p.u and 1.24 + 12.4i p.u. The results also illustrated that the buses' voltage did not violate their limits (0.9-1.1 p.u). Figure 10 points to the voltage at each bus for scenario number 24 as an example.

TEP
3-5 (1), 2-6 (3), 4-6 (2), 5-6 (1) 2-3 (1), 1-6 (2),3-4 (1), 3-6 (2), 4-6 (2), 5-6 (2) The results showed that the use of ESSs was essential to maintain power continuity and avoid cascading failures. The results indicated that NaS's capacities required for the events L2 and L1 were 198.35 and 192.79, respectively. The best locations for the ESSs, to achieve L1 security were at nodes 2 and 5. Whilst the locations of ESSs needed to achieve L2 security were at buses 2, 5, and 6. Figure 9 shows the short circuit current at scenario number24 as an example to verify that it did not exceed the maximum limit (6.5 p.u) and that the short-circuit constraints were achieved in each scenario. It is worth mentioning that the total size of FCLs required ranged between 0.59 + 5.9i p.u and 1.24 + 12.4i p.u. The results also illustrated that the buses' voltage did not violate their limits (0.9-1.1 p.u). Figure 10 points to the voltage at each bus for scenario number 24 as an example.  The previous results imply that the triggering event L2 was the worst, and the probability of occurring complete blackouts was high. The cost of projects needed for fulfilling L2 security exceeded 17.41 % compared to the cost of L1 projects.

Testing Results of the IEEE 24-Bus System
Similar to the Garver network, simulations on the 24-bus system confirmed the gravity of the L2 event compared to the L1 event. This is clearly shown in the size of ESSs required and the planning cost. There was an increase of about 108% in the size of ESSs to meet L2 security. The results pointed to the size of NaS batteries required for ensuring L1 security was MWh 29.43 at bus 6. Whilst 60.123 and 1.172 MWh of NaS batteries at The previous results imply that the triggering event L2 was the worst, and the probability of occurring complete blackouts was high. The cost of projects needed for fulfilling L2 security exceeded 17.41 % compared to the cost of L1 projects.

Testing Results of the IEEE 24-Bus System
Similar to the Garver network, simulations on the 24-bus system confirmed the gravity of the L2 event compared to the L1 event. This is clearly shown in the size of ESSs required and the planning cost. There was an increase of about 108% in the size of ESSs to meet L2 security. The results pointed to the size of NaS batteries required for ensuring L1 security was MWh 29.43 at bus 6. Whilst 60.123 and 1.172 MWh of NaS batteries at buses 6 and 21 were needed to fulfil L2 security.
According to Table 7, new circuits were installed at the routes 1-2, 1-3, 2-4, 2-6, 7-2, 7-5, and 7-8 to guarantee L1 and L2 securities. It is worth emphasizing that the reliance on the location of new circuits was not enough to achieve the short circuit-current constraints. The use of FCLs was necessary. The short-circuit current at scenario number 24, as an example, is depicted in Figure 11. The results indicated that the planning model was able to maintain the buses' voltage within their limits. They did not violate the upper and lower bounds (0.9-1.1 p.u), as shown in Figure 12.  Figure 11. Short-circuit current at scenario number 24 for the 24-bus system.

Figure 12.
Buses' voltage at scenario number 24 for the 24-bus system. Table 7. TEP's projects for the 24-bus system.   Based on all previous findings, it can be concluded that relying on a single strategy to avoid the cascading failures of the electrical networks does not guarantee their security and resilience. Two simple strategies were applied in this work, which significantly affected the size and number of projects required and hence the cost. For the Garver network, the impact was clearly shown in the number and the location of new circuits. While for the 24-bus system, the size of ESSs was affected by the type of cascading failures The TEP problem is complex, and the hybrid of two meta-heuristic algorithms effectively solved TEP models compared to individual algorithms. The two hybridization frameworks were superior in terms of standard deviation and quality of solutions. Further, the hybrid algorithm was more accurate in forecasting the electrical loads in the long term than some methods introduced in the literature.

Conclusions
In this work, a resilience TEP model was suggested to ensure networks' reliability and supply the future electrical demands. The short-circuit current and the cascading failure constraints were considered to achieve the resilience requirements. NaSs and FCLs' planning models were incorporated into the proposed TEP model to achieve these targets. Two simple scanning matrices, L1 and L2, were modeled as triggering events for cascading failures. L1 simulated a power system under a single circuit's outage between two nodes. L2 focused on the outage of all circuits between two nodes as the triggering event for cascading failures. The results showed that simulating a single strategy for cascading failures of power systems does not guarantee their security and resilience. The triggering event L2 was the worst. The probability of occurring complete blackouts was high. The impact was clearly shown in the number and the location of new circuits for the Garver network, while it affected the size of ESSs installed for the IEEE 24-bus system.
The proposed TEP problem was a large-scale non-linear optimization problem, and a hybrid of two meta-heuristic techniques, namely SO and SCA, was adopted. Two hybridization approaches were applied to enhance the exploration and exploitation stages. SO and SCA operated parallel in the first approach to calculate new variables in each iteration. While in the second approach, SO and SCA operated in a series mechanism. The results revealed that the second approach effectively solved small-scale systems, whereas the first approach efficiently solved larger systems.
Further, the results proved that both approaches were superior in terms of standard deviation and quality of solutions compared to individual algorithms. A long-term load forecasting technique based on SO-SCA was proposed. It was shown that the suggested method gave accurate results compared to PSO, SO, SCA, and LES.
This work is only concerned with investigating two simple cascading failure events, but more events that pose a greater threat to the network should be considered to get a secure network. More studies are required to evaluate the impact of including the transient stability constraints on enhancing power systems' securities and resilience. In addition, the planning model needs to be extended to study the effects of the optimal mix of transmission lines, distributed generators, FCLs, ESSs, reactive source compensators, demand response, and thyristor-controlled series compensators on improving power systems' securities and avoiding cascading failures. The performance of the hybrid SO-SCA should also be tested on larger systems with an extensive search space.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to their large size. 1443). The Author also acknowledges the support provided by King Abdullah City for Atomic and Renewable Energy (K.A.CARE) under K.A.CARE-King Abdulaziz University Collaboration Program.

Conflicts of Interest:
The author declares no conflict of interest.  Size of FCL installed in the route i-j