SciELO - Scientific Electronic Library Online

 
vol.26 issue71Monitoring Algorithm for Electrical Variables Implementing Blockchain and PythonUsing Environmental Education and Healthy Housing as Strategies for Household-Level Prevention of COVID-19 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.26 no.71 Bogotá Mar./Mar. 2022  Epub Apr 05, 2022

https://doi.org/10.14483/22487638.18635 

Investigación

Branch Optimal Power Flow Model for DC Networks with Radial Structure: A Conic Relaxation

Flujo de potencia óptimo de ramas para redes DC con estructura radial: una relajación cónica

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

Andrés Arias-Londoño2 
http://orcid.org/0000-0002-6210-3449

Alexander Molina-Cabrera3 
http://orcid.org/0000-0003-2448-2174

1 PhD in Engineering, Master’s degree in Electrical Engineering, electrical engineer. Assistant professor of the Department of Engineering, Universidad Distrital Francisco José de Caldas. Bogotá D. C., Colombia. Email: odmontoyag@udistrital.edu.co

2 PhD in Engineering, Master’s degree in Electrical Engineering, electrical engineer. Assistant professor of the Department of Engineering, Institución Universitaria Pascual Bravo, Medellín, Colombia. Email: andres.arias366@pascualbravo.edu.co

3 PhD in Engineering, Master’s degree in Electrical Engineering, electrical engineer. Associate professor of the Department of Engineering, Universidad Tecnológica de Pereira, Pereira, Colombia. Email: almo@utp.edu.co


Abstract

Objective:

This work involves a convex-based mathematical reformulation of the optimal power flow problem in DC networks. The objective of the proposed optimization model corresponds to the minimization of the power losses throughout all the network branches considering a convex conic model that ensures finding the global optimum solution.

Methodology:

This work is split into three stages. The first stage presents the mathematical model of optimal power flow for DC networks and all the geometric features that make it non-convex. The second stage presents the convex reformulation from a second-order conic relaxation. The third stage shows the main characteristics of the DC system under study, as well as the optimal solution of the power flow problem and its comparisons with some methods reported in the specialized literature.

Results:

The numerical validations demonstrate that the proposed convex optimal power flow model obtains the same solution as the exact model of the problem with an efficiency of 100%, which is in contrast with the variability of the results that are presented by the metaheuristic techniques reported as comparison methodologies.

Conclusions:

The proposed second-order conic relaxation ensured the convexity of the solution space and, therefore, the finding of the optimal solution at each execution, in addition to demonstrating that, for optimal power flow problems in DC networks, the numerical performance is better than most of the comparative metaheuristic methods and that the solution provided by the proposed relaxation is equivalent to that provided by the exact model.

Keywords: direct current networks; second-order conic relaxation; non-linear programming model; convex optimization

Resumen

Objetivo:

Este trabajo plantea una reformulación matemática de naturaleza convexa del problema de flujo de potencia óptimo en redes de corriente continua (DC). El objetivo del modelo de optimización propuesto corresponde a la minimización de las pérdidas de potencia en todas las ramas de la red considerando un modelo cónico convexo que garantice el hallazgo de la solución óptima global.

Metodología:

Este trabajo está dividido en tres etapas: la primera presenta el modelo matemático de flujo de potencia óptimo para redes DC y todas las características geométricas que lo hacen no convexo; la segunda presenta la reformulación convexa a partir de una relajación cónica de segundo orden; la tercera etapa presenta las principales características del sistema DC bajo estudio, además de la solución óptima del problema de flujo de potencia y sus comparaciones con algunos métodos reportados en la literatura especializada.

Resultados:

Las validaciones numéricas demuestran que el modelo de flujo de potencia óptimo convexo propuesto encuentra la misma solución el modelo exacto del problema y tiene una eficiencia del 100%, lo cual contrasta con la variabilidad de resultados que presentan las técnicas metaheurísticas reportadas como métodos de comparación.

Conclusiones:

La relajación cónica de segundo orden propuesta garantizó la convexidad del espacio de soluciones, y, por tanto, el hallazgo de la solución óptima en cada ejecución. Además, demostró que, para problemas de flujo de potencia óptimo en redes DC, tiene el mejor desempeño numérico que la mayoría de los métodos metaheurísticos comparativos; y que la solución provista por la relajación propuesta es equivalente a la proveída por el modelo exacto.

Palabras clave: redes de corriente continua; relajación cónica de segundo orden; modelo de programación no lineal; optimización convexa

INTRODUCTION

For decades, distribution networks have been traditionally designed under the alternating current (AC) paradigm (Starke et al., 2008; Gelani et al., 2019;Garcés et al., 2014). However, in recent years, the accelerated advance in power electronics, renewable generation, and energy storage systems has steered these distribution networks towards the direct current (DC) paradigm (Montoya et al., 2020a; Serra & De-Angelo, 2021). This shift is mainly supported by the following advantages: (i) multiple generation and energy storage technologies work directly with DC technologies (photovoltaic sources, batteries, supercapacitors, and superconducting coils), which implies that the number of power electronic interfaces can be reduced when these are connected to DC networks instead of AC networks (Lotfi & Khodaei, 2017; Hidalgo-Mora et al., 2014); (ii) DC distribution feeders are more efficient in terms of voltage profiles and power losses since reactive power and frequency are non-existing concepts within these grids (Garcés, 2018; Gil-González, et al., 2020; Grisales-Noreña, et al. 2020); and (iii) DC networks are easily controllable since the main goal is to control the voltage profile in all the buses of the grid without the need for frequency synchronization (Parhizi, et al., 2015).

To analyze DC distribution networks, the literature has proposed multiple approaches regarding power flow and optimal power flow analysis. Some of them are nonlinear analysis of DC grids with constant power loads (Simpson-Porco et al., 2015) and convergence analysis of the Newton-Raphson and Gauss-Seidel methods (Garcés, 2017, 2018). The most common approaches in the field are based on second-order, semidefinite programming and interior point methods under nodal representation of the DC network (Bahrami, et al., 2017; Li, et al., 2018; Montoya & Gil-González, 2021) and metaheuristic approaches based on particle swarm and genetic algorithm methods (Grisales-Noreña, et al., 2019). It is worth mentioning that, in the case of optimal power flow analysis, not all methods guarantee optimum global finding, as is the case of nonlinear programming models and metaheuristics (except semidefinite and second-order cone programming). Therefore, alternative methods are required (convex optimization approaches) which find the optimum global by transforming the non-convex solution space of the original model into a convex equivalent model via conic representation. For this reason, in this research, we proposed an alternative branch optimal power flow formulation for DC networks that has previously proposed for AC networks in (Farivar & Low, 2012), with the advantage that, for a DC grid with n nodes and l lines, the number of required variables is 2(n+l), whereas existing convex methods require (n 2 +n) variables (Li et al., 2018). This contribution may significantly reduce the efforts in terms of the processing time.

The remainder of this document is organized as follows: the Branch power flow formulation section presents the exact formulation of the branch optimal power flow model for DC distribution networks; Conic relaxation describes the second-order cone programming model with the relaxed branch optimal power flow; Computational validation presents the test system characteristics and the computational validation of the proposed convex model regarding metaheuristics and exact methods; and, finally, Conclusions presents the main concluding remarks of this research, as well as guidelines for future work.

BRANCH POWER FLOW FORMULATION

The power flow problem and its optimization variants (i.e., optimal power flow models) are classically formulated via nodal representation (Garcés, 2017). In these formulations, it is intended not to use current and branch power variables, which implies that these must be calculated once all the voltage variables are known. However, the power flow problem can be reformulated using branch and nodal variables at the same time by means of the branch-based power flow formulation (Farivar & Low, 2012). Consider the second Tellegen theorem applied to each node of the network except the slack node, which provides the following set of expressions:

where 𝑝𝑗𝑘 (𝑝𝑘𝑚) is the power flow through the line that connects nodes j and k (k and m, respectively), 𝑖𝑗𝑘 is the current flowing through the line that connects nodes j and k, 𝑝𝑘 is the power injection at node k defined as demand minus generation, and 𝑅𝑗𝑘 represents the resistive effect in line j-k. Note that Ε is the set that contains all the branches of the network. Ohm's law applied at each branch results in

with 𝑣𝑗 and 𝑣𝑘 being the voltage values at nodes j and k, respectively.

The power in the DC networks for each branch is defined as

To formulate the optimal power flow model, the scientific literature widely uses the minimization of power losses in all branches, which generates the following nonlinear programming optimization model:

Objective function:

Set of constraints:

Remark 1. The optimization model defined in (4) is known in the literate as the optimal branch power flow model (Farivar & Low, 2012), which is nonlinear and non-convex due to the square of the currents in (4b) and the product between voltages and currents in (4d).

Note that an equivalent optimization model can be obtained from (4) by making some algebraic manipulations. To this effect, let us pre-multiply (4c) by 𝑣 j , which produces

Now, if we raise (4c) to the square in both sides, we have:

Note that, if Equation (6b) is substituted into Equation (5b) and some algebraic manipulations are made, it yields

To obtain an equivalent model, let us define two auxiliary variables, 𝑙𝑗𝑘=𝑖𝑗𝑘 2 and 𝑢𝑗=𝑣𝑗 2.. With these new variables, the optimal branch power flow model (4) can be rewritten as follows:

Objective function:

Set of constraints:

Remark 2. The optimal branch power flow model defined in (8) is still non-convex due to the presence of the product between auxiliary variables 𝑢𝑗 and 𝑙𝑗𝑘 in equality constraint (8d). However, this complication can be dealt with by using a conic relaxation through its hyperbolic representation, as presented in the next section.

CONIC RELAXATION

Conic optimization is a subfield of convex optimization that allows relaxing some class of optimization problems using second-order cone constraints (Farivar & Low, 2012), which is especially attractive when faced with products among continuous variables. Even if the cones are nonlinear inequality constraints, they are convex since these constraints are in the interior space of the cone (Benson & Sauglam, 2013). Then, the conic relaxation of the branch optimal power flow model consists of rewriting (8) with its hyperbolic equivalent (Farivar & Low, 2012). To do so, let us consider the following relation:

Now, if we substitute (9) into (8d), then, we have

Observe that (10) can be rewritten using the Euclidean norm as follows:

Remark 3. The structure of Equation (11) is still non-convex due to the equality sign. However, as suggested byFarivar and Low (2012), this can be relaxed with a lower-equal symbol, which allows (11) to become a second order conic constraint.

With the conic relaxation of (12), the branch optimal power flow model defined in (8) takes the convex structure of (13).

Objective function:

Set of constraints:

Remark 4. The most important characteristic of the second-order cone programming branch optimal power flow model defined in (13) is that it guarantees finding the global optimum with interior point methods under well-defined voltage and demand conditions (Li et al., 2018). In addition, this convex relaxation creates a lower number of variables than the classical convex models reported in the literature, such as semidefinite programming (Bahrami et al., 2017) and SOCP models that only work with voltage variables (Li et al., 2018).

It is worth mentioning that this research does not present the methodology for a solution, since the main contribution of this research is indeed the convexification of the exact nonlinear programming model (8) using second-order cone programming as presented in (13). This implies that, due to the convexity of the solution space and of the objective function with any convex optimizer (e.g., CVX for MATLAB or CVXPY for Python), it is possible to reach the global optimum of the studied problem with a 100% of repeatability properties.

COMPUTATIONAL VALIDATION

To demonstrate the effectiveness and robustness of the proposed SOCP model presented in (13), we considered a large-scale medium-voltage distribution network composed of 69 nodes and 68 lines (radial configuration) with the configuration depicted in Figure 1.

Source: Montoya et al. (2019a)

Figure 1 Schematic connection among nodes for the 69-node test feeder 

The 69-node test feeder operates at 12.660 V. All the numerical information of this test feeder can be consulted in the work by Montoya et al. (2019a). For this test feeder, we considered the likely dispatch of three distributed generators located at nodes 21, 61, and 64, with maximum generation capabilities of 12 pu each (Montoya & Gil-González, 2021). In addition, to compare the performance of the proposed SOCP optimization model, we considered different metaheuristic approaches, typically reported in the literature for optimal power flow analysis in DC networks, namely black-hole optimization (BHO) (Velasquez et al., 2019), elephant swarm water search algorithms (ESWSA) (Montoya et al., 2019b), continuous genetic algorithms (GGA) (Montoya & Gil-González, 2021), sine-cosine algorithms (Giraldo, et al., 2019), vortex search algorithms (Montoya, et al., 2020b), and an exact approach based on interior points available in the GAMS software (Montoya & Gil-González, 2021).

Table 1 shows the numerical performance of all the comparative methodologies for optimal power flow analysis in DC grids. It is worth mentioning that all the metaheuristic methods were evaluated 100 times to determine the best possible solution reached by each one of them.

Table 1 Numerical results in the OPF problem with different comparative approaches 

Method Generation (kW) Power losses (kW)
BHO [460,21, 1170,28, 639,88] 5,025771
ESWSA [495,80, 1049,55, 699,62] 5,005033
CGA [401,31, 1191,55, 584,05] 4,982861
SCA [499,86, 1199,90, 564,26] 4,557555
VSA [455,47, 1200,00, 584,85] 4,454562
GAMS [453,21, 1200,00, 585,16] 4,454342
SOCP [453,22, 1200,00, 585,17] 4,454342

Source: Authors

The numerical results in Table 1 reveal that:

  • The best metaheuristic approach for dealing with optimal power problems in DC distribution networks is to the VSA approach. Its variations regarding the optimal solution are in the order of milliwatts, which constitutes an exact methodology from the numerical point of view.

  • The worst behavior regarding metaheuristics for optimal power flow solution is the result reported by the BHO. This is due to the fact that it is a simplification of the classical particle swarm optimizer, which is easily stuck in local optimal solutions (Gupta et al., 2016).

  • The proposed SOCP model for optimal power flow analysis in DC distribution networks allows reaching the global optimum solution for this problem since its solution matches with the interior point method available in GAMS. This is to be expected, as the solution of the SOCP models has been elaborated with interior point methods in the scientific literature (Benson & Sauglam, 2013).

Figure 2 presents the percentage of power losses minimization, considering that the base case without distributed generation has an initial power loss of about 153,847557 kW.

Source: Authors

Figure 2 Improvement of the power losses reached by the proposed SOCP reformulation and the comparative methods 

Note that the results in Figure 1 confirm that VSA, GAMS, and the proposed SOCP approach are the best methodologies regarding power losses minimization via optimal power flow analysis with distributed generation. Also, the most important implication of these results is that, in the case of DC networks, the total power losses can be minimized by up to 90%, which is not the case, with most of the results reported for the 69-nodes test feeder with distributed generation around 60% (Kaur, 2014). This difference between both technologies can be attributed the reactive effects inherent to AC networks, which do not arise in the case of DC distribution.

CONCLUSION

A second-order cone programming model for optimal power flow analysis in DC distribution networks that combines nodal and branch variables has been presented in this paper. The exact nonlinear programming model was convexified via the relaxation of the power sent from nodej to node k, i.e., 𝑝𝑗𝑘 = 𝑣𝑗 𝑖 𝑗𝑘, with its conic equivalent of hyperbolic representation. The numerical results in the 69-node test feeder demonstrate that the proposed SOCP model allows reaching the global optimal solution for the optimal power flow problem in DC distribution networks with distributed generators, given that its results are better than metaheuristic methods such as BHO, ESWSA, CGA, and SCA. In addition, numerically speaking, the only metaheuristic method that can reach a near-optimal solution is the VSA approach, which was compared with the proposed SOPC model and the interior point methods available in the GAMS software.

References

Bahrami, S., Therrien, F., Wong, V. W., & Jatskevich, J. (2017). Semidefinite Relaxation of Optimal Power Flow for AC-DC Grids. IEEE Transactions on Power Systems, 32(1), 289-304. https://doi.org/10.1109/TPWRS.2016.2543726Links ]

Benson, H. Y., & Saglam, U. (2013). Mixed-Integer Second-Order Cone Programming: A Survey. INFORMS TutORials in Operations Research, 13-36. https://doi.org/10.1287/educ.2013.0115Links ]

Farivar, M., & Low, S. H. (2012, December 10-13). Branch flow model: Relaxations and convexification [Conference presentation]. IEEE 51st IEEE Conference on Decision and Control (CDC), Maui, HI, USA. https://doi.org/10.1109/CDC.2012.6425870Links ]

Gupta, H., Gupta, A., Gupta, S. K., Nayak P., & Shrivastava, T. (2016). How effective is Black Hole Algorithm? [Conference presentation] 2nd International Conference on Contemporary Computing and Informatics (IC3I), Greater Noida, India. https://doi.org/10.1109/IC3I.2016.7918011Links ]

Gelani, H., Dastgeer F., Siraj, K., Nasir, M., Niazi, K., & Yang Y. (2019). Efficiency Comparison of AC and DC Distribution Networks for Modern Residential Localities. Applied Sciences, 9(3) 582. https://doi.org/10.3390/app9030582Links ]

Garcés, A., Correa, C. A., & Sánchez, A. (2014). Operación económica de dispositivos almacenadores de energía para disminución de pérdidas. Revista Tecnura, 18(42), 53-64. https://doi.org/10.14483/udistrital.jour.tecnura.2014.4.a04Links ]

Garcés, A. (2018). On the Convergence of Newtons Method in Power Flow Studies for DC Microgrids. IEEE Transactions on Power Systems, 33(5) 5770-5777. https://doi.org/10.1109/TPWRS.2018.2820430Links ]

Garcés, A. (2017). Uniqueness of the power flow solutions in low voltage direct current grids. Electric Power Systems Research, 151, 149-153. https://doi.org/10.1016/j.epsr.2017.05.031Links ]

Gil-González, W., Montoya, O. D., Grisales-Noreña, L. F., Ramírez-Vanegas, C. A., & Molina-Cabrera, A. (2020). Hybrid Optimization Strategy for Optimal Location and Sizing of DG in Distribution Networks. Revista Tecnura, 24(66), 47-61. https://doi.org/10.14483/22487638.16606Links ]

Giraldo, J. A., Montoya, O. D., Grisales-Noreña, L. F., Gil-González W., & Holguín, M. (2019). Optimal power flow solution in direct current grids using Sine-Cosine algorithm. J. Journal of Physics: Conference Series, 1403, 012009. https://doi.org/10.1088/1742-6596/1403/1/012009Links ]

Grisales-Noreña, L. F., Garzón-Rivera, O. D., Montoya, O. D., & Ramos-Paja, C. A. (2019). Hybrid Metaheuristic Optimization Methods for Optimal Location and Sizing DGs in DC Networks. In J. Figueroa-García, M. Duarte-González, S. Jaramillo-Isaza, A. Orjuela-Cañón, & Y. Díaz-Gutiérrez (Eds.) Communications in Computer and Information Science (vol. 1052, pp. 214-225). Springer. https://doi.org/10.1007/978-3-030-31019-6_19Links ]

Grisales-Noreña, L., Garzón Rivera, O., Ocampo Toro, J., Ramos Paja, C., & Rodríguez Cabal, M. (2020). Metaheuristic Optimization Methods for Optimal Power Flow Analysis in DC Distribution Networks. Transactions on Energy Systems and Engineering Applications, 1(1), 13-31. https://doi.org/10.32397/tesea.vol1.n1.2Links ]

Hidalgo-Mora, C. A., Díaz-Aldana, N. L., & Trujillo-Rodríguez, C. L. (2014). Diseño e implementación de controladores lineales para regulación del bus DC en convertidores VSC para HVDC. Revista Tecnura, 18(40), 48-61. https://doi.org/10.14483/udistrital.jour.tecnura.2014.2.a04Links ]

Kaur, S., Kumbhar, G., & Sharma, J. (2014). A MINLP technique for optimal placement of multiple DG units in distribution systems. International Journal of Electrical Power & Energy Systems, 63, 609-617. https://doi.org/10.1016/j.ijepes.2014.06.023Links ]

Li, J., Liu, F., Wang, Z., Low, S. H., & Mei, S. (2018) Optimal Power Flow in Stand-Alone DC Microgrids. IEEE Transactions on Power Systems, 33(5), 5496-5506. https://doi.org/10.1109/TPWRS.2018.2801280Links ]

Lotfi, H., & Khodaei, A. (2017). AC Versus DC Microgrid Planning. IEEE Transactions on Smart Grid, 8(1), 296-304. https://doi.org/10.1109/TSG.2015.2457910Links ]

Montoya, O. D., & Gil-González, W. (2021). A MIQP model for optimal location and sizing of dispatchable DGs in DC networks. Energy Systems, 12, 181-202. https://doi.org/10.1007/s12667-020-00403-xLinks ]

Montoya, O. D., Gil-González, W., & Garcés, A. (2019a). Power flow approximation for DC networks with constant power loads via logarithmic transform of voltage magnitudes. Electric Power Systems Research, 175, 105887. https://doi.org/10.1016/j.epsr.2019.105887Links ]

Montoya, O. D., Gil-González, W., & Holguín, M. (2019b). Optimal power flow studies in direct current grids: An application of the bio-inspired elephant swarm water search algorithm. Journal of Physics: Conference Series, 1403, 012010. https://doi.org/10.1088/1742-6596/1403/1/012010Links ]

Montoya, O. D., Serra, F. M., & de Angelo, C. H. (2020a). On the Efficiency in Electrical Networks with AC and DC Operation Technologies: A Comparative Study at the Distribution Stage. Electronics, 9(9), 1352. https://doi.org/10.3390/electronics9091352Links ]

Montoya, O. D., Gil-González, W., & Grisales-Noreña, L. F. (2020b). Vortex Search Algorithm for Optimal Power Flow Analysis in DC Resistive Networks with CPLs. IEEE Transactions on Circuits and Systems II, 67(8), 1439-1443. https://doi.org/10.1109/TCSII.2019.2938530Links ]

Parhizi, S., Lotfi, H., Khodaei, A., & Bahramirad, S. (2015). State of the Art in Research on Microgrids: A Review. IEEE Access, 3, 890-925. https://doi.org/10.1109/ACCESS.2015.2443119Links ]

Serra, F. M., & de Angelo, C. H. (2021). Control of a battery charger for electric vehicles with unity power factor. Transactions on Energy Systems and Engineering Applications, 2(1), 32-44. https://doi.org/10.32397/tesea.vol2.n1.3Links ]

Simpson-Porco, J. W., Dorfler, F., & Bullo, F. (2015). On Resistive Networks of Constant-Power Devices. IEEE Transactions on Circuits and Systems II, 62(8), 811-815. https://doi.org/10.1109/TCSII.2015.2433537Links ]

Starke, M., Tolbert, L. M., & Ozpineci, B. (2008, April 21-24). AC vs. DC distribution: A loss comparison [Conference presentation]. IEEE/PES Transmission and Distribution Conference and Exposition, Chicago, IL, USA. https://doi.org/10.1109/TDC.2008.4517256Links ]

Velásquez, O. S., Montoya, O. D., Garrido, V. M., & Grisales-Noreña, L. F. (2019). Optimal Power Flow in Direct-Current Power Grids via Black Hole Optimization. Advances in Electrical and Electronic Engineering, 17(1), 3069. https://doi.org/10.15598/aeee.v17i1.3069 Links ]

Cómo citar: Montoya., O.D. Arias-Londoño., A. y Molina-Cabrera., A. (2022). Branch Optimal Power Flow Model for DC Networks with Radial Structure: A Conic Relaxation. Tecnura, 26(71). https://doi.org/10.14483/22487638.18635

Future work As future work, it will be possible to embed the proposed SOCP model for the optimal power flow method in a master-slave optimization algorithm guided by a discrete metaheuristic method (master) to determine the optimal location and sizing of the distributed generators in DC distribution networks, where the SOCP (slave) proposed model is entrusted with determining the optimal sizes of these distributed generators for each possible location provided by the master algorithm.

Received: October 05, 2021; Accepted: October 20, 2021

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