Optimal Integration of Flexible Alternating Current Transmission Systems in Electrical Distribution Grids Using a Mixed-Integer Convex Model

: This research addresses the efﬁcient integration and sizing of ﬂexible alternating current transmission systems (FACTS) in electrical distribution networks via a convex optimization approach. The exact mixed-integer nonlinear programming (MINLP) model associated with FACTS siting and sizing aims for the minimization of the expected annual operating costs of the network (i


Introduction 1.General Context
Electrical distribution networks in medium-and low-voltage applications provide electricity to all end-users while observing quality, security, reliability, and efficiency criteria [1,2].These grids are typically operated at voltage levels lower or equal to 25 kV, mostly using alternating current (AC) technologies.Most of them are built with a tree structure (i.e., radial configuration) that helps to reduce investment costs and the complexity of coordinating protective devices [3,4].However, these radial configurations typically increase energy losses when compared to meshed topologies [5], and they deteriorate the voltage profiles (low voltage regulation) for end-users located far from the central substation nodes [6].Distribution companies typically use shunt compensation methodologies or grid topology modifications to address these technical challenges.In the case of shunt compensation, active and reactive power injection devices are commonly used (dispersed generation [7], batteries [8], capacitors [9], or flexible AC transmission systems [10,11], i.e., FACTS).In contrast, in the case of topology modifications, phase balancing or grid reconfiguration are common strategies for distribution companies.In the case of power loss reduction or voltage profile improvement, shunt reactive power compensation has demonstrated an adequate balance between saving operating costs regarding energy losses compared to the investment costs associated with compensation systems [9].
On the other hand, given the transition of electrical distribution networks from passive to active operation schemes [12], the use of fixed compensation strategies [13], such as capacitor banks or grid topology variations, results in low economic savings in comparison with dynamic compensation strategies using FACTS [10].
The use of FACTS compensation devices in electrical distribution networks has shown positive effects on the technical characteristics of electrical grids, especially in the case of energy losses minimization, with improvements of about 12.63% and 13.97% in the conventional IEEE 33-and 69-bus systems [10].In addition, for distribution networks in the Colombian context, energy losses can oscillate between 14.5% and 21.8% of the total energy input [10].This implies that the development of compensation technologies for reducing these losses is essential for distribution companies.In this regard, one of the advantages of reactive power compensation systems is that they are reliable [14] and economical [15], in addition to the fact that they have long useful lives [16].

Motivation
The efficient integration (nodal location and nominal size) and operation (on an hourly basis) of FACTS in electrical distribution networks pose fundamental challenges regarding the nature of this optimization problem as it belongs to the family of mixedinteger nonlinear programming (MINLP) [17].Given the complexity of the exact model, this research aims to propose an efficient solution methodology to locate and size FACTS in medium-voltage distribution networks by reformulating the MINLP model as a mixedinteger convex approximation [18].The main advantage of convex optimization is that it allows representing optimization problems in engineering and science while ensuring a global optimal solution [19].This research aims to provide the literature and the industry with an efficient and reliable solution methodology to address the studied problem, with superior performance compared to MINLP solvers in commercial software or metaheuristic algorithms [18].
It is important to mention that studying the problem regarding the optimal location and sizing of FACTS in electrical distribution networks is necessary since: i.
Dynamic reactive power compensation is an area of continuous development due to the quality impositions of regulatory offices on distribution companies, which aim to make their medium-voltage distribution grids efficient, reliable, and secure [20].To this effect, effective optimization algorithms must be proposed; ii.Most existing solution methodologies focus on metaheuristic optimization algorithms to deal with the exact MINLP formulation via decoupling-based approaches that work under a master-slave optimization strategy [21].However, even though these master-slave approaches are efficient and easily implemented in multiple programming languages, they do not allow ensuring a global optimum, given the random nature of heuristic-based optimizers [22].Therefore, this research takes advantage of mixed-integer convex programming to propose an efficient solution methodology to locate and size FACTS in electrical networks while allowing us to find the global optimum [23].

Literature Review
The problem regarding the optimal placement and sizing of FACTS in medium-and low-voltage distribution networks has been widely explored in the specialized literature.This section summarizes the most recent approaches in this research area.
The authors of [11] applied the hunter-prey algorithm to determine the optimal location and sizing of photovoltaic-SVC systems in electrical distribution networks, with the aim to minimize grid power losses and improve system voltage profiles.The main characteristic of a photovoltaic-SVC system is that it can use the converter interfacing with the photovoltaic (PV) system to add reactive power injection capabilities via control design [24].The numerical results presented in [11] for the IEEE 33-and 69-bus grids demonstrated the effectiveness of the proposed hunter-prey optimizer when compared to different combinatorial methods, such as the differential evolution algorithm, particle swarm optimization, the artificial rabbits algorithm, and the golden search optimizer.
The work by [25] presented a complete comparative analysis of multiple combinatorial optimizers in locating and sizing TCSCs and SVCs in power systems while aiming to minimize the total grid operating costs associated with energy losses.The voltage profiles were kept within an acceptable range.Numerical validations in the IEEE 30-and 57-bus systems were used to test the effectiveness of the whale optimization algorithm in comparison with multiple combinatorial optimizers.
In [10], the black widow optimization algorithm was proposed as a solution method to locate and size FACTS in medium-voltage distribution networks.The exact MINLP model was solved using a leader-follower optimization strategy, combining a discrete-continuous codification in the leader stage with a black widow optimizer using the successive approximations of power flow in the follower stage.Numerical results in the IEEE 33-, 69-, and 85-bus grids demonstrated the effectiveness of the leader-follower solution compared to that of the vortex search algorithm reported in [26].All this, while considering the minimization of investment and operating costs as objective functions.
The authors of [27] presented a complete review of the optimal location and sizing of SVCs in electrical distribution networks while considering the technical and economic aspects of the objective function.Five approaches regarding solution technologies were reviewed, which include artificial neural network techniques, analytical methods, combinatorial methods, and sensitivity approaches.This review demonstrated that SVCs are suitable FACTS for improving grid voltage profiles, enhancing stability indices, and reducing grid power losses with reasonable investments by utility companies.Thus, SVCs are a feasible economic solution for compensating reactive power in these networks.
The study by [28] presented a solution methodology to locate and size SVCs in distribution networks while considering load variations via a nodal sensitivity approach.Efficient numerical validation in practical distribution networks composed of 38 nodes demonstrated important improvements with regard to different objective functions, including energy losses cost reduction, voltage profile improvements, and stability margin enhancement.
Other solution methods for locating and sizing FACTS in distribution networks are the tabu search algorithm [29,30], the fractional levy flight bat algorithm [31], the particle swarm optimizer [32], the krill herd algorithm [33], and the barnacles mating optimization [34], among others.
All of the works reviewed above yielded good solutions.However, none of them can guarantee the optimal solution to the problem, and some require tuning parameters, which indicates that, if the test systems are changed, they may not reach reasonable solutions and thus require adjustments.

Contribution and Scope
In light of the above, it is worth noting that, in order to analyze it, the problem under study is standardized via the application of leader-follower optimization techniques (combinatorial methods to decide on the location of the FACTS and optimal power flows to determine their sizes).This provides a research opportunity in the area of convex optimization.In this sense, the main contributions of this research are presented below: i.
A reformulation of the exact MINLP model regarding the efficient location and sizing of FACTS in distribution networks via a mixed-integer second-order cone programming equivalent.This reformulation can guarantee the global optimum of the problem.Thus, it finds the best solution to the optimization problem according to its fitness function.ii.The evaluation of multiple operation scenarios where the FACTS can be operated as fixed or variable compensators.These scenarios demonstrated that variable compensation is the best option to reduce the annual grid operating costs.
Note that, within the scope of this research, the FACTS devices analyzed correspond to the shunt reactive power compensation elements studied by [10].These devices are the unified power flow controller (UPFC), the thyristor-controlled shunt compensator (TCSC), and the static var compensator (SVC).The investment costs of installing these devices are modeled as a cubic function based on the recommendations of [35].In addition, the active and reactive power curves are inputs provided by the distribution company at the terminals of the substation bus, representing the electrical network's daily behavior.Here, these curves are assumed as constants without noise (uncertainties).However, future research will be required to include the stochastic behavior of these curves in distribution systems analysis, operation, and control.

Document Structure
The remainder of this document is structured as follows.Section 2 presents the general MINLP model regarding the optimal location and sizing of FACTS in electrical distribution networks, with the aim to minimize the annual grid operating costs while including the investment costs of FACTS.Section 3 shows the second-order cone approximation based on the hyperbolic relation of the product between two variables, which generates a linear objective function with a set of linear and conic constraints that belongs to the family of convex optimization problems.Section 4 outlines the main characteristics of the IEEE test feeders under analysis.These systems are composed of 33, 69, and 85 nodes with medium-voltage profiles and typical radial structures.Section 5 describes the main numerical simulations carried out, a comparative study with literature reports and exact MINLP solvers, and a comparative analysis between all the studied FACTS devices, considering fixed and variable reactive power injections.Section 6 describes this article's main concluding remarks and possible future improvements.

Optimization Model
The problem regarding the optimal integration of FACTS devices in electrical distribution grids aims to minimize the annual costs related to energy losses and the investments made in FACTS installation.This problem generates an optimization model with the structure of a mixed-integer nonlinear programming (MINLP) model.This is because the optimization model contains binary/integer variables related to the location of the FACTS devices.At the same time, the model includes the continuous variables associated with power flows, nodal voltages, and the size of the FACTS devices.Figure 1 illustrates an example of a branch in an electrical network with FACTS devices installed.
In this paper, the following notation is used.The sets H, B, and N encompass all the analysis periods, network branches, and network nodes, respectively.R and C denote the sets of real and complex numbers, respectively.The subscripts k or m denote the nodes of the test system, the subscript l represents any branch connected between nodes k and m, and the subscript h represents a specific period under analysis.The superscripts s and r denote the sending and receiving power flows of the distribution line connected between nodes k and m, while the superscripts g and d represent variables related to generation and demand, respectively.The variable v denotes the nodal voltage, and s corresponds to apparent power in a node or distribution line.The operator • denotes the Euclidean norm, while the operators real(•) and imag(•) correspond to the real and imaginary part of the complex number, respectively.Finally, (•) * represents the conjugate of the complex number.

Objective Function
The main aim of this study is to incorporate FACTS devices into electrical distribution networks while aiming to minimize the annual equivalent operating costs represented by f .These costs consider the reduction of energy losses and the installation expenses associated with FACTS.To this effect, the following objective function is used: where f 1 is related to the annual energy loss costs; f 2 is associated with the FACTS investment costs; C and T are the average costs of energy losses and the number of days in a year, respectively; p s lh and p r lh represent the sending and receiving active power flows of a distribution line, respectively; ∆h represents the time interval analyzed on a daily basis (0.5 h); q FACT k is the nominal size of the FACTS; ω 1 , ω 2 , and ω 3 represent the polynomial coefficients of the objective function f 2 ; and k 1 > 0 and k 2 > 0 are constants denoting the annual investment costs for a planning horizon of ten years [18].

Set of Constraints
The constraints related to the optimal integration of FACTS in electrical distribution grids comprise several aspects, such as node power balance constraints for both active and reactive power, the limits regarding the maximum and minimum power flowing through the distribution lines, voltage regulation bounds, requirements regarding the number of FACTS devices to be installed, and their capacity for injecting or observing reactive power, among others.

Power Balance Equation
The nodal power balance equations for node k with regard to active and reactive power require calculating the sums of power generated or demanded, i.e., for active and reactive power.These sums are then equalized to the flow of injection power, be it active or reactive.The equations can be expressed as follows: where p g kh and p d kh are the active power generated and demanded; q g kh and q d kh are the reactive power generated and demanded; p s lh and q s lh are the sending active and reactive power flows of the distribution line l; p r lh and q r lh are the receiving active and reactive power flows of the distribution line l; q FACT kh is the reactive power injected by the FACTS at node k and time h; and A + and A − are the positive and negative values of the node-to-branch incidence matrix A, i.e., A = A + + A − .

Power Flow Equation
The active and reactive power flows transmitted through a specific distribution line, along with their corresponding maximum capacities, can be represented as follows: where v kh and v mh are the complex voltages, y l is the admittance of branch l, and s max l is the maximum apparent power flowing through the line.

Operating Regulations
All node voltage values in an electrical network must meet the limits set by regulatory policies, which are intended for an adequate network operation.These limits are where v 0h represents the nominal voltage at the slack node (substation) at time h; v nom is its nominal value, which usually has a per-unit value equal to 1.0; and v min and v max denote the minimum and maximum voltage values allowed in the electrical network.

Integration of FACTS
The optimal integration of the FACTS is divided into two parts, which include their proper location in the electrical grid as well as their size.The following constraints are defined to address these challenges: where q FACT k is the variable used to define the size of the FACTS to be installed, and q FACT max is its maximum capacity.This study's capacity goes from 0 Mvar to 2 Mvar at the distribution level, a value usually employed for electrical networks [18].z is a vector of the binary variable that denotes a FACTS location at node k, i.e., if z = 1 in position k, this indicates that a device will be located there (otherwise, z = 0).Finally, η is the maximum number of FACTS to be installed.

Interpretation of the Mathematical Model
The mixed-integer nonlinear optimization model described in (1)-( 16) aims to optimize FACTS integration in electrical distribution grids.This model incorporates both binary and continuous variables.The placement of these devices is determined by the former.The latter is associated with several variables, including the active and reactive power of the generators, power demands, the power flowing through the transmission lines, complex nodal voltages, and the size of the FACTS.
The following can be noted regarding this mixed-integer nonlinear optimization model.The objective function is defined by two terms.The first term, which implies an annual cost of f 1 , represents the energy losses of the electrical distribution system.In contrast, the second term calculates the investment costs of FACTS ( f 2 ).Expressions ( 2) and (3) represent the active and reactive power balance, respectively, for each node and time period under analysis.Equations ( 4) and ( 5) represent the active power flow that is transmitted through each branch of the transmission lines at each studied time interval.Equations ( 6) and ( 7) represent the reactive power flowing through the same transmission lines.Inequalities (8) and ( 9) limit the maximum apparent power that can flow through each branch of the transmission lines in any given period.Equation (10) sets the voltage at the substation, while inequalities ( 11) and ( 12) constrain the nodal voltages to their respective minimum and maximum values at each time step.Inequality (13) restricts the maximum value that the FACTS may reach.Inequality ( 14) sets the limit for the maximum reactive power that the FACTS may deliver or absorb at each node.Finally, inequality (15) determines the maximum number of FACTS that can be installed.

Convex Reformulation
The optimization model described in Equations ( 1)-( 16) is a mixed-integer nonlinear one, which is challenging to solve and falls into the category of problems with high computational complexity.Therefore, this type of problem is typically solved using a metaheuristic algorithm [18,36].However, metaheuristic algorithms cannot ensure the global optimum of the problem.Additionally, many of these algorithms require parameter tuning, which indicates that their performance is not always consistent.Another potential solution to this problem is introducing some relaxations that can transform the mixedinteger nonlinear model into a mixed-convex one, thus being able to guarantee the global optimum of the problem.Nevertheless, this also requires parameter tuning.

Approximation of the Objective Function z 2 to Linear Function
The objective function f 2 presented in (1) has a cubic form, which makes it a nonconvex function.Therefore, it is impossible to guarantee the global optimum of the problem.However, this objective function can only work with a linear coefficient ω 3 for a range of FACTS with values lower than or equal to 2 Mvar [18].The objective function z 2 is expressed as follows:

Convex Representation of the Active and Reactive Power Flow Equations
The products of the voltages presented in the active and reactive power flow equations, as described in (4)- (7), are equality constraints.Hence, they are non-convex constraints.However, it is possible to convert these constraints into convex ones by defining two auxiliary variables [37], as follows: where u kh ∈ R corresponds to the squared voltage at node k and time h, and w lh ∈ C denotes the product of the voltages in branch l at time h.Now, the auxiliary variables u kh and w lh can be substituted into the active and reactive power flow Equations ( 4)- (7) in order to obtain: The equations for the active and reactive power flows, as presented in ( 20)-( 23), depend on the auxiliary variables in (18) and (19).These equations are non-convex, so it is necessary to relax them via a hyperbolic shape, as shown below:

Proposed Mixed-Integer Convex Model
The optimization model described in (1)-( 16) can be transformed into a mixed-integer convex model by applying the abovementioned relaxation.The process is as follows: This mixed-integer convex model can reach the global optimum of the exact optimization model described in Equations ( 1)-( 16).This is possible only if the hyperbolic constraints are well-defined conditions, as demonstrated in [38].

Test System
This section presents the test systems employed to validate the proposed optimization model for the optimal integration of FACTS in the IEEE 33-, 69-, and 85-bus systems.Figure 2 depicts the topologies of three test systems, whose main features are as follows: i.
The IEEE 33-bus test system has 33 buses and 32 transmission lines in its radial configuration, as shown in Figure 2a.It has a substation at node 1 that works with 12.66 kV, as well as peak active and reactive demands of 3715 + j2300 kVA, respectively.These operating conditions generate active and reactive power losses of 210.9876 kW and 143.1283 kvar.The test system's peak demand, resistance, and reactance values are listed in Table 1.These values were taken from [39].ii.The IEEE 69-bus test system is fitted with 69 buses and 68 transmission lines, as depicted in Figure 2b.Its substation is located at node 1, which works at 12.66 kV, and its peak active and reactive demand is 3890.7 + j2693.6 kVA.Under these operating conditions, the active and reactive power losses generated are 210.9876kW and 143.1283 kvar, respectively.Table 2 lists this system's demand, resistance, and reactance values.These values were taken from [39].iii.The IEEE 85-bus test system has 85 buses and 84 transmission lines, as illustrated in Figure 2c.It has a substation at node 1 that works with 11 kV, as well as a peak active and reactive demand of 2570.28 + j2622.20 kVA.The test system's peak demand, resistance, and reactance values are listed in Table 3.These values were taken from [39].Typically, electrical distribution systems exhibit a daily load variation.The demand curves for the active and reactive power are considered as illustrated in Figure 3.These curves represent the typical behavior of an electrical distribution system in Colombia [40].The parameter values of the objective function z (1) are shown in Table 4.The costs associated with the FACTS were taken from [39].

Numerical Implementation
The proposed optimization model was implemented in the Yalmip toolbox (R20230622 version) [41], using the Gurobi 9.5.1 solver [42] in the MATLAB 2021a software.A Dell Inspiron 15 7000 Series (Intel Quad-Core i7-7700HQ @2.80 GHz) PC (Dell Inc., Round Rock, TX, USA; Intel Corporation, Santa Clara, CA, USA) with 16 GB RAM and 64-bit Windows 10 Home Single Language (Microsoft Corporation, Redmond, WA, USA) was used to carry out the simulations.Furthermore, the exact optimization model ( 1)-( 16) was also implemented in the GAMS software (23.5.1 version).
The four scenarios shown below were proposed to evaluate the performance of the proposed mixed-integer convex model ( 25)- (42).All scenarios consider a maximum of three FACTS.S1: The proposed model was compared against GAMS, the black widow optimization (BWO) algorithm proposed in [10], and the vortex search algorithm (VSA) presented by [26] for the IEEE 33-bus test system.
S2: The proposed convex model was evaluated in the IEEE 69-bus test system and compacted to the BWO and VSA, with the aim of installing the SVC devices with a fixed variable operation.
S3: The model was evaluated in a large test system, namely the IEEE 85-bus test system.Furthermore, it was compared to the BWO algorithm regarding the installation of SVCs with a fixed and variable operation.
S4: A comparison regarding the installation of different FACTS technologies was analyzed (i.e., SVC, TCSC, and UPFC devices).

Analysis of Scenario 1 (S1)
This scenario evaluated and compared the proposed optimization model against the BWO and VSA, and the GAMS software with three different solvers.For this comparison, the installation of only SVC devices was considered.Table 5 presents the results obtained for the objective function in the radial IEEE 33-bus test system.This scenario was analyzed by considering whether the FACTS-delivered/injected power was fixed or variable during daily operations.From the results shown in Table 5, it can be stated that: i.
The proposed convex model, as well as the BWO and VSA, and the BONMIM solver, reached the best configuration for the SVC devices.This configuration is the global optimum of the problem, as the convex model ensures it.According to the results, node 30 is the most sensitive point of the test system, given that the largest SVC is located in it.This occurs for both the fixed and the variable operation of SVC devices.ii.The variable operation of SCV devices reduces the value of the objective function by 14.24%.In contrast, the fixed operation reports a value of 12.63%.This indicates that it is better for the test system to implement a variable operation, as SVCs can inject or absorb reactive power according to the system requirements at each hour of the day.Furthermore, even though the size of the SVCs increases for a variable operation, the objective function value is lower than that of the fixed operation.The total size of the SVC devices for the fixed and variable operation is 0.6262 and 1.0102 Mvar, respectively, increasing the investment costs by 61.17%.

Analysis of Scenario 2 (S2)
This scenario studied the performance of the proposed optimization model and compares it against that of the BWO and VSA.These algorithms were only used for a fixed operation of the FACTS.In this scenario, no solver of the GAMS reaches convergence.As in the previous scenario, the installation of only SVC devices was considered, given that [10] only analyzes SVCs for the IEEE 69-bus system.The results of the objective function for the IEEE 69-bus test system are shown in Table 6, considering the fixed or variable operation of SVC devices during the day.From the results obtained in Table 6, note that: i.The solutions found by the BWO and VSA in [10,26] are the global optimum of the problem as the proposed convex model achieves the same configuration.Still, it is essential to note that the proposed convex model will always find the same values, while the BOW and VSA cannot guarantee these results.Furthermore, these algorithms require tuning many parameters, which can affect their performance.ii.According to the results obtained for the IEEE 69-bus test system, the node with the highest sensitivity is node 61, as it was selected to install the SVC device with the highest capacity.This behavior is the same for the fixed and variable SVC operations.iii.The optimal integration of SVC devices reduces the annual energy losses costs.This reduction is greater for the variable operation.These devices report reductions of 13.97% and 15.79% in the objective functions for the fixed and variable operation.The latter saves USD 3778.4/yearmore than the former.Despite this, the total size of the SVCs increases by 24.35% for the variable operation.

Analysis of Scenario 3 (S3)
This scenario evaluated the performance of the proposed optimization model in a large electrical distribution network, such as the IEEE 85-bus test system.Furthermore, the model was compared to the BWO algorithm, which was designed exclusively for the fixed operation of FACTS.Only SVC devices were considered for installation (see the literature reference [10]).The results regarding their location and size, as well as the objective function values and their reduction, are listed in Table 7 for this test system.This behavior is different in the other two test systems.However, these locations are expected given the characteristics of the IEEE 85-bus grid, i.e., its many circuit branches, distribution of loads, and large size.ii.The SVC device with the highest capacity was installed at node 61.However, unlike the other two test systems, this one did not report a node with heightened sensitivity, since there was no significant difference between the two largest SVCs installed.iii.The annual energy losses costs are reduced with the installation of SVCs.This reduction is significant for the variable operation.The reductions in the objective functions for the fixed and variable operation of the devices are 26.53% and 30.31%, respectively.This implies that the variable operation saves USD 9251.81/yearmore than the other.

Analysis of Scenario 4 (S4)
This scenario analyzes the impact of installing different FACTS technologies on the test systems, namely the SVC, TCSC, and UPFC devices.This analysis and comparison only consider a variable operation.Table 8 presents the results obtained for the objective function as well as the location and size of the FACTS for all test systems.These results only pertain to the proposed convex model, as it guarantees the global optimum of the problem.energy losses by injecting or absorbing reactive power.SVC devices maintain the same behavior in all test systems, reducing the costs of energy losses to a greater extent.ii.For the IEEE 33-bus test system, the annual operating costs were minimized to values ranging from USD 101,078.70 (for installed UPFC devices) to USD 96,676.76(for installed SVCs).The reductions in the objective function were 10.34%, 12.42%, and 14.24% for the UPFC, TCSC, and SVC devices, respectively.iii.For the IEEE 69-bus test system, the annual equivalent operating costs were reduced to values ranging from USD 105,316.40 (for installed UPFC devices) to USD 100,806.50 (for installed SVCs).The objective function was reduced by 12.03%, 14.07%, and 15.79% with respect to the benchmark case (i.e., the case without FACTS integration) for the UPFC, TCSC, and SVC devices, respectively.iv.For the IEEE 85-bus test system, the annual energy loss costs decreased to values ranging from USD 114,920.20 (for installed UPFC devices) to USD 107,777.80 (for installed SVCs).The annual equivalent operating costs were reduced by 25.69%, 28.25%, and 30.31% with respect to the benchmark case for the UPFC, TCSC, and SVC devices, respectively.

Conclusions and Future Works
This paper addressed the optimal integration of FACTS in electrical distribution grids, considering the minimization of the annual costs related to energy losses and the investments made in installing these devices.This integration was solved using a mixedinteger convex model, which was obtained by transforming the hyperbolic constraints of the MINLP model into second-order conic constraints.The effectiveness of our proposal was evaluated in three IEEE test systems, and it was compared to the BWO and VSA and some GAMS software solvers.The results showed that the proposed convex model found the global optimum of the studied problem.Fixed and variable operations were considered for the FACTS.For the fixed operation, the annual operating costs were reduced by 12.63%, 13.97%, and 26.53% for the IEEE 33-, 69-, and 85-bus test systems, respectively.These results were achieved by the proposed convex model and the BWO algorithm; the other solvers sometimes failed to converge or reached a worse solution.For the variable operation of SVC devices, the annual equivalent operating costs were reduced by 14.24%, 15.79%, and 30.31% for the IEEE 33-, 69-, and 85-bus test systems, respectively.These results were also achieved by the proposed convex model and the BWO algorithm.
The variable operation of SVC devices allowed for a more significant reduction in the annual energy losses costs, although the total size of the SVCs increased by 61.17%, 24.35%, and 27.71% for the IEEE 33-, 69-, and 85-bus test systems, respectively.However, this saved USD 6693.3, 3778.4,and 3778.4/yearmore than the fixed operation.
The installation of SVC devices was the best option to reduce the final objective function value in all test systems.The second best option was TCSC devices.The total size of the SVCs was greater in all test systems than the other FACTS, thus enabling a greater reserve of reactive power as demand grows.

Figure 1 .
Figure 1.Generic branch connection with FACTs devices in an electrical network. slack

Figure 3 .
Figure 3. in the considered active and reactive power curves.

Table 1 .
Electrical parameters of the IEEE 33-test system.

Table 2 .
Electrical parameters of the IEEE 69-test system.

Table 3 .
Electrical parameters of the IEEE 85-bus network.

Table 4 .
Parameter values for the objective functions.

Table 7 .
Cont.The proposed convex model and the BWO algorithm reach the best configuration with regard to the SVC devices, which constitutes the global solution to the problem.All SVC devices in the IEEE 85-bus test system are located in nodes with bifurcations.