The tests were conducted mainly in two balanced systems (33- and 69-bus systems) and two unbalanced systems (19- and 25-bus system). These four cases were first carried out with a fixed power demand and then with a variable demand. The voltage limits constraint was established between 0.93 p.u. and 1.05 p.u. The voltage unbalance index (VUI) and current unbalance index (CUI) for unbalanced systems were established as 3% [
41] and 30% [
38], respectively, with a specific analysis for the influence of these indexes conducted. To assess the results found using the SBAT, a comparison with the traditional SPSO [
46] and with a version of the selective harmony search (SHS) algorithm presented in [
48] was carried out. Tests considered 100 runs for each system, with the exception of SHS for the 69-bus system with variable demand. In this case, only 25 tests were performed due to the high computational time required for each run. In addition, a brief analysis of DNR considering VUI and CUI on a modified 123-bus system was conducted and presented. All simulations were run on a notebook equipped with an Intel-Core i7 (California, US, Intel Corporation,) processor running on a Microsoft-Windows 10 (Washington, US, Microsoft Corporation) environment.
5.1. Algorithm Parameters
One of the crucial points for extracting better performance of bio-inspired metaheuristics is the setting of parameters. The SBAT presents a considerable amount of parameters to be adjusted, which may require some time to be defined properly. In addition to the proper algorithm parameters are those related to iterations and population, also needing a correct adjustment. Also, the variable demand cases must present the characteristics of all loads considered in the system, together with the hourly costs.
Regarding the SBAT, the parameters set were:
,
,
,
,
and
. The parameters
and
were empirically set through tests with values varying from zero to one, setting them to the best values found, i.e., 0.8. Regarding the values of
,
,
and
, these were set as those usually adjusted for most applications, as defined in [
45] and cited in
Section 3. The SPSO parameters, learning factor 1
, global learning factor 2
, maximum inertia weight
, minimum inertia weight
and maximum velocity
V, were set as the default parameters of the PSO algorithm [
49]. As for SHS parameters, harmony memory considering rate
, pitch adjusting rate
and bandwidth
, they were set according to the values presented in [
50]. The only exception is the number of improvisations
, defined in the same way as the empirical test for the SBAT. These values were used in all simulations performed and are all described in
Table 1.
Population size, or harmony memory size (HMS) in the case of HS, and number of iterations for all algorithms were considered variable according to
, described in
Section 3. These values were defined to achieve a balance between convergence and computational time, aiming to extract a better performance from both points of view.
Table 2 presents these values for the four main systems tested.
As for the variable demand characteristics, the definitions found in [
22,
23] are followed and presented in
Table 3. Residential, commercial and industrial loads were respectively identified as 1, 2 and 3. For each load type defined in
Table 3, a cost (USD/kW) was associated with the respective hour together with a load factor. Each load located at a bus was characterized with a type described in
Table 3. The selection for all systems, except the 33-bus system, was made through the Python function
random.choice(), weighted with probabilities of 60%, 25% and 15% for choosing residential (Type 1), commercial (Type 2) or industrial (Type 3), respectively. The 33-bus system load distribution was taken from [
51] for a comparative basis. For the unbalanced systems, the load types per phase in each bus were the same, e.g., bus 1 loads at phases A, B and C were residential. Finally, to compute the total losses costs when dealing with variable demand, simulations were performed considering 24 h for each configuration found, where the real losses costs in the current topology were determined at each specific hour of the day. The total losses costs in one day are represented through the sum of each of these costs in all hours. This procedure results in a final topology for one day, thus avoiding a constant state-changing of the system switches and preventing their wear out. With these definitions, it is possible to simulate DNR for the proposed systems.
5.2. Balanced Systems
Firstly, a 33-bus and 37-branch balanced system was considered, presenting a base voltage of 12.66 kV. All system data are found in [
47].
Initially considering a fixed demand and characteristics provided in Reference [
47], the algorithm was tested in 100 runs and results are presented in
Table 4.
Figure 2 shows the voltage levels at each bus before and after the reconfiguration.
Results in
Table 4 show a 31.08% reduction in real losses, comparing initial and final distribution network configuration. From
Figure 2, it is possible to assume that voltage levels were improved compared to the initial configuration, going from a minimum of 0.912 p.u. at bus 18 to 0.936 p.u. at bus 32. All results obtained are compatible with the technical literature [
8,
15,
22,
52].
In
Table 4, a comparison between the SBAT, SPSO and SHS is also presented, showing which of the three found the best result. The SBAT overcame SPSO and SHS with a convergence rate of 96% and average losses of 139.74 kW.
The results provide a basis to consider the variable demand with the premises established in [
22,
23]. It is necessary to define the types of loads used in daily simulations. For the 33-bus case, this definition was (type of load/bus allocated): residential (2 to 4, 8 to 9, 11, 14 to 21, 23 to 25, 29, 31, 33), commercial (6 to 7, 10, 13, 22, 26 to 28, 30) and industrial (5, 12, 32).
Table 5 describes results for variable demand.
Figure 3 indicates minimum voltage levels for each hour along the day.
It is possible to verify from results in
Table 5 that the new configuration reduced the total losses costs by 31.57% in comparison with the initial topology. In addition, it is possible to notice that the result is different from the one found considering fixed demand, pointing to opening switches 9-7-14-28-32 instead of 9-7-14-32-37, showing the importance of considering a variable load demand during simulations when searching for an optimum DNR. The minimum voltage levels also improved, going from initial 0.926 p.u. minimum at bus 18 (20 h) and 0.981 p.u. maximum-minimum at bus 18 (4 h) to final 0.951 p.u. minimum at bus 33 (20 h) and 0.987 p.u. maximum-minimum at bus 32 (3 h). All results are compatible with the results presented in [
22,
23] for the same system.
A comparison presented in
Table 5 shows that the best solutions were also found by SPSO and SHS, but they had a worse performance in comparison with the SBAT. Selective PSO reached the best solution in 62% of all runs and an average total losses cost of 131.60 USD/kW. Selective HS reached the best solution in 72% of all runs and an average total losses cost of 130.02 USD/kW.
Figure 4 and
Figure 5 shows the final configuration and the open switches for fixed demand and variable demand respectively.
The second balanced system presents 69 buses and 73 branches and is widely used in the literature. It was firstly presented in [
53]. The base voltage of the system is 12.66 kV.
Table 6 describes the results obtained for the SBAT.
Figure 6 illustrates the respective initial and final voltage levels for each bus.
All results found for this system (losses and voltage levels) in
Table 6 and
Figure 6 point to an improvement in the operative indexes after the reconfiguration. The losses reduced by 56.19% and the minimum voltage level increased from 0.907 p.u. at bus 69 to 0.948 p.u. at bus 65. Again, the results found are compatible with those in the specialized literature [
54].
This case presents a particularity where different configurations, such as opening switches 55-61-70-69-14 or 56-61-70-69-14, show approximately the same percentage reduction and voltage levels as switches 58-61-70-69-14.
Table 6 also shows that the SBAT, SPSO and SHS reached the same solution, although the SBAT presented a better convergence to the best result (96%) and average final losses of 98.96 kW considering all runs. Selective PSO reached the best result in 31% of runs and an average final losses value of 108.07 kW. Selective HS reached the best result in only 8% of runs, although with a better average final losses value than SPSO (100.60 kW).
For the same system, now considering variable demand, the loads were defined as follows (load type/bus allocated): residential (9 to 10, 13 to 14, 16, 18, 20, 24, 26 to 29, 34, 39, 43, 45-46, 48, 50 to 51, 61 to 62, 64, 66 to 67, 69), commercial (7 to 8, 11, 17, 21, 33, 35, 37, 40, 49, 52, 55, 65, 68) and industrial (6, 12, 22, 41, 53 to 54, 59). As it is possible to notice, some buses are not connected to any type of load. Results obtained using the SBAT for the 69-bus balanced system with variable demand are presented in
Table 7 and
Figure 7.
In the same way as the 33-bus test system, results show improvements in several indexes, such as total losses costs and minimum voltage profile. The first presented a decrease of 54.56% and the second an improvement, going from initial 0.915 p.u. minimum at bus 65 (20 h) and 0.981 p.u. maximum-minimum at bus 65 (4 h) to final 0.951 p.u. minimum at bus 61 (20 h) and 0.989 p.u. maximum-minimum at bus 61 (4 h). For this system, the set of switches defined as the best by the SBAT was the same as the one presented for the fixed demand case (58-61-70-69-14). There is no comparative basis in the technical literature in order to compare the obtained results to other results, taking into account the same premises in this research. The final configuration both for fixed and variable demand with its open switches is shown in
Figure 8.
As presented for fixed demand,
Table 7 shows also the results found by SPSO and SHS considering the same condition. Again, SPSO was able to find the best result, but with a convergence rate to this result of 27% and average total losses costs of 93.23 USD/kW, lower than the results found by the SBAT, which were respectively 95% and 85.41 USD/kW. Selective HS also found the best result, with a lower convergence rate to the best (8% for only 25 runs), but with a higher average total losses costs (86.34 USD/kW) than SPSO.
5.3. Unbalanced Systems
The study of unbalanced systems with variable demand considering unbalance indexes aims to better represent the characteristics of a real distribution system. The first unbalanced case studied was a 19-bus and 20-branch system with different load levels in each phase characterizing the imbalance [
55]. The voltage base was 4.16 kV.
Initially, a fixed power demand was considered.
Table 8 describes results for 100 runs,
Table 9 shows the losses per phase before and after reconfiguration, and
Figure 9 shows the voltage levels for each bus and each phase before and after reconfiguration.
The results found by the SBAT for this first test indicate opening switches 10-11, with an active losses reduction of 38.98% and a reduction in losses per phase. Voltage levels presented in
Figure 9 show an improvement from initial to final configuration, from 0.950 p.u. minimum at bus 19 phase B and 0.999 p.u. maximum at bus 1 phase B to 0.969 p.u. minimum at bus 13 phase C and 0.999 p.u. maximum at bus 1 phase B. These results are compatible with those obtained in the literature [
55]. For comparison, SPSO and SHS found the same results.
The solution found attends all limits established, i.e., voltage limits, VUI and CUI. The maximum values obtained for VUI and CUI through the solution were 0.42% and 25.31%, respectively.
If the CUI value of 25.31% is considered high, the minimum CUI limit that presents a feasible solution is 21%. Now, the new configuration with minimum losses that attends all constraints indicates opening switches 9-20, with a minimum and maximum voltage of 0.948 p.u. at bus 19 phase C and 0.999 p.u. at bus 1 phase B, maximum VUI of 0.24% and maximum CUI of 20.63%, reducing the total losses only by 8.85% with a final value of 12.05 kW.
For a variable power demand, the loads were considered as the same type for all phases in each bus and as the following (type of load/bus allocated): residential (2 to 4, 7 to 9, 10, 15 to 16), commercial (5 to 6, 12 to 13, 17) and industrial (1, 11, 14, 18).
Results are presented in
Table 10, with the total losses costs (USD), open switches, percentage costs reduction (%), average costs (USD) and convergence rate (%) for 100 runs.
Table 11 shows the total losses costs per phase.
Figure 10 illustrates the minimum voltage profile for 24 h.
It is possible to infer from
Table 10,
Table 11 and
Figure 10 that there was an improvement comparing the final with the initial configuration. The total costs reduced by 37.49%, indicating opening switches 13-11, and the total costs per phase reduced from 3.69 USD/kW, 3.69 USD/kW and 3.85 USD/kW in phases A, B and C, respectively, to 2.4 USD/kW, 2.14 USD/kW and 2.48 USD/kW. Similar to the 33-bus balanced case, simulations considering variable demand demonstrated different solutions in comparison with the fixed demand case, which indicated opening switches 10–11.
As for the minimum voltage levels, the 19-bus system’s minimum and minimum-maximum rose respectively from 0.963 p.u. at bus 19 phase C (12 h) and 0.991 p.u. at bus 13 phase C (3 h) p.u. to 0.976 p.u. and 0.994 p.u. The final minimum and minimum-maximum values were located at bus 16 phase B (19 h) and at bus 13 phase C (3 h). SPSO and SHS found the same results, as shown in
Table 10.
The maximum VUI and CUI values obtained when opening switches 13 and 11 were respectively 0.21% and 26.02%, attending the established limits. Again, the value of VUI is low, but CUI may be considered high. The minimum value of CUI to obtain a feasible solution for this system with variable demand is 22%, which indicates opening switches 9–20, obtaining a maximum absolute VUI of 0.21%, a maximum absolute CUI of 21.18%, a minimum voltage of 0.960 p.u. at bus 19 phase C (12 h) and a maximum-minimum voltage of 0.991 p.u. at bus 16 phase C (3 h), with total losses costs of 10.22 USD/kW for a 24 h period.
Figure 11 and
Figure 12 shows respectively the final configuration for the fixed and variable demand.
The second unbalanced system tested is composed of 25 buses and 27 branches, as presented in [
55,
56,
57], with a base voltage of 4.16 kV and imbalance between phase loading. Results are presented in
Table 12,
Table 13 and
Figure 13.
Results for all algorithms (SBAT, SPSO and SHS) show the decrease in loss levels (10.79%) when comparing the initial open switches, 25-26-27, with those determined after the DNR, 22-15-17, as well as a decrease in losses per phase. In addition, there was an improvement in the general voltage levels, with the minimum value increasing from 0.927 p.u. at bus 12 phase B to 0.937 p.u. at bus 13 phase A. The maximum voltage level was 0.998 p.u. at bus 1 phase A in the initial and final configurations. The results are compatible with those found in [
55,
56,
57].
The maximum VUI and CUI values for the solution 22-15-17 are respectively 0.65% and 12.84%, attending the established limits of 3% and 30%. In the same way as the 19-bus system, if the CUI is considered high, a new limit must be set. The minimum limit value that presents a feasible solution is 12.83%, indicating opening switches 22-26-15, reducing losses to 141.85 kW and presenting minimum voltage of 0.937 p.u. at bus 13 phase A, maximum voltage of 0.998 p.u. at bus 1 phase A and maximum VUI and CUI of 0.65% and 12.828%, respectively. The VUI and CUI values were slightly lower than the original solution.
For the variable demand case, the loads were defined as the same type for all phases in each bus as follows (type of load/bus allocated): residential (1, 5, 7, 9 to 11, 16 to 18, 21, 23, 25), commercial (15, 20, 22, 24) and industrial (3 to 4, 6, 8, 12 to 13, 19).
Table 14 and
Table 15 and
Figure 14 present the results, total losses costs per phase and minimum voltage levels for 24 h, respectively.
Table 14 shows that the results for the three algorithms (SBAT, SPSO and SHS) indicate the same set of opened switches in the fixed demand case (25-15-17). The total losses costs diminished from 131.58 USD to 117.51 USD, representing a reduction of 10.71%.
Table 15 shows an improvement in the total losses costs per phase for a 24 h period. The minimum voltage levels in
Figure 14 also improved, from 0.947 p.u. minimum at bus 12 phase B (12 h) and 0.986 p.u. minimum-maximum at bus 15 phase A (3 h) to 0.951 p.u. minimum and 0.987 p.u. minimum-maximum. The final minimum and minimum-maximum values were located at bus 13 phase A (20 h) and at bus 13 phase A (3 h).
Figure 15 brings the final configuration for the 25-bus system considering both fixed and variable demand.
The maximum VUI and CUI presented for the solution 25-15-17 with variable demand were 0.46% and 13.41%, attending the limits established. For CUI, if the limit is considered high, the minimum limit value that presents a feasible solution reducing losses is 13%, which indicates opening switches 22-17-14, with a maximum VUI and CUI of 0.49% and 12.95%, respectively, minimum voltage of 0.945 p.u. at bus 13 phase B (12 h) and maximum-minimum voltage of 0.986 p.u. at bus 13 phase A (3 h), reducing losses costs to 119.39 USD/kW.
5.4. Additional Test
Additionally, a modified 123-bus system, presented in
Figure 16, was tested. This system was originally composed of voltage regulators, capacitor banks, single and three-phase sections, different types of cables and underground sections, among other characteristics described in [
58]. To assess the system only from the reconfiguration point of view, the system was modified, excluding all equipment.
However, when considering both VUI and CUI limits, this system did not present a feasible solution, as, although VUI presented values under 3%, CUI surpassed the value of 30% by a larger margin, as some sections of the system presented currents of low value (almost zero). For example, consider branch 22 of the modified system, which presented currents of magnitude 54.775 A in phase A, 0.004 A in phase B and 58.477 A in phase C. This happens also in this same branch of the system in its original state (considering all equipment) [
58], with currents of magnitude 55.853 A in phase A, 0.005 A in phase B, and 55.404 A in phase C, which points out a highly unbalanced characteristic from the point of view of currents.
A solution for both fixed and variable demand is only found when disregarding CUI constraint in the same way as the majority of papers dealing with DNR and unbalanced systems present. In addition, the minimum voltage level allowed, which is established as 0.93 p.u. in this paper for all systems, must be set to 0.9 p.u. for this system in fixed demand to find a minimum losses solution.
For fixed demand, the algorithm indicates opening switches 118-93 with losses of 105.68 kW, as opposed to the initial open switches 125-126 with losses of 109.12 kW. The voltage levels also improved, going from 0.912 p.u. minimum at bus 114 phase A to 0.918 p.u minimum at bus 114 phase A and maintaining 0.999 p.u as maximum at bus 150 phase B.
As for variable demand, the loads were divided as follows: residential (1, 2, 5, 6 to 9, 12, 19, 20, 28 to 30, 32 to 35, 37, 38, 42, 43, 45 to 50, 53, 55, 58, 59, 62, 65, 66, 73, 76, 83, 85 to 88, 90, 92, 94, 96, 100, 103, 104, 106, 107, 109, 111 to 113), commercial (4, 7, 17, 24, 52, 60, 63, 64, 70, 74, 79, 95, 98, 99) and industrial (16, 22, 31, 39, 41, 51, 56, 68, 69, 71, 75, 80, 82, 84, 102, 114). The algorithm indicates opening switches 125-93 with total losses costs of 95.6 USD as opposed to the initial set of open switches 125-126 with costs of 96.5 USD. The voltage levels, minimum and minimum–maximum, also improved, going from 0.928 p.u. to 0.931 p.u. at bus 114 phase A (20 h) and from 0.980 p.u. to 0.981 p.u. at bus 150 phase B (3 h).
This example helps to emphasize the influence of the unbalance indexes, especially CUI, when finding a solution for the DNR problem in unbalanced systems, as these constraints, when considered, can directly determine a feasible or unfeasible solution to the problem.
5.5. VUI and CUI Influence Analysis
The influence of VUI and CUI can be summarized in
Table 16 for the two main unbalanced systems studied (considering fixed and variable demand).
The results consolidated in
Table 16 show the direct influence of VUI and CUI in DNR results for these two systems. When considering the limits established for all systems, i.e., 3% for VUI and 30% for CUI, the results are the same as those found when disregarding these two constraints.
However, when slightly decreasing the CUI limits allowed in all systems to the minimum values that present a feasible solution, new topologies were found in comparison with those for the limit of 30%, hence reducing the current unbalance. These different solutions were seen for both demand cases in some systems, for example, the 25-bus unbalanced system, which presented three different solutions depending on the value of CUI considered and the demand level studied. In some cases, the value of VUI also decreased alongside the new value of CUI established, as in the case of 19-bus unbalanced system with fixed demand. In addition, the feasible solution with the minimum CUI increased the values of losses and losses costs in the majority of cases.
The presented analysis indicated that the DNR in unbalanced systems is indeed sensible to these indexes, particularly CUI, as it is the index that presented higher values in all systems, greater affecting the results found. This is clearly seen for the additional 123-bus system presented, as it is difficult to find a feasible solution when CUI is considered in the problem. Although all systems analysis, including that of the 123-bus system, could be easily carried out only in the VUI context or without these two constraints, as results show, the assessment of CUI may be important in some cases, as high values of CUI can cause the earlier stated difficulties related to protection [
37] and to power and energy losses [
38], emphasizing the importance of the assessment of this index in DNR problems.