Improved Cat Swarm Optimization for Simultaneous Allocation of DSTATCOM and DGs in Distribution Systems

This paper addresses a new methodology for the simultaneous optimal allocation of DSTATCOM and DG in radial distribution systems to maximize power loss reduction while maintaining better node voltage profiles under multilevel load profile. Cat Swarm Optimization (CSO) is one of the recently developed powerful swarm intelligence-based optimization techniques that mimics the natural behavior of cats but usually suffers from poor convergence and accuracy while subjected to large dimension problem. Therefore, an Improved CSO (ICSO) technique is proposed to efficiently solve the problem where the seeking mode of CSO is modified to enhance its exploitation potential. In addition, the problem search space is virtually squeezed by suggesting an intelligent search approach which smartly scans the problem search space. Further, the effect of network reconfiguration has also been investigated after optimally placing DSTATCOMs and DGs in the distribution network.The suggested measures enhance the convergence and accuracy of the algorithmwithout loss of diversity.Theproposedmethod is investigated on 69-bus test distribution system and the application results are very promising for the operation of smart distribution systems.


Introduction
The electric power industries have witnessed many reforms in recent years.The existing distribution systems are moving towards smart distribution systems to achieve larger socioeconomic and other nontangible benefits.The rise of smart grid is a boon not only to society as a whole but also to all who are involved in the electric power industry, its customers, and its stakeholders [1].Building of such distribution systems requires local generation of reactive and active power using distributed energy resources (DERs) such as Distribution Static Compensators (DSTATCOMs) and Distributed Generations (DGs).DSTATCOM is a power electronic-based synchronous voltage generator capable of providing rapid and uninterrupted capacitive and inductive reactive power supply [2].Various renewable and nonrenewable DG technologies are available on the market today, such as microturbines, fuel cells, combustion gas turbines, photovoltaic, wind turbines, and combined heat and power [3].The integration of renewable DG technologies such as photovoltaic and wind turbines is becoming more popular in distribution systems on account of smart grid initiatives and strict environmental laws.These components allow increased efficiency, more reliability, and better quality of electric service.Moreover, they also facilitate effective utilization and life extension of existing distribution system infrastructure [1].However, optimal placement and sizing of these components are the important issues to extract maximum possible benefits.
The optimal allocation of DSTATCOM and DGs in distribution systems is a highly nonlinear complex combinatorial problem which has to satisfy various equality and inequality constraints.In the recent past, several populationbased metaheuristic techniques such as Genetic Algorithm (GA), Ant Colony Optimization (ACO), Immune Algorithm (IA), Differential Evolution Algorithm (DEA), Firefly Algorithm (FA), Particle Swarm Optimization (PSO), Teaching-Learning Based Optimization (TLBO), Artificial Bee Colony (ABC), Harmony Search Algorithm (HSA), and Cuckoo Search Algorithm (CSA) have shown their potential to solve optimal DSTATCOM placement problem [4][5][6][7] or optimal DG placement problem [8][9][10][11][12].However, the simultaneous placement strategy can independently set and control the real and reactive power flow in distribution networks [13].
A lot of research work has been carried out to successfully optimize the siting and sizing problems of active and reactive components when allocated separately.However, only a few researchers have attempted simultaneous placement strategy.References [14][15][16][17] have shown mutual impact of these components on the performance of distribution networks.Abu-Mouti and El-Hawary [14] employed an ABC algorithm to determine the optimal size of DGs, power factor, and location to minimize power losses.A heuristic approach is suggested by Naik et al. [15] where a node sensitivity analysis is used to identify the candidate DER sites, and their optimal capacities are determined by suggesting heuristic curve fitting technique.Moradi et al. [17] proposed a combined imperialist competitive algorithm-(ICA-) GA method to solve this multiobjective optimization problem.In this method, first ICA is used to find the sites and sizing of DERs and then the operators of GA are employed to further refine these solutions.
The smart grid requires integrated solutions for available distributed resources that reflect their coexistence to achieve higher efficiency through loss minimization and good quality power supply.Distribution networks are reconfigured frequently with changing operating conditions, and it is one of the effective means to improve their performance.The network reconfiguration is a process that alters feeder topological structure by managing the open/close status of sectionalizing and tie-switches under contingencies or normal operating conditions [18].Changing network topology by reallocating loads from one feeder to another may balance loads among the feeders and decrease the real power losses [19].Therefore, this is another resource that can be utilized in conjunction with simultaneous placement of DSTATCOM and DG.This approach has possibly not been attempted till date.
Cat Swarm Optimization (CSO) is one of the recently established high performance computational techniques introduced by Chu and Tsai [20].CSO is inspired by the natural behavior of cats where two major behaviors of the cats are modeled into two submodels: seeking mode and the tracing mode.In the seeking mode, the cat looks around and seeks the next position to move to, whereas, in the tracing mode, the cat tracks some targets [21].The important property of CSO is that it provides local as well as global search capability simultaneously [22].It converges better and shows a better performance in finding the global best solution [21].It has been successfully applied to solve diverse engineering optimization problems such as linear antenna array synthesis [23], deployment of wireless sensors [21], IIR system identification [24], clustering [25], and linear phase FIR filter design [26].However, the exploration potential of CSO needs to be enhanced while subjected to large dimension problems by reviewing its seeking mode.Further, convergence and accuracy of the algorithm can be improved by suitably placing tentative solutions in the problem search space during the iterative process.
Several researchers [10,12,15,27], and many others, have squeezed the problem search space by restricting the number of candidate locations for placing these devices.They generated a node priority list using certain node sensitivitybased approach and then selecting top few nodes from it as the candidate sites to allocate these devices.This approach drastically reduces the problem search space and also the CPU time incurred.However, the sensitivities are normally calculated for the base case conditions, where no such devices are installed [28].Furthermore, when selecting only top nodes as the sensitive components, it did not give the true picture of the entire distribution network [29].Therefore, such approaches are unreliable and thus lead the algorithm to suboptimal solution.
In light of the above discussion, a new Improved CSO-(ICSO-) based method is proposed for the simultaneous allocation of DSTATCOM and DGs in radial distribution networks.The objective is to maximize power loss reduction while maintaining a better node voltage profile.The distribution network is reconfigured after the optimal placement of these devices to extract maximum possible benefits.The seeking mode of CSO is modified to enhance exploitation potential of the algorithm.In addition, an intelligent search is proposed to enhance the overall performance of the optimizing tool.
The remainder of the paper is organized as follows: The problem is formulated in Section 2. The description of the standard and proposed CSO algorithm is presented in Sections 3 and 4, respectively.Section 5 deals with simulation results and the analysis of results is discussed in Section 6.Finally, the conclusions drawn from this work are presented in Section 7.

Problem Formulation
The node voltage profile of distribution systems can be improved by installation of DERs as DGs and DSTATCOMs, network reconfiguration, tap changing transformers, and so forth.The proposed algorithm is installing these components, and then the distribution network is reconfigured.Therefore, a soft voltage constraint is treated as the part of objective function while installing DSTATCOM and DGs, and the solutions are accepted by imposing penalty so long as the voltage constraint violates within prespecified limits.However, a hard voltage constraint is necessary for the system operation.Thus it is employed while reconfiguring the distribution network after optimally placing these components.The amount of voltage profile improvement and the capacities of these components employed are not linearly related.Therefore, the advantage of soft voltage constraint employed in placing DSTATCOMs and DGs is that it results in lesser capacity allocation in the optimal solution.In other words, the whole burden of voltage profile improvement should not be imposed over these distributed components, as the network reconfiguration can successfully improve voltage profiles.Therefore, the objective function   is formulated to maximize power loss reduction while maintaining a better node voltage profile by proposing a voltage penalty factor approach as defined below: where PF  is the node voltage deviation penalty factor which is given by where subject to the following operational constraints.
(a) Power Flow Equations.The sum of the power purchased from utility grid and the total power generated by the different sources in the distribution system must be balanced by the local load demand and the power loss in the lines.For a radial network, a set of recursive equations are used to model the power flow in the network as shown by (4).A sample twobus system including DG and DSTATCOM units is shown in Figure 1.Consider (b) Branch Current Limit.The current flow in each branch must satisfy the rated ampacity of each branch: (c) Active and Reactive Compensation Limit at a Node.The active and reactive power injected by DG and DSTATCOM at each node must be within their permissible ranges: (d) System Compensation Limit.The sum of active and reactive power injected by DGs and DSTATCOMs at all candidate nodes should be less than system nominal active and reactive power demand, respectively: Equations ( 8) and ( 9) ensure nonrepetition of candidate nodes for DSTATCOM and DG allocation, respectively: DG, ̸ =  DG, ; , ∈ .
The distribution network is reconfigured after optimally placing these devices.The reconfiguration problem is solved to minimize real power loss while satisfying various network operational constraints.The mathematical formulation of reconfiguration problem for loss minimization is formulated below: subject to constraints defined by ( 4)-( 5) along with radial topology and node voltage constraints as given below.
(a) Radiality Constraint.The reconfigured network topology must be radial, that is, with no closed path: A hard voltage constraint is employed during the network reconfiguration as it is one of the important network operation strategies.All node voltages of the system must be maintained within acceptable operating limits during the optimization process: In the present work the codification proposed in [18] is used to solve the network reconfiguration problem.This is a rule-based codification to check and correct infeasible radial topologies while attempting this problem.

A Brief Overview of CSO
Cats initiate their move very slowly and cautiously after sensing the presence of a prey and finally chase it very quickly.By observing these features, two modes of operation for CSO are simulated, that is, seeking and tracing modes.These modes have been mathematically modeled for solving optimization problems and are combined together by defining a mixture ratio (MR).The position of cats represents the set of tentative solutions in the problem search space.Every cat has its own position composed of  dimensions, velocity for each dimension, a fitness value, which represents the accommodation of the cat to the fitness function, and a flag to identify whether the cat is in seeking mode or tracing mode [20].These two modes of operation can be described briefly as given below.

Seeking Mode.
In this mode, random mutations within narrow range are employed at predefined dimensions on the cats selected according to MR.Some terms related to this mode can be defined as follows: Seeking Memory Pool (SMP) is the number of copies generated for each cat; Seeking Range of Selected Dimension (SRD) is the predefined range of each dimension being selected for mutation; CDC is the Counts of Dimension to Change, that is, the number of dimensions to be mutated.According to MR, a definite number of cats are selected for this mode.Making SMP copies of each cat: randomly select CDC dimensions for each copy.Each of these selected dimensions is varied in the range [−SRD, +SRD].In this way each cat is mutated to facilitate local search.However, higher SRD (>2) causes unnecessarily larger search space; higher CDC (>0.6) causes much higher diversity of the individuals by mutation, thus delaying the convergence; lower MR (<0.2) causes the cats to spend most of their time resting and observing (seeking mode), thus avoiding local entrapment to suboptimal solutions [26].The higher the value of SMP (>4) is, the better the exploration will be, but at the cost of more CPU time.Thus in seeking mode, the best values of SRD, CDC, MR, and SMP should be decided by the usual trade-off between accuracy and CPU time of the algorithm.

Tracing
Mode.This mode corresponds to a local search technique for the optimization problem.Cat traces the target while spending high energy in this mode.The rapid chase of the cat is mathematically modeled as a large change in its position.Define position and velocity of th cat in the dimensional search space as   = ( ,1 ,  ,2 , . . .,  , ) and   = ( ,1 ,  ,2 , . . .,  , ).The global best position of the cat swarm is represented as  best = ( best,1 ,  best,2 , . . .,  best, ) [22].In this mode, the cat moves according to its own velocity among all its dimensions so that it can trace the prey, that is, the best fit cat.The velocity  , and position  , updates of the th cat are governed by the following relations: where  is a constant, usually taken as 2, and (⋅) is a random number in the range [0, 1].Whether a cat is in seeking or tracing mode, it survives, if its fitness is improved.At the end of each iteration, cats from both seeking and tracing modes merge together before initiating the next iteration.

Proposed ICSO
In CSO, the seeking mode provides local search whereas the tracing mode searches globally.However, it will be better if the current best cat is allowed to search locally.If it happens, the current best cat may upgrade its fitness, and later on this will positively influence the movement of all the cats going through the tracing mode.Moreover, it can also avoid possible local trappings.Therefore, seeking mode of the standard CSO is modified in the perception of the above facts as described below.

Modified Seeking Mode.
Like other swarm optimization techniques, the philosophy of CSO is "to follow the leader."If the fitness of the current best cat is improved by some means, the convergence of CSO would be improved.It is therefore suggested that the current best cat is mandatorily selected for the seeking mode.Further, the way of generating cats around the current best cat is also different in the proposed seeking mode, as described below.
The operators of CSO inherently generate continuous decision variables.However, the decision variables are strictly integers in this problem.Therefore, while employing local search around the current best cat, two cats are generated from this cat by employing ceiling and flooring of its decision variables.All possible combinations of cats are generated from these two cats.An illustration of modified seeking mode is shown in Figure 2 as a flowchart.For simplicity, it is shown only for two dimensions.The fitness of all these cats is evaluated and the best fit cat is updated, if a better cat is found.However, other cats of the seeking mode are updated as in the standard CSO.The local search around the best cat not only avoids local trappings but also facilitates the exploration of new search points in the problem search space during the tracing mode.

Intelligent Search.
While initializing or otherwise, it will be always better if all the tentative solutions spread in such a way that most of them lie near the promising region.But this is a difficult task.Nevertheless, a sufficient diversity is essential to explore new solution points in the problem search space.Several researchers, as mentioned earlier, have applied  perturbation-based node sensitivity approach to get the node priority list for optimal allocation of distributed resources and then selecting top few nodes from it to redefine the problem search space.However, none of the sensitivity-based approaches is foolproof so it is possible that the optimal node may not lie in the redefined problem search space.Therefore, such approaches are unreliable and may cause erroneous results.In the present work, candidate nodes are selected by proposing an intelligent search approach.The proposed approach is different than the conventional one in the sense that the node sensitivity is observed in terms of the change in objective function to be optimized, instead of the change in power loss.Moreover, the system nodes are arranged in the decreasing order of change in the objective function and then the candidate nodes are selected from this list using roulette wheel selection.In this way, the nodes are selected according to their probability of priority.Thus more chances are available to those nodes which are good for the allocation of these devices.However, this approach reserves the right of each system node to remain in the problem search space during the computational process.Therefore, using intelligent search the problem search space is virtually squeezed without loss of diversity.An illustration of intelligent search is shown in Figure 3.

Individual's Encoding.
The structure of the individuals for the proposed methods is shown in Figure 4 which is composed of candidate nodes and sizing for the respective candidate DSTATCOM and DGs.The candidate nodes are allocated using intelligent search approach, whereas the sizing of distributed resources is selected randomly within their respective predefined bounds as described by ( 6)- (7).

Termination Criterion.
Elitism is not required due to the intrinsic nature of CSO.Therefore, the termination criterion is taken as follows: "when either the maximum iteration count is exhausted or all cats acquire the same fitness, the evolutionary process stops."This criterion is selected to reduce the CPU time of the algorithm.The flowchart of the proposed ICSO method is presented in Figure 5.

Simulation Results
The proposed CSO method is investigated on 69-bus test distribution system taken from [30].The annual load profile is assumed to be piecewise segmented in three different load levels, that is, light, nominal, and peak, which are 50%, 100%, and 160% of the nominal system load, respectively, and the corresponding load durations are taken 2000, 5260, and 1500 hours as in [27].It has been assumed that DSTATCOM exchanges only reactive power with the network and DG is dispatchable and operated at unity power factor.The initial data of this system are presented in Table 1 and various design parameters selected are given in Table 2.The Newton-Raphson power flow method is used for load flow of the distribution system.The optimal solution obtained, for each load level, after 100 trials of proposed ICSO is presented in  respectively.The table also shows that the optimal network configuration is affected by the presence of these devices.The performance of the distribution network obtained using this optimal solution is presented in Table 4.
The table depicts power loss reduction after optimal allocation of these devices for each load level.It can be observed from the table that annual energy losses for this system are reduced significantly by about 88%.This shows that simultaneous placement strategy is very useful for distribution systems.The feeder power losses are further reduced when distribution network is optimally reconfigured after optimally placing DSTATCOM and DG in the distribution system.This causes a net annual energy loss reduction of about 94% from the base case network which is substantial.Thus, the reconfiguration of distribution network can be fruitfully utilized in order to enhance the energy efficiency of distribution systems.It is interesting to note that the minimum node voltage obtained by optimally placing DSTATCOM and DG is found to be 0.9493 p.u. at peak load level.However, this voltage is enhanced to 0.9751 p.u. by network reconfiguration.This shows the utility of the proposed penalty factor approach where a soft voltage constraint is imposed to permit the solutions having minimum voltage below 0.95 p.u. during evolutionary process.
The voltage profile obtained using the proposed method during all load conditions is presented and compared with base case in Figure 6.It can be observed from the figures that the voltage profile is improved by simultaneous placement of DSTATCOM and DG which is then further enhanced by network reconfiguration.This is true for all load levels considered.All voltages are found to be within permissible limits while integrated approach is employed.Thus the integrated solution provided by the proposed ICSO successfully achieved the desired objectives.
Not much literature is available for the validation of the proposed method.A comparison result with PSO [2] is presented in Table 5.The only scenario is considered where the number of sites for locating each of these components is being restricted to one.The table shows that the power loss reduction and minimum voltage provided by the solution using the proposed method are much better than those obtained with PSO [2].It happens because the penetration limits of DSTATCOMs and DGs considered in [2] are equal to the nominal load at the nodes where they have to be installed.This penetration limit is 0.888 MVAr and 1.244 MW, respectively, for node 61.However, PSO [2] explores only 0.1223 MW capacities for optimal DGs.Thus, the proposed ICSO is capable of generating a better solution than the existing PSO technique.

Discussion
The proposed ICSO has shown its potential to efficiently solve one of the large dimension problems of power system.It happens because of the proposed suggestions in the standard CSO.The comparison of the convergence characteristics of CSO and ICSO is presented in Figure 7.The figure also compares the convergence of PSO when applied for the problem with the same population size and maximum iterations.The figure shows a marked improvement in the convergence of ICSO, whereas the CSO and PSO seem to be getting trapped in local optima.It happens due to the poor exploitation of the search space.In ICSO, the intelligent search dispersed all individuals in the near vicinity of the promising region, without scarifying diversity.This causes better chances for good individuals to surf the problem search space amicably.Further, the modified seeking mode enhances

Conclusions
An integrated approach for the optimal placement of distributed resources and network reconfiguration using improved variant of CSO technique is presented.The objectives of maximum power loss reduction and voltage profile enhancements have been successfully achieved using the proposed modeling.It has been observed that the proposed integrated approach of simultaneous placing of DSTATCOMs and DGs in the distribution network and then optimally reconfiguring it can enhance the performance of the systems by a good margin.The suggested seeking mode in ICSO provides local random walk to the best fit cat and the proposed intelligent search scans the problem search space efficiently without loss of diversity.This causes dispersion of all tentative solutions in the close vicinity of the promising region which helps to obtain the global or near-global optima.
The comparison result validates that the proposed method is better than existing PSO and CSO methods.It has been observed that the proposed ICSO has better convergence and is also capable of generating fairly good quality solutions.

𝐶:
Acceleration coefficient : N u m b e ro fd e c i s i o nv a r i a b l e s : Each dimension of selected cat (1 ≤  ≤ ) : Total number of branches in the system   : Feeder

Figure 1 :
Figure 1: Single-line diagram of a two-bus system.

Figure 6 :
Figure 6: Comparison of voltage profiles for base case, after optimal allocation of devices and reconfiguration after optimal allocation of devices at (a) light load, (b) nominal load, and (c) peak load.

Table 1 :
Initial data of 69-bus system.

Table 3 :
Optimal solution for DSTATCOM and DG allocation.

Table 3 .
The table shows the capacity of DSTATCOM and DG in kVAr and kW,

Table 4 :
Distribution network performance using the proposed method.

Table 6 .
The table shows best, mean, and worst fitness of the sampled solutions.It can be observed from the table that the solution quality obtained using ICSO

Table 6 :
Comparison of solution qualities.
current at the th load level (p.u.)  loss, : Power loss for uncompensated system at the th load level (kW)  loss, : Power loss for compensated system at the th load level (kW)   : Real power for sending end of the th branch at the th load level (kW)  DSTAT , : Reactive power injection by DSTATCOM at the th node for the th load level (kVAr)  DSTAT ,min / DSTAT ,max : Minimum/maximum reactive compensation provided by DSTATCOMs at the th node (kVAr)   : Nominal reactive power demand of the system (kVAr)   : Reactive power for sending end of the th branch at the th load level (kVAr)   : L i n er e s i s t a n c eo ft h eth branch (Ω) (⋅): Random number in the range of [0, 1]  max / min : Maximum/minimum permissible node voltage (p.u.)  min : Minimum specified node voltage (p.u.)   : Voltage of the th node at the th load level (p.u.)   / +1 : V e l oc i t yo fth eth/( + 1)th cat Δ  : Maximum node voltage deviation of the th node at the th load level (p.u.)  best : Position of the cat that has the best fitness value   / +1 : Positionoftheth/( + 1)th cat   : L i n er e a c t a n c eo ft h eth branch (Ω) PF: Node voltage deviation penalty factor B  : C l o s e dl o o pa tt h eth load level.