SciELO - Scientific Electronic Library Online

 
vol.29 issue83Kinematics analysis of a hand prosthesis prototype for various types of gripsOccurrence of iron ores hosted in Los Tábanos Rhyodacite of San Juan del Cesar (La Guajira), Sierra Nevada of Santa Marta massif: petrologic significance author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Tecnura

Print version ISSN 0123-921X

Tecnura vol.29 no.83 Bogotá Jan./Mar. 2025  Epub Feb 02, 2026

https://doi.org/10.14483/22487638.20582 

Investigación

Application of the arithmetic optimization algorithm to the phase -balancing problem in three-phase asymmetric distribution networks

Aplicación del algoritmo de optimización aritmética al problema de balance de fases en sistemas de distribución asimétricos

Nelson Omero Tique Tique1  * 
http://orcid.org/0000-0001-7877-5525

Juan Camilo Castillo Sáenz2 
http://orcid.org/0000-0002-7993-1726

Oscar Danilo Montoya Giraldo3 
http://orcid.org/0000-0001-6051-4925

1 Electrical Engineering Student, Universidad Distrital Francisco José de Caldas. Bogotá DC, Colombia.Email: notiquet@correo.udistrital.edu.co

2 Electrical Engineering Student, Universidad Distrital Francisco José de Caldas. Bogotá DC, Colombia.Email: juccastillos@correo.udistrital.edu.co

3 Electrical Engineer, Master’s in Electrical Engineering, PhD in Engineering. Associate Professor, Department of Engineering, Universidad Distrital Francisco José de Caldas. Bogotá DC, Colombia.Email: odmontoyag@udistrital.edu.co


Abstract

Objective:

This article presents the application of the arithmetic optimization algorithm (AOA) to the phase balancing problem in three-phase asymmetric distribution systems. The algorithm is implemented in MATLAB and its performance is compared with the losses from the results of the networks presented in specialized literature using other optimization algorithms. The objective is to determine the optimal load connections at each node to reduce the active power losses (the objective function).

Methodology:

Based on the specialized literature, the IEEE standard test networks with 25 and 37 nodes are selected. A master-slave strategy is proposed, where the iterative swept power flow is programmed in its three-phase version (slave stage). Power losses are characterized under a standard phase connection condition on the load side. The arithmetic op-timization algorithm (AOA) (master stage) is then programmed and the most optimal connections to be implemented are defined. These are subsequently compared and analyzed with other methodologies and with the base case.

Results:

Using the selected IEEE test systems of 25 and 37 nodes from the specialized literature, the proposed master-slave methodology, incorporating AOA and the iterative sweep power flow in its three-phase version, effectively identifies an optimal solution. It determines the connections for the phases that balance the load in the nodes, reducing the losses of these systems by 4.155% and 19.249%, respectively.

Conclusions:

The study found that the AOA is effective in minimizing power losses in asymmetrical distribution systems, sligthly outperforming the CSA and CBGA optimization methods. The algorithm demonstrated robustness in solution accuracy and proved useful in validating results using standardized IEEE test systems. The findings highlight the importance of implementing efficient optimization techniques. Future work is suggested to explore the application of AOA to other electrical problems and to develop variants that enhance its parameters and reduce processing times.

Keywords: phase-balancing; three-phase power flow; load connection; arithmetic optimization algorithm; power los

Resumen

Objetivo:

Este artículo presenta la aplicación del algoritmo de optimización aritmética (AOA) al problema de balance de fases observado en los sistemas de distribución. Se utiliza para su implementación el software MATLAB. Los resultados se comparan con las pérdidas reportadas en las redes presentadas en la literatura especializada sobre estudios que aplicaron otros algoritmos de optimización. El objetivo es encontrar la conexión óptima de las cargas en cada nodo para reducir las pérdidas de potencia activa (Función objetivo).

Metodología:

Con base en la literatura especializada, se seleccionan las redes de prueba estándar IEEE de 25 y 37 nodos. Se propone seguir una estrategia maestro-esclavo, en la cual se programa el flujo de potencia barrido iterativo en su versión trifásica (etapa esclavo). Se caracterizan las pérdidas de potencia para una condición estándar de conexiones de fases en la carga y posteriormente se programa el algoritmo de optimización aritmética (AOA) (etapa maestro). Producto de ello, se definen las conexiones más óptimas a implementar, las cuales son luego comparadas y analizadas con otras metodologías y con el caso base.

Resultados:

Usando los sistemas IEEE de prueba seleccionados de 25 y 37 nodos a partir de la literatura especializada, se comprueba que la metodología propuesta maestro-esclavo, la cual emplea el algoritmo de optimización aritmética (AOA) y el flujo de potencia de barrido iterativo en su versión trifásica, logra encontrar una solución óptima. Asimismo, logra definir conexiones para las fases que balancean la carga en los nodos, reduciendo así las pérdidas de dichos sistemas en un 4.155% y 19.249% respectivamente.

Conclusiones:

El estudio encontró que el AOA es eficaz para minimizar la pérdida de potencia en sistemas de distribución de energía asimétricos, con un rendimiento ligeramente mejor que los métodos de optimización CSA y CBGA comparados. El algoritmo resultó ser robusto en cuanto a la precisión de la solución y útil para validar los resultados utilizando sistemas de prueba IEEE estandarizados. El estudio destaca la importancia de aplicar técnicas de optimización eficientes y propone trabajos futuros sobre la aplicación del AOA a otros problemas eléctricos, así como el desarrollo de variantes para mejorar sus parámetros y reducir los tiempos de procesamiento.

Palabras clave: balance de fases; flujo de potencia trifásico; conexión de cargas; algoritmo de optimización aritmética; pérdidas de potencia

Introduction

Distribution systems are a fundamental part of the electric power system; however, they account for the highest percentage of losses Swapna & Udaykumar, 2016. The growing electric power consumption, driven by various human activities, has led to an increased demand for a better quality of service (Carrión González, 2011). This increase in load also generates more power losses and risks to the system, such as overloading, which results in higher costs to the end user (Zhu et al., 1998). An increase in unbalanced load causes unbalanced voltages, even when the source voltage is balanced, which can damage equipment and reduce its useful life (Tuppadung & Kurutach et al., 2006). Overall, these situations lead to a reduction in the efficiency of the distribution system (Zhu et al., 1998).

A direct and efficient solution to mitigate the adverse effects of this issue is phase balancing, which involves reconfiguring or interchanging the phases connected to the load (Swapna & Udayku-mar, 2016). Several methods have been used to determine this phase reconfiguration, including meta-heuristic algorithms aimed at reducing the total losses of the distribution system. To achieve this, the losses for multiple-phase connection configurations are calculated using an iterative process. The configuration with the lowest percentage of losses is considered the optimal solution for load reconfiguration in demand nodes. The approach to finding an optimal solution varies depending on the method and algorithm used. In this study, the arithmetic optimization algorithm is implemented to address the phase balancing problem (Granada Echeverri et al., 2012).

Motivation

An electric power system consists of several stages, with the primarily objective of transporting energy from generation centers to the end users, delivering it economically, safely, and with the required quality levels (Ramirez Castaño, S., 2019). One such stage, energy distribution, constitutes a significant part of the electrical power system. As a result, it has received considerable attention in terms of quality regulation, enabling network operators to manage service conditions in favor of continuous improvement. This focus has driven ongoing processes of evolution and improvement, reflected in the type of techniques, equipment, and tools used for this purpose (CREG 2008). In this context, the present research aims to apply a new technique to the PB problem, seeking to minimize the power losses caused by three-phase distribution systems. For this purpose, the phase-balancing problem is addressed using the arithmetic optimization algorithm (AOA), with the goal of determining the set of optimal connections for the system loads, thus reducing the phase unbalance and finding an optimal operating point.

State of the art

The phase-balancing problem (PBP) has been addressed through various methods in the current literature. Zhu et al., 1998 formulated the optimal phase-balancing problem using linear-constrained integer programming (approximated modeling), aiming to reduce the unbalance in line currents in distribution systems. The authors (Zhu et al., 1999) proposed a solution to PBP using a metaheuristic algorithm called Simulated Annealing, seeking to minimize total operating costs. Chen & Cherng formulated a solution for the PBP using a genetic algorithm optimization technique. Their objective was to minimize the phase imbalance in the whole system. (Gandomkar et al., 2004) solved the PBP problem using a genetic-based optimizer. They sought to reduce power losses in distribution systems. In the study conducted by (Lin et al., 2005), PBP was resolved by employing a heuristic search algorithm with backward recurrence, which minimized the phase unbalance and reduced the neutral current. (Garcés-Ruiz et al., 2005) developed a solution methodology for the PBP using an ant colony optimization algorithm, which helped them minimize energy losses over a 24-hour period in distribution systems.

(Tuppadung & Kurutach et al., 2006) proposed a particle swarm optimization algorithm to solve the PBP. They reduced power losses and improved voltage profiles. In (Lin et al., 2008), the PBP was formulated using an expert system based on heuristic rules, aimed at minimizing system operating costs. (Huang et al., 2008), in turn, solved the PBP using an immunological optimization algorithm in which they reduced the operating costs of the system over a 24-hour period.

(Kuo & Chao 2010) approached the PBP using automatic phase shifting by mapping the information of the system—taking into account technical characteristics and load flow—to determine the phase rearrangement that results in the lowest power losses. (Srinivasan et al., 2011) and (Sathiskumar et al., 2012) resolved the PBP by developing a fuzzy hybrid greedy heuristic algorithm and a self-adaptive hybrid differential evolution algorithm, respectively. The main goal was to minimize the phase unbalance of the system based on load flow equations, capacity, and voltage restrictions, thus balancing the current flowing through the phases of the system and reducing the current flowing through the neutral. (Hooshmand & Soltani et al., 2012) proposed a methodology to determine PBP based on bacterial foraging guided by a particle swarm optimization algorithm. The study aims at minimizing the operating costs of the system, whether radial or meshed. In addition, (Granada Echeverri et al., 2012) solved the PBP employing a genetic-based optimizer. With this, their main goal was to minimize energy losses in distribution systems. As a result, there was evidence of improved voltage profiles and reduced congestion on the lines of the system.

(Grigoras, & Gavrilas, et al., 2016) approached the PBP using a heuristic method based on phase exchange for load balancing, aiming to minimize the degree of system unbalancing. Swapna & Udaykumar formulated the PBP on the IEEE 13-node test system using a balancing algorithm that evaluates possible phase exchange combinations for each candidate node until the combination that balances the system is chosen. (Nicusor Toma et al., 2018), in turn, presented a methodology for solving the PBP using the particle swarm optimization algorithm, with the primary objective of minimizing voltage imbalance in the distribution network. Fadhela & Abdulmuttalib T., proposed minimizing total voltage unbalance in distribution systems through a heuristic search algorithm, which determines how to change phases using contactors located in the load terminals.

(Gangwar & Chakrabarti et al., 2019) adopted an analytical approach to the problem, ensuring that phase interchanges were kept to a minimum. The objective was to reduce unbalancing in distribution systems and power losses, as a result. (Ovidiu I et al., 2009) proposed a different solution by employing a genetic optimizer to reduce energy losses over a 24-hour period (daily operation scenario).

(Cortés-Caicedo et al., 2021a) solved the PBP problem through a metaheuristic optimization technique called the vortex search algorithm, which requires discrete codification. The objective of the function was to minimize the power loss of the distribution system.

Finally, (Sarmiento & Ávila et al., 2022) offered an alternative solution for the PBP. Using a discrete version of the whale optimization algorithm, they aimed to reduce active power losses in all branches of the distribution system.

It is important to note that the main characteristics of the previously reviewed literature are as follows:

  1. Most optimization methodologies are based on combinatorial optimization approaches. These works under a master-slave optimization concept, i.e., the metaheuristic optimizer defines the nodal connections. The slave stage determines the value of the objective function using a three-phase unbalanced power flow approach.

  2. The most common objective function corresponds to minimizing expected grid power losses. This objective is effective in modeling the problem of load unbalance in asymmetric distribution systems due to the nonlinear relationships between voltages and demanded currents.

In summary, the PBP encompasses multiple objectives and solution techniques that lead to the search for improved power quality delivered by the grid operator to the end users in an efficient and cost-effective manner. According to the previously mentioned, the arithmetic optimization algorithm (AOA) has not been applied as a solution alternative for the PBP in asymmetric distribution systems. In this study, we will propose to use it to select the best load connection combination as a contribution to the current literature. It is also worth noting that, within the scope of this work, and consistent with the literature review, we adopt the minimization of power losses under peak load conditions as the objective function. In addition, all the impedance parameters of the systems are well-known, and the peak load information and their connections correspond to input data for the optimization model.

The remainder of this research is structured as follows: Section Mathematical formulation presents the general mathematical formulation of the PBP for three-phase unbalanced distribution grids. The proposed solution methodology based on the application of the AOA is described in Section Proposed solution methodology; Section Test Systems shows the 25 and 37 node test systems used for the AOA implementation; In section Results and comparisons, the results are analyzed and compared with other optimization algorithms; and finally, the section Conclusions and future works presents the conclusion of this study as well as proposing essential points for future work derived from this study.

Mathematical formulation

A mixed-integer nonlinear programming (MINLP) model represents the PBP in asymmetric distribution networks. Binary variables are the decision variables, which correspond to the connections established for each load presented in the system (Arias et al., 2019). On the other hand, continuous variables are associated with the classical power flow formulation, which uses magnitude and voltage angle for each node (Montoya et al., 2020b). Finally, the nonlinear nature is the result of the trigonometric functions present in the power balance equations and the products between the magnitudes of the nodes (Gil-González et al., 2020).

The MINLP formulation that describes the PBP is formulated below:

Objective Function

In this research, the objective function for the PBP aims at minimizing the total active power loss of the asymmetrical distribution system in (1).

where z represents the value of the objective function; Ω p is the set associated with the phases of the system; Ω N is the set associated with the nodes of the system. V p i and V q j are the voltage magnitudes and the phase p and phase q in the nodes i and j, respectively. δ p i and δ p j are the angles of the voltage

at the phase p in the nodes i and j, respectively. Y pq , is the magnitude of the admittance associated with the line connected between the nodes i and j in the phases p and q with an angle denominated as θ pq ij .

Restrictions

The set of constraints corresponds to the operational limitations found in any asymmetric distribution system. These are shown in equations (2) to (11). Equality constraints (2) to (4) correspond to the active power equilibrium at each node and phase.

where P gA , P gB , and P gC are the active power injected by a conventional generator in phases A, B, and C at each node i, respectively. P dA , P dB , and P dC are the active power demanded in each type of connection h, in phases A, B, and C at node i, respectively. X ih represents the binary decision variable, which takes the value of 1 if the connection h is selected for node i and in any other case.

Observe that the binary variable X ih will have to select among h possible connections for the loads of demand nodes i present in the distribution system, as shown in Table 1.

Table 1 Possible types of connections for the loads presented in the system 

Figure 1 presents a possible change in connection from the original ABC connection to the BAC connection, which implies a sequence change.

Source: Authors

Fig. 1 Phase-swapping according to the type of connection 

It is important to mention that if there are no three-phase loads on some of the nodes, then for these phases a zero value is assigned in the connection (Granada Echeverri et al., 2012).

The set of equality constraints (5) to (7) shows the reactive power balance constraint per node and phase, respectively.

where Q gA i , Q gB i , and Q gC i are the reactive power injected by a conventional generator in the phases A, B, and C at the i node, respectively. Q dA ih, Q dB ih, and Q dC ih , are the reactive powers demanded for the type of connection h, in phases A, B, and C, at the i node, respectively.

The set of box-type constraints (8) to (10) corresponds to the voltage regulation constraints per phase at node i at each phase.

where V min and V max correspond to the minimum and maximum voltage regulation limits permitted by regulatory entities in medium-voltage levels.

Finally, the equality constraint (11) defines the unique nature of each load connection per phase, that is, only one of the load combinations in Table 1 can be assigned to a particular node i.

Remark 1 The MINLP model that represents the PBP in three-phase asymmetric distribution networks defined from (1) to (11) is only applicable to three-phase networks with solidly grounded connection in the neutral wire at each demand node being all the constant power load connected in star form. In the case of loads connected in triangle, additional equations must be added to the optimization model.

Proposed solution methodology

In this research, the PBP in three-phase asymmetric distribution networks is addressed through the application of a master-slave optimization approach, where AOA is implemented in the master stage to define nodal load connections. The three-phase backward/forward power flow method is used in the slave stage to determine the value of the power losses of the system. The main aspects of the proposed optimization approach are described below.

Master stage: Arithmetic optimization algorithm

The AOA is an optimization algorithm recently developed that explores and exploits the solution space using the distribution of the four main arithmetic operations (multiplication, division, subtraction. and addition) (Hraiz et al., 2020). This algorithm can be classified as an optimizer based on population that evolves from a solution set through arithmetic advance rules to a final solution. The initial population is represented in (12) where x ij is the value of the variable j for the i th individual in the population. Note that n v is the number of variables, and n i is the number of individuals.

To explore the solution space, the AOA uses the acceleration function 13, which depends on its maximum and minimum values (MOA min) and (MOA max) as well as the current iteration counter t and the maximum number of iterations t max.

Note that MOAt can be considered a variable optimization factor that enhances the exploration and exploitation characteristics of the AOA.

Exploration phase

During the exploration phase, the AOA uses properties of the arithmetic operator division and multiplication, exploring the search area randomly in several regions and getting closer to finding a solution and generating the position of the individual ith using the evolution rule defined in (14).

The MOA function conditions this search phase, i.e., (r1 > MOAt).

where x t+1 indicates the j th position of the i th solution in the next iteration, and x t is the j th

ij bestj

position of the best solution obtained at the moment; x max and x min are the upper and lower limit

j j

value of the j th position, respectively and MOP t is a variable factor that depends of the iteration number with the structure of (15):

where β is a sensitive factor that helps improve exploration accuracy, which is set as 0.2 based on the experiments made by the authors (Abualigah et al., 2021). Note that ϵ is a factor that avoids divisions by zero, but it is small enough not to affect the results and can be selected the same as the tolerance ε used to test the convergence in power flow; moreover, r 2 is control parameter set as 0.5 based on the recommendations of (Abualigah et al., 2021).

Exploitation phase

During the exploitation phase, the AOA uses the advantages of the arithmetic operators, subtraction, and addition, since these produce small variations around the current solution X t , which allows it to get closer quickly to the nearest optimal solution that can be deduced after several iterations. The rule in this exploitation phase is defined in (16).

The evolution rules in (16) show that the new individual gets its new values by adding or subtracting positions (depends on the value of the random number r 3 which varies between 0 and 1 with uniform distribution), which are considered smooth variations around the current position, which in heuristics means, the exploitation of the solution space.

The procedures in this phase are similar to the partitions in the previous phase. However, the subtraction and addition operators usually prevent the optimization algorithm from getting stuck in a local nature solution (55). This procedure helps exploratory search strategies to find the optimal solution and maintain the diversity of candidate solutions.

It should be noted that each new individual generated in the exploration and exploitation phases must be reviewed to ensure that each new position x t+1 is between its maximum and minimum limits, i.e., x max and x min to ensureiability of the solution space.

Moreover, the first positions of the new solution vector X t+1 , are rounded close to integer values since they correspond to the nodes where the set of connections for the loads on each node will be defined; i.e., they must always be integer values. The application of the described AOA can be performed following Algorithm 1

Source: Authors

Slave Stage: Three-phase iterative sweep power flow

An acknowledged method for solving the power flow problem in distribution networks is the backward/forward power flow approach (also known as the iterative sweep method). This is based on representing the topology of the network through an incidence matrix that relates the connections between the nodes of the system and also represents the branches present in the network (Cortés-Caicedo et al., 2021a). Using the laws of Kirchoff, the currents in the nodes of the system are calculated, making a sweep from the end nodes to the generation nodes. From these obtained currents, the voltage drops in each segment of the network are calculated, this time sweeping from the generation node to the terminals (Shen et al., 2018).

The following equation solves the three-phase power flow (17), and taking as reference a distribution system with n nodes:

where the matrix Y dd3 represents the system’s topology, i.e., the connections between the demand nodes of the system. This matrix is square in dimension: 3(n − 1) × 3(n − 1). The vector I d3 represents the demand currents in each node, calculated from the power in each node and the type of connection, star or delta, at those nodes. The dimension of this matrix is 3(n − 1) × 1. The rectangular matrix Y d03 represents the system topology considering the generation and demand terms. The dimension of this matrix is 3(n − 1) × 3. Note that V 03 is the vector that represents the voltages for each phase of the system generation node with a dimension of 3 × 1. The vector V d3 with dimensions 3(n − 1) × 1, represents the three-phase voltages demanded in each of the load nodes of the system.

The algorithm starts defining parameters, such as the number of iterations and the maximum acceptable error allowed for the algorithm to stop. For the maximum iterations, a value of m max = 1000 iterations is defined, while for the maximum error allowed, named ϵ, a value of 1 × 10−10 is defined.

Remark 2 The iterative sweep power flow method in (17) depends on the type of load per node (or Y connection), i.e., the calculation of the demand current I m at the iteration m. To study in detail its algorithmic implementation, see Ref. (Shen et al., 2018).

Test systems

In this section, in order to show the efficiency of the optimization method proposed in the methodology section, the information of two unbalanced test systems is used, present in the specialized literature. The test systems correspond to the IEEE 25 and 37 node radial distribution networks reported by (Cortés-Caicedo et al., 2021a) to study the phase balancing problem. The features of the test systems are shown below.

IEEE 25-node test system

This is a radial unbalanced distribution system with 25 nodes, 24 lines, and 22 loads. The source is located at node 1 and the nominal voltage is 4.16 kV. The electrical configuration of this system can be seen in Figure 2. The total active and reactive power consumption is 1073 kW, and 792 kvar for phase A, 1083.3 kW and 801 kvar for phase B, and 1083.3 kW and 800 kvar for phase C. Load consumption per phase and impedance matrices per type of conductor are shown in the work of (Cortés-Caicedo et al., 2021a). Note that the electrical parameters, i.e., impedance, line lengths, and load values can be consulted in (Cortés-Caicedo et al., 2021a).

Source: Authors

Fig. 2 Single-phase diagram of the IEEE 25-bus grid. 

IEEE 37-node test system

The IEEE 37-node test system is a real radial unbalanced distribution system located in California, totally composed of underground lines. It has 37 nodes, node 1 being the source, 35 lines, and 25 loads. It has a transformer that operates in an unbalanced form and a voltage regulator. The nominal voltage is 4.8 kV. The electric configuration of this system can be seen in Figure 3. The total active and reactive power consumption is 727 kW and 357 kvar for phase A, 639 kW and 314 kvar for phase B, and 1091 kW and 530 kvar for phase C. The load consumption per phase and the impedance matrix information per type of conductor can be found in (Cortés-Caicedo et al., 2021a).

Source: Authors

Fig. 3 Single-phase diagram of the IEEE 37-bus grid 

Results and comparisons

In this section, the proposed master-slave methodology is tested. For comparative purposes, the AOA with the iterative sweep power flow proposed in this research is compared with two combinatorial optimization methods available in the current literature for solving the phase-balancing problem, i.e., the Chu & Beasley genetic algorithm (CBGA) (Montoya et al., 2021), and the crow search algorithm (CSA) (Cortés-Caicedo et al., 2021b), respectively. The following sections present the main results for the benchmark case of each test feeder (the power flow solution) as well as the results of the proposed optimization approach with analysis and discussion. Observe that all the numerical validations are carried out in the MATLAB software in its 2020a version on a laptop computer, Intel(R) Core(TM) i5-6200U @2.30 GHz processor with 12.00 GB of RAM running Windows 10 Home Single Language operating system running at 64 bits.

Parametrization of the proposed master-slave optimization method

To implement the proposed master-slave optimization method based on the hybridization of the AOA approach with the iterative sweep power flow method and its application to the problem of optimal phase balancing, the list of parameters in Table 2 is considered.

Table 2 Parametrization of the proposed master-slave optimization method 

Solution of the three-phase power flow problem

Table 3 presented the results obtained by implementing that method. Note that the power losses, iterations, and processing time are the same as those reported by (Cortés-Caicedo et al., 2021a).

Table 3 Power flow solution in the IEEE 25- and 37-bus grids with the iterative sweep power flow approach. 

The main characteristic of the numerical results in Table 3 is that for both test feeders, the efficiency of the iterative sweep power flow method was compared with the Newton-Raphson power flow method for unbalanced distribution grids available in the DIgSILENT software.

Results for the 25-Node test system

Table 4 presents the results obtained by the proposed AOA methodology in the 25-bus grid. It shows that the AOA is 4.1556% more efficient than the standard connection present in the benchmark case, resulting in a reduction of 3.1340 kW. Note that the proposed AOA is only 7.46 × 10−3 more effective in minimizing power loss compared to the CBGA and 0.0595% compared to the CSA.

Table 4 Performance of the power losses after implementing the phase-swapping in the 25-Node Test System 

The small difference between the proposed and compared optimizers in this test feeder shows that, in general, the phase-balancing problem can be solvable with different metaheuristics efficiently with minimum differences in the final solution value with respect to total grid power losses; nevertheless, more research is required since due to the large dimension of the solution space, it is impossible to ensure the global optimum finding, which implies that the results in Table 4, are local optimal solutions.

On the other hand, Figure 4 compares before and after implementing the AOA phase swapping methodology in the 25-node test system. It is possible to observe that 11.38 kW increased the losses in phase B, representing a 76.99% more power loss in that phase. In addition, reductions of about 11.22 kW and 3.3 kW were found for phases A and C, respectively, representing an improvement of approximately 30. 423% and 13. 88% in each phase.

Source: Authors

Fig. 4 Effect of the optimal phase-swapping on power losses in the 25-node test system for each phase. 

Figure 5 shows a comparison between the methodologies presented in Table 4 in terms of the percentage of power loss reduction. Note that all methodologies have a 4% minimum reduction, and AOA obtains a 4. 155% reduction, which means that it is a robust methodology capable of obtaining efficient results.

Source: Authors

Fig. 5 Reduction of the power losses in the 25-node test system for each combinatorial optimization method 

To better understand the performance obtained with AOA, Figure 6 shows the voltage profile for the phases A, B, and C, respectively. Note that Figure 6 shows that the voltage per load node follows a similar behavior to the benchmark case, in general terms, the unbalance present in the benchmark case has been reduced except in phase B, where the voltage drops increase compared to two other phases A and C where the voltage drops are reduced in those phases, phase C drops close to the benchmark case.

Source: Authors

Fig. 6 Voltage phase profiles of the 25-node test system, comparison of the benchmark case and the proposed optimization method 

Furthermore, after 100 consecutive evaluations of the proposed methodology AOA for the 25-node test system, the best solution is 72.2865 kW, and the CBGA found the best solution after doing 550 iterations is 72.2919 kW. In terms of processing time, it is important to note that the AOA took 55 seconds to find the optimal solution, 87.13% longer than the CSA and 194. 86% longer than the CBGA, which, although in practice does not represent something significant, the difference is in seconds of execution which, in comparison to the solution space of size 1.3162 × 1017, shows considerable efficiency in time and provides an excellent opportunity to analyze better ways of exploiting and exploring the solution space using the AOA to improve the performance of the same since this is a relatively recently developed optimization algorithm.

Results for the 37-Node Test System

Table 5 presents the results obtained from the AOA for the 37-bus system. The results and comparisons show that AOA reduces the total power losses by 19.2498% from the benchmark case, which means a loss reduction of 3.1319 kW. In terms of CBGA, AOA is 0.1604% better than CBGA, where AOA reduces losses by 98.8 W more than CBGA. With the CSA, the AOA is 0.2867% better than the CSA, where the AOA reduces losses by 176.8 W more than the CSA.

Table 5 Performance of the power losses after implementing the phase-swapping in the 37-node test system 

Figure 7 compares the power loss reduction per phase before and after the phase-swapping proposed by the AOA. Phases A and C present a power loss reduction, of 11.22 kW and 3.30 kW, respectively. Meanwhile, phase B presents a power loss increase of 11.38 kW. The total power losses are obtained by adding the power losses per phase, and this was previously presented in Table 5.

Source: Authors

Fig. 7 Effect of the optimal phase-swapping on power losses for the 37-node test system. 

Figure 8 shows the loss reduction in percentage terms. This compares the benchmark case with the other optimization algorithms and the AOA. The CBGA, CSA, and AOA have a loss reduction of 19.1201%, 19.0186%, and 19.2498% from the benchmark case, respectively. The optimization algorithms have similar behavior, all of them reduced the power losses by about 19%, which represents a reduction of more than 14 kW in total power losses in the system.

Source: Authors

Fig. 8 Reduction of the power losses in the 37-node test system. 

The voltage profiles per node of the system are presented in Figure 9, this complements the results presented in Figure 7. It shows the voltage behavior on each node compared with the benchmark case. The effects of phase-swap increase the voltage profiles on phases A and C. For phase C, all the node’s voltage profiles are better than the benchmark case, but in phase A, this does not occur, the voltage profile in some nodes decreases, as in nodes 25, 26, 27, 28, and 29. The opposite occurs in phase B, here, the voltage profiles decrease, which is caused by the power losses increment in this phase. Regarding processing times, the average time to processing was 178.75 s after 100 consecutive evaluations. The 37 node system has 2.843 × 1019 possible combinations for its 25 loads, which represents an enormous space solution. This means that the AOA works well for large systems.

Source: Authors

Fig. 9 Voltage phase profiles of the 37-node test system when compared the proposed approach with the benchmark case 

Conclusions and future work

The results obtained in this comparative study show that the AOA algorithm is highly effective in optimizing asymmetrical distribution power systems in terms of power loss minimization. The compared optimization methods, CSA and CBGA, also showed significant improvements in system efficiency, but the AOA shows a slightly better performance in solving power loss minimization problems. In particular, the AOA algorithm was found to be robust in terms of solution accuracy, suggesting that it can be a valuable tool for optimizing asymmetric distribution power systems in practice. Furthermore, the IEEE 25-node and 37-node test systems proved to be useful in validating the results obtained in simulations of real distribution power systems, suggesting that these models are a valuable tool for research and development in this field. In general, these findings highlight the importance of implementing efficient optimization techniques in asymmetric distribution power systems and suggest that the use of standardized test models can improve the operating characteristics of those systems.

For this research, the arithmetic optimization algorithm was presented to address the phase balance problem in radial distribution networks. With the AOA work, the best set of load connections was found for 25 and 37-node systems, respectively. The optimization methodology consisted of a master-slave solution approach. The slave stage was responsible for determining the operative state of the system as the voltage in all the system nodes, the power losses in each node, and the general operative characteristics of the distribution network under study. The master stage defined the initial conditions for the slave stage, i.e., the set of load connections applicable to each node.

The AOA shows that it can be compared with other optimization algorithms. Even AOA, a recently developed algorithm, showed better numerical results compared to CBGA and CSA. The results for the two test systems, that is, the IEEE 25 and 37 node grids, presented a reduction in total power loss of approximately 4.155 % and 19.249 %, respectively. These results confirmed that the AOA is an efficient optimization algorithm that can be extended to new optimization applications in electrical engineering.

Technological advances come with new challenges. Likewise, tools must be developed to address these problems. The AOA is presented as a new tool for solving problems like the one studied in this research. Once the effectiveness has been evaluated, showing that satisfactory results are obtained with AOA in terms of total loss reduction, efficiency should be mentioned. An average time of 55.0470 s was presented for the IEEE 25-node system, while the average time for the IEEE 37-node grid was 178.7548 s because the processing time for the algorithm in the master and slave stages is higher compared to CBGA and CSA. The algorithm explains this due to the two phases it composes, the exploration and exploitation phases, in addition to the parameters. This means more computational load and, therefore, a longer processing time.

Previously, it was mentioned that AOA is recent, and there is still more development on this algorithm to discover its potential. So, in future work, it is proposed to apply the AOA to electrical problems such as optimal location and sizing of compensators, batteries, and STATCOMS, and optimal selection of conductors or distribution networks. In addition, it is proposed to apply the developed variants of the AOA for its application phase-balancing problem in three-phase asymmetric distribution networks. Likewise, develop variants for AOA to improve the algorithm’s parameters to reduce the computational load and the processing times.

References

Abualigah L, Diabat A, Mirjalili S, Abd Elaziz M, Gandomi AH. The Arithmetic Optimization Algorithm. Computer Methods in Applied Mechanics and Engineering. 2021;376. https://doi.org/10.1016/j.cma.2020.113609Links ]

Arias J, Calle M, Turizo D, Guerrero J, Candelo-Becerra JE. Historical load balance in distribution systems using the branch and bound algorithm. Energies. 2019;127. https://doi.org/10.3390/en12071219Links ]

Carrión González JE. Análisis de métodos heurísticos de reconfiguración de sistemas de distribución para la reducción de pérdidas de potencia. Ingeniería Energética. 2011;3239-17. https://www.redalyc.org/articulo.oa?id=329127748002Links ]

Chen T-H, Cherng J-T. Optimal Phase Arrangement of Distribution Transformers Connected to a Primary Feeder for System Unbalance Improvement and Loss Reduction Using a Genetic Algorithm. Proceedings of the 21st International Conference on Power Industry Computer Applications. Connecting Utilities. PICA 99. To the Millennium and Beyond (Cat. No. 99CH36351). 1999 145-151. https://doi.org/10.1109/PICA.1999.779397Links ]

Correa CA, Bolaños RA, Ruiz Garcés A. Development and material properties of new hybrid medium density fibreboard from empty fruit bunch and rubberwood. Scientia et Technica. 2007;Año XIII37. https://hdl.handle.net/11059/270Links ]

Correa M OI. Estudio de reconfiguración y optimización de los alimentadores de la subestación Machala perteneciente a la Corporación Nacional de Electricidad SA-Regional El Oro. Bachelor’s thesis. Universidad Politécnica Salesiana. 2010. http://dspace.ups.edu.ec/handle/123456789/311Links ]

Cortés-Caicedo B, Avellaneda-Gómez LS, Montoya OD, Alvarado-Barrios L, Chamorro HR. Application of the vortex search algorithm to the phase-balancing problem in distribution systems. Energies. 2021;1451282. https://doi.org/10.3390/en14051282Links ]

Cortés-Caicedo B, Avellaneda-Gómez LS, Montoya OD, Alvarado-Barrios L, Álvarez-Arroyo C. An Improved Crow Search Algorithm Applied to the Phase Swapping Problem in Asymmetric Distribution Systems. Symmetry. 2021;1381329. https://doi.org/10.3390/sym13081329Links ]

Comisión de Regulación de Energía y Gas. Cartilla Distribución de Energía Eléctrica. Bogotá D.C, Colombia: Nuevas Ediciones S.A; 2008. https://www.creg.gov.co/sites/default/files/distribucin%2520de%2520energa%2520elctrica.pdfLinks ]

Handhal FK, Rashid AT. Load Balancing in Distribution System Using Heuristic Search Algorithm. 2018 International Conference on Advance of Sustainable Engineering and its Application (ICASEA). 2018 48-53. https://doi.org/10.1109/ICASEA.2018.8370954Links ]

Fahim M, Hassan E, Najjar BE. Single Phase Load Balancing in a Three Phase System at Distribution and Unit Level. IEEE International Conference on Industrial Technology (ICIT). 2018 1297-1301. https://doi.org/10.1109/ICIT.2018.8352365Links ]

Gandomkar M. Phase balancing using genetic algorithm. 39th International Universities Power Engineering Conference. 2004 1377-379. https://doi.org/10.1109/TPWRS.2008.926472Links ]

Gangwar P, Singh SN, Chakrabarti S. An Analytical Approach for Phase Balancing Considering Customer Load Profile. Proceedings of 2019 IEEE PES Innovative Smart Grid Technologies Europe, ISGT-Europe 2019. . 2019 https://doi.org/10.1109/ISGTEurope.2019.8905481Links ]

Garces A. A Linear Three-Phase Load Flow for Power Distribution Systems. IEEE Transactions on Power Systems. 2016;311827-828. https://doi.org/10.1109/TPWRS.2015.2394296Links ]

Garcés-Ruiz A, Granada-Echeverri M, Gallego-R RA. Balance de fases usando colonia de hormigas. Ingeniería y Competitividad. 2005;7243-52. http://hdl.handle.net/10893/1594Links ]

Gil-González W, Montoya OD, Rajagopalan A, Grisales-Noreña LF, Hernández JC. Optimal selection and location of fixed-step capacitor banks in distribution networks using a discrete version of the vortex search algorithm. Energies. 2020;13184914. https://doi.org/10.3390/en13184914Links ]

Gonzáles G EF, Morante B JA, Vicuña T WF. Estudio de un sistema de distribución y acometidas en baja tensión. tesis de grado. Universidad Politécnica Salesiana. Guayaquil. 2015. http://dspace.ups.edu.ec/handle/123456789/10236Links ]

Granada Echeverri M, Gallego Rendón RA, López Lezama JM. Optimal Phase Balancing Planning for Loss Reduction in Distribution Systems using a Specialized Genetic Algorithm. Ingeniería y Ciencia. 2012;815121-140. https://www.redalyc.org/articulo.oa?id=83524069006Links ]

Grigoras G, Gavrilas M. Phase Swapping of Lateral Branches from Low-Voltage Distribution Networks for Load Balancing. 2016 International Conference and Exposition on Electrical and Power Engineering (EPE 2016). 2016 715-718. http://doi.org/10.1109/ICEPE.2016.7781432Links ]

Grisales Noreña LF, Garzón Rivera OD, Montoya OD, Ramos Paja CA. Hybrid Metaheuristic Optimization Methods for Optimal Location and Sizing DGs in DC Networks. Communications in Computer and Information Science. 2019;1052214-226. https://hdl.handle.net/20.500.12585/9184Links ]

Hooshmand RA, Soltani S. Fuzzy optimal phase balancing of radial and meshed distribution networks using BF-PSO algorithm. IEEE Transactions on Power Systems. 2012;27147-57. https://doi.org/10.1109/TPWRS.2011.2167991Links ]

Hraiz MD, Martín García JA, Jiménez Castañeda R, Muhsen H. Optimal PV Size and Location to Reduce Active Power Losses While Achieving Very High Penetration Level With Improvement in Voltage Profile Using Modified Jaya Algorithm. IEEE Journal of Photovoltaics. 2020;1041166-1174. https://doi.org/10.1109/JPHOTOV.2020.2995580Links ]

Huang MY, Chen CS, Lin CH, Kang MS, Chuang HJ, Huang CW. Three-phase balancing of distribution feeders using immune algorithm. IET Generation, Transmission and Distribution. 2008;23383-392. https://doi.org/10.1049/iet-gtd:20070206Links ]

Kaveh A, Biabani Hamedani K. Improved arithmetic optimization algorithm and its application to discrete structural optimization. Structures. 2022;35748-764. https://doi.org/10.1016/j.istruc.2021.11.012Links ]

Kharrich M, Abualigah L, Kamel S, AbdEl-Sattar H, Tostado-Véliz M. An Improved Arithmetic Optimization Algorithm for design of a microgrid with energy storage system: Case study of El Kharga Oasis, Egypt. Journal of Energy Storage. 2022;51. https://doi.org/10.1016/j.est.2022.104343Links ]

Kuo CC, Chao YT. Energy management based on AM/FM/GIS for phase balancing application on distribution systems. Energy Conversion and Management. 2010;513485-492. https://doi.org/10.1016/j.enconman.2009.10.011Links ]

Laconico KCC, Aguirre Jr RA. Optimal Load Balancing and Capacitor Sizing and Siting of an Unbalanced Radial Distribution Network. IEEE PES GTD Grand International Conference and Exposition Asia (GTD Asia) 2019. 2019 939-944. https://doi.org/10.1109/GTDAsia.2019.8715965Links ]

Lin CH, Chen CS, Chuang HJ, Ho CY. Heuristic rule-based phase balancing of distribution systems by considering customer load patterns. IEEE Transactions on Power Systems. 2005;202709-716. https://doi.org/10.1109/TPWRS.2005.846077Links ]

Lin CH, Chen CS, Chuang HJ, Huang MY, Huang CW. An expert system for three-phase balancing of distribution feeders. IEEE Transactions on Power Systems. 2008;2331488-1496. https://doi.org/10.1109/TPWRS.2008.926472Links ]

Montano J, Garzón OD, Rosales Muñoz AA, Grisales Noreña LF, Montoya OD. Application of the arithmetic optimization algorithm to solve the optimal power flow problem in direct current networks. Results in Engineering. 2022;100654. https://doi.org/10.1016/j.rineng.2022.100654Links ]

Montoya OD, Gil-González W. On the numerical analysis based on successive approximations for power flow problems in AC distribution systems. Electric Power Systems Research. 2020;187. https://doi.org/10.1016/j.epsr.2020.106454Links ]

Montoya OD, Gil-González W, Giral DA. On the matricial formulation of iterative sweep power flow for radial and meshed distribution networks with guarantee of convergence. Applied Sciences (Switzerland). 2020;1017. https://doi.org/10.3390/app10175802Links ]

Montoya OD, Gil-González W, Grisales-Noreña LF. An exact MINLP model for optimal location and sizing of DGs in distribution networks: A general algebraic modeling system approach. Ain Shams Engineering Journal. 2020;112409-418. https://doi.org/10.1016/j.asej.2019.08.011Links ]

Montoya OD, Gil-González W, Orozco-Henao C. Vortex search and Chu-Beasley genetic algorithms for optimal location and sizing of distributed generators in distribution networks: A novel hybrid approach. Engineering Science and Technology, an International Journal. 2020;2361351-1363. https://doi.org/10.1016/j.jestch.2020.08.002Links ]

Montoya OD, Molina-Cabrera A, Grisales-Noreña LF, Hincapié RA, Granada M. Improved Genetic Algorithm for Phase-Balancing in Three-Phase Distribution Networks: A Master-Slave Optimization Approach. Computation. 2021;9667. https://doi.org/10.3390/computation9060067Links ]

Montoya OD, Gil-González WJ, Garcés-Ruíz A. Flujo de potencia óptimo para redes radiales y enmalladas empleando programación semidefinida Optimal Power Flow for radial and mesh grids using semidefinite programming. Tecnológicas. 2017;204029-42. [ Links ]

Toma N, Ivanov O, Neagu B, Gavrila M. A PSO Algorithm for Phase Load Balancing in Low Voltage Distribution Networks. 2018 International Conference and Exposition on Electrical And Power Engineering (EPE). 2018 0857-0862. https://doi.org/10.1109/ICEPE.2018.8559805Links ]

Ivanov O, Neagu BC, Gavrilas M, Grigoras G, Sfintes CV. Phase Load Balancing in Low Voltage Distribution Networks Using Metaheuristic Algorithms. International Conference on Electromechanical and Energy Systems (SIELMEN). 2019 1-6. https://doi.org/10.1109/SIELMEN.2019.8905900Links ]

Quille Pinto FS. Optimización de flujo de potencia en el sistema eléctrico ecuatoriano con programación no lineal bajo MATLAB. tesis de grado. Universidad Politécnica Salesiana. Quito. 2015. https://dspace.ups.edu.ec/handle/123456789/8066Links ]

Ramirez Castaño S. Redes de distribución de energía. Repositorio Universidad Nacional de Colombia sede Manizales. 2019. https://repositorio.unal.edu.co/handle/unal/7095Links ]

Rendón G RA, Granada E M, Escobar Z A. TECHNICAL LOSS REDUCTION USING CORRECTIVES STAGES FITTING - PHASE I. Revista Colombiana de Tecnologías de Avanzada. 2006;1. https://www.unipamplona.edu.co/unipamplona/portalIG/home_40/recursos/02_v07_12/revista_07/16112011/v07_14.pdfLinks ]

Rios MA, Castaño JC, Garcés A, Molina-Cabrera A. Phase Balancing in Power Distribution Systems: A heuristic approach based on group-theory. 2019 IEEE Milan PowerTech. 2019 1-6. https://doi.org/10.1109/PTC.2019.8810723Links ]

Ruiz AG, Manso JCG, Gallego Rendón RA. Solución al problema de balance de fases y reconfiguración de alimentadores primarios bajo un modelamiento trifásico usando Simulated Annealing. Sciencia et Technica. 2006;Año XII37-42. https://www.redalyc.org/articulo.oa?id=84920491059Links ]

Sarmiento Vargas GA, Ávila Furque EJ. Aplicación del algoritmo de caza de ballenas al problema de balance de fases en sistemas de distribución desbalanceados. tesis de grado. Universidad Distrital Francisco José de Caldas. 2022. https://hdl.handle.net/11349/30225Links ]

Sathiskumar M, Nirmal Kumar A, Lakshminarasimman L, Thiruvenkadam S. A self adaptive hybrid differential evolution algorithm for phase balancing of unbalanced distribution system. International Journal of Electrical Power and Energy Systems. 2012;42191-97. https://doi.org/10.1016/j.ijepes.2012.03.029Links ]

Schneider KP, Mather BA, Pal BC, Ten CW, Shirek GJ, Zhu H, Kersting W. Analytic Considerations and Design Basis for the IEEE Distribution Test Feeders. IEEE Transactions on Power Systems. 2018;3333181-3188. https://doi.org/10.1109/TPWRS.2017.2760011Links ]

Shen T, Li Y, Xiang J. A graph-based power flow method for balanced distribution systems. Energies. 2018;113511. https://doi.org/10.3390/en11030511Links ]

Singh B, Payasi RP, Shukla V. A taxonomical review on impact assessment of optimally placed DGs and FACTS controllers in power systems. Energy Reports. 2017;394-108. https://doi.org/10.1016/j.egyr.2017.07.001Links ]

Sirjani R, Mohamed A, Shareef H. Optimal placement and sizing of Static Var Compensators in power systems using Improved Harmony Search Algorithm. Przegląd Elektrotechniczny (Electrical Review). 2011;214-218. https://doi.org/10.1109/ICI.2011.71Links ]

Solera Losada NA, Villalba Jaramillo JP, Montoya Giraldo OD. Integración Óptima de Generadores Fotovoltaicos en Sistemas de Distribución DC a través de la Aplicación del Algoritmo de Optimización Aritmética Modificado. TecnoLógicas. 2022;2555. https://doi.org/10.22430/22565337.2418Links ]

Phase Balancing of Unbalanced Distribution Network through Hybrid Greedy-Fuzzy Algorithm. International Journal of Computer Applications. 34. https://www.researchgate.net/publication/235657105Links ]

Swapna M, Udaykumar RY. An Algorithm for Optimal Phase Balancing of Secondary Distribution Systems at Each Node. 2016 IEEE PES 13th International Conference on Transmission & Distribution Construction, Operation & Live-Line Maintenance (ESMO). 2016 1-5. https://doi.org/10.1109/TDCLLM.2016.8013241Links ]

Mansani S, Yaragatti RU. Backward Sweep Technique Based Phase Balancing Algorithm for Secondary Distribution System. 2017 IEEE PES Asia-Pacific Power and Energy Engineering Conference (APPEEC). 2017 1-6. https://doi.org/APPEEC.2017.8308968Links ]

Tuppadung Y, Kurutach W. The Modified Particle Swarm Optimization for Phase Balancing. TENCON 2006 - 2006 IEEE Region 10 Conference. 2006 1-4. https://doi.org/10.1109/TENCON.2006.344014Links ]

Zeng R, Abualigah L, Liu Q, Wang S. An Improved arithmetic optimization algorithm with forced switching mechanism for global optimization problems. Mathematical Biosciences and Engineering. 2022;191473-512. https://doi.org/10.3934/mbe.2022023Links ]

Zhu J, Bilbro G, Chow M-Y. Phase Balancing using Simulated Annealing. IEEE Transactions on Power Systems. 1999;1441508-1513. https://doi.org/10.1109/59.801943Links ]

Zhu J, Chow M-Y, Zhang F. Phase Balancing using Mixed-Integer Programming. IEEE Transactions on Power Systems. 1998;1341487-1492. https://doi.org/10.1109/59.736295Links ]

To cite this article: N. O. Tique Tique, J. C. Castillo Sáenz, y O. D. Montoya Giraldo, Application of the arithmetic optimization algorithm to the phase-balancing problem in three-phase asymmetric distribution networks,Tecnura, vol. 29, n.º 83, pp. 37-64, mar. 2025.. https://doi.org/10.14483/22487638.20582

Received: March 08, 2023; Accepted: October 11, 2024

*Corresponding author: notiquet@correo.udistrital.edu.co

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License