Robust Allocation of FACTS Devices in Coordinated Transmission and Generation Expansion Planning considering Renewable Resources and Demand Response Programs

Nowadays, the use of demand response programs (DRPs) in a variety of long-term and short-term planning problems has been explored. In this paper, a generation and transmission expansion planning (GTEP) model along with FACTS device allocation is presented. Furthermore, demand response programs are taken into account for more load ﬂexibility. The proposed model is presented as a multi-objective minimizing problem considering emission, cost, and voltage security index. Furthermore, the conventional Pareto optimization is adopted using fuzzy weighted sum method (FWSM) to achieve a single-objective model. The ﬁnal problem is constrained by equations of alternative current power ﬂow, operation and voltage security limits, planning model of shunt FACTS devices, and operation of the DRP. Adaptive robust optimization (ARO) is used to reach suitable models for the active power of renewable resources and power consumption. As a main search algorithm, a hybrid combination of water cycle algorithm (WCA) and ant lion optimization (ALO) is proposed to ﬁnd the optimum solution with a small standard deviation. The problem is tested on diﬀerent standard IEEE systems, and the results validate the operation and network security improvement due to optimal location of FACTS devices. According to the results, the economic and environmental status of the network has also improved.


Introduction
Among the development and expansion solutions to power systems is generation and transmission expansion planning (GTEP), which helps optimize the sitting and sizing of sources and lines, thus supplying the future energy demands [1].
is can be realized considering various technical, economic, and environmental objectives [2]. To overcome the high cost of operation and considerably high pollution levels and reduce them to the least amounts as much as possible, the electricity industry tries to utilize renewable energy sources (RESs) [3], more specifically wind turbines (WTs) [4]. Since the output power of renewable sources is uncertain and intermittent, the operation will vary every day and hour; hence, there will be unbalanced electricity throughout the system, and the cost of RESs will rise because of unbalanced operation [5]. is discussion is mostly interrelated with the reduced flexibility of the system, which can be realized by adopting flexible sources such as DRP and battery storage banks [6]. e GTEP should be done to meet the high energy demand and electricity consumption surges at peak-load circumstances such as holidays when voltage collapses are highly likely. One solution is to advise and take action beforehand [7]. By utilizing suitable models for GTEP, the abovementioned objectives can be satisfied.
Expansion planning methods for the power grids have widely been introduced in the literature. An adaptive robust model for GTEP is used in [8] in a system that includes wind turbines. Since the system energy demand is uncertain, storage systems, dynamic thermal rating systems, and optimum line switching approaches have been used to make the power system's response flexible. e model is subject to technical and economic constraints on uncertainty parameters. In an attempt to reduce the capital and operating cost, a stochastic co-optimization planning strategy is used [9], in which a linearized model-based AC power flow maintains the bus voltages, reactive power, and active power loss of the system. To compensate for reactive power, some equipment, including power generation sources, transmission lines, static var compensators (SVCs), and capacitor banks, was used in the same study. Moreover, the suggested model considers the loss of load expectation (LOLE). e authors in [10] discuss the GTEP problem subject to the vulnerability of the system in case of terrorist attacks and earthquakes. By doing this, the system will be reliable and safe for future operations.
One generation expansion planning-unit commitment method was presented to upgrade the available models [11]. Investment decisions made on RESs, storage devices, and thermal equipment are optimized. e model also considers real-time flexibility. Expansion planning is also discussed based on a mixed-integer linear robust multiobjective strategy [12]. Taking uncertainty effect into account when dealing with the demand and price of system elements, the researchers also introduce a structure based on information-gap decision theory (IGDT). Reference [13] considers the impact of reliability parameters on expansion planning, in which value-at-risk (VaR) and conditional VaR (CVaR) are also adopted.
is study uses Bender's decomposition to divide the planning problem to one investment problem and two subproblems. e aim is to assess the status of operating cost and reliability index.
yristor-controlled series compensators (TCSCs) are used in [14] in the transmission expansion model to increase the line capacity. In addition, short-circuit level is controlled by using superconducting fault current limiters (SFCLs). TCSCs and SFCLs together help reach a suitable optimal expansion plan. Reference [15] enhances the transmission network operation using a constant series capacitor (CSC) along with DRP for the GTEP. e conventional formulae of GTEP are mixed-integer nonlinear programming (MINLP). But [13] adopts a linearized approximating method (LAM) based on direct current (DC) power flow equations. Linearized AC power flow equations based on LAM have also been used in expansion planning [14]. Nonhybrid evolutionary algorithms (NHEAs) are also used to deal with expansion planning [16].
Various uncertainties include the amount of load, the energy price, RES output power, and impact of power sources planning. As a robust model can find an optimal solution when dealing with uncertainty parameters, it has been discussed in many works. In a min-max framework for the problem, ARO was used for modeling the demand and RES output power uncertainties [17]. e max term specifies the worst scenario, while the min term finds its optimal solution. A max term for the suggested technique is found according to the dual theory [18]. e energy demand and price of energy, and of EVs as uncertainties can be expressed by the boundary uncertainty robust optimization (BURO) [19], where uncertainties are represented by integer values between the upper and lower boundaries of the uncertainty. e robust model based on a linear problem has also been incorporated in [19], while the ARO model works using a nonlinear problem [20]. Table 1 summarizes the studies in this specific area.
Some of the points worth taking into account in the planning of power systems are as follows: (i) Security and pollution level are among the indices that need more focus in the GTEP of power systems. e optimality of an index does not necessarily guarantee the optimal status of another index. For instance, energy cost reduction equals higher power injection by power sources, hence occurrence of overvoltage and escalated loss. us, various effective indices need to be taken into account at the same time. (ii) e GTEP is MINLP, and a LAM is used to solve it.
Yet, LAM leads to considerable calculation errors. is is more sensible in cases where the GTEP is based on DC power flow. Moreover, this model cannot appropriately analyze voltage security, and implementation of reactive power planning is not possible. According to the linearized AC power flow, the calculation errors made in LAM for modeling the GTEP cannot be neglected [4]. Besides, the problem is solved using NHEAs, in which unique response conditions are not provided. To overcome this challenge, hybrid evolutionary algorithms (HEAs) are recommended, where decision variables are updated and unique approximate response conditions are met. (iii) In general, stochastic and probabilistic programming has been used to model uncertainties. But to ensure the optimal solution is obtained, robust planning is a suitable method. is is because it determines an optimal point that is robust to uncertainties. In this topic, various research works have been presented in the literature. However, they generally used a robust model such as ARO and BURO that can only be presented on linear problems. Nevertheless, based on the second research gap, using the LAM for GTEP bears a significant computational error. While ARO is applied to a nonlinear problem in [20], it has complexities, including dual gap and complementarity (equilibrium) constraints. To compensate for this, a combination of ARO and HEA is expected to provide a suitable capability.

2
International Transactions on Electrical Energy Systems is paper attempts to ll the abovementioned gaps, Figure 1. e placement of shunt FACTs devices when dealing with the GTEP problem is also addressed to predict the optimal operation situation, voltage security, and economic and environmental indices.
is strategy can be structured like three-objective optimization minimizing the cost of planning (which is the total capital cost of generation units and shunts FACTS devices, in addition to transmission lines together with generation units' operating cost), the pollution level released by generation units, and the index L max (concerning voltage security). e Pareto optimization method based on the sum of weighted functions method is employed to structure the objective function. e function relies on operation limits, equations of AC power ow, security of voltage, DRP operation model, operation of shunt FACTS devices, and planning model. e best solution compromised among objective functions is speci ed using the fuzzy decision method. Moreover, the ARO is adopted to model the energy demand and active power of RESs suitably. e HEA based on hybrid WCA-ALO solves the problem and nds an optimal solution. Notable innovations presented in this paper include the following: (i) Presenting a GTEP model to plan reactive power sources to supply reliable, clean energy.
(ii) Modeling operation, security, economic status, and environmental status indices at the same time. (iii) Robustly modeling the demand and output power of renewables using the MINLP model of the GTEP by adopting the HEA-based ARO. (iv) Adopting the hybrid WCA-ALO to nd the optimum solution in the nal response that has small standard deviation in special approximate response situations. e rest of this paper is organized as follows. e GTEP formulation will be described in Section 2. e uncertainty modeling and the problem resolving procedure by the HEAbased ARO are presented in Section 3. Section 4 illustrates the obtained results. Finally, Section 5 summarizes the conclusion.

Multi-Objective GTEP Strategy
e GTEP model in this section explains how to improve technical and economic indices for loads in the transmission network. A three-objective optimization problem is formulated for achieving the objectives of minimizing planning costs, greenhouse gases, and L max index of the voltage t∈T P DR dt y � 0, ∀d, y, subject to Equation (1) specifies the objective functions of the suggested method. e first part of (1) includes minimizing planning cost of transmission lines in addition to the active/ reactive resources. Costs for transmission lines, parallel FACTS, and GUs in the first and third terms are considered [1]. GUs' operating costs (fuel) are applied in the fourth term [4]. e transmission lines among n-b and b-n buses are similar, so a 1/2 factor is observed in the first section [1]. Equation (1) also addresses the environmental emission minimization (EM) resulting from the fossil fuel usage of nonrenewable GUs. us, it is not far-fetched that ecofriendly resources would take a larger portion of generation [21]. Equation (1), in its third part, tries to improve the situation of the network voltage security [22]. Here, L max is used for grid security analysis, a numerical value that occurs between zero and one. No-load operation is indicated by a value of zero, and the voltage drop is indicated by a value of one. As a result, the third part of (1) requires the minimization of index L max in order to ensure optimal voltage security.
e Pareto optimization method is used to formulate the objective function according to the summation of weighted functions [23]. Hence, various points are achieved for the values of ω VSI, ω EM, ω Cost and parameters, for VSI, EM, and cost functions. e Pareto front of the suggested scheme is plotted in 3-dimensional coordination. Furthermore, the function has the lowest and highest values that could be computed through three cases provided by ω Cost � 1, ω EM � 1, and ω VSI � 1. To obtain the optimal compromise answer among the aforementioned functions, the fuzzy decision-making strategy can be applied. Detailed information is described as a fuzzy decision algorithm (Algorithm 1). e planning model, including transmission lines, GUs, and parallel FACTSs, is considered in (2)-(5). Equations (2)-(4) show the investment budget equations for the GUs construction, parallel FACTSs, and transmission line, International Transactions on Electrical Energy Systems respectively. Equation (5) indicates a logical limit which pursues the structure of transmission lines between n-b and b-n buses since both lines present identical transmission lines. e grid constraints of AC-PF are provided in (6)-(9) [1,2], representing the balance between active and reactive power in the buses, along with availability of active/reactive power of transmission lines, respectively. Considering (8) and (9), when the binary variable of a transmission line presence (y TL ) is 1, it means that the line is joined to the grid; if not, it should be considered as disconnected.
In constraints (10)- (13), the planning and operation modeling of GUs and parallel FACTSs are given based on [1,24], respectively. Equations (10) and (11) denote the GU generator ability curve that shows the real power production limitation (with positive value) in addition to GU controllable reactive power (which can have inductive and capacitive operation modes). It is worthy of note that (10) could be applied for nonrenewable GUs since renewable ones, like wind farms, because of small emissions and operation costs, typically (12), supply real power of the utmost capacity commensurate with weather condition to the grid [20]. For parallel FACTS operation or planning modeling, (13) shows the controllable reactive power limitation of the parallel FACTSs in both modes of inductive and capacitive operation. When the binary variables y GY and y PaF are equal to 1, GUs or parallel FACTSs will be connected to the network as shown in (10)-(13); if not, they will not be considered in the grid. Considering (14) and (15), the responsive load operation modeling is indicated [20]. Such loads are based on the incentive-based demand response programs (IDRP) [25]. As a result, these loads decrease their usage in peak times when the electricity price is high. In peak times, because of great usage of grid, GUs that have higher fuel costs can be estimated for playing part in power provision; hence, the local marginal price (LMP) is higher [21]. In addition, in off-peak times, cheap GUs provide electricity, consequently, the electricity costs get lower, and also the responsive loads consume lower electricity in peak times. Hence, based on this rule, (14) shows the load controllable power limitation of DRP plan; (15) states that the whole responsive load's decreased electricity in peak times can be obtained by grid [20].
Moreover, in (16) and (17) [21] and (18) and (19) [22], the transmission network technical constraints, namely, security and operation, are shown, respectively. e constraint related to the network operation, as shown in (16) and (17), respectively, refers to the apparent power limitation that flows in transmission lines and the boundary of bus voltage [20]. In (18) and (19), the network voltage security model is shown; it computes the L max real and imaginary components in bus bars (PQ buses), respectively. It is worth mentioning that voltage-controlled bus (PV bus) has an assured distance from voltage drop point because of its capability of setting voltage. erefore, L max can be computed for PQ buses.

Uncertainty Modeling.
e suggested design in (1)-(19) includes uncertain variables like real and reactive loads, P D and Q D , and maximum WF real power, GU/P. In this study, in order to achieve a robust design for GTEP against the intermittent variables, the ARO is used for modeling the uncertainty noted in the suggested design.
is model needs a scenario known as "worst-case scenario" [26]. e uncertainties in this scenario create the worst conditions for the problem, i.e., reducing the feasibility region. erefore, the ARO is a suitable technique for modeling uncertainty and requires low computational efforts compared to probabilistic or stochastic programming [27]. Furthermore, WF would be uncertain because of the intermittent power generation. It is also predicted to be a proper resource to improve the operation, environmental emission, and voltage security indices of the transmission networks. Hence, robust modeling of its uncertainties obtains the robust ability of the WF to reach a network that consists of a proper technical situation. To put it differently, the response to the question "Is the WF able to improve power system technical situation in the worst-case scenario?" is given.
In the following, the matrix related to uncertainty parameters needed by ARO (ϑ) is introduced as (20), initially. According to this equation, the number of the rows should be equal to operational time number (n t ), and column number (n c ) is equal to n c � n y (2n d + n u ), where n y , n u , and n d denote the numbers of planning years, WFs, and load blocks, respectively. In the ARO, the uncertainty variable (ϑ) is the true value of uncertainties. Hence, the uncertainty matrix variables for the matrix ϑ are represented as (21). Afterward, according to [26], we need to determine the uncertainty set in the ARO method for the uncertainty variables, indicated for the t-th row of the matrix ϑ as (22).
6 International Transactions on Electrical Energy Systems In (22), ϑ indicates the uncertainty deviation, ϑ is the predicted value of uncertainty, and Δϑ denotes the uncertainty budget that is between zero and one. e zero value denotes the deterministic state of uncertainty ϑ, while the one value shows the worst-case scenario obtained by all the matrix elements (ϑ) [18]. Furthermore, the term indicates the changing amplitude of an uncertainty matrix element. Eventually, the uncertainties' set will be equal to unions Λ ϑ t as given in (50).
e worst-case scenario in the ARO method and its optimal value are obtained simultaneously [28]. In this method, it was observed that the worst state of the objective function belongs to the worst-case scenario compared with other scenarios. For instance, in minimizing the objective function of the main problem, a higher value would be achieved for the worst-case scenario compared with other scenarios. To mathematically model the problem, the term max is applied to the set of uncertainties in the main problem's objective function; hence, the true value of the uncertainty variable (ϑ) could be obtained. us, the uncertainty parameter (ϑ) could be replaced by the variable ϑ. Moreover, limitations of the robust problem based on the ARO would be identical to the constraints of the main problem and (23) [17]. As a result, the ARO-based robust model is described for the problem (1)- (19), as follows: constrained by equations (2) equations (6) − (19), where the uncertainty variables of the matrix ϑ are replaced by uncertainty parameters of (20), (26) equation (23).

International Transactions on Electrical Energy Systems
It is worthy of mention that based on problem (1)- (19), the binary variables y TL , y PaF , and y GU refer to the system construction status that is typically uncertainty independent and is frequent for the "wait and see" variable [17]. In contrast, "here and now" variables of the problem depend on uncertainty. erefore, the term ω Cost (1/2. n,b∈N C TL nb y TL nb + p∈PaF C PaF p y PaF (1) and (2)- (5) is independent of uncertainty variables. Hence, the robust model's objective function (24) involves two sections: e first one includes merely min term since this part is independent of ϑ. e second part of (25), because of its dependency on uncertainty variables, contains min-max term. is part of the problem contains variable z that presents the vector of the continuous variables (1)- (19), namely, L max real and imaginary components, voltage magnitude and angle of the network buses, and active/reactive power of different c. e term min in the second part of (24) is implemented for achieving the optimum value of z according to the main problem, (1)- (19), and the term max in this section achieves the uncertainty parameter values regarding the worst-case scenario. e equations of the problem in the robust model, (25)/(26), are similar to constraints (2)-(5)/(6)- (19); however, the variable ϑ is replaced with the parameter ϑ in (26). e uncertainty set is the other constraint in ARO formulation (27), which refers to the permitted changing values of the variable ϑ.

Problem Solving according to the HEA.
e suggested problem indicated in (24)-(27) is an MINLP problem. e dual model of this problem should be extracted for solving the problem using mathematical solvers until the min--max term changes to a max or min term. Nevertheless, the main formulation is hard to achieve and complicates the solution procedure because finding the dual model of an MINLP problem requires observing the duality gap equations and the equilibrium (complementarity) constraints [20]. erefore, a linear anticipation formulation is initially achieved for the robust model of the optimization problem in most research works. Nevertheless, according to [4], this technique is accompanied by some computational errors, and computational errors for power losses in the linear AC-OPF problem are greater than ten percent in comparison with the nonlinear modeling of the distribution networks. erefore, the solution achieved using this method consists of a small coefficient of reliability. To deal with the problem, this paper applies the HEA according to the hybrid algorithm based on WCA [29] and ALO [28] to obtain the optimum answer with a special response condition. e solver finds the optimum answer having the smallest error with the ultimate response because decision variables can be updated by both WCA and ALO procedures. Determination of the optimal compromised solution is given in subsection 4.1.b. e variables are divided into two groups for solving the suggested problem by the HEA algorithm: dependent and decision variables. e latter for the suggested scheme consists of y GU , y PaF , y TL , and P GU for nonrenewable GUs, Q GU , P DR , Q PaF , and ϑ, and the values that can be designated by the hybrid ALO + WCA algorithm are commensurate with the following equations: y PaF P GU uty ∈ Equation (10), ∀u ∈ GU − WF, t, y, Q PaF pty ∈ Equation (13) ∀p, t, y, e values of dependent variables, including L r , L im , P TL , Q TL , V, α, and P GU for renewable GUs, are defined according to (6)-(9), (12), and (18)- (19) and the amount of the decision variables. In the current section, firstly, P GU for renewable GUs is computed according to (12). Afterward, the variables P TL , Q TL , V, and α are achieved through the power flow constraints (6)- (9). e Newton-Raphson method is used in this study for solving the power flow equations. Eventually, (18) and (19) are applied to, respectively, compute L r and L im . In addition, a penalty function is applied for anticipating GU investment budget limitations, parallel FACTSs, and transmission lines, (2)-(4); IDRP operation constraint, (15); and system operation constraint, (16) and (17). In this method, the penalty function of a = b is denoted by λ.(b − a) and that of a ≤ b is described as μ×max(0, a − b) [30]. μ ≥ 0 and λ ∈ (− ∞, +∞) show Lagrange coefficients, where determinated and updated by ALO + WCA. Afterward, the fitness function (FF) equals the objective function aggregation of the original problem according to this technique and is identical to (36), (24), and the penalty function. Finally, the optimal amount of FF is reached based on decision values and dependent variables provided in the following equation: It is important to mention that as long as the convergence is provided, solutions of problem will be obtained; if not, according to the optimal FF value along with ALO + WCA solver rule, the decision variables will be updated. Here, convergence conditions are affordable when iteration number iter max achieves its maximum. Hence, the process of HEA problem solving can be as follows: (i) Generating initial population (initialization step): the hybrid ALO + WCA algorithm randomly produces N (size of the population) various values for Lagrange multipliers and decision variables commensurate with their allowable range.

Numerical Results
Here, the suggested scheme is implemented on modi ed IEEE 6 and 89-bus grid. MATLAB software is used for simulations. Here, the uncertainty deviation ( ϑ) is assumed to be the same amount for all intermittent parameters, which is considered as r.Δϑ, where r shows the uncertainty level. Next, the abilities of the noted strategies are implemented based on the numerical result achieved from various cases. Figure 2 depicts the single-line illustration of this grid [1]. e initial voltage and power values are 230 kV and 100 MVA, respectively. e allowable voltage range is [0.5-1.05] p.u. According to the grid assumptions, it is possible to construct a line parallel to transmission lines. In [20], all current and candidate lines' speci cations are provided, including construction cost, capacity, resistance, and reactance. It is possible to construct four renewable GUs of WF type and 13 nonrenewable GUs along with the presented GUs in the aforementioned grid. A1-A5 and B1-B8 show nonrenewable GUs and are, respectively, equipped in buses 1, 1, 2, 3, 6, 1, 2, 3, 6, 1, 1, 2, and 3. Moreover, W1 to W4 indicate wind farms that are, respectively, equipped in buses 6, 1, 1, and 2. All the existing and candidate GU speci cations such as b factor in cost of the operation, construction cost, and emission coe cients are represented in [20]. In addition, [31] reported the maximum active power generation. 30% of GU P is presumed for the least and most GUs' reactive power, and the minimum output of active power is supposed to be 0. Furthermore, in the operation cost of GUs, the coe cients a and c are, respectively, presumed to be 10×b and 10 − 4 ×b. e daily real power pro le of WFs can be obtained by multiplying their capacity by their daily power rate curve; this curve is represented in [1]. Furthermore, the consumption loads in buses 3-5 are 40%, 30%, and 30% of the total load of the grid. e planning horizon is considered for six years. For this horizon, the peak load of the network is 29-34 MW, with a coincidence coefcient of 70% and a power factor of 0.9. e hourly network load pro le for a day is also equal to the load factor curve and generation of the annual peak load shown in [1]. en, three parallel SVCs with the capacity of 5 MVA could be equipped in the network, which imposes $10/kVAr/year construction cost [24]. Moreover, the investment budget of parallel FACTC, GUs, and transmission lines can be, respectively, presumed to be $10 M, $100 M, and $20 M. Users would also have 30% contribution to DR program in this method. Table 2 presents the Pareto front of the suggested deterministic GTEP design using the fuzzy decision-making technique and the model denoted in (1)- (19). Table 2 lists the results of 0, 0.33, 0.5, and 1 of weighting coe cients, i.e., ωCost, ωEM, and ωVSI. Hence, based on  Table 2, the VSI, EM, and objective cost functions are not in the same directions. For instance, a decrease in VSI is proportional to a rise in planning cost since several GUs in the system are needed to enhance the voltage security of the network. It is equal to the increase in costs of planning and operation. Table 3 reports the optimal solution compromises between the aforementioned objective functions achieved with ALO + WCA, ALO, WCA, crow search algorithm (CSA) [32], and grey wolf optimization (GWO) [33]. According to Table 3, the maximum convergence iteration number (iter max ) and the population size (N) for the aforementioned algorithms are considered 5000 and 80, respectively. e quality of the responses for the mentioned algorithms is de ned by computing statistical indices like standard deviation (STD). In this regard, the solvers repeated the problem solution for 30 iterations. Table 3 reports the results for an uncertainty budget of 1 and various uncertainty levels. Eventually, the ALO + WCA algorithm gives the optimal solution compared with NHEAs, namely, CSA, GWO, WCA, and ALO. e optimal solution consists of the least VSI, EM, and cost values. For the computational time (CT) of 6.2 minutes in deterministic model (Δϑ = 1, r = 0), 711 was recorded for the least convergence iteration (CI) number. Additionally, compared with other algorithms, it has the minimum STD for this model (0.96%), which shows its low dispersion of the suggested GTEP problem ultimate solution.

Determination of the Optimal Compromise Solution.
e VSI, EM, and cost values will be increased by raising (reducing) the load (WF production power) in robust modeling by increment in r of the worst scenario considering the scenario utilized in the deterministic model. Furthermore, because of the decrease in the feasibility region of the suggested GTEP, the values of CT and CI increase in this condition. While SD value is almost the same with case of Δϑ = 1, r = 0, this is not true in NHEAs. Accordingly, the solution deviation reached by ALO + WCA could be indented from the feasibility region of the problem. Hence, the ALO + WCA solver approximates unique solution conditions, which aligns with the paper's second contribution noted in Section 1. In addition, the distance between the functions and the corresponding least values is lower in the optimal solution. erefore, in the compromising point, the cost function obtains 36% ((95.9-91.8)/(103.2-91.8)) far from its lowest value, that is, $91.8 M, in the deterministic model. In case of Δϑ = 1, r = 0, the values achieved are about 18.5% and 16% for VSI and EM functions, respectively.

Evaluating the Results of Network Expansion
Planning. e results of the mentioned grid expansion planning related to the suggested robust and deterministic GTEP strategies by an uncertainty budget of 1 are reported in Table 4. Based on  Table 4, three lines, T2, T6, and T7, are connected to the grid of Figure 2 for the two mentioned models. ese lines are placed between buses 2 and 3, 5 and 6, and 3 and 6, respectively. In addition, for minimizing the environmental emission, planning cost, and L max of security of voltage, every wind farm, W1-W4, in addition to units A1, A4, B5, B7, and B8, should be linked to the grid for various uncertainty levels provided in Table 4. Equation (1) presents the selection corresponding to the objective function. en, at the bus bars of the network, three SVCs are installed at buses 3-5, for all uncertainty values, because installing SVCs in these buses helps to reduce the consumers' reactive consumption, and the voltage safety margin because of high X/R ratio of transmission grid is notably associated with the reactive power consumption. Hence, it is predicted that SCVs have amending effect on the security of voltage. ey can also decrease losses of energy via line flow reduction and finally reduction in operation costs and emissions produced by GUs.
It is noteworthy that transmission line, SVC, and GU planning results are the same for both fixed and robust models (in the budget of uncertainty of one and different uncertainty levels). is is due to the construction status variable of the mentioned elements according to model (24)- (27), which does not depend on the decision variables; therefore, the same planning results are obtained for different levels of uncertainty. However, increasing the level of uncertainty compared to the definite model (Δϑ �1, r � 0) results in GU operation cost increment and ultimately increases the total planning cost of the network. According to Table 4, the amount of active and reactive load in the worst scenario rises due to the increase in the uncertainty value e best Pareto front compromise answer; Optimum solution of Pareto and the decision-maker favors; Step 1: Fuzzy membership function calculation Values′ computation of the linear fuzzy membership function (F i ) of any Pareto optimum front member: it is assumed as 0; end end Step 2: Obtain φm ∈ 1, 2, ..., n m , where m is a member of the Pareto front, n m is total members.
Step 3: e best compromise answer e best compromise answer is found via calculating max m φ m ALGORITHM 1: Pseudocode of the fuzzy decision support system. compared to the definite model-based scenario, but the level of active power generated by WFs decreases. erefore, the demand for energy from GUs has increased in these conditions, leading to increment in costs of GU operation and total planning of the network. In addition, it can be observed from Table 4 that an increase in the level of uncertainty to 0.21 results in changes in the active and reactive loads and WF real power output in the worst-case scenario. However, the obtained results for r > 0.21 are the same as the results for r � 0.21. is value of r is known as the maximum level of uncertainty (r max ) in this paper. is also means that the developed network is able to cover up to 21% of the prediction error of load uncertainties and active power generation of WFs. Suppose the prediction error of the mentioned uncertainties is more than 21%. In that case, it is expected that the planning results obtained in Table 4 for the mentioned network could not be successful in operation.

Examining the Abilities of the Suggested Method.
To test the suggested method, five case studies can be proposed as follows: (i) Analysis of load flow.
(v) GTEP Method by Both DRP and SVC Attribution. Table 5 shows the environmental status, economic status, voltage security, and operation results for the mentioned grid with an uncertainty budget of 1 and various uncertainty values. e operation indices include total energy losses (TEL) over the planning period and maximum voltage deviation (MVD). According to this table, the cost is minimal in Case I, but it merely indicates the grid power flow results. erefore, the cost function equals the AE1-AE3 and BE1 operation costs. Moreover, the grid and resources' technical limits are not taken into account in Case I; hence, it cannot be guaranteed that the mentioned sources will ensure consumption growth in the upcoming years. e results in Case II can be examined to prove this statement. In Case II, by considering the technical parameters, namely, power source, economic status, environmental status, voltage security, and operation limitations, it was seen that various resources and lines should be applied to the grid for fulfilling the consumption increase in the years ahead. With regard to SVC attribution and conducting the DR program in GTEP (Case V), the cost of planning is decreased according to Table 5by almost 16.8% ((115.3-95.9)/115.3) compared with Case II in the deterministic model (r = 0). In addition, Case V could decrease environmental emissions and voltage security by approximately 65.1% and 51.6%, respectively, compared with Case I in r = 0. In Case I, the highest TEL and MVD occur. Nevertheless, these parameters are improved in Case V by approximately 17.6% and 30.5%, in comparison with Case I in r = 0. Based on Table 5, the DR program strategy owns a considerable capability to decrease VSI, EM, and cost function in comparison with Case II. Furthermore, the DR program can also considerably influence the enhancement of MVD and TEL more than Case II. Hence, managing or storing the consumption would improve the environmental and economic conditions, grid voltage security, and operation indices. e SVC assignment can also potentially improve the security of voltage compared with Case II, and it has a significant impact on the improvement of the MVD status; however, the decrease of planning costs, TEL, and emissions is not considerable since these functions are less reactive power dependent. Nevertheless, the system voltage and VSI highly depend on reactive power because of the transmission grid high X/R ratio [22]. e results of the mentioned indicators for the maximum level of uncertainty and for different study cases are examined in Table 5. According to the table, an increase in the level of uncertainty in the worst-case scenario, in comparison with the strategy definite model-based scenario (r = 0), leads to an increase in cost, EM, VSI, TEL, and MVD parameters and finally results in the deterioration of economic status, environmental status, security, and network performance indicators of the model. ese conditions are because of an increment of load and decrease of active generation capacity of WFs in the worst-case scenario, according to Table 4. In addition, it should be noted that the proposed design for V is able to cover the high prediction error of wind and load and production capacity compared to other studies.
erefore, the suggested design can obtain more optimum values for the various mentioned indicators compared to cases I-IV, as well as in the case of forecast error or uncertainty of load and production capacity of WFs compared to cases. Other studies are more robust.

IEEE 89-Bus Grid.
is grid has 35 bus bars, 12 nonrenewable GUs, and 210 transmission lines; their characteristics are reported in [34]. Maximum and minimum value of voltage are considered to be 1.05 and 0.95 p.u., respectively [35]. Eventually, other characteristics like power factor of network, daily power rate curve of WF, CP, and other components, according to the information provided in Section 4.1, are assumed. Table 6 and Figure 3 represent results of this section. Table 6 indicates the planning results of understudied grid for various uncertainty levels in Δϑ = 1. According to it, 5 WFs, 4 nonrenewable GUs, 16 transmission lines, and 12 STATCOMSs ought to be considered in the planning horizon for fulfilling the consumption increase in years ahead. e placement of STATCOMs and WFs is represented in the prior paragraph; however, the constructed transmission line and nonrenewable GU reports are presented in Table 6, and their information can be found in [17]. e planning consists of costs of construction of transmission lines, GUs, and STATCOMs, being $197 M, $1,045 M, and $61 M, which are fixed for various uncertainty levels (r) since the binary variables corresponding to the constructed state of these elements are not dependent on uncertainty parameters. Nevertheless, the expected operation cost of GUs rises with the increase in r because the value of load (WFs active power) increases (reduces) compared to the deterministic model (r = 0), as shown in Table 4. On the other hand, according to the results of Table 6, the proposed design has a maximum uncertainty of 0.17; hence, this scheme in the network can have robust operation against the error of load uncertainty prediction and an active power of WFs up to 17%. Based on Table 6, the ALO + WCS algorithm can find the optimum solution of GTEP in the IEEE 89-bus network in convergence iteration (computational time) less than 1860 (18 minutes), with a final solution deviation similar to the results of the previous section, around 0.96%. In Figure 3, the environmental emission, security, and operation (TEL and MVD) indices in Case I and Case V are reported for uncertainty budget one and different uncertainty levels. According to this figure, the EM and VSI functions are improved by the suggested GTEP method (Case V) considering the power flow analysis (Case I) in the different uncertainty levels by about 42.7% (12526726-71812)/125267) and 51.7%, respectively, in the deterministic model (r � 0). Operation indices like TEL and MDV are likewise enhanced by almost 33.4% and 29.8% in the aforementioned situations (r � 0), respectively.  According to Figure 3, as uncertainty value increases, theses indices increase compared to the case r 0 because it raises (decreases) the load (active power of WFs) compared to the de nite model according to Table 4.

Conclusion
tA GTEP model with parallel FACTS attribution is presented here for transmission grid and responsive loads. e proposed problem is three-objective optimization whose functions are minimizing planning costs, reducing nonrenewable GU emissions, and minimizing the L max of voltage security; model of parallel FACTS planning, security of voltage, operation, AC power ow relations, and operation model of DRP are constraints. e formulation is done according to the Pareto optimization method, which is based on the sum of the weighted functions method, to obtain a single-objective model for the suggested design.
Moreover, the ARO is employed for modeling load uncertainties and real power from renewable resources. Next, to obtain the optimum solution, the hybrid ALO + WCA algorithm is applied. Eventually, based on presented results, this algorithm is able to gain the best compromise answer rather than NHEAs with the least computational time and convergence iteration number and with a standard deviation of approximately 0.96%. In other words, it is able to achieve unique solution's approximate conditions. With optimum placement of transmission lines, GUs, and paralleled FACTSs, the suggested GTEP strategy improved environmental status, security of voltage, losses of energy, and pro les of voltage by more than 42%, 51%, 18%, and 30%, respectively, based on optimum planning paid cost in comparison with power ow analysis. In addition, the uncertainty level increase of the solid model increases the load and decreases the active generation power of WFs in the worst-case scenario in comparison with the   T11, T17, T18, T31, T41, T55, T84, T98, T122, T139, T156,  T169, T178, T179, T200,  is also increases voltage drop, energy losses, operation costs, environmental emissions, and voltage safety index. But considering these cases, the proposed plan is able to find the optimum solution in the worst-case scenario so that it is able to cover the prediction error of uncertainty parameters between 17% and 21% for a specific planning result and obtain a robust solution.

Indices and Sets b:
Bus index u, p, d, n, t, y: Generation Voltage magnitude's lower and upper limits (p.u.) β: Participation rate of load block in the DRP χ, δ: Coefficients of emission equation for GU, respectively, in ton/MW and ton ω Cost , ω EM , ω VSI : Weighed coefficients in objective function.

Data Availability
Data sharing not applicable. No new data were created or analyzed in this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.