4.2. Implementation of JS-MFO for Solving OPF Problems
Two test systems were used to investigate the aptitude of the proposed JS-MFO algorithm to reach the near-global OPF solution for different cases of mono-objective optimization problems: the IEEE 30-bus, a real part of the USA network, and the RIM 27-bus. The modeling of the tested networks was carried out in the MATLAB 2009b environment using a personal computer running Windows XP Professional, Pentium P-IV CPU 3 GHz processor, and RAM of 1 GB.
A. IEEE 30-bus test system
The IEEE 30-bus model has 6 generators, 41 branches (37 lines and 4 transformers with off-nominal tap ratios), and 24 load buses as shown in
Figure 2. Shunt Volt Ampere Reactive (VAR) compensators are installed in buses 10, 12, 15, 17, 20, 21, 23, 24, and 29, where the reactive power injection is controlled between 0 and 5 MVAR as the lower and upper limits, respectively [
11]. The total system demand was (2.834 + j1.262) p.u. for the apparent power at the 100 MVA base. Bus 1 was taken as the slack bus. The upper and lower limits of active power generation, the reactive power limits, and the generator cost coefficients are taken from [
34].
In the base case, the system presents the following data: Fc (quadratic fuel cost) = 901.935 USD/h, F
Em = 0.23908 ton/h, F
Plosses = 5.8332 MW, F
VD = 1.1663 p.u., and L
maxr = 0.1731. The parameters of the JS, MFO. and JS-MFO algorithms are as follows: the population number (population size) N
POP is 50, the maximum number of iterations
MaxIt = 100, and the number of control variables is
D = 25. For the present test system, six cases (six objective functions) were studied to evaluate the performance of the proposed JS-MFO compared to that of JS and MFO in order to reach the optimal solution. In order to validate the proposed JS-MFO in dealing with and solving complex OPF problems, 30 independent trial runs were performed for each case study by implementing JS, MFO, and JS-MFO. Three complex cases, particularly 1, 3, and 4, were considered to provide the minimum, average, maximum, and standard deviation (SD) for each objective function as shown in
Table 3. A comparative report of each index between the JS, MFO, and JS-MFO algorithms, shown in the previous table, can clearly prove the superiority of JS-MFO over JS and MFO in providing the best solution quality and stability over the trial runs. The optimal control variable settings and the corresponding objective function value are depicted for each case in
Table 4 and
Table 5. Executing 30 independent runs for each case led to extracting the solution for the run that provided the best solution.
Case 1: Minimization of the quadratic total fuel cost function
The objective function in Equation (12) was evaluated for the optimal settings of the control variables by running the JS, MFO, and JS-MFO algorithms, and the simulation results are indicated in
Table 4 for case 1. The optimal total fuel costs obtained were 799.085 USD/h, 801.02 USD/h, and 800.79 USD/h using JS-MFO, JS, and MFO, respectively. The best solution of the total fuel cost for JS-MFO was compared with that of the JS and MFO algorithms and other optimization techniques in the literature shown in
Table 6, such as the Novel Improved Social Spider Optimization algorithm (NISSO) [
35], Firefly Algorithm (FFA) [
15], Moth Swarm Algorithm (MSA) [
15], Efficient Sine Cosine Algorithm (ESCA) [
36], Symbiotic Organism Search algorithm (SOS) [
35], Harris Hawks Optimization (HHO) [
16], Dragonfly Algorithm (DA) [
13], Adaptive Gaussian Teaching–Learning-based Optimization (AGLBO) [
14], PSO-PS [
21], DA-PSO [
23], HF-PSO [
24], Artificial Bee Colony based on hybrid fruit fly (HF-ABC) [
25] and HFA-JAYA [
27]. It is very clear from
Table 6 that the JS-MFO approach is able to reduce the fuel cost to the lowest value of 799.085 USD/h. The percent cost savings with respect to the base cases of JS, MFO and JS-MFO are 11.188%, 11.2137%, and 11.4028 %, respectively, which lead to annual cost savings for JS-MFO of 16,951 USD/year compared to JS and 14,936 USD/year compared to MFO. The evolution of the total fuel cost during the simulation is indicated in
Figure 3. It can be seen that the optimal total fuel cost converges towards the best optimal solution using JS-MFO.
Case 2: Minimization of Total Gas Emissions
The total gas emission function in (13) was minimized using the proposed JS-MFO, MFO, and JS algorithms. The gas emission coefficients are given in [
34]. The optimal solutions achieved using JS-MFO, MFO, and JS are 0.2047, 0.2048, and 0.2049 ton/h, respectively, based on the optimal setting of the control variables shown in
Table 4 for case 2.
Table 7 gives a statistical comparison report considering other optimization methods applied to the same test system and the same number of control variables as the Harris Hawks optimization algorithm (HHO) [
16], Turbulent Flow Of Water-Based Optimization algorithm (TFWO) [
17], hybrid DA-PSO [
23], Adaptive Gaussian AGTLBO [
14], and Social Spider Optimization algorithm (SSO) [
37]. By examining the previous table, the minimum total gas emissions obtained by the proposed JS-MFO algorithm was shown to be better than that of the JS and MFO algorithms and the comparison optimization techniques reported in the literature. The percent gas emission reductions were 14.30%, 14.34%, and 14.38% for JS, MFO, and JS-MFO, respectively, referring to the base case. It can be seen that the proposed hybrid technique JS-MFO achieves total gas emission reductions of 1.752 tons/year and 0.876 tons/year compared to JS and MFO, respectively. The convergence of the total gas emission function with respect to the number of iterations is depicted in
Figure 4, showing that the proposed JS-MFO algorithm gives faster convergence to the optimal solution with a better solution quality than the JS and MFO algorithms.
Case 3: Minimization of total fuel cost considering valve point effect
The objective function given by (14) is provided for the non-smooth curve of the total fuel cost function and that for the two generators in buses 1 and 2, where their cost coefficients are taken from [
10]. The fuel cost curves of the remaining generators keep the same characteristics as in case 1. The minimum total fuel costs obtained using the proposed JS-MFO, JS, and MFO algorithms are 918.073 USD/h, 922.48 USD/h, and 924.02 USD/h, respectively, considering the valve point effect. The optimal solution reflects the optimal control variables presented in
Table 4 for case 3. The percent cost savings considering the valve point effect for JS, MFO, and JS-MFO are 15.01%, 14.87%, and 15.42%, respectively, compared to the base case. It is clearly shown that JS-MFO gives the highest cost saving. The annual cost savings of JS-MFO correspond to 38,632 USD/year and 52,122 USD/year with respect to those of JS and MFO, respectively. The JS-MFO algorithm shows a better solution than that of the JS and MFO algorithms, and the other methods displayed in
Table 8 such as the Improved Salp Swarm Algorithm (ISSA) [
38], Salp Swarm Algorithm (SSA) [
38], Teaching–Learning-based Optimization Algorithm (TLBO) [
39], Gbest-guided ABC (Gbest-ABC) [
40], and ABC-based Grenade Explosion Method (GABC) [
41]. The progress of the non-quadratic total fuel cost with respect to the iteration evolutions is given in
Figure 5 for the JS, MFO, and JS-MFO algorithms. These curves confirm the superiority of JS-MFO over JS and MFO in attaining the best global solution.
Case 4: Minimization of Total Active Losses
Active transmission losses in Equation (15) are adopted as an objective function in this case and are minimized by carrying out the proposed JS-MFO, JS, and MFO techniques. The optimal total active losses resulting from the simulation are given in
Table 5 for case 4, along with their corresponding control and state variables. It is reported that JS-MFO is able to reduce the active transmission losses to the lowest value of 2.8810 MW, while JS and MFO lead to 3.232 MW and 3.049 MW, respectively. The percent loss savings of the presented techniques JS, MFO, and JS-MFO are, respectively, 44.59%, 47.72%, and 50.60% with reference to the base case.
For comparison purposes,
Table 9 depicts the optimal active transmission losses for JS-MFO, JS, MFO, and other techniques inspired from previous works in the literature such as AGTLBO [
14], FFA [
15], MSA [
15], DE-HS [
22], DA [
23], DA-PSO [
23], HFA-JAYA [
27], SS0 [
35], NISSO [
35], ESCA [
36], and Adaptive Multiple-Team Perturbation-Guiding JAYA [
42]. From this table, it can be observed that the results obtained using JS-MFO are better than those of the other methods presented in the current literature.
Figure 6 shows the evolution of total active losses with respect to the progression of iterations during the simulation, depicting the convergence characteristics of the JS and JS-MFO algorithms for this case.
Case 7: Minimization of the total fuel cost along with the deviation of the total load bus voltage
This case study is devoted to the simultaneous minimization of the quadratic total fuel cost and the total load bus voltage deviation as the objective function depicted in Equation (21). The problem is solved by running the JS, MFO, and JS-MFO algorithms separately. The simulation results for case 7 are listed in
Table 5 by showing the numerical optimal settings of the control and state variables.
Table 10 is provided to compare the optimal total fuel cost and the total load bus voltage deviation of the JS-MFO method to those of the JS and MFO algorithms and the other techniques given by recently published works in the literature. It is clearly shown that JS-MFO dominates all other comparison techniques in terms of the optimal fuel cost with a slightly higher optimal voltage deviation compared to AGTLBO [
14], HFA-JAYA [
27], ESCA [
36], PSO-SSA [
37], MVO [
43] and ECHT-DE [
44]. The ESCA is reported to present the best optimal voltage deviation among all other techniques, but it provides the highest optimal voltage deviation.
It is observed from the results that JS-MFO provides simultaneous reductions in quadratic total fuel cost and total voltage deviation of 10.90% and 91.48%, respectively, with reference to the base case, while the other existing techniques, JS, and MFO show less important results. By examining the convergence curves of the JS and JS-MFO algorithms given by
Figure 7, it can be seen that the JS-MFO algorithm gives the best final global optimal solution compared to the JS and MFO algorithms.
Case 8: Minimization of Quadratic Total Fuel Cost Along With Vulnerability Stability Enhancement
The present case study examines the minimization of the total fuel cost along with the enhancement of voltage stability by optimizing the objective function mentioned in (22). The optimal settings obtained for the control and state variables after the JS, MFO, and JS-MFO algorithm running processes are illustrated in
Table 5 for case 8. The optimal solution presented by JS-MFO gives the best simultaneous reductions in total fuel cost and voltage stability index than that of JS and MFO. By referring to the base case, JS-MFO presents, simultaneously, a percent cost saving of 11.37% and an improvement of the stability margin of 34.08%. A comparison in
Table 11 with other meta-heuristics-based optimization techniques such as ESCA [
36], PSO-SSA [
37], AMTPG-JAYA [
42], MVO [
43], and ECHT-DE [
44] shows that the proposed JS-MFO presents the best global solution among others.
B. Mauritanian electric power system RIM 27-bus:
B.1. Simulation of various cases for RIM 27-bus
The RIM 27-bus power grid is controlled by the SOMELEC society, involving three voltage levels: 225 kV, 90 kV, and 33 kV. This system model consists of six generation buses: 1 (slack bus), 8, 10, 23, 24, and 25, feeding a total load of (1.786 + j0.736) p.u. at a base power of 100 MVA, and three shunt compensators at buses 12, 14, and 16. Two solar power plants are installed on buses 8 and 23 with active power capacities of 50 MW and 15 MW, respectively, and a wind power plant on bus 10 with a capacity of 30 MW, as shown in
Figure 8. Also, the RIM 27-bus system has 21 load buses and 36 branches. The voltage magnitude of the generation buses is controlled between 0.92 p.u. and 1.08 p.u. The number of control variables for this system is 15. In the base case, the RIM 27-bus system presents the following data: F
Plosses = 6.6015 MW, F
VD = 1.6358 p.u., and L
maxr = 0. 522.
Case 4: Minimization of Total Active Losses for RIM 27-bus system:
In this case study, the aim is to minimize the transmission active losses of the RIM 27-bus system, adopting (15) as the objective function, and the simulation results are tabulated in
Table 12. The optimal active losses affected by JS, MFO, and JS-MFO are 3.35 MW, 3.56 MW, and 3.25 MW, which represent loss percent savings of 49.19%, 46.04%, and 50.67% compared to the base case, respectively. The best solution quality was found with the JS-MFO technique, which yielded the lowest total active power losses and the highest loss percent saving. The optimal total active transmission losses attributed to JS-MFO are 2.94% and 8.58% lower than those obtained by the JS and MFO algorithms. It is clearly illustrated in
Figure 9 that the JS-MFO algorithm dominates the JS and MFO algorithms in terms of solution quality and convergence characteristics.
Case 5: Minimizing of the total voltage deviation for the RIM 27-bus system:
The objective function presented in (16) is minimized in case 5 and depicts an optimal solution for the control and state variables shown in
Table 12. The RIM 27-bus power system security is enhanced via the minimization of the total voltage deviation F
VD using JS, MFO, and JS-MFO, signaling reductions of 59.97%, 55.25%, and 62.44%, respectively, compared to the base case. The JS-MFO pointed out the best optimal F
VD compared to that of JS and MFO. The optimal total voltage deviation achieved by JS-MFO is 0.6144 p.u., which is better than that of the JS and MFO algorithms, which correspond to 0.6547 p.u. and 0.732 p.u., respectively. The best performance of the JS-MFO algorithm is verified in
Figure 10, illustrating the convergence proprieties of JS-MFO overcoming those of the JS and MFO algorithms. Case 6: Voltage stability enhancement for the RIM 27-bus system:
The FLmax voltage stability index is optimized for case 6 in order to enhance voltage stability using the JS, MFO, and JS-MFO algorithms. The stability margin of the RIM 27-bus system is improved using JS, MFO, and JS-MFO, where FLmax reaches the optimal values of 0.3861, 0.4201, and 0.3389, respectively. The JS-MFO approach denoted a percent reduction of FLmax of 12.22% compared to JSO and 19.33% compared to MFO, which proves the capacity of JS-MFO to overcome the JS and MFO algorithms.
It appears from
Table 12 that the JS-MFO algorithm outperforms the JS and MFOs in solving the OPF problem. The JS-MFO simulation result provides 12.22% and 19.33% improvements in the stability margin compared to the JS and MFO algorithms, respectively.
B.2. Sizing of reactive power compensators for renewable sources
Table 13 points out the reactive power generation for each optimal state of cases 4, 5, and 6. The reactive power compensator sizing for each renewable energy source in generation buses 8, 10, and 23 is, respectively, −17.5/17.5 MVAR, −16/16 MVAR, and −15/15 MVAR.
Figure 11 indicates the voltage magnitude profile for each optimal state for cases 4–6, describing that the voltage magnitude is within the predefined limits between 0.92 p.u. and 1.08 p.u.