Multi‐objective economic operation of modern power system considering weather variability using adaptive cuckoo search algorithm

Introduction The energy requirement has been increased significantly across the world. For an instant, as per ‘executive summary on power sector-March 2019 by Central Electricity Authority (CEA)’, energy requirement across India is projected as 2.11% raise for March 2019 (108,665 MU) as compared with March 2018 (106,420 MU). Under this scenario, integration of renewable energy (RE) in power system becomes one of the counter measures Abstract


Page 2 of 29
Kavuturu and Narasimham Journal of Electrical Systems and Inf Technol (2020) 7:11 for both increasing electricity demand as well as global warming across the world [1,2]. According to the International Renewable Energy Agency (IRENA) statics report released in March 2019, the total RE integration in 2018 has been reached to 2364.4 GW including off-grid across the world [3]. Among various types of RE sources [4], solar photovoltaic (PV) has become the most attractive technology even at the end-user side. As compared to the fossil fuel-based conventional energy (CE), the RE sources have various advantages like (1) reduced active power losses, (2) improved voltage profile, (3) increased overall energy efficiency, (4) congestion relief across the system elements, (5) potential increase in service quality to the end-customers, etc., however, their intermittency nature of RE sources needs to be addressed potentially for avoiding the possible operational and controlling issues [5]. On the other side, the non-expanding transmission systems and various uncertainties have been resulted various power system blackouts across the world [6,7]. In order to overcome these potential problems, the flexible ac transmission system (FACTS) devices in power systems have been introduced and become one of the promising solutions for stability enhancement, loadability enhancement, loss minimization and economic operation, etc., as well as to prevent various likelihood uncertainties in real-time operation of many power systems across the world [8]. A comprehensive survey on different types of FACTS devices proposed for stability enhancement can be found in [9] and for transmission system security management in deregulation era can be found in [10,11]. In flip, the raising scenario of global warming becomes one of the major concerns in planning studies of almost all the sectors and particularly reducing the carbon emission in electric generation become more important in power system operation. According to 'climate change scenarios for India, ' the surface temperature is projected to increase as much as 3-4 °C toward the end of the twentyfirst century. Similarly, as per India Meteorological Department (IMD) reports, January 2019, annual mean temperature during 1901-2018 showed an increasing trend of 0.6 °C per 100 years. Under this scenario, it is essential to consider the impact of ambient temperature on power system operating condition, which has been neglected almost in all the works in the literature. Hence, consideration of RES uncertainty, change of ambient temperature and FACTS controls in economic operation of modern power systems cannot be negligible. Economic load dispatch (ELD) problem is one of the fundamental optimization problems toward minimum generation cost in power system operation. Traditionally, optimal power flow (OPF) is in practice for determining the solution of ELD problem without compromising in various operational constraints [12]. Apart from the conventional approaches like lambda-iterative method, the base point and participation factors method and the gradient methods, many authors have been focused on various meta-heuristic algorithm-based approaches in recent times for solving the OPF considering FACTS devices and/or uncertainty in power generation by RE generation sources. In [13], a hybrid algorithm moth swarm algorithm (MSA) and gravitational search algorithm (MSA-GSA) is proposed for OPF considering wind power. In [14], a modified version of the moth swarm algorithm (MMSA)-based OPF is proposed for handling the stochastic nature of wind generation. In [15], OPF based problem is solved by using a modified hybrid particle swarm optimization and gravitational search algorithm (PSO-GSA) with chaotic maps approach considering stochastic nature of wind generation and FACTS devices such as thyristor-controlled series capacitor (TCSC) and thyristor-controlled phase shifter (TCPS). In [16], distributed-FACTS (D-FACTS) and FACTS devices are proposed for attaining the substantial cost savings by assuring the congestion relief in transmission system as well as minimum RE power curtailment in the competitive electricity market environment. In [17], a novel tree-seed algorithm (TSA) is proposed for solving the OPF problem considering various continuous and discrete control variables. In [18], modified JAYA (MJAYA) algorithm is proposed for multi-objective OPF problem considering operating cost, emission, real power loss and voltage profile in the presence of different RE sources in the network. In [19], a hybrid algorithm using dragonfly algorithm (DA) and aging particle swarm optimization (APSO) is proposed for OPF considering wind power. In [20], novel improved social spider optimization (NISSO) algorithm is proposed for solving the OPF with various single-objective functions like operating cost, real power loss, total emission, voltage deviation and stability index. A modified sinecosine algorithm (MSCA) with Levy flights is proposed for solving single-objective OPF problem [21]. Considering uncertainties due to wind power and various factors in the grid, evidence theory (ET) and extended affine arithmetic (AA) are employed for OPF problem [22]. In [23], a three-stage real-time OPF (RTOPF) strategy is proposed to handle the error in forecasting the RE and demand loads. In [24], a modified moth swarm algorithm with an arithmetic crossover (MSA-AC) is proposed for solving the constrained optimization and OPF problem. In [25], optimal location of distribution generation (DG) is handled simultaneously while solving the conventional OPF problem over the 24 h of the day using sun flower optimization (SFO) algorithm. In [26], particle swarm optimization with an aging leader and challengers algorithm (ALC-PSO) is proposed for solving the different single-objective OPF problems considering thyristor-controlled series capacitor (TCSC) and thyristor-controlled phase shifter (TCPS). In [27], minimization of congestion cost under different contingency conditions is solved using symbiotic organisms search (SOS) algorithm. In [28], various single and multi-objective functions with different combinations of single-objective functions are solved using backtracking search optimization algorithm (BSA). In [29], a tri-population based DE-PSO-DE optimization algorithm using differential evolution (DE) and particle swarm optimization (PSO) is proposed for solving the economic load dispatch problem considering operating cost, real power loss, ramp rate and spinning reserve constraints. In [30], power system security enhancement via fuzzy harmony search algorithm (FHSA)-based OPF is proposed considering TCSC at pre-determined optimal location, which determined using line overload sensitivity index (LOSI). In [31], symbiotic organisms search (SOS) algorithm is proposed for solving the OPF considering thyristor-controlled series capacitor (TCSC) and thyristor-controlled phase shifter (TCPS). The objective functions include operating cost, valve-point loading effect, real power loss and combined economic and environmental cost. In [32], the OPF problem is solved using modified BAT optimization problem in the presence of interline power flow controller (IPFC) considering dynamic loading conditions. In [33], the design aspects of AC-DC smart micro-grid considering availability and cost of equipments is presented. The authors proposed the hybrid micro-grid with wind, solar and battery bank systems and optimized the economics of modern grid using multi-objective particle swarm optimization (MOPSO). In [34], the power quality issues with variable RESs have been mitigated using FACTS devices considering technical and economical aspects of power system.
From the above reviewed literature, it can be said that the OPF framework is used to solve a wide range of single-and/or multi-objectives (namely cost, voltage profile, voltage stability, loadability, security or congestion, multi-fuels, valve-point loading and effect emission) considering different types of FACTS devices and uncertainty in RES generation, load, contingencies and different operational equal and unequal constraints. Also, the OPF problem is solved with different heuristic algorithms effectively than conventional approaches like nonlinear quadratic programming [35], Newton [36] and interior-point [37] methods.
In recent times, consideration of environmental aspects in the economic operations and planning studies of power system engineering is getting attraction considerably. In [38,39], the impact of weather changes on distribution system load profile and losses are addressed. In [40], ambient temperature effect on load and consequently distribution system performance is analyzed and improved via optimal allocation of capacitor banks. In [41], weather changes on transmission system parameters are highlighted through transient stability analysis. In [42,43], the conventional NR power flow method is reformulated considering the impact of temperature on transmission system parameters. In [44], output characteristics of solar PV panels are analyzed under different weather conditions. In [45], the experimental investigations in Brighton in the southeast of the UK revealed that the solar PV system efficiency is reduced by 60% due to dusty or polluted climate. From these works, it can be concluded that the consideration of weather changes is also necessary in planning studies and for effective operation and control of power systems in the current as well as future scenario of increasing global warming and changing weather conditions [46]. Most importantly, it is also observed that the impact of weather changes on transmission system parameters, load profile, solar photovoltaic generation and consequently power system techno-economic operating condition is not addressed w.r.t. FACTS controls. Since OPF problem is nonlinear and complex in nature, it is easy to solve using heuristic approaches than non-heuristic approaches. But the heuristic approaches used in literature may not escape local optima while solving OPF with multi-objectives and many search variables. Hence, there is a need to reevaluate the performance of heuristic approaches by modifying the search process between local optima and global optima.

Methods of work
In light of the above introduction, ambient temperature effect-based load growth, photovoltaic (PV) generation and transmission system parameters are modeled and considered in the OPF framework in this paper. Also, optimal unified power flow controller (OUPFC) is allocated in the system and tuned its controlling parameters optimally for minimizing the multi-objective function with total generation cost of conventional energy (CE), real power losses, voltage deviation and proximity of voltage instability under different weather seasons. The multi-objective function is optimized using adaptive cuckoo search algorithm (ACSA). In the proposed ACSA, a dynamically increasing switching parameter in a power of three is adopted for adjusting the random walk between local optima and global optima [47]. The superiority of the proposed ACSA in solving the multi-objective, nonlinear complex optimization problem over basic cuckoo search algorithm (CSA), particle swarm optimization (PSO), chicken swarm optimization (CSO), flower pollination algorithm (FPA), grasshopper optimization algorithm (GHO), whale optimization algorithm (WOA) and ant lion optimization (ALO) is presented by illustrating various case studies on standard IEEE 14, 30, and 118-bus test systems.

Related concepts and modeling
In this section, the power injection modeling (PIM) of OUPFC device and voltage collapse point indicator (VCPI)-based strategy for determining its optimal location is explained. Also, the ambient temperature impact on PV system generation, system loading conditions and transmission line parameters are explained.

Optimal unified power flow controller
Among the FACTS devices, the second-generation FACTS devices like unified power flow controller (UPFC) and interline power flow controller (IPFC) have wide range of applications in power system operation and control. UPFC is a versatile device and can be able to control the power flow parameters of a line (i.e., bus voltage, phase angle and line reactance) individually or in combination. IPFC is able to control power flows in multiple transmission lines at a substation simultaneously or consecutively. The advanced version of UPFC by integrating conventional phase shift transformer (PST) is proposed as optimal unified power flow controller (OUPFC) with its operating features and steady-state mathematical modeling [49]. The basic difference between UPFC and OUPFC can be understood by comparing Fig. 1a, b, respectively. Generally, UPFC consists of shunt converter coupled with excitation transformer and series converter coupled with injecting transformer as shown in Fig. 1a. A similar configuration can be found in Fig. 1b for OUPFC except the transformers with triple winding. The secondary windings of these two transformers are connected by a phase shift transformer (PST) which can be controlled by static/mechanical switches to inject a voltage with fixed phase into the transmission line. On the other side, the tertiary windings of these two transformers are used for conventional UPFC configuration. According to PIM of OUPFC, the injection powers at bus i and bus j are given as follows: where V i and V j are the magnitudes and θ i and θ j are the voltage angles at bus-i and busj, respectively. Here, b se = 1/x se , where x se is the total reactance of the circuit includes transmission line and injecting transformer reactance, r is the radius of the UPFC operating region and γ upfc is the phase angle of the UPFC; a is the ration between injection voltage of PST and exciting transformer voltage and α pst is the phase angle of PST.

Optimal location of OUPFC
All the FACTS devices have their own mode of operation and designed to control some specific attributes, the location should be optimized to get better results. From the literature, the solution approaches for optimal location of FACTS devices problem can be classified into three categories like direct, strategic and optimization problem-based approaches. In direct approaches, the impact on a specific parameter is analyzed directly by integrating FACTS devices at a selected location. Based on the significant impact on a specified parameter, the locations are further ranked to finalize optimal location. The computational effort involved in direct approaches can overcome by using strategic approaches. These approaches are basically based on sensitivity analysis with control variables in load flow study. By analyzing the impact of change in control variable on a specific attribute, the location of FACTS devices is finalized. Similarly, various heuristic approaches which have been widely used for FACTS location and sizing can be found in [50,51].
Based on maximum power transfer, voltage collapse point indicator (VCPI) is proposed for assessing static voltage stability of the power system [52]. Consider a transmission line connected between bus-s (sending end) and bus-r (receiving end) having impedance Z sr = R sr + jX sr . For a connected load P r + jQ r at bus-r, VCPI is defined as 'the ratio between real power connected at receiving end (P r ) and maximum power (P r,max ) that can be loadable at critically stable operating point' and is given mathematically in (5).
(1) where V s is the sending end bus voltage; Ø r is the load power factor; Ø z is the impedance angle. According to VCPI, the line which has highest value can be treated as critical line due to lowest maximum transfer capability. For a specified operating condition, the lines are ranked as VCPI values and top ranked line is chosen for OUPFC integration.

Modeling of solar PV system
The power output of a solar PV module can be found as a function of global solar irradiance and ambient temperature using Eq. (8) [53].
where Correspondingly, the output of a PV system is modeled in proportional to the module output as given in Eq. (10).
Here, it is assumed that the PV systems are integrated to the grid via grid-compatible inverter, by which voltage and frequency regulation can be done effectively at the point of common coupling (PCC).
The impact of PV generation at a load bus if finally realized in terms of reduced load as given by: Correspondingly, there is a change in power factor ( cos φ 1 ) of that load bus and consequently, tan(Ø 1 ) is given by: In order to maintain power factor at the desire value (cos(Ø 1 )), the required reactive power to be supplied by PV inverter is given by: The new reactive loading condition at PV source bus-j is finally given by: where P max,m is the module maximum power under standard testing conditions (STC), V mmp and I mmp are the corresponding voltage and currents, respectively; V oc and I sc are the open-circuit voltage and short-circuit currents, respectively; G ref and T ref are the global irradiance and temperature at STC, respectively; G t and T a are the actual global irradiance and ambient temperature, respectively; IC PV is the installed capacity of the PV system in MW; P PV(ref ) and P PV(t) are the module output at reference and ambient temperature, respectively; P g,j is the real power generation of a PV system at ambient temperature in MW; S in is the inverter MVA rating; and Q in,j is the MVAr generation by the PV inverter; P d,j and Q d,j are the real and reactive power loads at bus-j, respectively.

Temperature dependency load modeling
The increased load demand due to raise in ambient temperature can be modeled as proposed [40] by: where P d(base) and P d(Ta) are the system load in MW at base case and ambient temperature, respectively; T a is the ambient temperature in degree Celsius; m is the slope of increasing power consumption per one degree Celsius.

Thermal modeling of transmission line resistance
The change in transmission line resistance due to change in ambient temperature can be modeled as follows [42]: where R(T ref ) and R(T a ) are the transmission line resistances at reference (T ref ) and ambient temperature (T a ), respectively; T c is the temperature constant depends on type of conductor material.

Problem formulation
The OPF problem can be formulated for single or multiple objectives by having a set of equal and inequal operating constraints.

Overall objective function
The overall objective function is formulated with multi-objectives to minimize total operating cost at CE sources (f 1 ), real power losses (f 2 ), average voltage collapse point indicator index (f 3 ) and average voltage deviation (f 4 ) simultaneously. The objective function expressed as, The individual objective functions are expressed in Eqs. (16)- (19) for real power loss, average line stability index and average voltage deviation index, respectively.
where a i , b i and c i are the cost coefficients of CE sources; P gi is power generation at CE bus0.00i; |V i,ref | and |V i | are the reference voltage magnitude and actual voltage magnitude for a particular operating condition, respectively; nCE, nb and nl are the number of CE generator buses, total buses in the system and number of lines, respectively; k is an index for the line connects between buses i and j; r k and I k are the resistance of line k and current through it, respectively; Y ii and θ ij are the shunt admittance at bus i and its angle, respectively; V i and V j and δ i & δ j are the voltage magnitudes and their angles at buses i and j, respectively.
Power flow equations corresponding to both real and reactive power balance equations are the equality constraints that can be written, for all the buses expect OUPFC incident buses, as Similarly, for the OUPFC incident buses, the real and reactive power balance equations can be written as where P inj,i and Q inj,i are the real and reactive power injections by OUPFC as given in (1)-(4).
The inequality constraints mainly real and reactive power generations at CE sources, bus voltages, load angles, tap-changer settings, shunt VAr injections and line flow limits are considered as variables in the optimization.

Cuckoo search optimization
This section presents the optimization procedure for scheduling of CE sources, OUPFC and system parameters under different ambient temperature conditions by sing basic cuckoo search algorithm (CSA) and its proposed improved version.

Overview of cuckoo search algorithm
Cuckoo search algorithm (CSA) is one of the recent nature-inspired algorithms introduced by Yang in 2009 which inspired by brood reproductive strategy of cuckoo birds to increase their population, and it is more effective than other same family algorithms such as bat-inspired algorithm (BIA), differential evolution (DE), simulated annealing (24) (SA), particle swarm optimization (PSO) and artificial bee colony (ABC) algorithms [55,54]. In fact, these algorithms are special cases of CSA and hence, CSA can outperform than these algorithms. The major difference between CSA to other similar algorithms is the balance between local random walk to global random walk by its switching parameter and hence effective in terms of convergence speed to reach global optima. The switching parameter pa ϵ [0, 1] used to control the balance between local and global random walk and are related mathematically as defined in Eq. (33).
where v t i and v t k are the current position selected by random permutation; α is positive step size scaling factor; s is step size; H is heavy-side function; Pa is switching parameter between local and global random walk; ε is random number from uniform distribution; ⊗ is entry-wise product of two vectors; L(s, ) is Lévy distribution used to define the step size of random walk.
Lévy flights essentially provide a random walk while their random steps are drawn from a Lévy distribution for large steps. One of the effective methodologies to generate step size is using Mantegnas equations using gamma distribution function, described by; where Γ () is the gamma distribution function.

Adaptive cuckoo search algorithm
Precisely, the improvements for CSA toward better efficiency have been addressed by various researcher using different probability distributions for defining the step size of random walk and dynamically adjustment for switching parameter between local and global random walks. In [57], Gauss distribution-based CSA has proposed and results have shown the Gauss-CS outperformed than originally Levy-CSA in terms of higher convergence rate with reduced average generation. In [58], Gaussian and Cauchy distributions have been proposed and proved their superiority than Levy-CSA by applying to the web documents clustering problem. Similarly, Cauchy distribution-based CSA also in [59] and has applied for solving the economic emission load dispatch problem with multiple fuels options, but the results have shown Levy-CSA is better than other Gaussian-CSA and Cauchy-CSA by having lesser computational time. In [60,61], Gamma distribution-based CSA has been proposed and proved the Gamma-CSA superiority than Levy-CSA in terms of accuracy and average time. On other side, lineally decreasing switching parameter instead of constant switching parameter in original CSA in [62] and sorting function instead of permutation [63] have been proposed for better efficiency. In (33) Page 12 of 29 Kavuturu and Narasimham Journal of Electrical Systems and Inf Technol (2020) 7:11 addition, the reader can found some other significant literature on CSA improvements in [64]. Instead of constant switching parameter, a dynamically increasing switching parameter in a power of three is adopted in the proposed ACSA for adjusting the random walk between local optimal and global optima [47]. Mathematically, (37)

Yes
Discovery alien-Cuckoo eggs by random biased walks using switching parameter as given in (36) and create new nests.
Read test system bus data, line data, ambient temperature, slop of increasing power consumption per one degree Celsius and CSA controlling parameters such as number of nests (nnest), probability rate (pa) and maximum iterations (itmax).

Start
Modify load profile at all the buses, resistance of transmission lines transformers as per ambient temperature and output power of PVs.  Solution approach using CSA In this section, the control variable in the multi-objective problem and the overall procedure for solving the proposed fitness function is described here.
Vector of control variables This paper is aimed to minimize the voltage deviation w.r.t. reference voltage, real power loss and voltage collapse point indicator by considering variations in RE generation and loading conditions with OUPFC optimal controls. Here, the vector of control variables consists of generation at conventional energy sources, generator bus voltage magnitudes, tap-changer settings, shunt MVAr injection and controlling variables of OUPFC device.
Fitness function In each iteration, the system bus data and line data are updated with new population comprises generator bus voltage magnitudes, tap-changer ratios, shunt MVAr injections, controlling variables of OUPFC device (UPFC voltage source magnitude and its angle as well as PST voltage regulation and its angle) and correspondingly power injections at their incident buses. By having NR power flow solution, total real power loss, average VCPI and average voltage deviation index are computed to update the overall objective function expressed in Eq. (16). The solution that minimizes the overall objective function is considered the best solution.
The overall procedure of the proposed ACSA approach for solving OPF problem is given as flowchart in Fig. 2.

Results and discussion
The effectiveness of the proposed ACSA approach for solving the OPF problem including OUPFC is illustrated using standard IEEE 14, 30 and 118-bus systems. The bus data, line data and cost coefficients of generator buses for all the test systems are taken from [65]. In general, the PV systems are not dispatchable and hence their minimum and maximum limits are considered equal to the actual generation at the specified weather conditions and not considered as controlling variable in the OPF problem. The MVA rating of PV inverter is assumed as equal to the installed plant capacity of PV system and their operating power factor is assumed commonly as 0.95 leading, and hence, they can also able to support reactive power as modeled in Eq. (13).
In all the case studies, it is considered as reference temperature (T ref ) = 25 °C, temperature constant (T c ) = 228.1 °C for hard-drawn aluminum conductors and 225 °C for aluminum transformer [42,54]. The operating constraints are considered as mentioned in section Problem formulation, step 6. Also, the parameters for ACSA are considered as follows: number of nests (n nest ) = 25, probability rate (p a ) = 0.25, and maximum iterations (it max ) = 50.

IEEE 14-bus system
This test system consists of 5 generator buses, 9 load buses in which bus-9 is having shunt VAr injection. The buses are interconnected by 20 transmission lines, in which 3 are having tap-changing transformers. The line MVA ratings are given in Appendix. It has totally 259 MW real power and 73.5 MVAr reactive power loads, respectively. It is assumed that the existing all generator buses in the standard test system as CE sources. Scenario 1 Since scenario (1) is similar to the normal test system without PV sources and OUPFC, the loading on the system remains equal to 259 MW real and 73.5 MVAr reactive loads, respectively. In case (a), the test system is assumed for weather impact only on system loading conditions and transmission system parameters. The optimal schedule of CE sources and different objective function values are given in Table 1 as scenario (1) under case (a). The operating cost (F 1 ) for optimal schedule is equal to 8054.41 $/h. For this schedule, the system has 8.819 MW real power losses (F 2 ), average voltage stability index (F 3 ) is 0.006 and average voltage deviation index (F 4 ) is equal to 0.0006.
In order to simulate case (b), buses-10, 11 and 12 are treated as solar PV generation sources with 25 MW installed capacity of each. TP250 series of TATA Power Solar module [48] is considered for all the PV systems, and its specifications are as follows: Also, the following strategy is followed for identifying OUPFC location. By excluding the transmission lines which are incident to generator buses, shunt VAr injection buses and the lines having tap-changers, the remaining lines are ranked based on their VCPI values. Among all, the line connected between bus-13 and bus-14 is stood in first place by having highest value of VCPI = 0.0099 and considered as optimal location for the integration of OUPFC.
In case (b), the actual loading condition is still remaining same due to load growth m = 0, but the CE sources are optimally scheduled for only effective loading (of 192.08 MW and 23.46 MVAr real and reactive power, respectively), which can obtain by reducing the total PV generation of 75.2 MW and 27.8 MVAr real and reactive powers, respectively. The corresponding operating cost (F 1 ) at CE sources is reduced to 5197.37 $/h from 8054.41 $/h. By having this optimal schedule at CE sources, PV generations and OUPFC controls in the network, the real power losses (F 2 ) have reduced to 8.2745 from 8.819 MW. But notably, the average voltage stability index (F 3 ) is increased to 0.0088 from 0.006 and average voltage deviation index (F 4 ) is increased to 0.0027 from 0.006. Scenario 2 In scenario 2, the weather conditions are considered little bit higher than scenario 1. Due to increased slop of load growth, the real power loading conditions are increased to 269.5 MW from 259 MW. As compared to case (a) in scenario 1, the operating cost is increased 8482.53 $/h from 8054.41 $/h, losses are increased to 8.9987 MW from 8.819 MW, stability index increased to 0.0067 from 0.006 and voltage deviation index is increased to 0.0007 from 0.0006, respectively. This operating condition is clearly indicating the negative impact on system operating condition due to rise in ambient temperature and load growth.
In case (b), the PV generation is also decreased to 51.65 MW from 75.2 MW and 21.57 MVAr from 27.8 MVAr due to rise in ambient temperature. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls  in the network are given in Table 2 as scenario 2 under case (b). It is observed that the operating cost of CE is decreased to 6448.67 $/h from 8482.53 $/h, whereas losses, stability index and voltage deviation indexes are increased considerably. Scenario 3 Under this scenario, the test system loading is increased to 282.8 MW. The generation is increased to 292.00 at CE sources including losses. Correspondingly, the system has F 1 = 9024.33 $/h, F 2 = 9.1974 MW, F 3 = 0.0076 and F 4 = 0.0008. In case (b), the test system has OUPFC controls and PV generation equal to 47.45 MW and 22.11 MVAr; consequently, the CE sources are scheduled optimally for 244.31 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 7128.36 $/h, losses (F 2 ) = 8.8649 MW, stability index (F 3 ) = 0.0073 and voltage deviation index (F 4 ) = 0.0007. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2  For this optimal schedule, the system has an operating cost (F 1 ) = 7895.83 $/h, losses (F 2 ) = 8.9808 MW, stability index (F 3 ) = 0.0071 and voltage deviation index (F 4 ) = 0.0007. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2 as scenario 4 under case (b).
Scenario 5 Under this scenario, the test system loading is increased to 317.8 MW. The generation is increased to 327.65 MW at CE sources including losses. Correspondingly, the system has F 1 = 10,460.34 $/h, F 2 = 9.8233 MW, F 3 = 0.0112 and F 4 = 0.0019. In case (b), the test system has OUPFC controls and PV generation equal to 41.54 MW and 25.2 MVAr; consequently, the CE sources are scheduled optimally for 285.47 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 8766.25 $/h, losses (F 2 ) = 9.2046 MW, stability index (F 3 ) = 0.0079 and voltage deviation index (F 4 ) = 0.0008. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2 as scenario 5 under case (b).

IEEE 30-bus system
The test system has 6 generators buses (assumed as CE sources), 24 load buses of which 2 buses are having shunt MVAr injections and connected by 41 transmissions in which 4 lines are having tap0.00changers. It has totally 283.4 MW and 126.2 MVAr real and reactive power loads, respectively. In this test system, 3 locations namely buses 12, 17, and 21 are assumed as PV sources of 25 MW capacity each. For the optimal location of OUPFC, the methodology explained in IEEE 14-bus system is followed here also. As per highest value of VCPI = 0.0168, line 38  including losses. Correspondingly, the system has F 1 = 8870.97 $/h, F 2 = 11.1527 MW, F 3 = 0.0058 and F 4 = 0.0013. In case (b), the test system has OUPFC controls and PV generation equal to 75.2 MW and 36.82 MVAr; consequently, the CE sources are scheduled optimally for 218.32 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 5956.04 $/h, losses (F 2 ) = 10.1207 MW, stability index (F 3 ) = 0.0067 and voltage deviation index (F 4 ) = 0.0018. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2 as scenario 1 under case (b). Scenario 2 Under this scenario, the test system loading is increased to 305.9 MW from 283.4 MW. The total generation at CE sources is 317.65 MW at CE sources including losses. Correspondingly, the system has F 1 = 9797.44 $/h, F 2 = 11.7513 MW, F 3 = 0.0082 and F 4 = 0.0024. In case (b), the test system has OUPFC controls and PV generation equal to 51.65 MW and 30.59 MVAr; consequently, the CE sources are scheduled optimally for 265.77 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 7732.51 $/h, losses (F 2 ) = 11.5181 MW, stability index (F 3 ) = 0.0063 and voltage deviation index (F 4 ) = 0.0012. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2 as scenario 2 under case (b).
Scenario 3 Under this scenario, the test system loading is increased to 334.4 MW from 283.4 MW. The total generation at CE sources is 347.17 MW at CE sources including losses. Correspondingly, the system has F 1 = 10,984.77 $/h, F 2 = 12.7711 MW, F 3 = 0.0117 and F 4 = 0.0046. In case (b), the test system has OUPFC controls and PV generation equal to 47.45 MW and 31.13 MVAr; consequently, the CE sources are scheduled optimally for 299.04 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 9055.7 $/h, losses (F 2 ) = 12.0917 MW, stability index (F 3 ) = 0.0086 and voltage deviation index (F 4 ) = 0.0024. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2 as scenario 3 under case (b).
Scenario 4 Under this scenario, the test system loading is increased to 368.9 MW from 283.4 MW. The total generation at CE sources is 374.89 MW at CE sources including losses. Correspondingly, the system has F 1 = 13,329.5 $/h, F 2 = 5.9863 MW, F 3 = 0.0162 and F 4 = 0.0058. In case (b), the test system has OUPFC controls and PV generation equal to 44.171 MW and 32.37 MVAr; consequently, the CE sources are scheduled optimally for 333.76 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 10,876.59 $/h, losses (F 2 ) = 9.0299 MW, stability index (F 3 ) = 0.0143 and voltage deviation index (F 4 ) = 0.0059. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2 as scenario 4 under case (b).
Scenario 5 Under this scenario, the test system loading is increased to 409.4 MW from 283.4 MW. The total generation at CE sources is 422.04 MW at CE sources including losses. Correspondingly, the system has F 1 = 16,415.97 $/h, F 2 = 12.6449 MW, F 3 = 0.0222 and F 4 = 0.0131. In case (b), the test system has OUPFC controls and PV generation equal to 41.54 MW and 34.22 MVAr; consequently, the CE sources are scheduled optimally for 379.08 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 14,693.45 $/h, losses (F 2 ) = 11.2169 MW, stability index (F 3 ) = 0.0177      and voltage deviation index (F 4 ) = 0.0086. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 2 as scenario 5 under case (b).

IEEE 118-bus system
This test system consists of 54 generator buses and 64 load buses in which 10 buses are with shunt VAr injections and 2 buses are with shunt reactors. All the buses are interconnected by 186 interconnected transmission lines, in which 9 are with tap-changing transformers. It has totally 4242 MW real and 1438 MVAr reactive loads, respectively. It is assumed that the existing all generator buses in the system as CE sources and 5 load buses, i.e., 53, 54, 70, 71 and 78 are installed with 50 MW solar PV systems. According to VCPI values, lines are ranked and line-54 with 0.0182 (connected between bus-30 to bus-38) and line-126 with 0.0112 (connected between bus-68 to bus-81) are top-ranked and chosen for two OUPFC devices integration simultaneously. In this test system, the MVA limits for transmission lines are not considered in the optimization problem. The test system loading details under different scenarios and correspondingly obtained objective functions are given in Table 3. Also, the optimal generation schedules of all CE sources are given in Table 4. There are certain generators (i.e., 4, 8, 24, 62, 77, 85, 90, 91, 92, 99 and 116) which are not picked up generation under any scenario and hence not listed in Table 4.
Scenario 1 Under this scenario, the test system loading remains equal to 4242 MW and 1438 MVAr. The total generation at CE sources is 4315.5234 MW at CE sources including losses. Correspondingly, the system has F 1 = 129,440.95 $/h, F 2 = 73.5234 MW, F 3 = 0.0038 and F 4 = 0.0005. In case (b), the test system has 2 OUPFC controls and PV generation at 5 locations equal to 250.66 MW and 81.66 MVAr; consequently, the CE sources are scheduled optimally for 4065.33 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 119,544.86 $/h, losses (F 2 ) = 73.9873 MW, stability index (F 3 ) = 0.004 and voltage deviation index (F 4 ) = 0.0006. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 3 as scenario 1 under case (b).
Scenario 2 Under this scenario, the test system loading is increased to 4330.5 MW from 4242 MW. The total generation at CE sources is 4404.8 MW at CE sources including losses. Correspondingly, the system has F 1 = 133,007.71 $/h, F 2 = 74.3031 MW, F 3 = 0.004 and F 4 = 0.0006. In case (b), the test system has OUPFC controls and PV generation equal to 172.15 MW and 58.37 MVAr; consequently, the CE sources are scheduled optimally for 4232.37 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 126,156.87 $/h, losses (F 2 ) = 74.0204 MW, stability index (F 3 ) = 0.0039 and voltage deviation index (F 4 ) = 0.0005. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 3  For this optimal schedule, the system has an operating cost (F 1 ) = 137,072.44 $/h, losses (F 2 ) = 75.324 MW, stability index (F 3 ) = 0.0045 and voltage deviation index (F 4 ) = 0.0006. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 3 as scenario 4 under case (b).
Scenario 5 Under this scenario, the test system loading is increased to 4737.6 MW from 4242 MW. The total generation at CE sources is 4814.26 MW at CE sources including losses. Correspondingly, the system has F 1 = 149,412.16 $/h, F 2 = 76.6637 MW, F 3 = 0.0051 and F 4 = 0.0007. In case (b), the test system has OUPFC controls and PV generation equal to 138.46 MW and 58.88 MVAr; consequently, the CE sources are scheduled optimally for 4675.14 MW. For this optimal schedule, the system has an operating cost (F 1 ) = 143,825.23 $/h, losses (F 2 ) = 75.997 MW, stability index (F 3 ) = 0.005 and voltage deviation index (F 4 ) = 0.007. The optimal schedule at CE sources and correspondingly objective functions with PV sources and OUPFC controls in the network are given in Table 3 as scenario 5 under case (b).
Case (a) and (b) of scenario 5 are simulated with basic CSA, flower pollination algorithm (FPA), chicken swarm optimization (CSO) and particle swarm optimization (PSO), grasshopper optimization algorithm (GOA), whale optimization algorithm (WOA) and ant lion optimization (ALO) algorithms and the results given in Table 5 for case (a) and in Table 6 for case (b), respectively. By observing the results, ACSA is outperformed than CSA and PSO by providing the lowest objective function values. Due to space constraints, the voltage profiles obtained only with ACSA for case (a) and case (b) in all the scenarios are given in Fig. 3 and Fig. 4, respectively. From the results, the test system has minimum voltage in case (a) is 1.0448 p.u. (see Fig. 3) and in case (b), it is 1.0478 p.u. (see Fig. 4). This indicates the improved voltage profile with RE and FACTS devices in the network and consequently enhanced stability, reduced losses. Comparison of different algorithms in terms of overall objective functions for scenario 5: case (a) and case (b) is given in Fig. 5 and Fig. 6, respectively. In both the cases, ACSA is outperformed than other algorithms with lowest overall objective function value and less computational time.

Conclusion
In this paper, the impact of changing weather conditions on transmission system parameters, PV generation and loading profile are modeled considering the effect of ambient temperature and analyzed w.r.t. power system economic operation. In order to overcome the transmission system dispatchable problems for economic schedule, optimal unified power flow controller (OUPFC) is proposed to integrate optimally in the system. Primarily the location of OUPFC is determined using voltage collapse point indicator (VCPI) and its controlling variables involved in power injection modeling and various system operating constraints are optimized. A multi-objective function is formulated to minimize operating cost, real power losses, average voltage collapse point indicator and average voltage deviation index and solved using adaptive cuckoo search algorithm