Next Article in Journal
Blue Carbon: Comparison of Chronosequences from Avicennia marina Plantation and Proteresia coarctata Dominated Mudflat, at the World’s Largest Mangrove Wetland
Next Article in Special Issue
Optimal Feeder Reconfiguration and Placement of Voltage Regulators in Electrical Distribution Networks Using a Linear Mathematical Model
Previous Article in Journal
Spatial-Temporal Evolution Characteristics of Agricultural Intensive Management and Its Influence on Agricultural Non-Point Source Pollution in China
Previous Article in Special Issue
An Energy Management System for PV Sources in Standalone and Connected DC Networks Considering Economic, Technical, and Environmental Indices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Phase-Balancing in Three-Phase Distribution Networks Considering Shunt Reactive Power Compensation with Fixed-Step Capacitor Banks

by
Daniel Federico A. Medina-Gaitán
1,
Ian Dwrley Rozo-Rodriguez
1 and
Oscar Danilo Montoya
1,2,*
1
Grupo de Compatibilidad e Interferencia Electromagnética (GCEM), Facultad de Ingeniería, Universidad Distrital Francisco José de Caldas, Bogotá 110231, Colombia
2
Laboratorio Inteligente de Energía, Facultad de Ingeniería, Universidad Tecnológica de Bolívar, Cartagena 131001, Colombia
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(1), 366; https://doi.org/10.3390/su15010366
Submission received: 24 November 2022 / Revised: 21 December 2022 / Accepted: 22 December 2022 / Published: 26 December 2022
(This article belongs to the Special Issue Smart Grid Technologies and Applications)

Abstract

:
The black hole optimization (BHO) method is applied in this research to solve the problem of the optimal reactive power compensation with fixed-step capacitor banks in three-phase networks considering the phase-balancing problem simultaneously. A master–slave optimization approach based on the BHO in the master stage considers a discrete codification and the successive approximation power flow method in the slave stage. Two different evaluations are proposed to measure the impact of the shunt reactive power compensation and the phase-balancing strategies. These evaluations include a cascade solution methodology (CSM) approach and a simultaneous solution methodology (SSM). The CSM approach solves the phase-balancing problem in the first stage. This solution is implemented in the distribution network to determine the fixed-step capacitor banks installed in the second stage. In the SSM, both problems are solved using a unique codification vector. Numerical results in the IEEE 8- and IEEE 27-bus systems demonstrate the effectiveness of the proposed solution methodology, where the SSM presents the better numerical results in both test feeders with reductions of about 32.27 % and 33.52 % , respectively, when compared with the CSM. To validate all the numerical achievements in the MATLAB programming environment, the DIgSILENT software was used for making cross-validations. Note that the selection of the DIgISLENT software is based on its wide recognition in the scientific literature and industry for making quasi-experimental validations as a previous stage to the physical implementation of any grid intervention in power and distribution networks.

1. Introduction

1.1. General Context and Motivation

In the last two decades, population growth and technological advances have caused a considerable increase in the use of electrical energy at all voltage levels for residential, industrial, and commercial applications [1,2]. Because they directly connect end users and the electricity service, electrical distribution networks are the systems with the most accelerated growth compared to large-scale power and generation systems [3,4]. Due to the operating values of medium-voltage distribution networks, typically between 1 and 25 kV, these networks exhibit higher power loss percentages than transmission and sub-transmission systems. These energy losses in distribution grids can be between 6 and 18% of the energy purchased at the terminals of the substation, whereas transmission systems show values between 1.5 and 2.5% [5].
On the other hand, given the policies implemented by the regulatory entities of the electrical sector, the values regarding energy losses in distribution networks must be continuously reduced by utilities to improve the quality, and the distribution efficiency of electrical energy [6,7]. To this effect, utilities need to design efficient maintenance and operation plans that allow for reaching the expected energy losses with minimal investment costs [8,9].
Due to the diversity of alternatives to improve the efficiency of the electrical service, distribution companies can select one or more of the following options to reduce their energy losses in distribution systems. The first alternative, valid for three-phase networks, involves using shunt reactive power compensators with fixed-capacitor banks or static var compensators [10,11]. The second option is using dispersed generation sources and battery energy storage systems. However, their costs are very high when compared to the first alternative, and the main application of distributed sources and batteries is related to active energy support to reduce energy purchasing costs for planning periods that oscillate between 5 and 20 years [12], not to reduce the energy losses during distribution activity. A detailed study regarding the optimal integration of renewable generation in distribution networks using multi-objective optimization was presented by the authors of [13], which constitutes an essential reference to understand the importance of having a multi-criteria decision algorithm to install dispersed generation sources in hybrid AC–DC distribution grids from economic and technical perspectives. The third option is an efficient grid reconfiguration, i.e., modifying the grid topology by using available tie-lines [14]. The final alternative corresponds to the optimal phase balancing, i.e., the redistribution of the load connections at all the network nodes to reduce excessive voltage drops in charged phases and the magnitude of the current in these phases [5].
This research proposes an efficient alternative to minimize the total grid power/energy losses based on the above-mentioned options. To this effect, the two cheapest options were selected, i.e., the optimal selection and location of fixed-step capacitor banks combined with optimal phase-swapping in all the network nodes.

1.2. State of the Art

In the literature, many approaches focus on integrating fixed-step capacitor banks into electrical distribution networks and solving the optimal phase-swapping problem in three-phase asymmetric networks. This subsection presents some of the most recent works in those areas.

1.2.1. Optimal Placement of Fixed-Step Capacitor Banks

The authors of [15] presented a solution methodology for the optimal selection and location of fixed-step capacitor banks in medium-voltage distribution grids, intending to reduce grid power losses and improve the voltage profiles. The proposed methodology was based on applying the crow search algorithm, a bio-inspired combinatorial optimization methodology. Numerical results were obtained in two test feeders composed of 9 and 33 buses, with better solutions when compared to the classical particle swarm optimization (PSO) method. The work by [16] presented an interesting case study associated with installing capacitor banks in a distribution network that provides energy service to the Tehran metro. The main characteristic of this approach is the use of the ETAP software and its optimization tool concerning the optimal placement of capacitor banks through the implementation of a specialized genetic algorithm, which minimizes the total investment costs in compensators for an expected analysis period. Numerical results showed that, for a planning period of 5 years, the cumulative net profit regarding the reduction of energy losses costs while considering the capacitors’ investment and operating costs amounts to more than 300,000 dollars. These results revealed the positive impact of using fixed-step capacitor banks to improve performance in electrical distribution networks with bad lagging power factors. The work by [10] proposed the application of the whale optimization algorithm to locate and select fixed-step capacitor banks in radial distribution grids. The optimization process considered two objective functions regarding operating cost reduction and power loss minimization. The IEEE 34-bus grid and the IEEE 85-bus grid were selected as test feeders to validate the effectiveness of the proposed optimization approach in comparison with the methods reported in the literature, such as the PSO method, the plant growth simulation algorithm, and the bacterial foraging optimization algorithm, among others. The authors of [17] presented the application of the flower pollination algorithm to locate and select fixed-step capacitor banks in radial distribution networks to minimize the annual costs of energy losses while including the investment costs of the capacitor banks. Numerical results in the IEEE 33-, 34-, 69-, and 85-node grids demonstrated the effectiveness of the proposed optimization algorithm when compared to an analytical method and an improved fuzzy-logic genetic algorithm. The study by [18] applied the cuckoo search algorithm to locate and size fixed-step capacitor banks in radial distribution networks. The IEEE 34- and 69-bus grids were employed in all numerical validations. The objective was to minimize the total grid power losses and improve the voltage profiles along the distribution feeder. A comparative analysis with the classical PSO approach and the plant growth simulation algorithm demonstrated the effectiveness of the cuckoo search algorithm in solving the studied problem. Additional optimization algorithms to locate and capacitor banks are the hybrid honey bee colony algorithm [19], the tabu search algorithm [20,21], the vortex search algorithm [22], and the gravitational search algorithm [23], among others. In the case of three-phase asymmetric distribution grids, the authors of [24] presented the application of the imperialist competitive algorithm to locate and size capacitor banks while considering three-phase networks with harmonic pollution. The work by [25] presented the effect of an unbalanced grid regarding the presence of single-, two-, and three-phase loads on the final location and size of fixed-step capacitor banks while considering the minimization of the total grid power losses and the reduction of voltage imbalances.
The works mentioned above have two main characteristics: (i) the use of combinatorial optimization techniques (i.e., metaheuristics) to deal with the nonlinearities and non-convexities of the exact mixed-integer nonlinear programming model that represents the studied problem; and (ii) the fact that most of the optimization approaches use single-phase equivalents for the distribution grid modeling. This means that more research is required regarding asymmetric distribution networks.

1.2.2. Optimal Phase-Balancing in Three-Phase Asymmetric Grids

The authors of [26] presented the application of a specialized Chu and Beasley genetic algorithm to solve the problem regarding optimal phase-balancing in three-phase asymmetric networks with multiple constant power loads. Numerical results in the IEEE 37-bus grid demonstrate the effectiveness of the proposed genetic algorithms. However, no comparisons with other optimization methods were reported. The work by [27] presented the application of a mixed-integer approximation to solve the optimal load balancing problem in three-phase distribution networks while using a current approximation method. The study by [5] proposed the application of the vortex search algorithm to solve the optimal phase-balancing problem in three-phase distribution networks while aiming to minimize the expected value of the grid power losses. The IEEE 8-, 25-, and 37-bus grids were used as test feeders. Numerical comparisons with the classical Chu and Beasley genetic algorithm demonstrated the effectiveness of the proposed vortex search algorithm in minimizing the objective function. The authors of [28] presented a heuristic optimization algorithm based on measuring the phase current at the point of load connection. With these measurements, each load is reconfigured to minimize the expected current imbalance. Numerical results and comparisons with the phase commitment algorithm and the modified leap frog optimization method confirmed the proposal’s effectiveness in test feeders composed of 8, 15, and 30 asymmetric loads. The study by [29] presented the solution of the optimal phase-balancing problem in three-phase asymmetric networks operated at low voltage levels (i.e., secondary distribution networks). The IEEE 13-node test feeder was selected to evaluate the heuristics-based solution proposal. However, no comparisons with literature reports were presented. The authors of [30] proposed applying artificial intelligence-based methods using artificial neural networks that consider information provided by smart meters. The main contribution of this paper was the application of neural networks to an actual distribution feeder belonging to the Irbid district electricity company. The work by [31] proposed applying the hurricane optimization algorithm to solve the optimal phase-swapping problem in three-phase asymmetric networks. The objective function corresponded to minimizing the total grid power losses under peak load conditions. The IEEE 8-, 25-, and 37-bus grids were employed for numerical validations. A comparative analysis with the Chu and Beasley genetic algorithm and the vortex search algorithm demonstrated the effectiveness of the proposed optimization method regarding the final power loss value.
Note that most of the optimization methods mentioned above for solving the problem under study are based on applying combinatorial optimization algorithms (metaheuristics) to obtain a high-quality solution, which implies that the application of recently developed combinatorial algorithms shows excellent promise for research in this field, as is the case of the black hole optimization (BHO) technique used in this work.

1.3. Contribution and Scope

Considering the above, no solution methodologies in the current literature simultaneously address the problems regarding phase-balancing and reactive power compensation. Thereupon, this work makes the following contributions:
  • A complete comparative analysis regarding reactive power compensation with fixed-step capacitor banks and the phase-balancing approach to reduce peak power losses in distribution grids by presenting cascade and simultaneous solution methodologies;
  • The application of the BHO method to determine the size and location of the capacitors, as well as the best set of load connections in three-phase nodes, using a discrete codification.
It is worth mentioning that this research only considers an objective function based on the minimization of a technical aspect of the distribution grid, i.e., the expected power losses under peak load conditions. In addition, the selection of test feeders for validating the proposed optimization method is based on the literature reports where the studied problems were solved separately. In addition, note that the selection of the BHO approach is based on simple evolution rules and is highly efficient in solving combinatorial optimization problems, added to the fact that, according to our exploration of the literature, this algorithm has not been applied to solve both studied optimization problems using cascade or simultaneous solution methodologies, which has been identified as a research gap to be filled by this work. On the other hand, it is essential to acknowledge the strong connection between sustainability and the issue under study, given that part of the electrical energy converted into heat in distribution lines is obtained from fossil fuels. Reducing this makes distribution networks more sustainable, as less damaged energy resources are required to provide electrical energy to all end users connected to those grids.

1.4. Document Structure

The remainder of this research document is structured as follows: Section 2 presents a general formulation of the load flow problem for three-phase asymmetric distribution networks, which is based on the three-phase version of the successive approximations power flow method; Section 3 presents the general characteristics of the mathematical models regarding the optimal location of the fixed-step capacitor banks and optimal phase-swapping in three-phase grids; Section 4 describes the general aspects of the BHO technique by presenting its general conception, mathematical formulation, and algorithmic implementation; Section 5 shows the main characteristics of the distribution systems under analysis, as well as the parametrization regarding the sizes of the fixed-step capacitor banks; Section 6 presents all the numerical validations in the 8- and 25-bus grids, with a complete analysis and discussion; and Section 7 lists the main conclusions derived from this work, as well as some possible future works.

2. Load Flow Analysis

This section presents the main aspects of the power flow problem for three-phase unbalanced networks. A detailed description of a recently developed power flow approach based on the successive approximations power flow method is provided, as it is at the heart of any combinatorial optimization method applied to improve the performance of electrical networks.
The load flow problem is one of the most classical problems in the field. It has been studied for six decades and is related to determining the steady-state conditions of an electrical network when it has at least one slack source and multiple constant power loads [32]. The solution of the power flow problem in three-phase distribution networks with asymmetric loads, as well as the case of single-phase systems, requires the application of numerical methods and a set of nonlinear equality constraints involved in this problem [33]. Once the load flow problem has been solved, it is possible to determine all the electrical characteristics of the electrical network under analysis, i.e., the total grid losses, the maximum voltage regulation, and the maximum loadability in the transmission/distribution lines, among others [34]. This research adopts the successive approximations method for solving the load flow in three-phase asymmetric networks, as reported by [5]. The main characteristics of this power flow method are as follows:
  • It is a derivative-free load flow method and corresponds to a generalization of the classical backward/forward load flow method;
  • The convergence of the successive approximations method is linear due to the absence of derivatives in its formulation;
  • By applying Banach’s fixed-point theorem, it is possible to ensure its convergence to the load flow solution if and only if the distribution grid operates far from the voltage collapse point [35].
The formulation of the successive approximations load flow method for general three-phase asymmetric distribution grids starts with the general definition of complex power (Tellegen’s second theorem) for a generic node k, as defined by (1):
S k = diag V k I k ,
where S k is the complex average power consumption at node k, which has a complex voltage V k , and a net current injection I k . Note that A represents the conjugate value of the complex variable/parameter A, and diag B is an operator that makes the vector B a diagonal matrix. If the conjugate operator is taken on both sides of (1), then the following equivalent equation is obtained:
S k = diag V k I k .
According to [36], the relationship between currents and voltages for a three-phase system is given by the following equations:
I k a b c = Y k a b c V k a b c ,
where V k a b c , I k a b c , and Y k a b c are defined as follows:
V k a b c = V k ( a ) V k ( b ) V k ( c ) ,
I k a b c = I k ( a ) I k ( b ) I k ( c ) ,
Y k a b c = Y k ( a a ) Y k ( a b ) Y k ( a c ) Y k ( b a ) Y k ( b b ) Y k ( b c ) Y k ( c a ) Y k ( c b ) Y k ( c c ) ,
where V k ( a ) , V k ( b ) , and V k ( c ) are the a-, b-, and c-phase voltages at node k; I k ( a ) , I k ( b ) , and I k ( c ) are the net injected currents in phases a, b, and c at node k; and Y k a b c represents the nodal admittance matrix that relates phases a, b, and c at node k.
Suppose that the definitions of the three-phase variables in (4) and (5) are replaced into the general complex power of node k in (2). In that case, the general expression for the three-phase complex power is obtained as presented below:
S k a b c = diag V k a b c I k a b c , .
In addition, if the nodal admittance matrix (6) is substituted into (7), the results of (8) and (9) yield the following:
S k a b c = diag V k a b c Y k a b c , V k a b c , ,
S k a b c , = diag ( V k a b c , ) Y k a b c V k a b c .
For the analysis of our load flow, the nodes are classified into generation nodes (g) and load or demand nodes (d), where the slack node is taken as the generating node and the rest as demand nodes. With this in mind, the following formulation can be made:
I g I d = Y g g Y g d Y d g Y d d V g V d
where I g is a complex variable that contains the net current injection in the generation node, I d is a complex vector that contains all the demanded currents of the network, V g is a complex variable that contains the voltage output of the generation source, and V d is a complex vector that contains all the demanded voltages in the consumption nodes. Now, by substituting the definition in (10) into (9), the following set of general equalities is found:
S g a b c , = diag V g a b c , Y g g a b c V g a b c + Y g d a b c V d a b c ,
S d a b c , = diag V d a b c , Y d g a b c V g a b c + Y d d a b c V d a b c ,
where S g a b c and | S d a b c | are the apparent powers of phases a, b, and c generated by the slack node, as well as the apparent power of phases a, b, and c consumed by the load nodes, respectively; | V g a b c | is the three-phase voltage output of the generation sources, and | V d a b c | is a complex vector with all the unknown voltages of the remaining nodes. It is worth noting that the negative sign in (12) is associated with the direction of the injected power, as it leaves the demand nodes and arrives at the generation nodes.
Now, if what is described in (12) is mathematically reorganized, it is possible to determine the unknown voltage for each load node:
V d a b c = Y d d a b c 1 diag V d a b c , 1 S d a b c , + Y d g a b c V g a b c
Due to the nonlinear structure of Equation (13), it is necessary to implement an iterative process to reach a numerical solution with an acceptable convergence ( ε ).
V d ( a b c , t + 1 ) = Y d d a b c 1 diag V d a b c , , t 1 S d a b c , + Y d g a b c V g a b c , |
where t is defined as the iterative counter, and the evaluation of the recursive formula reaches convergence when the criterion in (15) is met:
max V d a b c , t + 1 V d a b c , t ε
As recommended by [5], the value assigned for the ε -parameter is 1 × 10 10 .

3. Phase-Balancing and Fixed Capacitor Banks in Distribution Grids

This section describes the main aspects of the problems regarding phase-balancing and optimal reactive power injection with fixed-step capacitor banks in asymmetric three-phase distribution networks.

3.1. Phase-Balancing Problem

In three-phase asymmetric distribution networks, one of the most common problems concerns the current imbalances along the distribution feeder and the voltage imbalances in the demand nodes [37]. This is due to the presence of multiple single-, two-, and three-phase loads along the distribution grid, which are not adequately connected to the grid, i.e., without considering an efficient expansion/operation plan [38].
To deal with the problems caused by load imbalances in three-phase distribution networks, one of the most widely known and cheapest alternatives is optimal phase-balancing using mathematical optimization [5]. In three-phase systems, for a general three-phase load, there are six different possible connections H ( i ) between phases a, b, and c. These configurations are listed in Table 1.
Remark 1.
The main characteristic of the type of load connections in Table 1 is that the first three connections are defined in a positive sequence, and the other three in a negative sequence, which means that one must check whether there are rotating machines at the points of load connection. In this sense, an inversion of the input voltage sequence can damage these devices due to torque inversion [26].
The main idea of optimal phase-swapping in three-phase asymmetric networks is to minimize the expected power losses under particular load conditions. These power losses can be calculated via Equation (16):
min p loss = real S g a b c S d a b c ,
where p loss is the expected value of the power losses, calculated as the difference between the power generation and power consumption.
The minimization of the power losses defined in (16) is subject to the expected electrical behavior of the distribution network under analysis, which is determined by the active and reactive power balance at each node and phase of the distribution grid. These are defined generically in Equations (17) and (18). In addition, the voltage regulation limits are defined in inequality constraint (19):
P V k a b c , Θ k a b c , V m a b c , Θ m a b c , Y k a b c , Y m a b c , H = 0 ,
Q V k a b c , Θ k a b c , V m a b c , Θ m a b c , Y k a b c , Y m a b c , H = 0 ,
V k a b c , min V k a b c V k a b c , max .
Note that Equations (17) and (18) represent the general active and reactive power balance constraints at each node k at each one of the phases, which are functions of the voltage variables V k a b c and V m a b c , their angles Θ k a b c and Θ m a b c , and the parametric information of the distribution network (i.e., the admittance matrices Y k a b c and Y m a b c , and the load connections H defined in Table 1).
Remark 2.
The solution of the optimization model (16)–(19) requires the application of efficient optimization methods due to the intrinsic nonlinear, non-convex characteristics of the active and reactive power balance equations, as well as the inclusion of integer variables to represent each possible load connection at each node of the three-phase asymmetric distribution grid.
In this research, the set of load connections at each node (i.e., the values of the variables H i ) is defined by applying the black hole optimizer while using a discrete codification.

3.2. Capacitor Banks

The optimal placement of capacitor banks in electrical distribution networks, in conjunction with solving the optimal phase-balancing problem, is one of the most common and cheapest approaches to reducing the expected cost of the energy losses in distribution grids [10]. An additional advantage of including capacitor banks for reactive power compensation in distribution grids is that these devices also contribute to voltage profile improvements along the distribution feeder, especially in nodes far from the substation bus [15].
The main characteristic of installing capacitor banks in distribution networks is using fixed reactive power compensation values. This research employs the parametric information reported by [10], where capacitors with sizes of 150 kvar are implemented while considering the possibility of installing 0 to 12 per node in parallel (see the sizes of the capacitors in Table 2).
Equation (20) presents the final value of the set of capacitor banks located at node k, defined as V R C ( k ) :
V R C ( k ) = n Q ( c ) ,
where n is the number of capacitor banks connected in parallel, each of them with a size of Q ( c ) = 150 kvar.
Once the sizes of the set of capacitor banks that will be connected to the node k are known, they must be subtracted from the reactive power consumption of these nodes to quantify their effect in terms of power losses by using the three-phase power flow solution presented in Section 2. The new reactive power demand in the nodes where the capacitor banks must be placed is presented in Equation (21):
Q D k ( a ) Q D k ( b ) Q D k ( c ) = Q d k ( a ) V R C ( k ) Q d k ( b ) V R C ( k ) Q d k ( c ) V R C ( k ) ,
where Q d k ( a ) , Q d k ( b ) , and Q d k ( c ) correspond to the reactive power demand values for phases a, b, and c at node k, respectively; and Q D k ( a ) , Q D k ( b ) , and Q D k ( c ) correspond to the value of reactive power demanded by phases a, b, and c in node k after placing and connecting the fixed-step capacitor banks.
Remark 3.
Note that, to determine the best set of capacitor banks that must be installed at each node k, it is necessary to implement efficient optimization methods, given the discrete nature of the solution space and the nonlinear relation with the active and reactive power constraints in the power flow problem.
In this research, the set of capacitor banks (i.e., the values of the variables R C ( k ) ), is defined by applying the black hole optimizer while using a discrete codification.

4. Solution Methodology Based on the BHO Approach

Bio-inspired computing is a field of study based on natural phenomena that seek to solve complex optimization problems [39]. At present, the complexity of optimization problems is increasing almost exponentially, which makes it increasingly challenging to solve said problems using classical algorithms. This has motivated the research and development of bio-inspired optimization algorithms [40].
The BHO algorithm is a bio-inspired metaheuristic method based on the dynamic interaction that exists between the center of a black hole and the cosmic matter transiting around it, whose purpose is to solve large nonlinear optimization problems by using simple mathematical formulas [41,42]. BHO was introduced in 2013 by [43], who illustrated it as an iterative method based on population, where the best candidate solution is selected as the black hole and the rest of the candidates are the cosmic matter revolving around it, i.e., stars.
This research considers two optimization strategies to solve the problems regarding optimal phase-balancing and optimal reactive power compensation in three-phase asymmetric distribution networks: the cascade optimization method and the simultaneous optimization approach.
  • Cascade solution methodology (CSM): The cascade optimization approach consists of a sequential solution of the studied optimization problems. First, the BHO method is applied to define the best set of load connections in all the nodes of the network, i.e., to solve the optimal phase-swapping problem. Second, fixing the solution of the phase-balancing problem, BHO is applied once again to obtain the set of fixed-step capacitor banks that will be connected to the three-phase nodes. The final solution corresponds to the grid power losses after locating and sizing the set of capacitor banks.
  • Simultaneous solution methodology (SSM): The simultaneous solution strategy consists of using a unified codification to represent the problems regarding optimal phase-swapping and the optimal location of fixed-step capacitor banks with a single vector. Each candidate solution is evaluated with the three-phase power flow approach in order to determine the final value of the grid power losses. The final solution corresponds to the black hole’s location with the best load connections and capacitor sizes, i.e., the minimum value for the grid power losses.
Both of the proposed solution methodologies, which are based on applying the BHO method, are explained in Algorithms 1 and 2.
Algorithm 1 Cascade solution methodology
  1:
Stage 1
  2:
Data: BHO parameters and test feeder information.
  3:
for t 1 = 1: t m a x do
  4:
    if  t 1 = = 1 then
  5:
        Generate the initial population of stars for the phase-balancing problem;
  6:
        Evaluate each star with the three-phase power flow algorithm;
  7:
        Determine the expected power losses for each star (objective function value);
  8:
        Establish the position of the black hole (best current solution)
  9:
     else
10:
        Generate the descending population for the new set of load connections;
11:
        Evaluate each new star in the three-phase power flow algorithm;
12:
        Determine the expected power losses for each star (objective function value);
13:
        Update black hole position;
14:
        Replace the stars absorbed by the black hole;
15:
     end if
16:
end for
17:
Stage 2
18:
Data: Fix the solution of the phase-balancing problem and define the set of capacitor banks available
19:
for t 2 = 1: t m a x do
20:
    if  t 2 = = 1 then
21:
        Generate the initial population of stars for the capacitor sizing and location problem;
22:
        Evaluate each star with the three-phase power flow algorithm;
23:
        Determine the expected power losses for each star (objective function value);
24:
        Establish the position of the black hole (best current solution)
25:
    else
26:
        Generate the descending population for the new set of capacitor banks;
27:
        Evaluate each new star in the three-phase power flow algorithm;
28:
        Determine the expected power losses for each star (objective function value);
29:
        Update black hole position;
30:
        Replace the stars absorbed by the black hole;
31:
    end if
32:
end for
Algorithm 2 Simultaneous solution methodology
  1:
Data: BHO parameters and test feeder information.;
  2:
for t : 1 = 1: t m a x do
  3:
    if  t 1 = = 1 then
  4:
     Generate the initial population of stars for the phase-balancing and capacitor sizing and location problems;
  5:
     Evaluate each star with the three-phase power flow algorithm;
  6:
     Determine the expected power losses for each star (objective function value);
  7:
     Establish the position of the black hole (best current solution)
  8:
   else
  9:
     Generate the descending population for the new set of load connections and capacitor sizes and locations;
10:
       Evaluate each new star with the three-phase power flow algorithm;
11:
       Determine the expected power losses for each star (objective function value);
12:
       Update black hole position;
13:
       Replace the stars absorbed by the black hole;
14:
    end if
15:
end for
Remark 4.
The detailed aspects of the implementation of the BHO method (i.e., evolution rules, black hole updating, and stars replacement, among others) can be consulted in [43].

5. Test Feeders

This section presents the main characteristics of the test feeders used for validating the solution methodologies: IEEE 8- and 25-bus grids. Both networks are asymmetric and highly unbalanced. The main characteristics of these networks are described below.

5.1. IEEE 8-Bus Grid

The IEEE 8-bus grid is composed of a slack node located at node 1 and seven demand nodes, all connected by seven lines with a radial structure [5]. This test feeder operates with a medium voltage value of 11 kV and total per-phase active and reactive power values of 1005 kW and 485 kvar, 785 kW and 381 kvar, and 1696 kW and 821 kvar. The electrical configuration of this test feeder is depicted in Figure 1.
The parametric information regarding three-phase loads and impedances for the IEEE 8-bus grid is presented in Table 3 and Table 4.

5.2. IEEE 25-Bus Grid

The IEEE 25-bus grid is a radial medium-voltage distribution network that operates at a voltage of 4.16 kV in terminals of the substation [44]. It has 25 nodes, 22 nodes with unbalanced loads, two step nodes, and one substation bus located at node 1. The electrical configuration of this test feeder is presented in Figure 2.
The parametric information regarding three-phase loads and impedances for the IEEE 25-bus grid is presented in Table 5 and Table 6.

6. Numerical Results

This section shows the results obtained when the BHO algorithm was implemented in the 8- and 25-node IEEE test systems for the cases described above, i.e., the required new configurations and locations of the capacitor banks.
To define the best parameters for the proposed BHO in the CSM and SSM cases, multiple evaluations in both test feeders were conducted, modifying the number of individuals in the population and the number of iterations. In addition, two different computers were used to corroborate the effectiveness of the optimization method. After these evaluations, the set of parameters selected for the BHO approach was 1000 iterations per evaluation and a dimension of the initial population equivalent to 30 stars.

6.1. Benchmark Case Evaluation

To corroborate that the three-phase power flow solution presented in Section 2 is efficient in solving the power flow problem, a comparative analysis was performed with benchmark cases in both test feeders, i.e., without optimal phase-balancing and capacitor connection, while using the DIgSILENT software with the Newton–Raphson method. Note that the selection of the DIgSILENT software for validating the power results is based on its wide acceptance in the current literature and the industry for power system studies, as it is considered to be a quasi-experimental validation [45,46]. Table 7 presents the power flow solution obtained by using the MATLAB programming environment and the DIgSILENT software.
The numerical results in Table 7 allow for confirming that the successive approximations power flow method presented in Section 2 is efficient in solving power flow problems in three-phase distribution grids. In comparison with the Newton–Raphson method available in the DIgSILENT software, less than 0.011 % was obtained in both test feeders. This difference can be attributed to the convergence error that can be set for the Newton–Raphson in the DIgSILENT software, which is about 1 × 10 06 , whereas the MATLAB power flow was set with a convergence error of about 1 × 10 10 .

6.2. Results for the IEEE 8-Bus Grid

In this test feeder, the first evaluation concerned the effect of the CSM on the final grid power losses of the network. The first stage, as presented in Figure 1, is the solution to the optimal phase-swapping problem. Table 8 presents the solution to the phase-balancing problem using the BHO method and its evaluations in the DIgSILENT software.
Note that the comparison between MATLAB and DIgSILENT confirms the effectiveness of the proposed optimization method to redefine the load connections, reaching final power losses values of about 10.5868 kW, i.e., a reduction of 24.34 % with respect to the benchmark case. It is worth mentioning that the solution reached by the BHO approach is equivalent to the solution reported in [5] with the application of the vortex search algorithm.
To verify the effect of integrating capacitor banks in three-phase asymmetric distribution grids, the BHO approach was implemented in order to locate these devices on the benchmark test system, i.e., without employing the phase-balancing approach. These results are listed in Table 9, including the comparison between the MATLAB power flow and the DIgSILENT one.
The numerical results in Table 9 show that, with respect to the benchmark case, the optimal placement of fixed-step capacitor banks allowed a power losses reduction of 13 % , i.e., 1.8193 kW, which is lower in comparison with the phase-swapping solution reported in Table 8. It also has a positive impact on the final grid power losses, and the combination of both approaches for the CSM and the SSM will perhaps yield better results in terms of power loss minimization since both approaches independently allow for reducing this indicator with respect to the benchmark case.
Table 10 presents the numerical comparison between the proposed CSM and SSM presented in Figure 1 and Figure 2.
The numerical results presented in Table 10 allow for noting that:
  • Both solution methodologies found the same objective function value regarding the final grid power losses, i.e., 8.7780 kW. This implies a reduction of about 37.27 % with respect to the benchmark case;
  • For the IEEE 8-bus grid, the CSM and the SSM are equivalent in terms of the final objective function value. However, due to the nonlinearities and non-convexities of the solution space for the studied problems, the final load connections in both solution methods differ from each other.
It is worth mentioning that the SSM and CSM found the same objective function values due to the small size of the solution space for the studied problem in the IEEE 8-bus grid. As demonstrated by the authors of [5] for the optimal phase balancing problem, the solution of the IEEE 8-bus grid is the global optimum, and, given the results obtained by this system, perhaps the size and location of the fixed-step capacitor banks presented in this paper is also optimal.
To demonstrate that the solutions reached with both methodologies in the MATLAB programming environment are equivalent to those of the three-phase power flow approach in the DIgSILENT software, Table 11 lists the voltage magnitudes per phase.
In addition, to verify the effect of optimal phase-balancing while considering capacitor banks in three-phase asymmetric networks, the behavior of the voltage profiles per phase is presented in Figure 3. This figure shows the voltage profile for the benchmark case, as well as for the CSM and the SSM.
The behavior of the voltage profile for the IEEE 8-bus grid (Figure 3) allows for making the following remarks. (i) The voltage profile in phase c for the benchmark case (Figure 3a) is the worst profile when compared to phases a and b. However, this is an expected behavior, as c is the most charged phase in the IEEE-8 bus grid. (ii) The plots in Figure 3b,c are identical because the CSM and the SSM found the same objective function value (see Table 10). However, the main result is the fact the all voltages are uniform, i.e., the load redistribution in all the nodes and the reactive power compensation allowed for reaching a more balanced operation at bus 8 for phase b with a minimum voltage of 0.9560 pu.

6.3. Results for the IEEE 25-Bus Grid

As previously shown with the IEEE 8-bus grid, the first stage of the CSM reorganizes the system loads in order to reduce losses. Table 12 presents the solution to the phase balance problem using the BHO method and its verification in DIgSILENT.
The results in Table 12 allow for stating that the evaluation of the optimal solution provided by the BHO algorithm in the IEEE 25-bus grid shows similarities between the MATLAB and the DIgSILENT software, where the minimum value was reached when the three-phase power flow approach was implemented in MATLAB, i.e., in comparison with the Newton–Raphson technique available in the DIgSILENT software. This difference is explained by the convergence errors used in each algorithm, given that the successive approximations method uses the difference between two consecutive voltages, with a maximum error of 1 × 10 10 , whereas the Newton–Raphson approach uses the mismatch between calculated and specified powers, with a precision of 1 × 10 06 .
On the other hand, in order to observe the effect of using fixed-step capacitor banks for shunt reactive power compensation in the IEEE 25-bus grid, Table 13 presents a comparative analysis between MATLAB and DIgSILENT once the BHO was implemented.
The results in Table 13 show that the expected power losses in the IEEE 25-bus grid after installing capacitor banks constitute a reduction of about 28.29% with respect to the benchmark case (i.e., 75.4207 kW), which implies that, in this test feeder, shunt reactive power compensation allowed for higher power loss reductions when compared to the phase-balancing case, unlike those reported by the IEEE 8-bus grid. These differences between both test feeders are explained by the nonlinearities and non-convexities of the optimization models that represent both studied problems.
To show the effect of using the CSM in the IEEE 25-bus grid, Table 14 presents the solutions obtained by this approach, which is based on the algorithm depicted in Figure 1, as well as a comparison with the SSM presented in Figure 2.
The results in Table 14 allow for making the following remarks. (i) The CSM allows for finding a final objective function value of about 51.2485 kW, i.e., a reduction of about 32.05% with respect to the benchmark case, while the SSM has a final objective function value of about 50.1396 kW, i.e., a 33.52% reduction in the total grid power losses. (ii) The CSM and the SSM differ in the final set of nodal connections for loads and in the size and location of the fixed-step capacitor banks. However, this is an expected behavior, as the CSM approach focuses on minimizing each problem separately, i.e., the first stage solves the phase-swapping problem, and the second stage solves the optimal reactive power compensation problem. Nevertheless, the solution to the second problem is completely conditioned by that obtained in the first stage. In contrast, the SSM method addresses both problems simultaneously, which implies no conditioning by the solution of an initial stage, thus allowing for a better exploration and exploitation of the solution space.
To confirm that the solutions reached with both methodologies in the MATLAB programming environment are equivalent to the solution of the three-phase power flow approach obtained via the DIgSILENT software, Table 15 lists voltage magnitudes per phase.
Furthermore, in order to verify the effect of optimal phase-balancing while considering capacitor banks in three-phase asymmetric networks, the per-phase behavior of the voltage profiles is presented in Figure 4. This figure shows the voltage profile of the benchmark case, as well as that of the SSM, since, as observed in Table 15, it yields better objective function values in comparison with the CSM.
The behavior of the voltage profile for the IEEE 25-bus grid in Figure 4 allows for observing that: (i) the minimum voltage of the benchmark case is defined by phase at node 13, with a magnitude of 0.9584 pu, i.e., a voltage regulation of about 4.16 % (see Figure 4a); and (ii) after implementing the solution provided by the SSM, all the voltage profiles in phases b and c become similar, which is caused by the simultaneous effect of phase-swapping and shunt reactive power compensation (see Figure 4b,c). However, the worst voltage profile now corresponds to phase b at node 13, with a magnitude of about 0.9675 pu, which means a regulation of about 3.25%.

6.4. Complementary Analysis

To summarize, the main results regarding the reductions in the total grid power losses are presented in Table 16 for the cases in which the phase-balancing problem and the optimal location of capacitor banks are implemented separately, as well as when the cascade and solution methodologies are employed.
The results in Table 16 reveal the following: (i) When the phase-balancing problem and the optimal selection and location of capacitor banks are solved independently, important variations occur in the expected energy loss reductions, and it is not possible to generalize a tendency. As explained in previous sections, this can be attributed to the nonlinearities and non-convexities of both problems, in addition to the particularities of each test feeder, i.e., the number of nodes, branch impedances, and longitudes, the total power load, and the per-phase unbalanced load. (ii) It is evident that the combination of the phase-balancing problem with the selection and location of shunt capacitors in both methodologies (i.e., CSM and SSM) allows for better improvements in the final objective function when compared to the separate use of each approach. (iii) There is a need to explore alternative solution methods with convex approximations or new combinatorial optimizers in order to validate the effectiveness of the BHO, which constitutes a research opportunity that remains open for future studies.

7. Conclusions

This paper addressed two classical engineering problems concerning distribution networks. The first problem involved optimal phase-swapping in three-phase asymmetric distribution networks, and the second was optimal shunt reactive power compensation using fixed-step capacitor banks. A master–slave solution methodology was implemented to solve both problems, based on applying the BHO approach combined with the successive approximations power flow method. The main contribution of this research is its proposal and application of two solution methodologies to solve both optimization problems. The first approach was the CSM, where the first stage solves the phase-balancing problem, and the second one solves the shunt reactive power compensation problem by fixing the set of load configurations provided by the solution of the first stage. The second approach was the SSM, which solves both studied problems simultaneously.
Numerical results in the IEEE 8- and 25-bus grids demonstrated the following:
  • For the IEEE 8-bus grid, both solution methodologies found the same value regarding the total grid power losses, with a reduction of about 32.27% concerning the benchmark case. The size of the solution space explains the coincidence in both solutions in this test feeder, which can be regarded as undersized for real distribution grids. This implies that an efficient optimization method can find the optimal global solution.
  • In the IEEE 25-bus grid, the solution of the SSM was better than that of the CSM by about 1.1089 kW. This result is expected since the solution space increased in size. The SSM is more capable of exploring and exploiting the solution space, as the solution of the phase-balancing problem does not condition it.
  • More research is required to include, within the proposed optimization model, the possibility of having shunt reactive power compensation per phase (unbalanced reactive power compensation) combined with the phase-balancing solution employing CSM and SSM. In addition, a complete comparative analysis with different combinatorial optimizers is needed in order to verify the effectiveness of the proposed solution methodology addressed in this research.
In future works, it will be possible to conduct the following studies: (i) applying new combinatorial optimization methods to deal with the studied problems while using the SSM; (ii) considering daily load profiles and different end-users (residential, industrial, and commercial) in the optimization model; (iii) allocating static distribution compensators in three-phase grids for hourly reactive power control according to the grid requirements; and (iv) evaluating additional objective functions while considering daily demand curves, as is the case of the total required investments in compensation devices and grid interventions and the expected yearly reduction costs regarding energy losses.

Author Contributions

Conceptualization, O.D.M.; Methodology, D.F.A.M.-G., I.D.R.-R. and O.D.M.; Software, D.F.A.M.-G. and I.D.R.-R.; Validation, D.F.A.M.-G. and I.D.R.-R.; Formal analysis, O.D.M.; Investigation, I.D.R.-R.; Resources, O.D.M.; Writing—original draft, D.F.A.M.-G., I.D.R.-R. and O.D.M.; Writing—review & editing, O.D.M.; Visualization, O.D.M.; Supervision, O.D.M.; Project administration, O.D.M.; Funding acquisition, O.D.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chowdhury, P.K.R.; Weaver, J.E.; Weber, E.M.; Lunga, D.; LeDoux, S.T.M.; Rose, A.N.; Bhaduri, B.L. Electricity consumption patterns within cities: Application of a data-driven settlement characterization method. Int. J. Digit. Earth 2019, 13, 119–135. [Google Scholar] [CrossRef] [Green Version]
  2. Mutumba, G.S.; Odongo, T.; Okurut, N.F.; Bagire, V. A survey of literature on energy consumption and economic growth. Energy Rep. 2021, 7, 9150–9239. [Google Scholar] [CrossRef]
  3. Nogueira, T.; Sousa, E.; Alves, G.R. Electric vehicles growth until 2030: Impact on the distribution network power. Energy Rep. 2022, 8, 145–152. [Google Scholar] [CrossRef]
  4. Arefi, A.; Shahnia, F.; Ledwich, G. (Eds.) Electric Distribution Network Management and Control; Springer: Singapore, 2018. [Google Scholar] [CrossRef]
  5. Cortés-Caicedo, B.; Avellaneda-Gómez, L.S.; Montoya, O.D.; Alvarado-Barrios, L.; Chamorro, H.R. Application of the Vortex Search Algorithm to the Phase-Balancing Problem in Distribution Systems. Energies 2021, 14, 1282. [Google Scholar] [CrossRef]
  6. Arefi, A.; Olamaei, J.; Yavartalab, A.; Keshtkar, H. Loss reduction experiences in electric power distribution companies of Iran. Energy Procedia 2012, 14, 1392–1397. [Google Scholar] [CrossRef]
  7. Abagiu, S.; Lepadat, I.; Helerea, E. Solutions for energy losses reduction in power networks with renewable energy sources. In Proceedings of the 2016 International Conference on Applied and Theoretical Electricity (ICATE), Craiova, Romania, 6–8 October 2016. [Google Scholar] [CrossRef]
  8. Hesaroor, K.; Das, D. Annual energy loss reduction of distribution network through reconfiguration and renewable energy sources. Int. Trans. Electr. Energy Syst. 2019, 29, e12099. [Google Scholar] [CrossRef]
  9. Chernykh, A.G.; Barykina, Y.N.; Morozevich, O.A. Development of methods for minimizing energy losses in electrical networks. IOP Conf. Ser. Earth Environ. Sci. 2022, 1070, 012006. [Google Scholar] [CrossRef]
  10. Prakash, D.; Lakshminarayana, C. Optimal siting of capacitors in radial distribution network using Whale Optimization Algorithm. Alex. Eng. J. 2017, 56, 499–509. [Google Scholar] [CrossRef]
  11. Sulaiman, M.H.; Mustaffa, Z. Optimal placement and sizing of FACTS devices for optimal power flow using metaheuristic optimizers. Results Control Optim. 2022, 8, 100145. [Google Scholar] [CrossRef]
  12. Valencia, A.; Hincapie, R.A.; Gallego, R.A. Optimal location, selection, and operation of battery energy storage systems and renewable distributed generation in medium–low voltage distribution networks. J. Energy Storage 2021, 34, 102158. [Google Scholar] [CrossRef]
  13. Ghiasi, M. Detailed study, multi-objective optimization, and design of an AC-DC smart microgrid with hybrid renewable energy resources. Energy 2019, 169, 496–507. [Google Scholar] [CrossRef]
  14. Salau, A.O.; Gebru, Y.W.; Bitew, D. Optimal network reconfiguration for power loss minimization and voltage profile enhancement in distribution systems. Heliyon 2020, 6, e04233. [Google Scholar] [CrossRef]
  15. Askarzadeh, A. Capacitor placement in distribution systems for power loss reduction and voltage improvement: A new methodology. IET Gener. Transm. Distrib. 2016, 10, 3631–3638. [Google Scholar] [CrossRef]
  16. Ghiasi, M.; Olamaei, J. Optimal capacitor placement to minimizing cost and power loss in Tehran metro power distribution system using ETAP (A case study). Complexity 2016, 21, 483–493. [Google Scholar] [CrossRef]
  17. Tamilselvan, V.; Jayabarathi, T.; Raghunathan, T.; Yang, X.S. Optimal capacitor placement in radial distribution systems using flower pollination algorithm. Alex. Eng. J. 2018, 57, 2775–2786. [Google Scholar] [CrossRef]
  18. Devabalaji, K.; Yuvaraj, T.; Ravi, K. An efficient method for solving the optimal sitting and sizing problem of capacitor banks based on cuckoo search algorithm. Ain Shams Eng. J. 2018, 9, 589–597. [Google Scholar] [CrossRef] [Green Version]
  19. Taher, S.A.; Bagherpour, R. A new approach for optimal capacitor placement and sizing in unbalanced distorted distribution systems using hybrid honey bee colony algorithm. Int. J. Electr. Power Energy Syst. 2013, 49, 430–448. [Google Scholar] [CrossRef]
  20. Ogita, Y.; Mori, H. Parallel Dual Tabu Search for Capacitor Placement in Smart Grids. Procedia Comput. Sci. 2012, 12, 307–313. [Google Scholar] [CrossRef]
  21. Faiz, M.; Chang, A.Q.; Memon, N.; Pathan, A.Z.U. Optimal Capacitor Placement Using Tabu Search Algorithm to Improve the Operational Efficiency in GEPCO Network. Int. J. Electr. Electron. Eng. 2021, 8, 1–5. [Google Scholar] [CrossRef]
  22. Gil-González, W.; Montoya, O.D.; Rajagopalan, A.; Grisales-Noreña, L.F.; Hernández, J.C. Optimal Selection and Location of Fixed-Step Capacitor Banks in Distribution Networks Using a Discrete Version of the Vortex Search Algorithm. Energies 2020, 13, 4914. [Google Scholar] [CrossRef]
  23. Khan, N.A.; Ghosh, S.; Ghoshal, S.P. Binary Gravitational Search based Algorithm for Optimum Siting and Sizing of DG and Shunt Capacitors in Radial Distribution Systems. Energy Power Eng. 2013, 05, 1005–1010. [Google Scholar] [CrossRef]
  24. Zadeh, A.Z.; Andami, H.; Talavat, V.; Ebrahimi, J. Optimal Capacitor Placement in the Unbalanced Distribution Networks Contaminated by Harmonic through Imperialist Competitive Algorithm. Res. J. Appl. Sci. Eng. Technol. 2014, 7, 1230–1235. [Google Scholar] [CrossRef]
  25. Murty, V.; Kumar, A. Capacitor Allocation in Unbalanced Distribution System under Unbalances and Loading Conditions. Energy Procedia 2014, 54, 47–74. [Google Scholar] [CrossRef] [Green Version]
  26. Echeverri, M.G.; Rendón, R.A.G.; Lezama, J.M.L. Optimal Phase Balancing Planning for Loss Reduction in Distribution Systems using a Specialized Genetic Algorithm. Ing. Cienc. 2012, 8, 121–140. [Google Scholar] [CrossRef] [Green Version]
  27. Khodr, H.; Zerpa, I.; de Jesu’s, P.D.O.; Matos, M. Optimal Phase Balancing in Distribution System Using Mixed-Integer Linear Programming. In Proceedings of the 2006 IEEE/PES Transmission & Distribution Conference and Exposition: Latin America, Caracas, Venezuela, 15–18 August 2006. [Google Scholar] [CrossRef]
  28. El Hassan, M.; Najjar, M.; Tohme, R. A Practical Way to Balance Single Phase Loads in a Three Phase System at Distribution and Unit Level. Renew. Energy Power Qual. J. 2022, 20, 173–177. [Google Scholar] [CrossRef]
  29. Mansani, S.; Udaykumay, R.Y. An optimal phase balancing technique for unbalanced three-phase secondary distribution systems. In Proceedings of the 2016 IEEE 7th Power India International Conference (PIICON), Bikaner, India, 25–27 November 2016. [Google Scholar] [CrossRef]
  30. Alhmoud, L.; Nawafleh, Q.; Merrji, W. Three-Phase Feeder Load Balancing Based Optimized Neural Network Using Smart Meters. Symmetry 2021, 13, 2195. [Google Scholar] [CrossRef]
  31. Cruz-Reyes, J.L.; Salcedo-Marcelo, S.S.; Montoya, O.D. Application of the Hurricane-Based Optimization Algorithm to the Phase-Balancing Problem in Three-Phase Asymmetric Networks. Computers 2022, 11, 43. [Google Scholar] [CrossRef]
  32. Pandey, A.; Jereminov, M.; Wagner, M.R.; Bromberg, D.M.; Hug, G.; Pileggi, L. Robust Power Flow and Three-Phase Power Flow Analyses. IEEE Trans. Power Syst. 2019, 34, 616–626. [Google Scholar] [CrossRef] [Green Version]
  33. Garces, A. A Linear Three-Phase Load Flow for Power Distribution Systems. IEEE Trans. Power Syst. 2016, 31, 827–828. [Google Scholar] [CrossRef]
  34. Marini, A.; Mortazavi, S.; Piegari, L.; Ghazizadeh, M.S. An efficient graph-based power flow algorithm for electrical distribution systems with a comprehensive modeling of distributed generations. Electr. Power Syst. Res. 2019, 170, 229–243. [Google Scholar] [CrossRef]
  35. Shen, T.; Li, Y.; Xiang, J. A Graph-Based Power Flow Method for Balanced Distribution Systems. Energies 2018, 11, 511. [Google Scholar] [CrossRef] [Green Version]
  36. Gnanambal, K.; Marimuthu, N.; Babulal, C. Three-phase power flow analysis in sequence component frame using Hybrid Particle Swarm Optimization. Appl. Soft Comput. 2011, 11, 1727–1734. [Google Scholar] [CrossRef]
  37. Singh, M.K.; Gupta, A.R.; Kumar, A. Analysis of Unbalanced Radial Distribution System with SVR and Impact of Phase Shifter Angle considering Time Varying Loads. Smart Sci. 2021, 1–16. [Google Scholar] [CrossRef]
  38. Swapna, M.; Udaykumar, R.Y. An algorithm for optimal phase balancing of secondary distribution systems at each node. In Proceedings of the 2016 IEEE PES 13th International Conference on Transmission & Distribution Construction, Operation & Live-Line Maintenance (ESMO), Columbus, OH, USA, 12–15 September 2016. [Google Scholar] [CrossRef]
  39. Kumar, S.; Datta, D.; Singh, S.K. Black Hole Algorithm and Its Applications. In Studies in Computational Intelligence; Springer International Publishing: Cham, Switzerland, 2014; pp. 147–170. [Google Scholar] [CrossRef]
  40. Abualigah, L.; Elaziz, M.A.; Sumari, P.; Khasawneh, A.M.; Alshinwan, M.; Mirjalili, S.; Shehab, M.; Abuaddous, H.Y.; Gandomi, A.H. Black hole algorithm: A comprehensive survey. Appl. Intell. 2022, 52, 11892–11915. [Google Scholar] [CrossRef]
  41. Velasquez, O.S.; Montoya, O.D.; Garrido, V.M.; Grisales, L.F. Optimal Power Flow in Direct-Current Power Grids via Black Hole Optimization. Adv. Electr. Electron. Eng. 2019, 17, 24–32. [Google Scholar] [CrossRef]
  42. Arenas-Acuña, C.A.; Rodriguez-Contreras, J.A.; Montoya, O.D.; Rivas-Trujillo, E. Black-Hole Optimization Applied to the Parametric Estimation in Distribution Transformers Considering Voltage and Current Measures. Computers 2021, 10, 124. [Google Scholar] [CrossRef]
  43. Hatamlou, A. Black hole: A new heuristic optimization approach for data clustering. Inf. Sci. 2013, 222, 175–184. [Google Scholar] [CrossRef]
  44. Ramana, T.; Ganesh, V.; Sivanagaraju, S. Distributed Generator Placement In addition, Sizing in Unbalanced Radial Distribution System. Cogener. Distrib. Gener. J. 2010, 25, 52–71. [Google Scholar] [CrossRef]
  45. Ganesh, S.; Perilla, A.; Torres, J.R.; Palensky, P.; van der Meijden, M. Validation of EMT Digital Twin Models for Dynamic Voltage Performance Assessment of 66 kV Offshore Transmission Network. Appl. Sci. 2020, 11, 244. [Google Scholar] [CrossRef]
  46. Bifaretti, S.; Bonaiuto, V.; Pipolo, S.; Terlizzi, C.; Zanchetta, P.; Gallinelli, F.; Alessandroni, S. Power Flow Management by Active Nodes: A Case Study in Real Operating Conditions. Energies 2021, 14, 4519. [Google Scholar] [CrossRef]
Figure 1. Single-phase diagram of the IEEE 8-bus system.
Figure 1. Single-phase diagram of the IEEE 8-bus system.
Sustainability 15 00366 g001
Figure 2. Single-phase diagram of the IEEE 25-bus system.
Figure 2. Single-phase diagram of the IEEE 25-bus system.
Sustainability 15 00366 g002
Figure 3. Behavior of the voltage profile in the IEEE 8-bus grid: (a) benchmark case; (b) CSM; and (c) SSM.
Figure 3. Behavior of the voltage profile in the IEEE 8-bus grid: (a) benchmark case; (b) CSM; and (c) SSM.
Sustainability 15 00366 g003
Figure 4. Behavior of the voltage profile in the IEEE 25-bus grid: (a) benchmark case; (b) MATLAB; and (c) DIgSILENT.
Figure 4. Behavior of the voltage profile in the IEEE 25-bus grid: (a) benchmark case; (b) MATLAB; and (c) DIgSILENT.
Sustainability 15 00366 g004
Table 1. Possible load connections in a particular node.
Table 1. Possible load connections in a particular node.
Index H i ConnectionSequence
1ABCNo change
2BCA
3CAB
4ACBChange
5CBA
6BAC
Table 2. Set of possible sizes for the capacitor banks.
Table 2. Set of possible sizes for the capacitor banks.
Capacitor Banks
Bank123456789101112
Size150300450600750900105012001350150016501800
Table 3. Load information and conductor types for the IEEE 8-bus system.
Table 3. Load information and conductor types for the IEEE 8-bus system.
LineNode iNode jConductorLenght (Ft) P jA (kW) Q jA (kvar) P jB (kW) Q jB (kvar) P jC (kW) Q jC (kvar)
11215280519250259126515250
2232528000259126486235
325352800000226109
427352804862350000
534452800000324157
638552800026712900
75665280000014570
Table 4. Impedance matrix per type of conductor in the IEEE 8-bus grid.
Table 4. Impedance matrix per type of conductor in the IEEE 8-bus grid.
Conductor ( Ω /mi)
10.093654 + 0.0429300i0.031218 + 0.0134310i0.031218 + 0.0134310i
0.031218 + 0.0134310i0.093654 + 0.0402930i0.031218 + 0.0134310i
0.031218 + 0.0134310i0.031218 + 0.0134310i0.093654 + 0.0402930i
20.156090 + 0.0671550i0.052030 + 0.0223850i0.052030 + 0.0223850i
0.052030 + 0.0223850i0.156090 + 0.0671550i0.052030 + 0.0223850i
0.052030 + 0.0223850i0.052030 + 0.0223850i0.156090 + 0.0671550i
30.046827 + 0.0201465i0.015609 + 0.0067155i0.015609 + 0.0067155i
0.015609 + 0.0067155i0.046827 + 0.0201465i0.015609 + 0.0067155i
0.015609 + 0.0067155i0.015609 + 0.0067155i0.046827 + 0.0201465i
50.062436 + 0.0268620i0.020812 + 0.0089540i0.020812 + 0.0089540i
0.020812 + 0.0089540i0.062436 + 0.0268620i0.020812 + 0.0089540i
0.020812 + 0.0089540i0.020812 + 0.0089540i0.062436 + 0.0268620i
60.078045 + 0.0335775i0.026015 + 0.0111925i0.026015 + 0.0111925i
0.026015 + 0.0111925i0.078045 + 0.0335775i0.026015 + 0.0111925i
0.026015 + 0.0111925i0.026015 + 0.0111925i0.078045 + 0.0335775i
Table 5. Load information and conductor types for the IEEE 25-bus system.
Table 5. Load information and conductor types for the IEEE 25-bus system.
LineNode iNode jConductorLength (ft) P jA (kW) Q jA (kvar) P jB (kW) Q jB (kvar) P jC (kW) Q jC (kvar)
11211000000000
22315003621.628.819.24226.4
326250043.228.833.6243030
434150057.643.24.83.44830
5318250057.643.238.428.84836
645250043.228.828.819.23624
742324008.664.84.83.86042
8672500000000
9682100043.228.828.819.23.62.4
107925007250.438.428.84830
11714250057.63638.428.86042
12716250057.64.33.828.84836
1391025003621.628.819.23226.4
141011230050.431.72414.43624
151112320057.6364833.64836
161113320064.821.633.621.13624
18141523007.24.34.82.963.6
171417330057.643.233.6245438.4
191820250050.43638.428.85438.4
20182134005.84.33.42.45.43.8
21201934008.66.54.83.464.8
22212234007250.457.643.26048
232324240050.43643.230.74.83.6
24242534008.66.54.82.964.2
Table 6. Impedance matrix per type of conductor in the IEEE 25-bus grid.
Table 6. Impedance matrix per type of conductor in the IEEE 25-bus grid.
Conductor ( Ω /mi)
10.3686 + 0.6852i0.0169 + 0.1515i0.0155 + 0.1098i
0.0169 + 0.1515i0.3757 + 0.6715i0.0188 + 0.2072i
0.0155 + 0.1098i0.0188 + 0.2072i0.3723 + 0.6782i
20.9775 + 0.8717i0.0167 + 0.1697i0.0152 + 0.1264i
0.0167 + 0.1697i0.9844 + 0.8654i0.0186 + 0.2275i
0.0152 + 0.1264i0.0186 + 0.2275i0.981 + 0.8648i
31.9280 + 1.4194i0.0161 + 0.1183i0.0161 + 0.1183i
0.0161 + 0.1183i1.9308 + 1.4215i0.0161 + 0.1183i
0.0161 + 0.1183i0.0161 + 0.1183i1.9337 + 1.4236i
Table 7. Comparison of power losses obtained using MATLAB and DIgSILENT in the benchmark cases for the IEEE 8- and IEEE 25-bus grids.
Table 7. Comparison of power losses obtained using MATLAB and DIgSILENT in the benchmark cases for the IEEE 8- and IEEE 25-bus grids.
IEEE 8-Bus Grid
SoftwarePhase a (kW)Phase b (kW)Phase c (kW)Total (kW)
MATLAB1.715792.330489.9462413.99252
DIgSILENT1.715792.329489.9486913.99396
IEEE 25-bus grid
MATLAB36.8800814.7859823.7545375.42059
DIgSILENT36.8843614.7834823.7553275.42315
Table 8. Power losses after solving the optimal phase-balancing problem in the IEEE 8-bus grid.
Table 8. Power losses after solving the optimal phase-balancing problem in the IEEE 8-bus grid.
SoftwarePhase a (kW)Phase b (kW)Phase c (kW)Total (kW)
MATLAB3.846352.741213.9993010.58686
DIgSILENT3.846462.741033.9994210.58691
Table 9. Power losses after solving the problem regarding optimal reactive power compensation in the IEEE 8-bus grid.
Table 9. Power losses after solving the problem regarding optimal reactive power compensation in the IEEE 8-bus grid.
SoftwarePhase a (kW)Phase b (kW)Phase c (kW)Total (kW)
MATLAB1.945361.545998.6818712.17322
DIgSILENT1.941831.547928.6897312.17948
Table 10. Comparison of solution methodologies for the IEEE 8-bus grid.
Table 10. Comparison of solution methodologies for the IEEE 8-bus grid.
CSMSSM
MATLABDIgSILENTMATLABDIgSILENT
Load connections3, 5, 1, 3, 2, 2, 53, 5, 1, 3, 2, 2, 53, 5, 1, 3, 4, 2, 53, 5, 1, 3, 4, 2, 5
Capacitor banks location2, 32, 32, 32, 3
Capacitor banks type3, 13, 13, 13, 1
Phase a (kW)2.263403.220203.218423.22020
Phase b (kW)3.218423.295533.296213.29553
Phase c (kW)3.296212.262712.263402.26271
Total (kW)8.778048.778438.778048.77843
Table 11. Voltage magnitude per phase in the IEEE 8-bus grid after implementing the CSM and/or the SSM.
Table 11. Voltage magnitude per phase in the IEEE 8-bus grid after implementing the CSM and/or the SSM.
MATLABDIgSILENT
Node V a V b V c V a V b V c
11.000001.000001.000000.999990.999991.00002
20.998010.998120.998530.998000.998120.99855
30.996580.996450.997960.996560.996440.99798
40.996620.996500.997660.996610.996490.99768
50.997730.997970.998620.997720.997960.99863
60.997880.997630.998740.997780.997620.99869
70.998120.998250.997850.998100.998240.99786
80.996670.995940.998040.996650.995940.99806
Table 12. Power losses after solving the optimal phase-balancing problem in the IEEE 25-bus grid.
Table 12. Power losses after solving the optimal phase-balancing problem in the IEEE 25-bus grid.
SoftwarePhase a (kW)Phase b (kW)Phase c (kW)Total (kW)
MATLAB25.2232726.4085620.710.1972.34203
DIgSILENT26.7360725.6826419.9581772.37688
Table 13. Power losses after solving the problem regarding optimal reactive power compensation in the IEEE 25-bus grid.
Table 13. Power losses after solving the problem regarding optimal reactive power compensation in the IEEE 25-bus grid.
SoftwarePhase a (kW)Phase b (kW)Phase c (kW)Total (kW)
MATLAB27.923269.4893316.6715154.08410
DIgSILENT27.829059.4135016.7869954.02954
Table 14. Comparison of solution methodologies for the IEEE 25-bus grid.
Table 14. Comparison of solution methodologies for the IEEE 25-bus grid.
CSMSSM
MATLABDIgSILENTMATLABDIgSILENT
Load connections3, 2, 4, 5, 5, 3, 3, 3, 1, 5, 5, 3, 4, 6, 5, 5, 2, 6, 5, 4, 5, 2, 4, 43, 3, 2, 4, 3, 3, 4, 4, 2, 2, 3, 3, 4, 1, 1, 3, 2, 1, 2, 4, 3, 4, 4, 2
Capacitor banks location7, 3, 47, 3, 43, 10, 73, 10, 7
Capacitor banks type2, 1, 12, 1, 12, 1, 12, 1, 1
Phase a (kW)17.4918317.4753618.6221018.62736
Phase b (kW)18.4030718.3719419.3387819.25855
Phase c (kW)15.3535515.3033712.1787312.11915
Total (kW)51.2484651.1506750.1396150.00506
Table 15. Voltage magnitude per phase in the IEEE 25-bus grid after implementing the CSM and/or the SSM.
Table 15. Voltage magnitude per phase in the IEEE 25-bus grid after implementing the CSM and/or the SSM.
CSMSSM
MATLABDIgSILENTMATLABDIgSILENT
Node V a V b V c V a V b V c V a V b V c V a V b V c
11.000001.000001.000001.000001.000000.999991.000001.000001.000000.999991.000001.00001
20.991130.990410.993320.990770.990880.993060.990110.990100.994800.989680.989700.99456
30.989230.988670.991840.988810.988270.991530.988190.987820.993830.987690.987380.99356
40.988470.988310.991530.988010.987880.991200.986100.985890.992420.985600.985450.99214
50.987610.987610.990520.987150.987180.990190.985010.985050.991800.984510.984610.99153
60.984570.983100.987090.984040.982610.986700.982460.983290.989280.981830.982740.98892
70.978960.977840.982410.978270.977200.981910.976800.977570.985260.975970.976840.98478
80.984420.980770.985980.983890.980280.985600.980340.983300.987750.979710.982750.98739
90.973790.972170.977690.973100.971530.977180.973240.973620.982500.972300.972790.98194
100.969850.968340.973720.969150.967690.973210.971570.970720.980620.970520.969790.98000
110.968000.966450.971830.967300.965800.971320.969900.968500.978930.968850.967570.97830
120.967070.965520.970770.966360.964870.970260.968950.967430.978030.967900.966500.97740
130.967300.965410.971240.966600.964760.970720.969200.967460.978340.968150.966530.97771
140.976060.975390.980010.975360.974740.979500.973740.974530.983640.972810.973740.98307
150.975990.975280.979930.975290.974640.979420.973640.974470.983560.972810.973740.98307
160.977800.977360.981430.977110.976710.980920.975950.977170.983930.975120.976440.98344
170.974530.974440.978290.973830.973800.977780.972160.972850.982710.971330.972120.98222
180.985100.984630.987390.984670.984230.987080.984260.983140.989860.983760.982690.98959
190.983440.983120.985830.983010.982720.985510.982720.981460.988350.982220.981020.98808
200.983620.983460.986060.983190.983060.985750.983060.981640.988590.982550.981190.98832
210.982510.982190.984530.982080.981790.984100.981640.980190.987540.981130.979740.98726
220.980150.979950.981790.979720.979550.981470.979240.977450.985340.988740.977000.98507
230.987190.987010.989890.986730.986580.989550.984110.984550.991490.983610.984100.99122
240.986020.986880.988880.985550.986450.988550.983030.984400.990390.982530.983960.99012
250.985680.986660.988710.985210.986230.988370.982870.984170.990060.982370.983730.98978
Table 16. Percentage analysis of the effect of using capacitors and phase-balancing in three-phase asymmetric networks.
Table 16. Percentage analysis of the effect of using capacitors and phase-balancing in three-phase asymmetric networks.
IEEE 8-Bus Grid
ApproachLosses (kW)Reduction (%)
Benchmark case13.992520.00
Phase-balancing only10.5868624.34
Capacitors only12.1732213.00
CSM8.7780437.27
SSM8.7780437.27
IEEE 25-bus grid
ApproachLosses (kW)Reduction (%)
Benchmark case75.420590.00
Phase-balancing only72.342034.08
Capacitors only54.0841028.29
CSM51.2484632.05
SSM50.1396133.52
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Medina-Gaitán, D.F.A.; Rozo-Rodriguez, I.D.; Montoya, O.D. Optimal Phase-Balancing in Three-Phase Distribution Networks Considering Shunt Reactive Power Compensation with Fixed-Step Capacitor Banks. Sustainability 2023, 15, 366. https://doi.org/10.3390/su15010366

AMA Style

Medina-Gaitán DFA, Rozo-Rodriguez ID, Montoya OD. Optimal Phase-Balancing in Three-Phase Distribution Networks Considering Shunt Reactive Power Compensation with Fixed-Step Capacitor Banks. Sustainability. 2023; 15(1):366. https://doi.org/10.3390/su15010366

Chicago/Turabian Style

Medina-Gaitán, Daniel Federico A., Ian Dwrley Rozo-Rodriguez, and Oscar Danilo Montoya. 2023. "Optimal Phase-Balancing in Three-Phase Distribution Networks Considering Shunt Reactive Power Compensation with Fixed-Step Capacitor Banks" Sustainability 15, no. 1: 366. https://doi.org/10.3390/su15010366

APA Style

Medina-Gaitán, D. F. A., Rozo-Rodriguez, I. D., & Montoya, O. D. (2023). Optimal Phase-Balancing in Three-Phase Distribution Networks Considering Shunt Reactive Power Compensation with Fixed-Step Capacitor Banks. Sustainability, 15(1), 366. https://doi.org/10.3390/su15010366

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop