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):
where
is the complex average power consumption at node
k, which has a complex voltage
, and a net current injection
. Note that
represents the conjugate value of the complex variable/parameter
A, and
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:
According to [
36], the relationship between currents and voltages for a three-phase system is given by the following equations:
where
,
, and
are defined as follows:
where
,
, and
are the a-, b-, and c-phase voltages at node
k;
,
, and
are the net injected currents in phases a, b, and c at node
k; and
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:
In addition, if the nodal admittance matrix (
6) is substituted into (
7), the results of (
8) and (
9) yield the following:
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:
where
is a complex variable that contains the net current injection in the generation node,
is a complex vector that contains all the demanded currents of the network,
is a complex variable that contains the voltage output of the generation source, and
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:
where
and
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;
is the three-phase voltage output of the generation sources, and
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:
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 (
).
where
t is defined as the iterative counter, and the evaluation of the recursive formula reaches convergence when the criterion in (
15) is met:
As recommended by [
5], the value assigned for the
-parameter is
.
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 = 1: do - 4:
if 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 = 1: do - 20:
if 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 = 1: do - 3:
if 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]. 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 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.