# Transmission Network Expansion Planning Considering Phase-Shifter Transformers

Authors: Celso T. Miasaki; Edgar M. C. Franco; Ruben A. Romero
Journal: Journal of Electrical and Computer Engineering
 (2012)

---

## Abstract

This paper presents a novel mathematical model for the transmission network expansion planning problem. Main idea is to consider phase-shifter (PS) transformers as a new element of the transmission system expansion together with other traditional components such as transmission lines and conventional transformers. In this way, PS are added in order to redistribute active power flows in the system and, consequently, to diminish the total investment costs due to new transmission lines. Proposed mathematical model presents the structure of a mixed-integer nonlinear programming (MINLP) problem and is based on the standard DC model. In this paper, there is also applied a specialized genetic algorithm aimed at optimizing the allocation of candidate components in the network. Results obtained from computational simulations carried out with IEEE-24 bus system show an outstanding performance of the proposed methodology and model, indicating the technical viability of using these nonconventional devices during the planning process.


---

## Body

## 1. Introduction

Transmission network planning begins with the establishment of power demand growth scenarios, in accordance with forecasts along the time. Given these scenarios, one can verify the eventual need to broaden and to strengthen the network. In case electric service conditions are not satisfied, there should be proposed a plan that has coherence among the power supply availability, demand, and installation of new equipments in the network. Integration of these new equipments in the network, aimed at maintaining suitable technical and operating conditions, requires planning of the allocation of such reinforcement.Main objective of the transmission expansion planning is to obtain theoptimal expansion plan, while fulfilling operating and economic constraints.Formulation of a mathematical representation for the transmission expansion planning problem begins with some assumptions, where accuracy and complexity are considered in the model construction. Regularly, the problem is represented by a Mixed-Integer Nonlinear Programming (MINLP) problem that presents many local optima solutions for real-life systems. This high number is due to the possible expansion plans that shows the association of the specified optimal operational mode. Therefore, a basic problem consists in defining the least-cost expansion alternative that satisfies all operating constraints.In static long-term transmission expansion planning (typically with a planning horizon of more than 5 years), all investments are carried out in a single-year planning horizon, whereas for the multistage it is divided into several stages.Static planning is aimed at searchingwhere and whattype of circuits should be constructed for the network to operate properly. This type of planning applies some simplifications; for instance, voltage magnitudes constraints are neglected, albeit in short-term planning are essential. Voltage and stability problems are not taken into account in this approach either. Even though considering these simplifications, the problem still remains complex with unknown optima solutions for various real-life systems.In the technical literature, DC and transportation models are static mathematical models often used to solve the transmission expansion planning problem. These models consider only the addition of transmission lines (TLs) and conventional transformers.Here a novel transmission expansion strategy is proposed. An improved model considers the inclusion of a new kind of device, in this case, a flexible alternating current transmission system (FACTS) device.The literature concerning the use of FACTS devices is wide. However, most papers treat only the operational improvement by using FACTS devices [1–3]. This topic is considered as the operation planning, in which FACTS devices are used, for example, to redistribute the active power flow in order to eliminate congestion problems and/or to optimize the power dispatch. In this way, the PS is considered a FACTS device with the characteristic of redistributing the active power flow.Feasibility of employing the PS as a candidate component in the long-term transmission expansion planning process is analyzed regarding a static and centralized planning model. Nevertheless, the present model can be extended to the multistage planning [4], to the competitive planning [5], or to the planning with security constraints [6].Instead of trying to consider all functions concerning the planning problem in a single model, this work is focused solely on the core of the network synthesis; for which the mathematical modeling and the solution technique is addressed.The proposed model is based on the DC model, which is the most employed one in planning problems; consequently, only the active power flow is considered. Other aspects such as performance analysis (reliability and stability analysis, reactive planning, AC power flow, and short-circuit calculation) relevant to transmission expansion are beyond the scope of this paper. Nonetheless, in general, after obtaining a basic solution, all those analyses can be carried out.


## 2. Mathematical Model

This section introduces the classical mathematical model and the proposed one.


### 2.1. Classical Model

The mathematical formulation of the DC model for transmission network expansion planning problem, when considering solely the installation of transmission lines and/or conventional transformers, assumes the following form:(1)Minimizev=∑(i,j)cijnijsubject to(2)S⋅f+g=d,(3)fij-γij⋅(nij0+nij)⋅(θi-θj)=0,(4)|fij|≤(nij0+nij)⋅f¯ij,(5)0≤gi≤gi¯,(6)0≤nij≤nij¯,(7)nijInteger,(8)fijUnbounded,(9)θjUnbounded,(10)(i,j)∈Ω.The objective function (1) represents the investment cost due to new transmission lines during the specified planning horizon.The first set of constraints (2) represents Kirchhoff’s Current Law (KCL) equations (one constraint per bus), and the second one (3) represents Kirchhoff’s Voltage Law (KVL) (one constraint per branch). In this model, the transmission lines or regular transformers are represented without distinction by decision variables nij.Set of constraints (4) refers to the capacity of transmission circuits (lines and/or transformers) in terms of the required absolute value since the power flow can flow in both ways.Constraint (5) represents the limits for generation buses and (6) represents the limits for the added circuits at each candidate path i-j. Remaining constraints refer to the self-characteristics of the variables.From the operational research standpoint, system (1)–(10) is a MINLP whose resolution is complex, especially for large-size electric systems. Main source of complexity in the problem is due to variables nij that need to be integer.For some types of algorithms utilized for the transmission network planning problem, it is more suitable to carry out alterations to the basic modeling for allowing the application of the solution techniques. An alteration commonly used is the insertion of new variables that represent the load shedding associated with all load buses of the system. This resource can also be seen as an artificial generation aimed at turning the problem always viable during computational implementations [7]. Note that this artificial generation, from the mathematical standpoint, is solely an increase in the number of variables of the original problem. Additionally, all variables corresponding to the artificial generators in the final solution should be equal to zero, so both problems (original and modified) will have equivalent solutions.


### 2.2. Model with Phase Shifters

When PS is considered in the transmission network planning problem, the DC model assumes the following form:(11)Minimizev=∑(i,j)cijnij+∑dij⋅(nij0+nij)⋅δij+α∑rksubject to(12)S⋅f+g+r=d,(13)fij-γij⋅(nij0+nij)⋅(θi-θj+φij⋅δij)=0,(14)|fij|≤(nij0+nij)⋅f¯ij,(15)0≤gi≤gi¯,(16)0≤nij≤nij¯,(17)0≤rkj≤dk,(18)nijInteger,(19)fijUnbounded,(20)θjUnbounded,(21)φijUnbounded,(22)δijUnbounded,(23)(i,j)∈Ω.The objective function (11) corresponds to the investment cost. First two terms refer to the costs due to the insertion of transmission lines and PS, respectively. Last term refers to a penalty for the configuration under analysis, due to the presence of artificial generation. The penalty parameter α should be properly weighted in order to turn those configurations less attractive with artificial generation.In constraint (12), corresponding to KCL, vectors g and r indicate existent and artificial generations, respectively.Installation of PS in one or more lines can be represented as the combination of the buses’ angles and the angle supplied by the equipment.Angle between terminal voltages of a transmission line can be modified by installing a PS. Therefore, power flow equations (KVL) are affected when these devices are inserted. Thus, the function of PS appears in the KVL, which redirects the active power flow.The PS is considered a component with negligible reactance that can be placed in series with a transmission line or a conventional transformer.In this work, the angular difference of the PS was considered an unbounded variable; however, a limit can be set without modifying the mathematical model significantly.Presence of variableδij∈{0,1} in (13) makes this relation even more nonlinear in relation to the DC model and, consequently, more complex since four types of variables (nij, θj, φij, and δij) are multiplied.Set of constraints (14), (15), (16), and (17) refers to the transmission capacity of the circuits, limits of generation buses, limits of circuits added in each candidate path i-j, and limits of artificial generation buses, respectively. Remaining constraints refer to the characteristics of variables.Therefore, the proposed mathematical model is more complex than the classical one, due to the characteristics of a PS.In traditional models, a transmission line or a conventional transformer in a pathk-m can be represented in Figure 1.Figure 1
Transmission line.The PS is considered with zero impedance connected in series with a transmission line or a conventional transformer in pathk-m as shown in Figure 2.Figure 2
Transmission line and phase-shifter transformer.It is also considered that when a pathk-m is selected, the PS is allocated to every transmission line present in path k-m. Therefore, the number of PS at each selected path is equal to nPSij=(nij0+nij)·δij. The PS can be allocated in existing transmission lines as well as in those created during the optimization process.Problem (11)–(23) represents a MINLP, which is more complex than the model (1)–(10). However, metaheuristic techniques employed for (1)–(10) can also be applied to (11)–(23) after performing some modifications.


### 2.3. Illustrative Example: 3-Bus System

An example consisting of a network with three buses is presented in order to illustrate the application of the PS. System data is shown in Tables1 and 2.Table 1
Lines data of the 3 buses system.
LineReactance (pu)Maximum flow (MW)1-20.333351–30.500402-30.50040Table 2
Generation and demand data of the 3-buse system.
BusGeneration (MW)Demand (MW)170020603010Results obtained by performing the linearized DC power flow are illustrated in Figure3.Figure 3
Three-bus system.Notice that there is a load shedding ofr2=3.75MW at bus 2. Therefore the system is not operating adequately and new transmission lines are required. On the other hand, lines 1–3 and 2-3 are operating below their capacity. In this case, there is no possibility of increasing the transmitted power through the lines without overloading line 1-2. An interesting solution for this problem is to install a PS at 1–3 (see Figure 4), so the system operates without load shedding.Figure 4
Load flow solution with a phase shifter.Thus, this example shows how a PS is able to redirect the active power flow. This property will be employed in the long-term transmission planning, in which basic components are transmission lines, conventional transformers, and PS.


### 2.4. Phase-Shifter Transformer Features

PSs have the ability to redirect active power flows in the network. This feature provides a dynamic operational mode since it makes increasing the utilization of existing circuits possible. Consequently, as it can be verified from relation (13), the PS acts directly on KVL.Another important aspect is the use of relaxed models. In general, optimal solutions for relaxed models are not feasible for more accurate or constrained ones. Thus, it is probable that the optimal solution obtained by the transportation model, where the KVL (13) is relaxed (dropped), presents the KVL constraints violated in the DC model. With the inclusion of PS in determined positions, a feasible solution obtained with transportation model becomes also feasible for the DC model since the KCL constraints are already satisfied for both models. However, an efficient inclusion should be carried out in order to satisfy the optimality conditions of the DC model, once there exist expansion proposals with lower costs. It is worthwhile to notice that if PS costs are equal to zero, then the optimal solution for the transportation model and the DC model presents the same added lines, that is, presents the same topology for both mathematical models. Only difference with the conventional approach is that some PS are included in order to satisfy the KVL.The primary objective of this work is to verify the operation of PS and the technical feasibility of considering such type of equipment in long-term transmission expansion planning. In case a reduced-cost PS is employed, the optimal solution will be the same of the transportation model with addition of PS. On the other hand, higher costs will inhibit the presence of PS in the optimal solution, tending to the solution given by classical DC model. Finally, if the costs are competitive to transmission lines, an intermediary solution will be provided.


## 2.1. Classical Model

The mathematical formulation of the DC model for transmission network expansion planning problem, when considering solely the installation of transmission lines and/or conventional transformers, assumes the following form:(1)Minimizev=∑(i,j)cijnijsubject to(2)S⋅f+g=d,(3)fij-γij⋅(nij0+nij)⋅(θi-θj)=0,(4)|fij|≤(nij0+nij)⋅f¯ij,(5)0≤gi≤gi¯,(6)0≤nij≤nij¯,(7)nijInteger,(8)fijUnbounded,(9)θjUnbounded,(10)(i,j)∈Ω.The objective function (1) represents the investment cost due to new transmission lines during the specified planning horizon.The first set of constraints (2) represents Kirchhoff’s Current Law (KCL) equations (one constraint per bus), and the second one (3) represents Kirchhoff’s Voltage Law (KVL) (one constraint per branch). In this model, the transmission lines or regular transformers are represented without distinction by decision variables nij.Set of constraints (4) refers to the capacity of transmission circuits (lines and/or transformers) in terms of the required absolute value since the power flow can flow in both ways.Constraint (5) represents the limits for generation buses and (6) represents the limits for the added circuits at each candidate path i-j. Remaining constraints refer to the self-characteristics of the variables.From the operational research standpoint, system (1)–(10) is a MINLP whose resolution is complex, especially for large-size electric systems. Main source of complexity in the problem is due to variables nij that need to be integer.For some types of algorithms utilized for the transmission network planning problem, it is more suitable to carry out alterations to the basic modeling for allowing the application of the solution techniques. An alteration commonly used is the insertion of new variables that represent the load shedding associated with all load buses of the system. This resource can also be seen as an artificial generation aimed at turning the problem always viable during computational implementations [7]. Note that this artificial generation, from the mathematical standpoint, is solely an increase in the number of variables of the original problem. Additionally, all variables corresponding to the artificial generators in the final solution should be equal to zero, so both problems (original and modified) will have equivalent solutions.


## 2.2. Model with Phase Shifters

When PS is considered in the transmission network planning problem, the DC model assumes the following form:(11)Minimizev=∑(i,j)cijnij+∑dij⋅(nij0+nij)⋅δij+α∑rksubject to(12)S⋅f+g+r=d,(13)fij-γij⋅(nij0+nij)⋅(θi-θj+φij⋅δij)=0,(14)|fij|≤(nij0+nij)⋅f¯ij,(15)0≤gi≤gi¯,(16)0≤nij≤nij¯,(17)0≤rkj≤dk,(18)nijInteger,(19)fijUnbounded,(20)θjUnbounded,(21)φijUnbounded,(22)δijUnbounded,(23)(i,j)∈Ω.The objective function (11) corresponds to the investment cost. First two terms refer to the costs due to the insertion of transmission lines and PS, respectively. Last term refers to a penalty for the configuration under analysis, due to the presence of artificial generation. The penalty parameter α should be properly weighted in order to turn those configurations less attractive with artificial generation.In constraint (12), corresponding to KCL, vectors g and r indicate existent and artificial generations, respectively.Installation of PS in one or more lines can be represented as the combination of the buses’ angles and the angle supplied by the equipment.Angle between terminal voltages of a transmission line can be modified by installing a PS. Therefore, power flow equations (KVL) are affected when these devices are inserted. Thus, the function of PS appears in the KVL, which redirects the active power flow.The PS is considered a component with negligible reactance that can be placed in series with a transmission line or a conventional transformer.In this work, the angular difference of the PS was considered an unbounded variable; however, a limit can be set without modifying the mathematical model significantly.Presence of variableδij∈{0,1} in (13) makes this relation even more nonlinear in relation to the DC model and, consequently, more complex since four types of variables (nij, θj, φij, and δij) are multiplied.Set of constraints (14), (15), (16), and (17) refers to the transmission capacity of the circuits, limits of generation buses, limits of circuits added in each candidate path i-j, and limits of artificial generation buses, respectively. Remaining constraints refer to the characteristics of variables.Therefore, the proposed mathematical model is more complex than the classical one, due to the characteristics of a PS.In traditional models, a transmission line or a conventional transformer in a pathk-m can be represented in Figure 1.Figure 1
Transmission line.The PS is considered with zero impedance connected in series with a transmission line or a conventional transformer in pathk-m as shown in Figure 2.Figure 2
Transmission line and phase-shifter transformer.It is also considered that when a pathk-m is selected, the PS is allocated to every transmission line present in path k-m. Therefore, the number of PS at each selected path is equal to nPSij=(nij0+nij)·δij. The PS can be allocated in existing transmission lines as well as in those created during the optimization process.Problem (11)–(23) represents a MINLP, which is more complex than the model (1)–(10). However, metaheuristic techniques employed for (1)–(10) can also be applied to (11)–(23) after performing some modifications.


## 2.3. Illustrative Example: 3-Bus System

An example consisting of a network with three buses is presented in order to illustrate the application of the PS. System data is shown in Tables1 and 2.Table 1
Lines data of the 3 buses system.
LineReactance (pu)Maximum flow (MW)1-20.333351–30.500402-30.50040Table 2
Generation and demand data of the 3-buse system.
BusGeneration (MW)Demand (MW)170020603010Results obtained by performing the linearized DC power flow are illustrated in Figure3.Figure 3
Three-bus system.Notice that there is a load shedding ofr2=3.75MW at bus 2. Therefore the system is not operating adequately and new transmission lines are required. On the other hand, lines 1–3 and 2-3 are operating below their capacity. In this case, there is no possibility of increasing the transmitted power through the lines without overloading line 1-2. An interesting solution for this problem is to install a PS at 1–3 (see Figure 4), so the system operates without load shedding.Figure 4
Load flow solution with a phase shifter.Thus, this example shows how a PS is able to redirect the active power flow. This property will be employed in the long-term transmission planning, in which basic components are transmission lines, conventional transformers, and PS.


## 2.4. Phase-Shifter Transformer Features

PSs have the ability to redirect active power flows in the network. This feature provides a dynamic operational mode since it makes increasing the utilization of existing circuits possible. Consequently, as it can be verified from relation (13), the PS acts directly on KVL.Another important aspect is the use of relaxed models. In general, optimal solutions for relaxed models are not feasible for more accurate or constrained ones. Thus, it is probable that the optimal solution obtained by the transportation model, where the KVL (13) is relaxed (dropped), presents the KVL constraints violated in the DC model. With the inclusion of PS in determined positions, a feasible solution obtained with transportation model becomes also feasible for the DC model since the KCL constraints are already satisfied for both models. However, an efficient inclusion should be carried out in order to satisfy the optimality conditions of the DC model, once there exist expansion proposals with lower costs. It is worthwhile to notice that if PS costs are equal to zero, then the optimal solution for the transportation model and the DC model presents the same added lines, that is, presents the same topology for both mathematical models. Only difference with the conventional approach is that some PS are included in order to satisfy the KVL.The primary objective of this work is to verify the operation of PS and the technical feasibility of considering such type of equipment in long-term transmission expansion planning. In case a reduced-cost PS is employed, the optimal solution will be the same of the transportation model with addition of PS. On the other hand, higher costs will inhibit the presence of PS in the optimal solution, tending to the solution given by classical DC model. Finally, if the costs are competitive to transmission lines, an intermediary solution will be provided.


## 3. Solution Mathodology

Metaheuristic algorithms are specially suited for problems that present large search space with many local optima, such as the transmission expansion planning problem. The nonlinearity of the problem concerning the KVL is higher than the conventional model, thus degenerating even more the performance of more accurate methods. For instance, simulated annealing, genetic algorithms, and tabu search represent efficient methods for solving such problems. This work employed a modified version of the genetic algorithm presented in [4, 8] for solving the mathematical optimization problem. Indeed, extra modifications become necessary in order to guarantee the acceptable performance of the genetic algorithm.


## 4. Genetic Algorithm

This section presents the genetic algorithm developed for the planning problem considering the addition of PS.


### 4.1. Encoding

Each individual in a population (chromosome) is a proposed solution for the problem. In this work, an individual is encoded considering only the integer and binary variables. Remaining variables (continuous variables) are obtained from the linear programming (LP) solution. Thus, transmission lines and transformers are represented by decimal encoding (variablenij), whereas the PSs are represented by binary variables (δij).An example of this chromosome of length2nl is shown in Figure 5. The first nl positions show the number of transmission lines added for each configuration and the last nl positions indicate whether a PS was added or not in each path. The existing transmission lines are not coded but they are taken into account when operational variables are calculated, as well as when PSs are added in series to the transmission lines at each path. It can be noticed that in path 2, two transmission lines have been added with the corresponding PS. In path 3, one transmission line and the necessary number of PSs have been added, and so on.Figure 5
Encoding proposal (chromosome).In [8], an excellent performance of the genetic algorithm was obtained with the coding above, which is also in accordance with the coding proposals suggested for genetic algorithms presented in [4, 9, 10].


### 4.2. Objective Function and Infeasibility

The objective function of any solution proposal is found by solving an LP problem. The LP determines the exact values of the operational variables, which makes verifying the operation feasibility of a determined investment proposal possible, that is, whether the system presents load shedding to the implemented expansion proposal. Considering that an investment proposals is characterized by variables nijs and δijs, the load shedding is obtained by solving the following LP problem: (24)Minimizews=∑k∈Γrksubject to(25)S⋅f+g+r=d,(26)fij-γij⋅(nij0+nijs)⋅(θi-θj+φij⋅δijs)=0,(27)|fij|≤(nij0+nijs)⋅f¯ij,(28)0≤gi≤gi¯,(29)0≤rk≤dk,(30)fijUnbounded,(31)θjUnbounded,(32)φijUnbounded,(33)(i,j)∈Ω.For each solution proposal, the objective function is calculated with the following expression:(34)vs=∑(i,j)∈Ωcijnijs+∑dij⋅(nij0+nijs)⋅δijs+α⋅ws.In the genetic algorithm, every solution proposal is considered, including the infeasible ones. The infeasible configurations (with load shedding) are eliminated gradually by selection process, since these configurations are penalized by parameterα in the objective function.


### 4.3. Selection Process

The selection is based on tournament withk=2, that is, a game where two topologies of the current population participate. Recent publications show that selection by tournament is the most efficient method as long as parameter k is set adequately.


### 4.4. Crossover

The single point crossover was employed in this work. The crossover point is chosen randomly and a descendant, which has a parcel of its parents from the crossover point, is created. The random point was generated from an interval of 1 to (nl-1), that is, considering only the first half of the chromosome. In the second part the same operation is executed as shown in Figure 6. The objective of the operation is to transmit the existing information, which is present in both sectors of the chromosome to only one descendant in order to avoid quality degeneration.Figure 6
Single point crossover.


### 4.5. Mutation

The mutation operator acts in the following form. Considering transmission lines, the application of mutation operation means the addition or removal of one transmission line added during the optimization process (nij+1 or nij-1). For PS, it means the change of an allele, that is, inclusion or exclusion of the PS in the path.The mutation operation should be executed respecting the following conditions:(1)
the maximum number of transmission lines in the path;(2)
before adding a PS in the selected path, the existence of a transmission line must be checked; if it is an empty path, a transmission line has also to be added;(3)
when a transmission line is removed, the corresponding PS is removed, in case it exists.Concluding, a PS can only be inserted to an existing transmission line, whereas the number of PS is equal to the number of transmission lines in a branch. In Figure7, two cases of mutation are presented. The first case shows the random addition of PS to an existing transmission line, and the second case, mutation is applied to the transmission line that is removed (in this case there is no transmission line in the base topology); then the PS is removed in order to keep the process coherent.Figure 7
Mutation.


### 4.6. Proposed Genetic Algorithm

The general structure of the implemented genetic algorithm is similar to that presented in [8], except for some modifications. The basic structure of the algorithm is the following.(1)
Set the control parameters and generate the initial population. Make the initial population the current population.(2)
Calculate the objective function of the current population by solving one LP for each element (topology) of the current population. Update the incumbent solution whenever possible.(3)
If the stopping criterion is satisfied, stop the process. Otherwise go to step 4.(4)
Execute selection by tournament withk={2;3}.(5)
Execute one point crossover.(6)
Implement specialized mutation.(7)
Form the current population and go to step 2.Some details of the algorithm are presented in the next section.


### 4.7. Details and Improvement of the Genetic Algorithm

We present, briefly, details of the algorithm and some improvements made to the genetic algorithm.


#### 4.7.1. Generation of the Initial Population

Generation of the initial population is made by a controlled random process. Basically, it defines the number of paths in which the transmission lines are added and the maximum number of PS. Regarding transmission lines are defined randomly number of branches where the lines are inserted, position, and number of transmission lines (subject to the limits of added lines). In the case of PS, the number of branches and the position are selected randomly. In general, experience in transmission planning indicates that the number of branches to be added should be small, whereas the number for PS should be even smaller.The performance coefficient (35) can be employed for ranking the most interesting transmission lines. Observe that lines with a larger capacity that is also cost-effective have the priority in the addition: (35)PCij=fij¯γijcij.After ranked, the circuits are separated into two groups with different size (75%, 25%), whereas the largest one presents the most interesting transmission lines in terms of capacity usage. The initial population is formed by 80% up to 100% of elements belonging to the largest group. Another option to the initial population generation is to employ constructive heuristic algorithms as in [11, 12].


#### 4.7.2. Control Parameters: Crossover and Mutation Rates

The employed crossover rate wasRC=0.8. Two mutation rates were employed: RM=0.1 or RM=0.6 (the value depends on the current population diversification level). The mutation rate is defined differently from the classical concept and it applies for an entire solution proposal (topology). For example, a rate of 0.1 means that there is 10% of probability for a topology to suffer mutation. The population size varied from 40 to 200 elements.


#### 4.7.3. Mutation Details

Mutation is executed in the following way: the power flow of each topology is stored in four matrices considering the load level of each circuit. They are separated in intervals of 25, 50, 75, and 100% of capacity. The load level is calculated by means of the relationship (36): (36)FCij=fij(nij0+nij)⋅fij¯.Mutation operation is executed based on the probability of 70% of circuit removal and 30% of circuit addition. In case that circuit addition operation is selected, a random selection is carried out over the branches from the matrix with most loaded branches (up to 75%). If the matrix is empty, the next load level matrix is searched. For the circuit removal operation, the search starts from the least loaded branches (matrices with branches on 25% and 50% of its capacity). With this strategy, transmission lines with low utilization are removed from the system, whereas in regions where higher utilization is observed lines are inserted.


#### 4.7.4. Stopping Criterion

The algorithm stops when a defined total number of iterations is reached or when there is no improvement of the incumbent solution after a specified number of iterations.


#### 4.7.5. Elitism and Diversification Control

The population diversity is controlled by changing the mutation rate. A measurement of diversity is given by (37):(37)RDiv=Pop-NRPop.The diversification rate is calculated after mutation. If this rate is below 50%, there is used a mutation rate of 0.6; otherwise a rate of 0.1 is considered. This mechanism was applied aimed at maintaining diversity and at exploring new search spaces.Another strategy used was the elitism, in which the parent topologies are compared with the descendants and then the two best topologies are preserved in the current population.


## 4.1. Encoding

Each individual in a population (chromosome) is a proposed solution for the problem. In this work, an individual is encoded considering only the integer and binary variables. Remaining variables (continuous variables) are obtained from the linear programming (LP) solution. Thus, transmission lines and transformers are represented by decimal encoding (variablenij), whereas the PSs are represented by binary variables (δij).An example of this chromosome of length2nl is shown in Figure 5. The first nl positions show the number of transmission lines added for each configuration and the last nl positions indicate whether a PS was added or not in each path. The existing transmission lines are not coded but they are taken into account when operational variables are calculated, as well as when PSs are added in series to the transmission lines at each path. It can be noticed that in path 2, two transmission lines have been added with the corresponding PS. In path 3, one transmission line and the necessary number of PSs have been added, and so on.Figure 5
Encoding proposal (chromosome).In [8], an excellent performance of the genetic algorithm was obtained with the coding above, which is also in accordance with the coding proposals suggested for genetic algorithms presented in [4, 9, 10].


## 4.2. Objective Function and Infeasibility

The objective function of any solution proposal is found by solving an LP problem. The LP determines the exact values of the operational variables, which makes verifying the operation feasibility of a determined investment proposal possible, that is, whether the system presents load shedding to the implemented expansion proposal. Considering that an investment proposals is characterized by variables nijs and δijs, the load shedding is obtained by solving the following LP problem: (24)Minimizews=∑k∈Γrksubject to(25)S⋅f+g+r=d,(26)fij-γij⋅(nij0+nijs)⋅(θi-θj+φij⋅δijs)=0,(27)|fij|≤(nij0+nijs)⋅f¯ij,(28)0≤gi≤gi¯,(29)0≤rk≤dk,(30)fijUnbounded,(31)θjUnbounded,(32)φijUnbounded,(33)(i,j)∈Ω.For each solution proposal, the objective function is calculated with the following expression:(34)vs=∑(i,j)∈Ωcijnijs+∑dij⋅(nij0+nijs)⋅δijs+α⋅ws.In the genetic algorithm, every solution proposal is considered, including the infeasible ones. The infeasible configurations (with load shedding) are eliminated gradually by selection process, since these configurations are penalized by parameterα in the objective function.


## 4.3. Selection Process

The selection is based on tournament withk=2, that is, a game where two topologies of the current population participate. Recent publications show that selection by tournament is the most efficient method as long as parameter k is set adequately.


## 4.4. Crossover

The single point crossover was employed in this work. The crossover point is chosen randomly and a descendant, which has a parcel of its parents from the crossover point, is created. The random point was generated from an interval of 1 to (nl-1), that is, considering only the first half of the chromosome. In the second part the same operation is executed as shown in Figure 6. The objective of the operation is to transmit the existing information, which is present in both sectors of the chromosome to only one descendant in order to avoid quality degeneration.Figure 6
Single point crossover.


## 4.5. Mutation

The mutation operator acts in the following form. Considering transmission lines, the application of mutation operation means the addition or removal of one transmission line added during the optimization process (nij+1 or nij-1). For PS, it means the change of an allele, that is, inclusion or exclusion of the PS in the path.The mutation operation should be executed respecting the following conditions:(1)
the maximum number of transmission lines in the path;(2)
before adding a PS in the selected path, the existence of a transmission line must be checked; if it is an empty path, a transmission line has also to be added;(3)
when a transmission line is removed, the corresponding PS is removed, in case it exists.Concluding, a PS can only be inserted to an existing transmission line, whereas the number of PS is equal to the number of transmission lines in a branch. In Figure7, two cases of mutation are presented. The first case shows the random addition of PS to an existing transmission line, and the second case, mutation is applied to the transmission line that is removed (in this case there is no transmission line in the base topology); then the PS is removed in order to keep the process coherent.Figure 7
Mutation.


## 4.6. Proposed Genetic Algorithm

The general structure of the implemented genetic algorithm is similar to that presented in [8], except for some modifications. The basic structure of the algorithm is the following.(1)
Set the control parameters and generate the initial population. Make the initial population the current population.(2)
Calculate the objective function of the current population by solving one LP for each element (topology) of the current population. Update the incumbent solution whenever possible.(3)
If the stopping criterion is satisfied, stop the process. Otherwise go to step 4.(4)
Execute selection by tournament withk={2;3}.(5)
Execute one point crossover.(6)
Implement specialized mutation.(7)
Form the current population and go to step 2.Some details of the algorithm are presented in the next section.


## 4.7. Details and Improvement of the Genetic Algorithm

We present, briefly, details of the algorithm and some improvements made to the genetic algorithm.


### 4.7.1. Generation of the Initial Population

Generation of the initial population is made by a controlled random process. Basically, it defines the number of paths in which the transmission lines are added and the maximum number of PS. Regarding transmission lines are defined randomly number of branches where the lines are inserted, position, and number of transmission lines (subject to the limits of added lines). In the case of PS, the number of branches and the position are selected randomly. In general, experience in transmission planning indicates that the number of branches to be added should be small, whereas the number for PS should be even smaller.The performance coefficient (35) can be employed for ranking the most interesting transmission lines. Observe that lines with a larger capacity that is also cost-effective have the priority in the addition: (35)PCij=fij¯γijcij.After ranked, the circuits are separated into two groups with different size (75%, 25%), whereas the largest one presents the most interesting transmission lines in terms of capacity usage. The initial population is formed by 80% up to 100% of elements belonging to the largest group. Another option to the initial population generation is to employ constructive heuristic algorithms as in [11, 12].


### 4.7.2. Control Parameters: Crossover and Mutation Rates

The employed crossover rate wasRC=0.8. Two mutation rates were employed: RM=0.1 or RM=0.6 (the value depends on the current population diversification level). The mutation rate is defined differently from the classical concept and it applies for an entire solution proposal (topology). For example, a rate of 0.1 means that there is 10% of probability for a topology to suffer mutation. The population size varied from 40 to 200 elements.


### 4.7.3. Mutation Details

Mutation is executed in the following way: the power flow of each topology is stored in four matrices considering the load level of each circuit. They are separated in intervals of 25, 50, 75, and 100% of capacity. The load level is calculated by means of the relationship (36): (36)FCij=fij(nij0+nij)⋅fij¯.Mutation operation is executed based on the probability of 70% of circuit removal and 30% of circuit addition. In case that circuit addition operation is selected, a random selection is carried out over the branches from the matrix with most loaded branches (up to 75%). If the matrix is empty, the next load level matrix is searched. For the circuit removal operation, the search starts from the least loaded branches (matrices with branches on 25% and 50% of its capacity). With this strategy, transmission lines with low utilization are removed from the system, whereas in regions where higher utilization is observed lines are inserted.


### 4.7.4. Stopping Criterion

The algorithm stops when a defined total number of iterations is reached or when there is no improvement of the incumbent solution after a specified number of iterations.


### 4.7.5. Elitism and Diversification Control

The population diversity is controlled by changing the mutation rate. A measurement of diversity is given by (37):(37)RDiv=Pop-NRPop.The diversification rate is calculated after mutation. If this rate is below 50%, there is used a mutation rate of 0.6; otherwise a rate of 0.1 is considered. This mechanism was applied aimed at maintaining diversity and at exploring new search spaces.Another strategy used was the elitism, in which the parent topologies are compared with the descendants and then the two best topologies are preserved in the current population.


## 4.7.1. Generation of the Initial Population

Generation of the initial population is made by a controlled random process. Basically, it defines the number of paths in which the transmission lines are added and the maximum number of PS. Regarding transmission lines are defined randomly number of branches where the lines are inserted, position, and number of transmission lines (subject to the limits of added lines). In the case of PS, the number of branches and the position are selected randomly. In general, experience in transmission planning indicates that the number of branches to be added should be small, whereas the number for PS should be even smaller.The performance coefficient (35) can be employed for ranking the most interesting transmission lines. Observe that lines with a larger capacity that is also cost-effective have the priority in the addition: (35)PCij=fij¯γijcij.After ranked, the circuits are separated into two groups with different size (75%, 25%), whereas the largest one presents the most interesting transmission lines in terms of capacity usage. The initial population is formed by 80% up to 100% of elements belonging to the largest group. Another option to the initial population generation is to employ constructive heuristic algorithms as in [11, 12].


## 4.7.2. Control Parameters: Crossover and Mutation Rates

The employed crossover rate wasRC=0.8. Two mutation rates were employed: RM=0.1 or RM=0.6 (the value depends on the current population diversification level). The mutation rate is defined differently from the classical concept and it applies for an entire solution proposal (topology). For example, a rate of 0.1 means that there is 10% of probability for a topology to suffer mutation. The population size varied from 40 to 200 elements.


## 4.7.3. Mutation Details

Mutation is executed in the following way: the power flow of each topology is stored in four matrices considering the load level of each circuit. They are separated in intervals of 25, 50, 75, and 100% of capacity. The load level is calculated by means of the relationship (36): (36)FCij=fij(nij0+nij)⋅fij¯.Mutation operation is executed based on the probability of 70% of circuit removal and 30% of circuit addition. In case that circuit addition operation is selected, a random selection is carried out over the branches from the matrix with most loaded branches (up to 75%). If the matrix is empty, the next load level matrix is searched. For the circuit removal operation, the search starts from the least loaded branches (matrices with branches on 25% and 50% of its capacity). With this strategy, transmission lines with low utilization are removed from the system, whereas in regions where higher utilization is observed lines are inserted.


## 4.7.4. Stopping Criterion

The algorithm stops when a defined total number of iterations is reached or when there is no improvement of the incumbent solution after a specified number of iterations.


## 4.7.5. Elitism and Diversification Control

The population diversity is controlled by changing the mutation rate. A measurement of diversity is given by (37):(37)RDiv=Pop-NRPop.The diversification rate is calculated after mutation. If this rate is below 50%, there is used a mutation rate of 0.6; otherwise a rate of 0.1 is considered. This mechanism was applied aimed at maintaining diversity and at exploring new search spaces.Another strategy used was the elitism, in which the parent topologies are compared with the descendants and then the two best topologies are preserved in the current population.


## 5. Tests and Results

In order to analyze the performance of the proposed genetic algorithm and to demonstrate viability of the mathematical model with PS, tests were carried out with the IEEE-24 bus system. This system has 41 circuit paths, 8550 MW of load, generation capacity of 10215 MW, and five different generation plans whose data is present in [5].For all simulations, a fixed cost for the PS was adopted, that varies among 2 and 120 million of dollars. In all the tests was used a value ofα=106 US$/MW.The simulations were divided in three stages. In the first stage, equipments of low value were added. In the following stage, equipments of high value were added. Finally, in the last stages, intermediate costs were adopted; that is, in each stage different values were fixed for the PS.


### 5.1. Low-Cost Phase Shifters

With the purpose of testing the algorithm, the first simulations set the cost of PS as low asdij=2×106 (US$). Table 3 shows the results obtained.Table 3
Five transmission expansion plans considering low-cost phase-shifter transformers.
CircuitsPlanP0PlanP1PlanP2PlanP3PlanP4TLPSTLPSTLPSTLPSTLPSn01-031n03-091n04-091n03-241n06-1011111n07-0822122n08-0911n08-101n10-1211n11-1411n12-2311n13-231n14-161111n15-161n16-172211n16-1911n17-1811n17-221n19-222n13-1411Partial cost102422663108110141884Total cost106232318124192A very interesting fact in the simulations with low-cost PS refers to transmission lines. For the five plans illustrated in Table3, it is noticed that the lines added by the proposed specialized genetic algorithm were the same of the optimal solution for the transportation model. This fact was confirmed by utilizing a branch-and-bound algorithm presented in [13].Results also confirmed that the optimal solution for the transportation model is not feasible for the DC model since KVL constraints are violated. Nonetheless, during the simulation, this problem was overcome by adding low-cost PS in strategic branches.


### 5.2. High-Cost Phase Shifters

Here there is considereddij=120×106 (US$) for high-cost PS. For the same five generation plans, obtained results are illustrated in Table 4.Table 4
Five transmission expansion plans considering high-cost phase-shifter transformers.
CircuitsPlanP0PlanP1PlanP2PlanP3PlanP4TLPSTLPSTLPSTLPSTLPSn01-0511n03-24111n06-1011111n07-0822122n09-111n10-121111n11-141n14-1611111n15-2411n16-172211n16-191n17-1822n20-231Partial cost15203900392021803460Total cost152390392218346Again, allocation of lines deserves importance. Simulations carried out with high-cost PS indicate that there was no addition of PS in none of the five plans of expansion of the system. The justification for the absence of PS, in all the topologies, is that they are too expensive. Thus, they did not take part of the optimal solutions. Consequently, there is a tendency of solely adding transmission lines.Another interesting fact that deserves emphasis is that the results are the same found to the conventional DC model.


### 5.3. Intermediary-Cost Phase Shifters

PSs with arbitrary intermediary costs were employed in order to produce expansion proposals with intermediary values when compared to extreme solutions obtained in previous simulations. The PS was considered with cost ofdij=55×106 (US$) for every branch in plan P1. Obtained results are shown in Table 5.Table 5
PlanP1 of the system IEEE-24bus considering intermediary-cost phase-shifter transformers.
CircuitsPlanP0PlanP1PlanP2TLPSTLPSTLPSn03-24111n06-10111n07-08222n08-101n09-121n10-1111n10-121n12-231n14-16111n16-17222n16-19111n17-18111Partial cost276110276110276110Total cost386386386The proposed algorithm found alternative optimal solutions for the planP1. Observe that circuit 10-11, which is not present in the previous cases (PS with low and high costs), appears now. Notice that in the three topologies found for the plan P1 were added two phase-shifter transformers in the network.It is important to mention that all tests have been carried out for PS with nonrealistic cost values in order to test, from the theoretical point of view, the feasibility of modeling such devices as expansion components for electrical systems. Additional network transmission planning bibliography can be found in [14]. It should be observed that this novel model could be solved by other metaheuristics as presented in [9, 15, 16] after small adaptations.


## 5.1. Low-Cost Phase Shifters

With the purpose of testing the algorithm, the first simulations set the cost of PS as low asdij=2×106 (US$). Table 3 shows the results obtained.Table 3
Five transmission expansion plans considering low-cost phase-shifter transformers.
CircuitsPlanP0PlanP1PlanP2PlanP3PlanP4TLPSTLPSTLPSTLPSTLPSn01-031n03-091n04-091n03-241n06-1011111n07-0822122n08-0911n08-101n10-1211n11-1411n12-2311n13-231n14-161111n15-161n16-172211n16-1911n17-1811n17-221n19-222n13-1411Partial cost102422663108110141884Total cost106232318124192A very interesting fact in the simulations with low-cost PS refers to transmission lines. For the five plans illustrated in Table3, it is noticed that the lines added by the proposed specialized genetic algorithm were the same of the optimal solution for the transportation model. This fact was confirmed by utilizing a branch-and-bound algorithm presented in [13].Results also confirmed that the optimal solution for the transportation model is not feasible for the DC model since KVL constraints are violated. Nonetheless, during the simulation, this problem was overcome by adding low-cost PS in strategic branches.


## 5.2. High-Cost Phase Shifters

Here there is considereddij=120×106 (US$) for high-cost PS. For the same five generation plans, obtained results are illustrated in Table 4.Table 4
Five transmission expansion plans considering high-cost phase-shifter transformers.
CircuitsPlanP0PlanP1PlanP2PlanP3PlanP4TLPSTLPSTLPSTLPSTLPSn01-0511n03-24111n06-1011111n07-0822122n09-111n10-121111n11-141n14-1611111n15-2411n16-172211n16-191n17-1822n20-231Partial cost15203900392021803460Total cost152390392218346Again, allocation of lines deserves importance. Simulations carried out with high-cost PS indicate that there was no addition of PS in none of the five plans of expansion of the system. The justification for the absence of PS, in all the topologies, is that they are too expensive. Thus, they did not take part of the optimal solutions. Consequently, there is a tendency of solely adding transmission lines.Another interesting fact that deserves emphasis is that the results are the same found to the conventional DC model.


## 5.3. Intermediary-Cost Phase Shifters

PSs with arbitrary intermediary costs were employed in order to produce expansion proposals with intermediary values when compared to extreme solutions obtained in previous simulations. The PS was considered with cost ofdij=55×106 (US$) for every branch in plan P1. Obtained results are shown in Table 5.Table 5
PlanP1 of the system IEEE-24bus considering intermediary-cost phase-shifter transformers.
CircuitsPlanP0PlanP1PlanP2TLPSTLPSTLPSn03-24111n06-10111n07-08222n08-101n09-121n10-1111n10-121n12-231n14-16111n16-17222n16-19111n17-18111Partial cost276110276110276110Total cost386386386The proposed algorithm found alternative optimal solutions for the planP1. Observe that circuit 10-11, which is not present in the previous cases (PS with low and high costs), appears now. Notice that in the three topologies found for the plan P1 were added two phase-shifter transformers in the network.It is important to mention that all tests have been carried out for PS with nonrealistic cost values in order to test, from the theoretical point of view, the feasibility of modeling such devices as expansion components for electrical systems. Additional network transmission planning bibliography can be found in [14]. It should be observed that this novel model could be solved by other metaheuristics as presented in [9, 15, 16] after small adaptations.


## 6. Conclusions

Nowadays, modern elements, such as FACTS devices, are playing an important role in transmission systems. In this way, inclusion of such devices jointly with classical components is of importance for the transmission expansion planning problem. Thus, this work was aimed at presenting the technical feasibility of considering phase-shifter transformers as components for the long-term transmission expansion planning, jointly with conventional transformers and transmission lines.A novel methodology was proposed for the inclusion of phase-shifter transformers in the mathematical model that represents the transmission planning problem. The proposed model is more complex than the model DC. However, the mathematical problem was solved adequately with genetic algorithms.Tests have shown the model consistency as well as the high performance of the algorithm.This contribution extends the utilization of classical components during the expansion-planning problem, to modern elements, such as the FACTS devices.



---

*Source: 527258-2012-03-22.xml*