Optimal reactive power dispatch by chaotic biogeography based optimization

This paper proposes a chaotic biogeography based optimization (CBBO) for solving optimal reactive power dispatch (ORPD) problem. Based on biogeography based optimization (BBO) theory proposed by Dan Simon in 2008, a new artificial intelligence with full models and equations have been used to achieve the best solution for objective function of ORPD such as total power loss, voltage deviation and voltage stability index while satisying various constraints of power balance, voltage limits, transformers tap changer limits and switchable capacitor bank limits. The BBO has been enhanced its search ability by adding chaotic theory. Therefore, the proposed CBBO can obtain better solutiong quality than BBO for optimization problems. The proposed method has been tested on the IEEE-30 and IEEE-118 bus systems and the obtained results have been verified with other methods. The result comparison has indicated that the CBBO can be a promise method for dealing the ORPD problem


INTRODUCTION
The main objective of optimal reactive power dispatch (ORPD) [1] in electrical power system is to minimize the objective function via the optimal adjustment of the power system control variables, while at the same time satisfying various equality and inequality constraints.Some objective functions in ORPD to evaluate the quality of power system is real power loss, voltage deviation at load buses [2], voltage stability index [3].The equality constraints are the power flow balance equations, while the inequality constraints are the limits on the control variables and the operating limits of the power system dependent variables.
The problem control variables include the generator bus voltages, the transformer tap settings, and the reactive power of shunt compensator, while the problem dependent variables include the load bus voltages, the generator reactive powers, and the power line flows.
There are various techniques ranging were introduced to solve ORPD, from conventional methods to artificial intellgence based methods.These conventional methods have been used for approaching the ORPD is linear programming (LP) [4], mixed-integer programming (MIP) [5], interior point method (IPM) [6], dynamic programming (DP) [7] and quadratic programming (QP) [8].The convetinal optimizations are easily to be carried out, the results is acceptable but can be trapped in local minima and this optimization can not act on the discrete variables.Recently, meta-heuristic search methods become more popular in doing with ORPD.Several methods, most of them are based on the biological model like evolutionary and behavior in species, were used such as evolutionary programming (EP) [9], genetic algorithm (GA) [10], differential evolution (DE) [11], ant colony optimization (ACO) [12] and particle swarm optimization (PSO) [13].These methods can improve the solutions for ORPD although it is more complex and slow in performance.
In this project, we discuss about an evolutionary algorithm that was found in 2008 by Dan Simon [14], called Biogeography Based Optimization.It is based on the migration and mutation of species in natural and the status of ecosystem in different time.By supplying the full theory and model of CBBO, we proved the useful of this algorithm by testing on IEEE-30 bus system and IEEE-118 bus system.The results is compared with the other paper to evaluate the advantage or disadvantage of this method.

FORMULATION OF ORPD
The ORPD problem is built based on the mathematics concepts 0 ) , ( 0 (1) where ( , ) F x u called the objective function whose output is the minimum value we want.

g( , )
x u is the equality constraints and h( , ) x u is the inequality constraints.
Applied the above to the ORPD problem, x is the containing vector of the controlled variables: the voltage and phase of load, reactive power of the generators and real power of slack bus.(P , ,..., , ,..., , ,..., ) u is the containing vector of the controlling variables: voltage of generators, tap-setting of transformers and the reactive power at compensator.
The objective function is depended on the target of optimization.Normally, there are three functions used: -The total active power loss in transmission: where k g is the conductance of branch k , i V is the voltage magnitude at bus i and ij  is the voltage angle different between bus i and j .
-Voltage deviation at loaded buses for voltage profile improvement: where sp i V is the standard value to evaluate the deviation, normally set at 1 p.u.
-Voltage stability index for voltage stability enhancement: max ( , ) max{ }; 1,..., where ( , )  g x u is the equality constraints, it follow the power conservation law: This equation can be spread: ( cos sin ) 0 ( sin cos ) 0 where gii, Bii are the transfer conductance and susceptance between bus i and bus j ; , b) The voltage limitations: d) The compensator capacitor limitations: e) The power flow limitations: where i S is the maximum power flow between bus i and bus j .( ) So with the penalty function, the objective function will be rewritten as:

Migration [14]
Mathematical models of biogeography describe how species migrate from one island to another, how new species arise, and how species become extinct.Geographical areas that are well suited as residences for biological species are said to have a high habitat suitability index (HSI) [14].The variables that characterize habitability are called suitability index variables (SIVs) [14].SIVs can be considered the independent variables of the habitat, and HSI can be considered the dependent variable.Habitats with a high HSI tend to have a large number of species, while those with a low HSI have a small number of species.Habitats with a high HSI have many species that emigrate to nearby habitats, simply by virtue of the large number of species that they host.Habitats with a high HSI have a low species immigration rate because they are already nearly saturated with species.Therefore, high HSI habitats are more static in their species distribution than low HSI habitats.By the same token, high HSI habitats have a high emigration rate; the large number of species on high HSI islands have many opportunities to emigrate to neighboring habitats.
The parameters below is used for BBO investigation: -Habitat suitability index (HSI): to evaluate the capability of the island for the creatures.
-Suitability index variables (SIVs): the independent variables such as rainfall, temperature, humidity… -Immigration rate:  -Emigration rate:  Considering the immigration curve.The maximum possible immigration rate to the habitat is which occurs when there are zero species in the habitat.As the number of species increases, the habitat becomes more crowded, fewer species are able to successfully survive immigration to the habitat, and the immigration rate decreases.The largest possible number of species that the habitat can support is at which point the immigration rate becomes zero.Now considering the emigration curve.If there are no species in the habitat then the emigration rate must be zero.As the number of species increases, the habitat becomes more crowded, more species are able to leave the habitat to explore other possible residences, and the emigration rate increases.The maximum emigration rate is which occurs when the habitat contains the largest number of species that it can support.
The equilibrium number of species is, at which point the immigration and emigration rates are equal.However, there may be occasional excursions from due to temporal effects.Positive excursions could be due to a sudden spurt of immigration, or a sudden burst of speciation.Negative excursions from could be due to disease, the introduction of an especially ravenous predator, or some other natural catastrophe.It can take a long time in nature for species counts to reach equilibrium after a major perturbation.
With the linear curves, the value of ,   whether there are s species in habitat can be written as: where: -E is the highest emigration rate -I is the highest immigration rate -n is the maximum species in habitat We use the emigration and immigration rates of each solution to probabilistically share information between habitats.If a given solution is selected to be modified, then we use its immigration rate  to probabilistically decide whether or not to modify each suitability index variable (SIV) in that solution.If a given SIV in a given solution i S is selected to be modified, then we use the emigration rates  of the other solutions to probabilistically decide which of the solutions should migrate a randomly selected SIV to solution i S .
The BBO migration strategy is similar to the global recombination approach of the breeder GA and evolutionary strategies in which many parents can contribute to a single off-spring, but it differs in at least one important aspect.In evolutionary strategies, global recombination is used to create new solutions, while BBO migration is used to change existing solutions.Global recombination in evolutionary strategy is a reproductive process, while migration in BBO is an adaptive process; it is used to modify existing islands.
As with other population-based optimization algorithms, we typically incorporate some sort of elitism in order to retain the best solutions in the population.This prevents the best solutions from being corrupted by immigration.

Mutation [14]
Now, consider the probability P s that the habitat contains exactly S species.By calculate the limit of the changing time of habitat, 0 t   , we have the probability equation: If a given solution S has a low probability P s , then it is surprising that it exists as a solution.It is likely to mutate to some other solution.This can be implemented as a mutation rate m that is inversely proportional to the solution probability: where max m is the user-defined parameters.This mutation scheme tends to increase diversity among the population.Without this modification, the highly probable solutions will tend to be more dominant in the population.This mutation approach makes low HSI solutions likely to mutate, which gives them a chance of improving.It also makes high HSI solutions likely to mutate, which gives them a chance of improving even more than they already have.

Application BBO to ORPD problem
Step 1: Set the initial value for the BBO variables.The i-th species in BBO is a vector of controlling variables: The starting value of id X is defined by: m in max min (X X ) Step 2: Set the value of BBO algorithm.
Step 3: Run the Power-flow by Newton-Raphson method and check the constraint of controlling variables.
Step 4: Calculate the fitness value and compute ,   .
Step 5: Do the migration step Step 6: Do the mutation step Step 7: Back to the step 3 for the next iteration.
If the variables after step 5 and 6 is not satisfied the constraints, we optimize them by set the threshold for the variables:

Chaos theory and application in BBO algorithm
In BBO algorithm, we used the random value to define whether migration, mutation or not.It is absolutely incidental process.Various researches before and my results have pointed that this process complied with Normal (Gaussian Distribution) [15].The solutions complied with this distribution have very high probability near average point, means that the solutions is concentrated at a specific value which is not the minimum value.(see the Figure 2) To demolish this disadvantage of BBO, chaos theory was used to supply the comparing value in migration or mutation step.Chaos theory is used to research about systems that seem to be chaotic but can be predicted.This is applied in dynamic systems that is sensitive with initial conditions and have unlimited dimensions.This is popular applied in Soil Mechanics, Solar-system, Liquid convection, Geography and Economics.
A chaotic map in this paper is a reflect: [0,1]→[0,1] by the recursive function: x n is the value of chaotic map at n-th iteration.The orbit of function can be easily predict by the characteristics of value and convergence.Each chaotic map has unique characteristics and with the different initial values, we have the different displays of the graph of function.The chaotic maps used in this project are list below [16].

RESULTS
We use the CBBO algorithm to apply in the IEEE-30 bus and IEEE-118 bus system to calculate and evaluate with the other recent project.With a chaotic map, we run with the initial value in {0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9}.The algorithm is simulated on MATLAB 2012 R2012b and the CPU: Intel core i5, 2.4 Ghz, 2.00 GB RAM.

IEEE-30 bus System
The IEEE-30 bus system is available in [17] with the data in the two following tables.In this paper, the power flow solutions for the systems are obatined from Matpower toolbox [18].In test system, the generators are located at buses 1, 2, 5, 8, 11, 13 and the available transformers are located on lines 6-9, 6-10, 4-12 and 27-28.The switchable capacitor banks will be installed at buses 10, 12,15,17,20,21,23,24 and 29 with the minimum and maximum values of 0 and 5 MVAr, respectively.The limits for controls variables are given in [20], generation active power in [21], and power flow transmisson lines in [22].The number of population is set to 10, the maximum iterations is 200 and the results were got by 50 independent runs.The comparion results were from [19].Two following figures shows the results of 50 independent runs of "random" BBO and a random CBBO in optimal total power loss, respectively, to clear the optimization of chaos theory to BBO.
Clearly, only 6% of solutions in BBO is far away the average point but the CBBO have high probability (18%) of values near the minimum value of computing.The minimum value in CBBO is better a lot than the BBO's.The results in CBBO is presented in below table with comparing results by three criteria: total power loss, voltage deviation and voltage stability index, respectively.

IEEE-118 bus System
The IEEE-118 bus system is available in [17] with the data in the two following table The limits of variables is similar with IV.a.The limits for controls variables are given in [20], generation active power in [21], and power flow transmisson lines in [22].The number of population is set to 30, the maximum iterations is 200 and the results were got by 50 independent runs.The comparion results were from [19].The results in CBBO is presented in below table with comparing results by three criteria: total power loss, voltage deviation and voltage stability index, respectively.

CONCLUSIONS
In this paper, a new artificial intelligence based method BBO has been presented with full overview and results.With the optimization by chaos theory, CBBO have high probability for searching and approach the minimum value of objective function of the ORD probem better than BBO algorithm.For the result comparison, the method is shown more useful with the large searching space with more variables althoungh the CBBO is not effective in searching voltage deviation and voltage stability index value.By testing on the IEEE-30 bus and IEEE-118 bus systems, the proposed method has shown that it is more effective for large scale systems.Therefore, the proposed CBBO is very favoable for solving the large-scale ORPD problem.

2 (
inequality constraint, we use the Static Square method.The objective function F not only have the output value but also adding the penalty function

Figure 1 .
Figure 1.Linear curve of fitness -migration

Figure 2
Figure 2 Values of min loss P

Figure 3
Figure 3 Values of min loss P )

Table 2 .
The structure of the experimented IEEE-30 bus system

Table 3 .
Basic values in IEEE-30 bus test system

Table 4 .
Result by CBBO methods for the IEEE-30 bus system with power loss objective and comparison

Table 5 .
Result by CBBO methods for the IEEE-30 bus system with voltage deviation objective and comparison

Table 6 .
Result by CBBO methods for the IEEE-30 bus system with voltage stability index objective and

Table 7 .
The structure of the experimented IEEE-118 bus system

Table 8 .
Basic values in IEEE-118 bus test system

Table 9 .
Result by CBBO methods for the IEEE-118 bus system with power loss objective and comparison

Table 10 .
Result by CBBO methods for the IEEE-118

Table 11 .
Result by CBBO methods for the IEEE-118