DOI:
https://doi.org/10.14483/23448393.21419Publicado:
2023-10-17Número:
Vol. 28 Núm. 3 (2023): Septiembre-diciembreSección:
EditorialThree-Phase Power Flow Tool for Electric Distribution Grids: A Julia Implementation for Electrical Engineering Students
Palabras clave:
Electrical Engineering, students (en).Descargas
Referencias
C. Yang, Y. Sun, Y. Zou, F. Zheng, S. Liu, B. Zhao, M. Wu, and H. Cui, "Optimal power flow in distribution network: A review on problem formulation and optimization methods," Energies, vol. 16, no. 16, p. 5974, Aug. 2023. https://doi.org/10.3390/en16165974 DOI: https://doi.org/10.3390/en16165974
A. Garcés, "A linear three-phase load flow for power distribution systems," IEEE Trans. Power Syst., vol. 31, no. 1, pp. 827-828, Jan. 2016.
https://doi.org/10.1109/TPWRS.2015.2394296 DOI: https://doi.org/10.1109/TPWRS.2015.2394296
B. Cortés-Caicedo, L. S. Avellaneda-Gómez, O. D. Montoya, L. Alvarado-Barrios, and H. R. Chamorro, "Application of the vortex search algorithm to the phase-balancing problem in distribution systems," Energies, vol. 14, no. 5, p. 1282, Feb. 2021.
https://doi.org/10.3390/en14051282 DOI: https://doi.org/10.3390/en14051282
O. D. Montoya, J. S. Giraldo, L. F. Grisales-Noreña, H. R. Chamorro, and L. Alvarado-Barrios, "Accurate and efficient derivative-free three-phase power flow method for unbalanced distribution networks," Computation, vol. 9, no. 6, p. 61, May 2021. https://doi.org/10.3390/computation9060061 DOI: https://doi.org/10.3390/computation9060061
J. M. Perkel, "Julia: Come for the syntax, stay for the speed," Nature, vol. 572, no. 7767, pp. 141-142, Jul. 2019. https://doi.org/10.1038/d41586-019-02310-3 DOI: https://doi.org/10.1038/d41586-019-02310-3
Cómo citar
APA
ACM
ACS
ABNT
Chicago
Harvard
IEEE
MLA
Turabian
Vancouver
Descargar cita
Most of the concepts involved in Electrical Engineering have been developed from the 20th century to the present day. One of the most popular studies in this field is the power flow problem for electrical systems composed of multiple nonlinear loads (i.e., loads with constant power consumption) 1. This problem is formulated by applying Kirchhoff’s voltage law to each node in the system, which results in a set of nonlinear algebraic equations that consider steady-state conditions 2. In this vein, this editorial note presents a tutorial for Electrical Engineering students regarding the general implementation of the three-phase power Flow problem in distribution networks via the successive approximations method 3.
General three-phase power flow formulation
The general formula for the power flow solution in three-phase networks is obtained after applying Kirchhoff’s voltage law to each network node, which produces the set of constraints defined by
where Yds3φ ∈ C3(n−1)×3 and Ydd3φ ∈ C3(n−1)×3(n−1) are complex sub-matrices obtained for the nodal three-phase admittance matrix Ybus3φ (note that this matrix can be obtained using the node-to-branch incidence and the admittance primitive matrices 3); Id3φ ∈ C3(n−1)×1 is the complex vector of demanded currents, which is a nonlinear function of the demanded voltages Vd3φ ∈ C3(n−1)×1 and the constant complex power loads; and Vs3φ ∈ C3×1 is the complex power voltage output at the terminals of the substation. Note that C denotes the set of complex numbers.
Remark 1. The calculation of the demanded current vector Id3φ depends on the load connection type, i.e., Wye (Y ) or Delta (Δ).
To illustrate the general calculation of the three-phase demanded currents, consider the schematic load connections presented in Figs. 1a and 1b.
Figure 1: Constant power load: a) Y connection and b) Δ connection
For these load connections, when Tellegen’s second theorem is applied (i.e., the relation between voltages, currents, and powers), the following compacted three-phase current formulas are obtained:
where the matricesMand H take the following form
and (·) is the complex conjugate operation applied to the argument.
Iterative power flow solution
Considering the general power flow formula in (1) and the definition of the demanded current in 4, an iterative power flow formula can be obtained (4).
where
for all k in the set of the demanded buses.
Remark 2. Eqs. (5) and (6) are iteratively solved until the desired error of convergence is reached, i.e.,
with ε being the convergence’s tolerance.
Julia implementation
Julia is an efficient programming environment that combines the advantages of Python, Matlab, and R with the efficiency of well-known programming languages such as C++ or Fortran 5. To illustrate the general implementation of the three-phase power flow problem in Julia, consider a small distribution grid composed of seven buses and six lines, which is unbalanced and operates with a line-to-ground voltage of 23 kV at the terminals of the substation. All the information on this system is directly presented in the Julia scripts. Fig. 2 illustrates the first part of the three-phase power flow problem’s implementation in distribution networks.
Figure 2: Parametric information associated with the three-phase 7-bus network
Fig. 2 presents the parametric information of the network, i.e., the distribution line impedance (in this example, the three-phase impedance matrix is assumed to be balanced with values only in its diagonal) and the constant power loads. Note that, if node_data.Type [k,1] == 1, the load associated with the k-node is Y -connected; otherwise, it is Δ−connected.
The next step is calculating the nodal admittance matrix. This process is shown in Fig. 3.
Once the nodal admittance matrix has been constructed, the iterative power flow formula defined in (5), the current calculation in (6), and the stopping criterion in (7) are implemented in Julia using a function called Successive3f(node_data,Ybus3f,error,tmax,Sb) (Fig. 4). This implementation is intuitive, and the student/researcher can easily follow all the procedures for solving the three-phase power flow problem via the successive approximations method.
Finally, to obtain the general power flow solution for the 7-bus grid used in this example, the power flow function calculates all the phase voltages, which allows obtaining the total grid apparent power losses and reporting the voltage magnitudes and angles per phase. This process is depicted in Fig. 5. It is worth mentioning that, this figure, the apparent power losses are calculated using Eq. (8).
Figure 3: General construction of the three-phase nodal admittance matrix
where V corresponds to the voltage output obtained by calling the power flow function implemented in Fig. 4 as V3f = Successive3f(node_data,Ybus3f,error,tmax,Sb).
The power flow solution for this numerical example is plotted in Fig. 6.
The student/researcher interested in the Julia-based three-phase power flow solution presented in this editorial note should thoroughly examine all the codifications and follow each one of the presented plots (see Figs. 2 to 6) to reach the same numerical solution presented in Fig. 6.
Conclusion
This editorial note presented an intuitive algorithm based on the successive approximations power flow method with the aim of solving the three-phase power flow problem in electric distribution networks via the Julia programming environment (version 1.9.2). The main idea of this tutorial is to provide a new power flow tool for Electrical Engineering students and researchers, an easily implementable open-source algorithm for conducting studies in unbalanced distribution networks.
Figure 4: Iterative three-phase power flow solution
Figure 5: Calling the three-phase power flow function and reporting the results
Figure 6: Three-phase power flow solution in the studied 7-bus grid
Acknowledgements
Acknowledgements
The first author wants to thank the Master’s Program in Electrical Engineering of Universidad Tecnológica de Pereira and the Unbalanced Distribution Grids Analysis course, wherein this material has been developed to teach graduate and undergraduate Electrical Engineering students during the second semester of 2023.
Referencias
Licencia
![Creative Commons License](http://i.creativecommons.org/l/by-nc-sa/4.0/88x31.png)
Esta obra está bajo una licencia internacional Creative Commons Atribución-NoComercial-CompartirIgual 4.0.
A partir de la edición del V23N3 del año 2018 hacia adelante, se cambia la Licencia Creative Commons “Atribución—No Comercial – Sin Obra Derivada” a la siguiente:
Atribución - No Comercial – Compartir igual: esta licencia permite a otros distribuir, remezclar, retocar, y crear a partir de tu obra de modo no comercial, siempre y cuando te den crédito y licencien sus nuevas creaciones bajo las mismas condiciones.