Multi-Stage Dynamic Transmission Network Expansion Planning Using LSHADE-SPACMA

: This paper introduces a multi-stage dynamic transmission network expansion planning (MSDTNEP) model considering the N-1 reliability constraint. The integrated planning problem of N-1 security and transmission expansion planning is essential because a single line outage could be a triggering event to rolling blackouts. Two suggested scenarios were developed to obtain the optimal conﬁguration of the Egyptian West Delta Network’s realistic transmission (WDN) to meet the demand of the potential load growth and ensure the system reliability up to the year 2040. The size of a blackout, based on the amount of expected energy not supplied, was calculated to evaluate both scenarios. The load forecasting (up to 2040) was obtained based on an adaptive neuro-fuzzy inference system because it gives excellent results compared to conventional methods. The linear population size reduction—Success-History-based Differential Evolution with semi-parameter adaptation (LSHADE-SPA) hybrid—covariance matrix adaptation evolution strategy (CMA-ES) algorithm (LSHADE-SPACMA)—is applied to solve the problem. The semi-adaptive nature of LSHADE-SPACMA and the hybridization between LSHADE and CMA-ES are able to solve complex optimization problems. The performance of LSHADE-SPACMA in solving the problem is compared to other well-established methods using three testing systems to validate its superiority. Then, the MSDTNEP of the Egyptian West Delta Network is presented, and the numerical results of the two scenarios are compared to obtain an economic plan and avoid a partial or total blackout.


Introduction
Transmission network expansion planning (TNEP) is a process to determine an optimal strategy for where, when, and how many transmission facilities are needed to extend the current power system transmission network. TNEP seeks to fulfill the demand of future load growth and extra generators while retaining the power system's reliability and safety performance [1][2][3].
TNEP usually considers the minimization of objective functions related to cost such as investment, operation, and reliability and subject to technical, financial, and service quality constraints. Technical constraints are synonymous with generator and branch this shortcoming and optimize the fitness function of the DPSO. Torres and Castro [22] implemented an improved local PSO (LPSO) algorithm to solve the STNEP problem.
Moreover, Huang and Dinavahi [23] presented the multi-group PSO (MGPSO) algorithm to solve the DSTNEP problem. The MGPSO is based on the DPSO framework with many improvements, including an initialization method, a multi-group co-evolution strategy, and a mutation mechanism. The integer-based particle swarm optimization (IBPSO) technique [24] and multiverse optimization (MVO) technique [25] were recently applied to solve TNEP; however, these techniques cannot guarantee the global optimum.
The present study employed a meta-heuristic algorithm called LSHADE-SPACMA to solve a MSDTNEP for the West Delta Network (WDN) up to year 2040 with an embedded N-1 security constraint. The semi-adaptive nature of LSHADE-SPACMA improved the exploration capability and exploitation tendency of the algorithm's ability to avoid local optima stagnation [26]. Moreover, the hybridization between LSHADE and CMA-ES is powerful in solving complex optimization problems [26].
The planning period (2016 to 2040) was divided into five stages, and two scenarios were suggested and solved to guarantee N-1 security and decrease the amount of energy not supplied based on the DCTNEP model. The load forecasting up to year 2040 was calculated based on adaptive neuro-fuzzy inference system (ANFIS), because it is efficient for long-term load forecasting based on a set of statistical tests [24]. The performance of LSHADE-SPACMA to solve the TNEP problem was compared to other well-established meta-heuristic methods such as IBPSO, MVO, and the HM. Three testing systems: the Garver 6-bus test system, Egyptian WDN, and 93-bus Colombian system were used to validate the proposed algorithm's capability to solve the TNEP problem. A comparison of two planning scenarios presented for WDN expansion was also conducted.
The following sections are organized as follows: Section 2 presents mathematical formulations of models adopted in this study; long-term load forecasting using ANFIS and LSHADE-SPACMA algorithms are introduced in Sections 3 and 4, respectively; the characteristics of the testing systems and results are discussed in Section 5; and conclusions are presented in Section 6.

Problem Formulation
The main objective of the MSDTNEP problem is the minimization of the total investment cost and expected energy not supplied subjected to a set of constraints over the entire planning horizon. The constraints are the power flow equations, line limit equations, and N-1 security.
Due to the TNEP problem's complexity, the DCPF model was used in the TNEP models applied in this study. The main assumptions in the traditional DCPF model are: network loss is negligible, line resistance (active power loss) is negligible or insignificant (i.e., R << X), and magnitude of bus voltage is set to 1.0 per unit (flat voltage profile), where R is the line resistance and x is the reactance. Based on these three assumptions, active power injection at bus i (P i ) is obtained as follows [4]: (1)

Deterministic Static TNEP (DSTNEP) Model without Security Constraint
The standard DSTNEP based on the lossless DC model can be written as follows: Which is subject to the following equality and inequality constraints: P Gi − P di = P ij , i, j = 1, . . . , N, i = j, . . . , N, where X ij is an integer. In the DSTNEP model, Equation (2) represents the objective function aimed to minimize total investment cost. Equation (3) represents the nodal balance, where the net power injection at bus i is equal to the difference between total generation and total loads connected to bus i. Equations (4) and (5) calculate the active power flow for existing lines (P ij ) as determined by the product of line susceptance (B ij ) and voltage phase angle difference θ i − θ j and should be less than or equal to the active power flow limit of the i-j right of way.

Multi-Stage Dynamic TNEP (MSDTNEP) Model
The formulation of the MSDTNEP model to minimize investment cost and amount of energy not supplied was computed using Equations (7)-(14) [4][5][6]. The objective function (Equation (7)) corresponds to the addition of the investment cost in each period adequately transferred to the initial period using an interest rate λ. To maintain the transmission system's reliability, North America Electrical Reliability Corporation (NERC) has published a series of standards in which all the balancing authorities within North American interconnection must comply [27]. N-1 security constraint is widely embedded in TNEP models to improve transmission capacity and security of power grids and is presented in Equation (14) [5,27].
which is subject to: . . , n, j = 1, 2, . . . , n, i = j (9) The N-1 security constraint was modeled using Equation (14): where the adequate values for C r in real networks can vary between 0.5 and 1 million USD/MW.

Applied Strategy
The proposed strategy to solve the TNEP problem ( Figure 1) starts with an initial load forecast for WDN up to year 2040 using ANFIS. Next, the LSHADE-SPACMA algorithm was implemented to solve MSDTNEP and DSTNEP problems as long as the current run is less than or equal to the maximum number of runs (max_run). If the current run exceeded the max_run, both the lower bound and stage number were updated. This process was repeated until the current stage exceeded the maximum number of stages (max_stage). The max_stage for MSDTNEP and DSTNEP were 5 and 1, respectively. process was repeated until the current stage exceeded the maximum number of stages (max_stage). The max_stage for MSDTNEP and DSTNEP were 5 and 1, respectively.

Load Forecasting Using ANFIS
The ANFIS was used to obtain long-term load forecasting of WDN up to year 2040. Artificial intelligence techniques such as the fuzzy logic controller (FLC) and the artificial neural network (ANN) are commonly used in long-term load forecasting [24][25][26][27][28][29].
Recognizing the non-linear relationship between the selected input and the given output is the FLC's key benefit in many applications in diverse fields. ANNs can "learn" because of the versatility of neuron-linking weights; thus, the solution can be adapted to boost efficiency. ANFIS combines ANN and FLC to enable this hybrid approach to adapt to various problems through machine learning using an iterative mechanism to detect non-linear relationships. ANFIS rapidly overcomes the complexities and robustness of the system by modeling its non-linear function with reasonable precision. The basic steps of ANFIS for long-term load forecasting are summarized as follows: Step 1: Define the input and output of the model. Input is the historical year and output is the actual peak load data.
Step 2: Collect all data from previous periods and normalize scales,

Load Forecasting Using ANFIS
The ANFIS was used to obtain long-term load forecasting of WDN up to year 2040. Artificial intelligence techniques such as the fuzzy logic controller (FLC) and the artificial neural network (ANN) are commonly used in long-term load forecasting [24][25][26][27][28][29].
Recognizing the non-linear relationship between the selected input and the given output is the FLC's key benefit in many applications in diverse fields. ANNs can "learn" because of the versatility of neuron-linking weights; thus, the solution can be adapted to boost efficiency. ANFIS combines ANN and FLC to enable this hybrid approach to adapt to various problems through machine learning using an iterative mechanism to detect non-linear relationships. ANFIS rapidly overcomes the complexities and robustness of the system by modeling its non-linear function with reasonable precision. The basic steps of ANFIS for long-term load forecasting are summarized as follows: Step 1: Define the input and output of the model. Input is the historical year and output is the actual peak load data.
Step 2: Collect all data from previous periods and normalize scales, Step 3: Divide data into two sets-training data and test data. Training data are about 70-90% of the data available.
Step 4: Run and estimate all reasonable ANFIS and determine the type of membership function and number of linguistic variables.
Step 5: Select the best ANFIS model.
Step 6: Project input variables using an autoregressive model in years (defined for the future).
Step 7: Predict total peak loads using selected ANFIS.

LSHADE with Semi-Parameter Adaptation Hybrid with CMA-ES
Differential Evolution (DE) is a stochastic population-based optimization algorithm, starting with randomly generated individuals that evolve through probabilistic operators such as recombination and mutation [30][31][32]. The efficiency of DE depends heavily on population size (Np), the chosen mutation and/or crossover strategy, the mutation-associated scale factor (F), and the recombination-associated crossover rate (Cr). Adaptive methods for online modification of control parameters during DE were explored to eliminate the need to tune the parameters. In the literature, novel mutation techniques were introduced.
In Mohamed et al. [26], a semi-parameter adaptation scheme was implemented to enhance LSHADE technique's performance. Besides, a hybridization strategy between LSHADE-SPA and CMA-ES was presented and considered successful on a wide range of optimization issues. In this section, the specifics of LSHADE-SPACMA are summarized, which is a recent improvement of LSHADE and is introduced in Mohamed et al. [26].

Initialization
To initialize LSHADE, a population of candidate solutions (decision vectors) is generated randomly within the specified upper and lower bounds. The mth component of the decision vector was formulated by the following Equation [26]: where rand (0, 1) returns a uniformly distributed random number in [0, 1] and superscript 'o' represents initialization. If 'd' is the dimension of decision vector, then k = 1, 2, . . . , Np and m = 1, 2, . . . , d.

Mutation
Different DE mutation strategies were considered to generate a mutant vector corresponding to each population member D G k . In Mohamed et al. [26], the 'current-to-pbest/1' strategy is used and given by: The p-value here is a control parameter that is supposed to enhance exploitation and exploration processes. Random indices (r 1 and r 2 ) were picked from the population's concatenation with an external archive that houses parent vectors that have generated successful vectors. The scale factor F is a positive control parameter used to scale the difference vector.

Crossover
For each generation G, the target vector is combined with mu G using Equation (17) to generate the trial vector u G k,m based on Cr [26]: k yields an as-good-as or better fitness function value than D G k , then u G k is set to D G+1 k . Otherwise, the old vector D G k is reserved. The selection scheme for a minimization problem follows [26]:

Linear Population Size Reduction (LPSR)
Linear Population Size Reduction (LPSR) is employed to boost the efficiency of LSHADE-SPA. The population size of the LPSR will be reduced by linear function as follows [22,26,27]: 4.1.6. Semi-Parameter Adaptation (SPA) of Scaling Factor (F) and Crossover Rate (Cr) Parameter adjustments had a significant effect on the efficiency of DE. Each problem has appropriate values for each parameter. SPA for F and Cr is proposed in [26] and is implemented in this work. It consists of two parts, in which the first part is activated during the first half of the search, and the second part is activated during the second half.
First Part of SPA SPA's idea is focused on adapting Cr using original LSHADE adaptation when the current number of iterations is less than half of the maximum number of iterations. While the F parameter is tuned using uniform distribution randomly within a specific limit as given in Equation (20) [26].
However, Cr k values are adapted based on the following Equation [26]: where Mcr k is a randomly chosen memory slot that stores successful means of preceding generations. The memory index k is randomly selected from the range [1, h]. h is the size of memory. All M cr initial values are 0.5, and one memory slot M cr is updated at the end of each generation using the arithmetic mean of Cr k values, that succeed to generate new individuals.

Second Part of SPA
In this part, the adaptation is concentrated on F, and LSHADE adaptation is used to adjust the F k parameter as follows [26]: where σ is the standard deviation for Cauchy distribution and equals 0.1, M F is a memory slot randomly chosen to store successful means of preceding generations. One memory slot M F is modified using the Lehmer mean of F k values at the end of each generation to successfully generate new individuals. The F k value of the previous five generations of the first part of the SPA was used to initialize the memory slot M F for the second part of SPA. The Cr parameter adaptation process remains as is during the second part. However, due to the nature of the LSHADE parameter adaptation, the Cr parameter is gradually frozen to the adapted values. According to the SHADE parameter adaptation, when successful individuals do not produce all Cr values in a generation, the corresponding memory slot is adjusted to the terminal value. Therefore, M Cr will not be updated until the search is over.

Covariance Matrix Adaptation Evolution Strategy (CMA-ES)
Among its many variants, CMA-ES can solve various types of optimization problems efficiently. In CMA-ES, a multivariate normal distribution is used to model the search space. New individuals are produced using Gaussian distribution that accounts for the population's path over successive generations. CMA-ES automatically adapts the mean vector m, covariance matrix C, and step size σ c . CMA-ES steps are summarized according to [26]: Generate an initial population and then calculate the fitness function.

2.
Gaussian distribution is used to produce new individuals, thus: 3.
Update m using the best µ individuals according to Update σ c and C.

5.
Steps 2 and 3 are repeated until a stop criterion is met.

LSHADE-SPACMA Hybridization Framework
In order to enhance LSHADE-SPA's performance, a hybridization framework between the recent version of CMA-ES and LSHADE-SPA is applied [26]. Each individual D in the population produces an individual u offspring using either LSHADE or CMA-ES according to the probability variable (FCP) class. FCP values are chosen randomly from memory slots that are set to M FCP . One memory slot is modified at the end of each generation based on the output of each algorithm. As a result, populations were gradually assigned to a better performance algorithm. Updates were produced using individuals that successfully created new individuals. The memory slot of the M FCP is then updated based on the following formula: where C is the learning rate, and ∆ AlG1 is the improvement rate for each algorithm calculated according to Equation (25): where prob min and prob max are the minimum and maximum probabilities set to each algorithm, respectively. Thus, to execute both algorithms together, FCP values must be maintained in the range of 0.2 to 0.8. The variable ω AlG1 is the summation of differences between the old and new fitness values for each individual belonging to the algorithm in Equation (26): where f is the fitness function. u old and u are the old and the offspring individuals, respectively. n is the number of individuals belonging to algorithm AlG1.

Pseudocode of LSHADE-SPACMA for Solving the TNEP Problem
The pseudocode of the LSHADE-SPACMA algorithm for solving the TNEP problem is presented in Table 1. It starts with setting the population size and dimension of the problem and defining the maximum number of problem evaluations, the lower bound, and Appl. Sci. 2021, 11, 2155 9 of 20 the upper bound of decision variables. The population of Np individuals was initialized according to Equation (15), and initial values of SHADE and CMA control parameters were introduced. Steps 1 through 10 are repeated until a stop criterion is reached.  • Initialize CMA parameters.
• Evaluate objective function and constraint violation using Equation (2) to Equation (6) or using Equation (7) to Equation (14) for each individual in POP. • Increase function evaluation counter NFE by N p .
Step 5 Update POP and Fitness function according to the evaluation of U.
Step 6 Store successful FCP, F, and Cr Step 7 Delete individuals archive randomly End Step 8 Update M Cr and M FCP If nfes < max_nfes/2 Update M F end Step 9 • Calculate N G+1 according to Equation (19).

End
Step 10 • Update CMA parameters.

Simulation Results and Analyses
An assessment of three standard test systems-Garver 6-bus test system, Egyptian WDN, and 93-bus Colombian system-were carried out to verify the proposed optimization technique's capability to solve the TNEP problem. Static plan results, using LSHADE-SPACMA, were compared to results obtained in previous studies [4,13,20,21,28]. Finally, MSDTNEP of WDN up to 2040 with embedded N-1 security constraints is employed based on two scenarios. The following case studies are conducted on MATLAB r2017a on a DELL PC, and its model name is 'OptiPlex7050', including an 'Intel ® Core™ i7' CPU at 2.6 GHz and 16 GB RAM. In the pioneering paper [4], the Garver network was discussed and has since been used by many other researchers to compare various TNEP approaches. Figure 2 shows that it comprises 6 buses and 6 lines. Bus 6 is not connected to the rest of the network in the initial configuration, and the current demand for buses 1-5 is 760 MW, while the installed generation capacity for these interconnected nodes is only 215 MW. As a result, to avoid power from not being supplied, the TNEP would have to facilitate interconnecting bus 6 to the rest of the network.


If NG<NG+1 -Sort individuals in P based on their fitness values and delete the NG-NG+1 member.

-
Resize archive (| |) according to the new POP.

End
Step 10  Update CMA parameters.

Simulation Results and Analyses
An assessment of three standard test systems-Garver 6-bus test system, Egyptian WDN, and 93-bus Colombian system-were carried out to verify the proposed optimization technique's capability to solve the TNEP problem. Static plan results, using LSHADE-SPACMA, were compared to results obtained in previous studies [4,13,20,21,28]. Finally, MSDTNEP of WDN up to 2040 with embedded N-1 security constraints is employed based on two scenarios. The following case studies are conducted on MATLAB r2017a on a DELL PC, and its model name is 'OptiPlex7050', including an 'Intel ® Core™ i7' CPU at 2.6 GHz and 16 GB RAM.

Garver 6-Bus Test System
In the pioneering paper [4], the Garver network was discussed and has since been used by many other researchers to compare various TNEP approaches. Figure 2 shows that it comprises 6 buses and 6 lines. Bus 6 is not connected to the rest of the network in the initial configuration, and the current demand for buses 1-5 is 760 MW, while the installed generation capacity for these interconnected nodes is only 215 MW. As a result, to avoid power from not being supplied, the TNEP would have to facilitate interconnecting bus 6 to the rest of the network. Regarding the standard Garver system, it is important to stress that the optimal identified solution is the same as the solution mentioned in the literature to validate the proposed algorithm's efficiency in solving the TNEP problem. The final configuration Regarding the standard Garver system, it is important to stress that the optimal identified solution is the same as the solution mentioned in the literature to validate the proposed algorithm's efficiency in solving the TNEP problem. The final configuration and the optimal solution of TNEP for the Garver 6-bus network system using the LSHADE-SPACMA technique are depicted in Table 2 and Figure 2b. This solution includes one new branch between nodes 3 and 5, two new branches between nodes 4 and 6, and four new branches between nodes 2 and 6. The proposed method gives 13 paths with costs of 200 million USD. It provides the same results obtained by Garver [4] using the linear programming (LP) method in 6.68 s, decreasing the computational burden.

93-Bus Colombian System
The Colombian test system had 93 buses and 155 circuits. It was designed to meet a 9750 MW load system, and 155 control variables were optimized to minimize the total investment cost. The initial configuration for the system can be seen in Escobar et al. [17].
The best topology for the Colombian system using LSHADE-SPACMA is depicted in Table 3, which produced the same results obtained by Escobar et al. [13]. Both methods added 6 circuits at an investment cost of 316.44 million USD, while LSHADE-SPACMA got the global optimum solution in 475.19 s. The results demonstrate that LSHADE-SPACMA is sufficient to solve large-scale systems and complex problems in one run.

West Delta Network (WDN) System Planning
WDN is a 66 kV transmission network and a section of the Unified Egyptian Network. The initial configuration of WDN with proposed routes is shown in Figure 3

West Delta Network (WDN) System Planning
WDN is a 66 kV transmission network and a section of the Unified Egyptian Network. The initial configuration of WDN with proposed routes is shown in Figure 3

Validation of LSHADE-SPACMA Technique to Solve the TNEP Problem for WDN
WDN was expected to meet the predicted peak load of the 2195.8 MW. A new site generator at bus 53 and 31 candidate routes have been proposed [24,25,33]. The new configuration of WDN is shown in Table 4. LSHADE-SPACMA adds seven new right of way routes and yields 17.28 million USD for an investment cost of the added transmission lines, and the consumed time of one run is 164.78 s. The results demonstrate, moreover, that LSHADE-SPACMA has the superiority compared to the HM [33], IBPSO [24], and MVO [25], because they acquired a high investment cost of 21.24, 22.19, and 20.645 million USD, respectively. The results presented in Table 5 explore that the power flow in the transmission lines is within its prescribed limits.

Validation of LSHADE-SPACMA Technique to Solve the TNEP Problem for WDN
WDN was expected to meet the predicted peak load of the 2195.8 MW. A new site generator at bus 53 and 31 candidate routes have been proposed [24,25,33]. The new configuration of WDN is shown in Table 4. LSHADE-SPACMA adds seven new right of way routes and yields 17.28 million USD for an investment cost of the added transmission lines, and the consumed time of one run is 164.78 s. The results demonstrate, moreover, that LSHADE-SPACMA has the superiority compared to the HM [33], IBPSO [24], and MVO [25], because they acquired a high investment cost of 21.24, 22.19, and 20.645 million USD, respectively. The results presented in Table 5 explore that the power flow in the transmission lines is within its prescribed limits.  In this section, MSDTNEP of WDN without N-1 security constraint is employed up to year 2040. The actual peak load data of WDN from the year 2008 to 2015 are used in the load forecasting process [24]. The load forecasting up to year 2040 using ANFIS, linear, parabolic and exponential trends is presented in Figure 4. The results show that ANFIS has the lowest mean absolute error (MAE) and the best results; therefore, it is applied in this work. The predicted loads from 2016 to 2040 using ANFIS are shown in Table 6. In this section, MSDTNEP of WDN without N-1 security constraint is employed up to year 2040. The actual peak load data of WDN from the year 2008 to 2015 are used in the load forecasting process [24]. The load forecasting up to year 2040 using ANFIS, linear, parabolic and exponential trends is presented in Figure 4. The results show that ANFIS has the lowest mean absolute error (MAE) and the best results; therefore, it is applied in this work. The predicted loads from 2016 to 2040 using ANFIS are shown in Table 6.

Year
Predicted Load Year Predicted Load Figure 4. WDN load forecasting using ANFIS, linear, parabola and exponential trends.

Year Predicted Load Year Predicted Load
MSDTNEP without including an N-1 security constraint for WDN from 2016 up to 2040 is shown in Table 7. The planning period is divided into five stages and each stage comprises five years. The considered interest rate is 0.2. The dynamic planning results in adding three circuits in the first and second stages, two circuits in the third stage, and four circuits in the fourth and fifth stages to guarantee that the power system operates economically and reliably. The results also show that the generation station's capacity at bus 53 should be increased to 1315.7 MW to meet the possible load demand in 2040. Dynamic planning of WDN in 2030 (third stage) shows that total circuits added from 2016 to 2030 are eight; however, the static plan in Section 5.2.1 gives seven circuits. As in dynamic planning, each stage depends on previous stages. In the second and third stages, six new routes are added to ensure reliable operation for the power system and meet each bus's load demand in these periods. Hence, these routes are imposed in the planning of the third stage. To study the impact of including N-1 security constraint on WDN planning, two scenarios were suggested to keep the power system reliable, secure, and low cost: Scenario 1: The system was planned with the same routes proposed in Fathy et al. [24] and load shedding was applied.
Scenario 2: The maximum number of circuits in each route was supposed to be four, as per Egyptian requirements. Table 8 and Figure 5 show that the planning of WDN with routes proposed by Fathy et al. [24] for scenario 1 cannot ensure system security without load shedding. Scenario 1 results in the addition of 15 circuits with load shedding of 162.68 MW in the first stage, 6 circuits with load shedding of 43.23 MW in the second stage, 1 circuit with load shedding of 127.4 MW in the third stage, 4 circuits with load shedding of 134.19 MW in the fourth stage, and 5 circuits with load shedding of 160.46 MW in the fifth stage. The total investment cost would be 287.46 million USD.  0  0  2  2  2  5-7  1  1  1  1  1  41-39  0  2  2  2  2  5-8  1  1  1  1  1  39-38  1  2  2  2  2  5-22  1  1  1  1  1  8-41  2  2  2  2  2  5-29  1  1  1  1  1  8-51  1  1  2  2   Scenario 2 proposes 25 circuits in the first stage, 8 circuits in the second stage, 20 circuits in the third stage, 7 circuits in the fourth stage, and 12 circuits in the fifth stage, as depicted in Table 8 and Figure 6. The total investment cost is 57.33 million USD, lower than scenario 1. Hence, scenario 2 is recommended.
It is clear that steady-state power system insecurity based on scenario 1 overloads transmission lines and causes successive interruptions leading to partial blackouts ( Figure 7 and Table 8). The results demonstrate that increasing the capacity at bus 53 cannot mitigate the rolling blackout problem, and the available solution is to install distribution generators at various locations, causing the total investment cost to increase. On the other hand, MS-DTNEP of WDN based on scenario 2 leads to a better plan for WDN and lower investment cost, because it accounts for several features such as annual load growth, prevention of rolling blackouts, and N-1 security criterion. Blackouts are avoided in scenario 2 by installing additional circuits of lower cost compared to the cost of installing new distribution generators (Figure 7).   Table 8 and Figure 6. The total investment cost is 57.33 million USD, lower than scenario 1. Hence, scenario 2 is recommended. It is clear that steady-state power system insecurity based on scenario 1 overloads transmission lines and causes successive interruptions leading to partial blackouts (Figure 7 and Table 8). The results demonstrate that increasing the capacity at bus 53 cannot mitigate the rolling blackout problem, and the available solution is to install distribution generators at various locations, causing the total investment cost to increase. On the other hand, MSDTNEP of WDN based on scenario 2 leads to a better plan for WDN and lower investment cost, because it accounts for several features such as annual load growth, prevention of rolling blackouts, and N-1 security criterion. Blackouts are avoided in scenario 2 by installing additional circuits of lower cost compared to the cost of installing new distribution generators (Figure 7).

Single Circuit
Two circuits Four circuits

Conclusions and Future Works
The implementation of the LSHADE-SPACMA algorithm was investigated to solve the TNEP problem. Promising results were obtained using three sample systems and confirmed the potential of this approach. The results show that the LSHADE-SPACMA algorithm produces a better (low-cost) solution than IBPSO, MVO, and HT for the WDN. Moreover, the LSHADE-SPACMA reached optimum solutions for both large-and small-scale systems such as the Colombian system and the Garver test system in a single run.
Ignoring N-1 reliability constraint in the TNEP model threatens the security of power systems. A suggested approach to tackle this problem was applied to the Egyptian WDN, and its performance was evaluated. The MSDTNEP up to year 2040 including an N-1 security constraint, and it was assessed based on two scenarios. Scenario 1 adopted plans based on the candidate routes with the implementation of load shedding. Scenario 2, the suggested solution, assumed that the maximum number of circuits in each line could be increased to 4.
ANFIS was employed to predict load demand growth until 2040. The input and the output of ANFIS are the historical year and the actual peak load data, respectively. The numerical results show that dynamic planning leads to a more accurate and realistic assessment of the network. Moreover, the N-1 security constraint has a significant impact on the final configuration of WDN. The results also demonstrate that the planning of WDN with the candidate lines proposed in scenario 1 cannot guarantee system security without imposing a load shedding or partial blackouts. It should be noted that scenario 1 results in the addition of 31 circuits with load shedding of 627.96 MW at the end of the planning horizon. The total investment cost is estimated to be 287.46 million USD. Scenario 2 does not cause blackouts and is recommended. This scenario requires the addition of 62 circuits with no load shedding and an investment cost of 57.33 million USD.
The results demonstrated that the suggested strategy is considered a promising solution to maintain N-1 security and avoid rolling blackouts. It is applicable for all power systems because it is not complex and has a low investment cost [34].
Egypt's vision is to achieve an energy sector that meets national sustainable development standards and optimizes the use of renewable resources to support economic growth and protect the environment's health and quality [35,36]. Hence, a new paradigm of transmission and generation planning for WDN is needed, which considers

Conclusions and Future Works
The implementation of the LSHADE-SPACMA algorithm was investigated to solve the TNEP problem. Promising results were obtained using three sample systems and confirmed the potential of this approach. The results show that the LSHADE-SPACMA algorithm produces a better (low-cost) solution than IBPSO, MVO, and HT for the WDN. Moreover, the LSHADE-SPACMA reached optimum solutions for both large-and small-scale systems such as the Colombian system and the Garver test system in a single run.
Ignoring N-1 reliability constraint in the TNEP model threatens the security of power systems. A suggested approach to tackle this problem was applied to the Egyptian WDN, and its performance was evaluated. The MSDTNEP up to year 2040 including an N-1 security constraint, and it was assessed based on two scenarios. Scenario 1 adopted plans based on the candidate routes with the implementation of load shedding. Scenario 2, the suggested solution, assumed that the maximum number of circuits in each line could be increased to 4.
ANFIS was employed to predict load demand growth until 2040. The input and the output of ANFIS are the historical year and the actual peak load data, respectively. The numerical results show that dynamic planning leads to a more accurate and realistic assessment of the network. Moreover, the N-1 security constraint has a significant impact on the final configuration of WDN. The results also demonstrate that the planning of WDN with the candidate lines proposed in scenario 1 cannot guarantee system security without imposing a load shedding or partial blackouts. It should be noted that scenario 1 results in the addition of 31 circuits with load shedding of 627.96 MW at the end of the planning horizon. The total investment cost is estimated to be 287.46 million USD. Scenario 2 does not cause blackouts and is recommended. This scenario requires the addition of 62 circuits with no load shedding and an investment cost of 57.33 million USD.
The results demonstrated that the suggested strategy is considered a promising solution to maintain N-1 security and avoid rolling blackouts. It is applicable for all power systems because it is not complex and has a low investment cost [34].
Egypt's vision is to achieve an energy sector that meets national sustainable development standards and optimizes the use of renewable resources to support economic growth and protect the environment's health and quality [35,36]. Hence, a new paradigm of transmission and generation planning for WDN is needed, which considers the system's increasingly decentralized and stochastic nature due to the predicted high penetration of renewable energy sources in the Egyptian network.

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.

Conflicts of Interest:
The authors declare no conflict of interest. Active power injection at bus i Bij Susceptance of route between bus i and j θ i , θ j Voltage angles at bus i and j, respectively P Gi , P di Active power generation source and the load demand (MW) at bus i, respectively C ij Cost of circuit between buses i and j X Line reactance X 0 ij , X ij_max , X ij Initial number of circuits, the maximum number of circuits, and the actual number of circuits between buses i and j, respectively P ij , P ij_max Active power flow and the active power flow limit in the i-j right of way (MW), respectively P sh,i Amount of load shedding C r Penalty parameter that penalizes, in the objective function, any system load shedding