Abstract

In this paper, stochastic fractal search method (SFS) is employed for solving the optimal reactive power flow (ORPF) problem with a target of optimizing total active power losses (TPL), voltage deviation (VD), and voltage stability index (VSI). SFS is an effective metaheuristic algorithm consisting of diffusion process and two update processes. ORPF is a complex problem giving challenges to applied algorithms by taking into account many complex constraints such as operating voltage from generators and loads, active and reactive power generation of generators, limit of capacitors, apparent power limit from branches, and tap setting of transformers. For verifying the performance, solutions of IEEE 30 and 118-bus system with TPL, VD, and VSI objectives are found by the SFS method with different control parameter settings. Result comparisons indicate that SFS is more favorable than other methods about finding effective solutions and having faster speed. As a result, it is suggested that SFS should be used for ORPF problem, and modifications performed on SFS are encouraged for better results.

1. Introduction

In the power system, optimal reactive power flow (ORPF) is not only one of the best famous optimization problems but also a very complex problem. In the ORPF problem, two variables need to be considered such as control variables and dependent variables. Control variables are voltage of generation buses, on load tap-changer setting of transformers and generated reactive power of capacitor banks, while dependent variables are voltage of load buses, apparent power flow of transmission lines, and reactive power of generators. So, the major objectives of such ORPF problem is to find control variable so that others have values falling into a permitted operating range [1, 2]. Traditionally, the ORPF problem concentrates on reducing three individual objectives such as power losses of transmission lines, voltage deviation, and voltage stability index. So, a power system economically and stably operates when these goals are fully achieved.

In the last decades of the 20th century, the ORPF problem has been successfully addressed by many conventional methodologies called deterministic methodologies such as the Newton method [3], linear programming [47], interior point method [8, 9], quadratic programming method [10, 11], and dynamic programming method [12]. With appearance of the mentioned methods, they proved their strong points in dealing with the ORPF problem having linear constraints and differentiable functions for application, but a large system or more complicated constraints and their applicability must be stopped to make rooms for new methods which have a promising ability.

Luckily, developing computer science supported researchers much in creating new population-based methods to handle drawbacks of conventional methods. These methods have been successfully and widely applied to solving the ORPF problem, consisting many original methods, improve/modified methods, or combined/hybrid methods. They have been constantly developed and have become a big method family such as particle swarm optimization (PSO) family [1317], differential evolution (DE) family [1821], and genetic algorithm (GA) family [2225], while many standard methods have been also applied in [2634]. Sahli et al. [16] presented a combination between particle swarm optimization and tabu search (PSO-TS) by incorporating the best search function of PSO and TS. It was capable of finding the global solution and avoiding to fall into local optimum. By the way of evaluation of results attained from the standard IEEE 30-bus system with objective of power loss minimization, PSO-TS has seen better solution quality than other methods as conventional DE, PSO, and TS. Furthermore, a modified version of PSO called modified pseudogradient search-particle swarm optimization (MPG-PSO) has been proposed in [17]. MPG-PSO has the most powerful ability in the PSO family due to applying pseudogradient theory for determining the best velocity direction. As a result, the method has overtaken other PSO methods involving PSO using the time-varying inertia weight factor (PSO-TVIW), PSO using time-varying acceleration coefficients (PSO-TVAC), self-organizing particle swarm optimization using time-varying acceleration coefficients (SPSO-TVAC), PSO using constriction factor (PSO-CF), pseudogradient-based PSO (PG-PSO), PSO using stochastic weight trade-off factor (SWT-PSO), and SWT-PSO using pseudogradient method (PGSWT-PSO). Differential evolution (DE) family has been offered for ORPF-like traditional differential evolution (DE), hybrid ant system and differential evolution method (HAS-DE) [20], and hybrid double differential evolution technique and modified teaching learning technique (DDET-MTLT) [21]. In [20], Huang and Huang have replaced the selection operation of DE by the ant system to enhance the global search capability and avoid falling into local minima and decrease computational time. DDET-MTLT [21] was a combination of double differential evolution technique (DDET) and modified teaching learning technique (MTLT). The obtained results of DDET-MTLT have been compared to some methods on IEEE-30 and IEEE 118-bus systems. Besides, variants of genetic algorithm have been applied for ORPF such as genetic algorithm (GA) [22], enhanced genetic algorithm (EGA) [23], modified NSGA-II (MNSGA-II) [24], self-adaptive real coded genetic algorithm (SARCGA) [25], and hybrid evolutionary programming technique HEP [25].

In addition to the three above method family, other standard methodologies have been also applied for solving ORPF problem such as gravitational search algorithm (GSA) [26], ant lion optimizer (ALO) [27], quasi-oppositional teaching learning based optimization (QOTLBO) [28], teaching learning based optimization (TLBO) [28], Pooled-neighbor swarm intelligence algorithm (PNSA) [29], hybrid Nelder–Mead simplex-based firefly algorithm (HFA-NMS) [30], chaotic krill herd algorithm (CKHA) [31], artificial bee colony algorithm (ABC) [32], exchange market algorithm (EMA) [33], backtracking search algorithm (BTSA) [34], and harmony search algorithm (HSA) [35]. In summary, all methods have demonstrated their qualification for addressing almost constraints of ORPF problem with acceptable solutions.

In this article, we present a standard stochastic fractal search (SFS) with the goal determining minimization of three different individual objectives of ORPF problem such as minimizing TPL, reducing of VD, and enhancing VSI. The standard stochastic fractal search [36] was developed in 2014 by Salimi. It has been applied for addressing twenty-three standard benchmark functions and has proven its proficiency in finding optimal solutions better than many methods available in this literature. In [37], Tran at al. have applied an improved SFS (ISFS) method and SFS method for dealing with the ORPF problem. Only the IEEE 30-bus system with three objectives consisting of total power losses, voltage deviation and L-index have been employed for comparing and evaluating the real performance of the SFS method and ISFS method. From the minimum, average, and maximum fitness functions obtained from the three objectives, ISFS has been considered to be more effective than the SFS method. However, Tran at al. [37] have not taken the setting of control parameters into account. In fact, SFS has three important control parameters consisting of walk factor, population size, and the maximum number of iterations. Among the three parameters, the walk factor has high impact on working performance of the diffusion process; meanwhile, population size has high contribution to the first update and the second update processes. So, in this paper, we focus on the setting of control parameters to overcome such mentioned shortcoming of the work in [37]. Furthermore, we also expand study cases by considering both IEEE 30-bus system and IEEE 118-bus system with the three mentioned objectives. As a result, the novelty and the main contribution of the paper are as follows:(i)Finding optimal solutions for IEEE 30-bus and 118-bus transmission power networks of ORPF problem by using the SFS method(ii)Testing the real performance of the SFS method with the change of population size and iterations(iii)Tuning the best walk factor for determining more appropriate equation for diffusion of the SFS method(iv)Illustrating the fluctuations of search process of SFS method with different settings;the simulation can support to evaluate the real performance of SFS and the impact of each parameter on the real performance of the SFS method(v)Demonstrating the effectiveness and robustness of the SFS method by comparing total power losses, load bus voltage, and enhancement of voltage profile

Apart from the introduction, other parts of the paper are as follows: single objective functions and constraints of electric components are mathematically formulated in Section 2. The structure of the SFS method consisting of diffusion and the first and the second update techniques are described in detail in Section 3. Computation steps of solving ORPF problem by using the SFS method are shown in Section 4. The results obtained by SFS and other methods from two standard IEEE transmission power networks with 30 buses and 118 buses are compared and discussed in Section 5. The whole work of the paper is summarized in and concluded in Section 6.

2. Formulation of ORPF Problem

The ORPF problem is constructed by considering minimization of total active power losses (Ploss), reduction of load bus voltage deviation (VD), and enhancement of voltage stability as objectives and taking into account operating voltage of generators and loads, apparent power limit of branches, reactive power generation of generators, and active and reactive power balance equations. The structure of ORPF problem is described in detail as follows.

2.1. Objectives of ORPF Problem

Three main objectives of the ORPF problem considered in the paper consist of minimization of total active power losses, minimization of voltage deviation at load buses, and minimization of L-index in which minimization of L-index is corresponding to the enhancement of voltage stability. The three objectives can be seen in formulas (1)–(3) as follows:where Gij is the conductance of conductor ij; βi and βj are the phases of voltage at buses i and j, respectively; Vloadi is the voltage of load bus i; Vref is expected to be the voltage equaling 1.0 pu; and Lj is called L-index of-bus j [16].

Basically, L-index is within the range from 0 to 1 in which 0 is the best value and 1 is the worst value. The power system is considered to be working stably when L-index is close to 0, and it is working unstably or it will be collapsed in some seconds if L-index is close to 1. So, the main task to keep the power system working stably and economically is to reduce L-index close to 0. However, it is hard to tune control parameter for obtaining 0 value for L-index.

2.2. Constraints of ORPF Problem

ORPF are constrained by equality and inequality constraints covering the whole transmission power network. The equality constraints are considering active and reactive power balance while the inequality constraints are upper and lower limitations of electricity components and working parameters of power network.

The two equality constraints are as follows:where and are real and unreal terms of admittance of conductor ij, respectively. The inequality constraints are limitations of reactive power output and voltage of generators, the reactive power output of capacitors, tap values of transformers, voltage of load buses, and apparent power of lines. All the inequality constraints are as follows:

3. Stochastic Fractal Search Algorithm

3.1. Diffusion Technique

SFS is an improved version of fractal search and was developed by Salimi in 2014 [36]. SFS is constructed by three main processes including diffusion and two different update processes. Consequently, SFS has three new solution generations in each iteration in which the diffusion process plays the most important role. The diffusion process uses Gaussian random walk for generating new solutions as follows:where Gbest is the best solution among the set of points. CI is the current iteration.

As seen from equations (12) and (13), the diffusion process can be accomplished by using either equations (12) or (13) meanwhile the two equations have a major difference. Equation (12) uses Gaussian random walk around the best solution and an updated step by using (Gbest − Xs) while equation (13) only employs random walk around Xi. Due to the difference, SFS must propose one more control parameter, called walk factor (WF) to control the use of either equations (12) or (13). The walk factor is compared to a random number and the result will lead to a decision of used equation. If walk factor is higher than the random number, equation (12) is employed for producing new solution for the sth solution. Otherwise, equation (13) is selected. The diffusion technique can be performed by using Algorithm 1.

 Set a value to WF
For s = 1 to Npop
  Randomly produce a random number λs for the sth solution
  If WF is higher than λs
    Using equation (12) for updating the sth new solution
  Else
    Using equation (13) for updating the sth new solution
  End
end
3.2. Second Update Technique

The first update is the second new solution generation and is performed by the following equation:where Xrandom1 and Xrandom2 are solutions chosen randomly in the current set of solutions and is the new solution of solution Xs.

3.3. Third Update Technique

After implementing the first update, the second update equivalent to the third new solutions is carried out for updating solutions of the second generation and is expressed bywhere Xrandom1, Xrandom2, and Xrandom3 are three random solutions chosen from the current set and RNs is a random number ranging from 0 to 1 for solution s.

4. Implementation of SFS for ORPF Problem

4.1. Initialization

The ORPF problem can be solved by using mathpower 4.1 programming after a set of control variables is predetermined by the SFS problem. In fact, the control variable set that has to be inserted into the program is the voltage of thermal generators (VG1, …, VGNG), tap values of transformers (T1, …, TNt), and reactive power output of capacitor banks (Qc1, …, QcNc). Thus, each point (corresponding to each solution) of the SFS method must represent all the variables and is randomly produced as follows:where Xmin and Xmax are the lowest and highest values of such control variables.

4.2. Fitness Function

As a result, the fitness function is constructed equaling the sum of objective function and penalty terms [17]. For three single objective functions including total power loss, voltage deviation, and L-index, three corresponding fitness functions are formulated as follows:where , VD and L-index are obtained by using equations (1)–(3), VPF is the penalty factor for the violation of the dependent variable, and , and are penalty terms corresponding to the violation of QGi in constraint (6), Vloadi in constraint (10), and Sl in constraint (11).

The penalty terms are determined by

4.3. Termination Condition

The three main processes of the SFS method will be terminated when current iteration (CI) is equal to the maximum iteration (MI), which is predetermined initially.

5. Numerical Results

In the section, we test the performance of the SFS method for the ORPF problem with two systems having 30 buses and 118 buses under considering three objectives such as power loss, voltage deviation, and L-index. The method is executed on Matlab program language and a computer with the processor of Core i7, 2.4 GHz, and 4 GB of RAM.

5.1. Results Obtained on IEEE 30-Bus System

In the section, we implement the SFS method for optimizing total power losses, voltage deviation, and L-index of the IEEE 30-bus system by setting different values to control parameters such as population size, maximum iteration, and walk factor. The IEEE 30-bus system consists of 6 generators, 24 loads, and 41 branches, 9 VAR compensators, and 4 transformers [38]. In order to indicate the impact of control parameters on results, we perform three testing cases as follows:(1)The first testing case: tune different values for population size while fixing the maximum iteration at value of 200(2)The second testing case: tune different values for the maximum iteration while fixing population size to the best value obtained by the first testing case(3)The third testing case: tune different values for walk factor while fixing population size and the maximum iteration at the best value obtained by the first testing case and the second testing case

The three testing cases are presented in Section 5.1.1, while the comparison of results from SFS method and other ones is shown in Section 5.1.2 and Section 5.2.

5.1.1. Analysis of Control Parameters on Results Obtained by SFS Method

The results in terms of the best, mean, worst power losses, standard deviation, and simulation time for each run (ST) are shown in Tables 13 corresponding to the first, second, and third testing cases. The results from Table 1 indicate that population size of 10 is high enough for finding the best solution while increasing the population size to 15 and 20 cannot find a better solution but average and maximum power losses can be improved whereas simulation time is increased. Figure 1 also gives the same evaluation since the curve in blue has high number of solutions with better fitness than other curves.

Table 2 show the results obtained by setting different values to the maximum iteration while keeping the population size at 10. The maximum iteration is increased from 50 to 250 with a step of 50. The best power loss implies that SFS method can find the best solution at different values of the maximum iteration such as 200 and 250; meanwhile, smaller number of iterations cannot reach the same best solution. The best optimal solution of MI = 200 and MI = 250 has the same quality; meanwhile, the mean fitness function value of 50 runs from MI = 250 is smaller. In spite of the advantage of MI = 250, SFS should adopt the most appropriate number of iterations is 200 because it can find as good as solution but it can reduce simulation time. Figure 2 shows that the curve in blue has many solutions with the same fitness as those from curve in black.

From Tables 1 and 2 as well as from Figures 1 and 2, it can result in the selection of population size and the maximum iteration in which the best value of the comer is 10 while the most appropriate value of the later is 200. By using the two values, walk factor’s impact on performance of the SFS method is tried by setting it to 0, 0.2, 0.4, 0.6, 0.8, and 1. Table 3 and Figure 3 show that the walk factor to be 1 is the best value for reaching the best solution, whereas other values cannot reach the same best solution. As observing from the best power loss, the walk factor with higher value can support the SFS method find superior solutions. In fact, power loss is the highest at WF = 0, and it decreases since the factor approaches to 1. So, it can conclude that diffusion technique becomes more effective if equation (12) is used to replace equation (13).

5.1.2. Result Comparisons for IEEE 30-Bus System

The results in terms of minimum, average, maximum, and standard deviation accompany with control parameters including MI and Npop from SFS and other methods are reported in Tables 46 for power loss, voltage deviation, and L-index, respectively. Besides, saving percentage (%) of the SFS method compared with each one is also calculated and reported in the tables for further comparisons. Saving percentage from these tables can see that SFS outperforms most methods excluding PSO-TS [16] and ISFS [37] for power loss objective and QOTLBO [28] for voltage deviation objective. Saving values show that these methods get improvement over SFS by 0.137%, 0.295%, and 2.45%; however, only ISFS [37] has found better solution than SFS, meanwhile PSO-TS [16] has not reported MI and Npop for comparison of convergence speed and recalculated minimum of QOTLBO is 0.1031, which is much higher than reported value of 0.0856. Clearly, QOTLBO is less effective than SFS. SFS can get improvement over other methods for power loss optimization from 0.7% to 8.365%, for voltage deviation optimization from 1.68% to 57.51%, and for L-index from 13.26% to 32.82%. Clearly, the improvement is significant and optimal solution of SFS is much better than those from other methods. Furthermore, convergence speed of SFS is also faster or approximate with other compared methods since SFS uses MI = 200, Npop = 10, while those from others are 100 and 50, 200 and 20, or MI = 30,000.

Optimal solutions of three optimization cases are given in Table 7.

5.2. Results Obtained on IEEE 118-Bus System

This section uses the IEEE 118-bus system consisting of 54 generator buses, 64 load buses, and 186 transmission lines, 14 VAR compensators, and 9 transformers [37] (Table 8). We implement the SFS method for optimizing total power losses of the system by setting different values to population size and the maximum iteration. In the first trial, we increase the population size from 10 to 100 with a change of 10 and fixing the maximum iteration to 150 and 200. In the second trial, population size is selected to be 50 and 75; meanwhile, the maximum iteration is increased from 50 to 400 with a step of 50. The optimal solution obtained by SFS method is given in Table 9. In addition, the results in terms of the best power loss, mean power loss, and success rate (SR) in percentage are, respectively, shown in Tables 9 and 10. From the two tables, it is clear that SFS can improve search ability since the population size and the maximum iteration are increased but simulation time is also increased. 113.021 MW in Table 9 indicates SFS can find the best solution with the setting of Npop = 60, 70, 80, 90, and 100, and MI = 150, and Npop = 50, 60, 70, 80, 90, and 100 and MI = 200. Clearly, Npop = 60 and MI = 150 are the best selection for the case of tuning population size and fixing the maximum iteration. SFS cannot improve result better if the population size is set to higher than 60 for MI = 150 and 70 for MI = 200. Similarly, 113.021 MW in Table 10 indicates that Npop = 50 and MI = 200 and Npop = 75 and MI = 150 are the best selection for finding the best performance of the SFS method. Improvement of the best solution fails if increasing MI to higher than 200 for Npop = 50 and to higher than 150 for Npop = 75. So, it can conclude that the best selection for population size and the maximum iterations are Npop = 60 and MI = 150. The results from the setting are reported for comparisons with other methods shown in Table 11. The saving percentage of SFS compared with other methods is from 0.09% to 9.1%, while SFS is only less effective than QOTLBO [28] by −0.66%. The comparison indication can conclude that optimal solution of SFS is the second best among all compared methods. Furthermore, search speed of SFS is in the fastest group once MI and Npop are, respectively, 150 and 60, while those from others are 200 and 40, 300 and 15, and 100 and 50. In summary, SFS can find optimal solutions with higher quality than most compared methods; however, compared with some best methods, SFS is also less potential. Consequently, we suggest SFS should be used as an optimization tool for the ORPF problem and it is more promising if SFS is improved by proposing modifications.

6. Conclusion

In the paper, we apply the SFS method for finding optimal solutions of the ORPF problem for different objectives consisting of power loss, voltage deviation, and voltage stability index. Two different power systems with 30 buses and 118 buses are employed for running the SFS method and results found by SFS together with control parameters are compared with those from other methods. As a result, SFS becomes one of the best methods searching the best optimal solutions for each case and its search speed is also faster than most methods. SFS can own the outstanding points, thanks to its construction consisting of three new solution generations, diffusion process, first update process, and the second update process. However, the performance of SFS still copes with constriction, leading to worse results than several methods about solution and speed. Thus, we suggest SFS should be used for finding solutions of the ORPF problem but modifications should be performed on the conventional SFS for improving the search ability.

Abbreviations

Nbus:Number of all buses in network
Vi, Vj:Voltage at bus i and bus j
Nload:Number of loads
Pdi, Qdi:Active power and reactive power of load i
Qci:Reactive power generation of capacitor i
NG:Number of available generators
QGi,min, QGi,max:Minimum and maximum limits of reactive power for generator i
VGi,min, VGi,max:Minimum and maximum voltage of generator i, respectively
Qci,min, Qci,max:Minimum and maximum reactive power output of capacitor at bus i
Nc:Number of capacitor buses
Ti,min, Ti,max:Minimum and maximum tap changer of transformer at bus i
Nt:Number of transformer buses
Vloadi,min, Vloadi,max:Minimum and maximum operating voltage of load i
Sl,max:Capacity of branch l.

Data Availability

The information of transmission lines and loads in IEEE 30-bus power transmission power network and IEEE 118-bus power transmission power network used to support the finding of this study have been taken from [38].

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was funded by Funds for the Science and Technology Development of the University of Danang, under project no. B2019-DN01-19.