Relation between the Branch-to-Node Incidence and Triangular Matrices in Radial Distribution Networks
Relación entre la matriz de incidencia rama-nodo y la matriz triangular superior en redes de distribución radial
O. D. Montoya, A. Molina-Cabrera, and J. C. Hernandez, "A comparative study on power flow methods applied to AC distribution networks with single-phase representation," Electronics, vol. 10, no. 21, p. 2573, oct 2021.
G. Li, X.-P. Zhang, and X.-F. Wang, "Accelerated newton-raphson power flow," Eur. T. on Electr. Power, vol. 22, no. 4, pp. 504-517, mar 2011.
A. Marini, S. Mortazavi, L. Piegari, and M.-S. Ghazizadeh, "An efficient graph-based power flow algorithm for electrical distribution systems with a comprehensive modeling of distributed generations," Electr. Power Syst, vol. 170, pp. 229-243, may 2019.
T. Shen, Y. Li, and J. Xiang, "A graph-based power flow method for balanced distribution systems," Energies, vol. 11, no. 3, p. 511, feb 2018.
A. Garcés-Ruiz, "Power Flow in Unbalanced Three-Phase PowerDistribution Networks Using Matlab: Theory, analysis, and quasi-dynamic simulation," Ing, vol. 27, no. 3, p. e19252, aug 2022.
M. C. Herrera-Briñez, O. D. Montoya, L. Alvarado-Barrios, and H. R. Chamorro, "The Equivalence between Successive Approximations and Matricial Load Flow Formulations," Applied Sciences, vol. 11, no. 7, p. 2905, mar 2021.
S. Ouali and A. Cherkaoui, "An Improved Backward/Forward Sweep Power Flow Method Based on a New Network Information Organization for Radial Distribution Systems," Journal of Electrical and Computer Engineering, vol. 2020, pp. 1-11, jan 2020.
How to Cite
The power flow solution is a classical problem in electrical engineering which has been studied for more than 60 years 1. One of the most widely used methods corresponds to the Newton-Raphson approach, which is currently employed for analyzing power systems with meshed configurations and multiple generation sources, i.e., it is typically employed for power systems at high-voltage levels 2. In the case of medium- and low-voltage levels, power flow analysis is performed with graph-based approaches that consider the radial topology of the network in order to propose derivative-free solution methods 3), (4. Said graph-based approaches reduce the processing times required to obtain the power flow solution because they do not use recursive matrix inversions during the iterative solution process, unlike derivative-based power flow methods 5.
Two of the most widely known power flow solution methodologies for radial distribution networks in medium-voltage levels are the backward/forward power flow 4 and the triangular-based power flow 3. Both methods are based on the tree structure that represents radial distribution networks. The main characteristic of both methods is that they are formulated using the branch-to-node (Ad) and the triangular (T) matrices. However, both solution methods are mathematically different. The power flow method based on the triangular matrix is only applicable to strictly radial electrical distribution networks. Still, it is possible to find an analytical relation between both matrices, as demonstrated in this editorial note.
Theorem 1. The branch-to-node incidence matrix A and the upper-triangular matrix T are related through a negative inverse, including a transposed operation, i.e.,
Proof. To demonstrate Equation (1), let us recur to the distribution grid presented in Fig. 1. Note that Jb represents the current flowing through line b, and Ik corresponds to the net demanded current at node k.
It is worth mentioning that the branch and nodal currents can be related using a matricial operation, as presented in Equation (2):
where is the vector that contains all the branch currents, and contains all the demanded currents, respectively.
Now, if the Kirchhoff’s first law is applied to each node of the network, except to the substation node, i.e., node 1, as illustrated in Figure 1, the relation between nodal injected currents and branch currents is reached as presented in (3).
For more details regarding the construction of the branch-to-node incidence matrix, please refer to 6. Now, in order to demonstrate that Equation (1) is true, if (3) is substituted into (2), then (4) is obtained:
which clearly shows that
where l is the number of distribution lines in the distribution grid under analysis, i.e., for Figure 1 l = 6.
In addition, if the operation in (5) is made for matrices T and A⊤ d in (2) and (3), the following result is reached.
which confirms that
thus completing the proof.
Remark 1. The result in Equation (1) is particularly important for power flow studies since it demonstrates that, for strictly radial distribution grids, the branch-to-node incidence matrix Ad and the upper-triangular matrix T are related. This is, it will be possible to demonstrate that the upper-triangular power flow method is indeed a particular case of the classical and well-known backward/forward power flow method 7.
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License.
From the edition of the V23N3 of year 2018 forward, the Creative Commons License "Attribution-Non-Commercial - No Derivative Works " is changed to the following:
Attribution - Non-Commercial - Share the same: this license allows others to distribute, remix, retouch, and create from your work in a non-commercial way, as long as they give you credit and license their new creations under the same conditions.