Mathematical Explanations of a Paradox Observed in a HVAC (High Voltage Alternating Current) Untransposed Overhead Line

: The constructive asymmetry of the untransposed overhead lines of a high voltage alternating current is the cause of a great number of difﬁculties in their operation and modeling. In order to model the operating regimes of such lines, the symmetrical component method, based on constructive symmetry and thus the symmetry of the equivalent phase parameters, is inappropriate, which is why many research papers have been dedicated to either setting up improved modeling methods or to returning to phase coordinate modeling. This paper intends to justify a paradox found on some untransposed overhead lines of a high voltage alternating current during the no-load operating conditions by performing phase coordinate modeling. In such a situation, the transmission or distribution operators measured signiﬁcant negative values for the active powers on one or two phases at the beginning of the lines. Considering the case of a real untransposed overhead line operating under no-load conditions, the paper starts from presenting the recorded electrical values. Then, the paper moves on to outlining the Carson’s simpliﬁed computing relations for calculating the series and shunt primitive equivalent parameters and Kron’s transformation relationships for calculating the phase equivalent parameters. After applying them to the real line, the calculation of the power ﬂow for the no-load operating conditions, which is applied to an equivalent scheme of the line consisting of nine identical octopoles, is performed. Both the untransposed line and its transposed variant are studied here. The values of the electrical amounts obtained by the calculation for the untransposed line are basically similar to those obtained by measuring on the real line, which gives a mathematical conﬁrmation of the so-called paradox. Its occurrence represents the effect of the asymmetry of the equivalent phase capacities, which causes a redistribution of the active powers between the phases of the network to which the overhead line operating in no-load conditions is connected.


Introduction
The symmetrical components method of the positive, negative, and zero sequence is the most widely used method for analyzing the normal operating conditions of the electrical power systems (EPS) [1]. It is, however, an approximate method based on a simplifying hypothesis with which the mathematical model is constructed, considering the symmetry of the equivalent, own, or mutual, longitudinal or transversal phase parameters. As it is known, the symmetrical components method uses three symmetric three-phase networks characterized by fictitious electrical amounts (impedances, unbalance is mainly caused by the existence of shunt currents corresponding to the capacitive couplings between phase and ground conductors respectively between phase conductors and grounding wires. It was also found that the magnetic couplings between the phases (which correspond to the mutual inductive impedances) cause relatively low unbalances.
A large portion of the existing studies have focused on the study of the effects of impedance asymmetries in three-phase systems, a particular form of the symmetrical component method. This method gives up the neglect of phase impedance asymmetry, so that the sequence impedance matrix has non-zero, non-diagonal elements. This means that the equivalent positive, negative, and zero sequence circuits are no longer independent, and they are coupled by impedances that cannot be neglected [15,16]. In this context, Bellan et al. [17] showed that in the case of the untransposed HVACTOHL, in the negative sequence equivalent circuit, there is a voltage known as the voltage asymmetry (unbalance) emission. Bellan et al. also presented a methodology for constructing an equivalent two-port network by which the equivalent positive and negative sequence circuits are coupled. This method is based on the equivalent scheme of an ideal transformer having a complex transformation ratio. In the case of untransposed HVACTOHL asymmetries having common paths, the application of the symmetric coordinate method results in equivalent impedances of coupling between the equivalent positive sequence scheme and the zero sequence scheme [18]. In this study, Bellan [18] presented an equivalent two-port network by which the coupling between the two circuits is modeled, namely the influence of the positive sequence scheme on the zero sequence scheme. This influence is modeled by a voltage source controlled by the positive sequence current, which is dependent on the asymmetry level of the phase parameters. In the same context, Youssef et al. [19] presented an original method of modeling the three-phase or six-phase circuits belonging to the real networks of alternative currents containing untransposed TOHL, which involves determining the coupling impedances between the fictitious sequences' networks. This method is applicable to the analysis of any type of faults, including those that are simultaneous and indifferent to the number of phases involved. The non-zero, non-diagonal elements of the sequence impedance matrix were also used in [20], this time to analyze the influence of phase transposition, in particular the influence of TOHL lengths, on the voltage asymmetry in the buses of the network.
In the same context, regarding the concern for the correct modeling of the electric line behavior, Wang et al. [21] presented the results of research performed to determine the influence of the conductor temperature on the equivalent parameters of a HVACTOHL and thus on the power flow in the system to which the line belongs.
The abandonment of symmetrical component modeling and the use of phase coordinates in the study of overhead and underground electric lines is also required when dealing with steady state or fault regimes of passive conductors (earth wires for overhead lines, metallic screens and armors for cables, and enclosures for gas insulated lines). For this purpose, the multiconductor cell analysis (MCA) method is used, which allows for a more accurate modeling of overhead or underground electric lines. An example of an efficient use of the MCA method is the direct calculation, in steady state and faulty regimes, of the current induced in the earth wires and ground returns alongside single or double circuit overhead lines (OHLs) [22].
The present article joins previous studies dedicated by the same authors to the numerical evaluation of the effects of impedance and admittance asymmetries on OHL operating regimes [23][24][25]. The main contributions of the article consist of the phase components modeling of a real 110 kV untransposed OHL with the goal of validating a "paradox" observed by the specialists of a distribution system operator during the no-load operating condition. For this, the equivalent longitudinal and shunt phase parameters were first calculated, and then, the calculation of the power flow for the studied regime was performed using an equivalent scheme of nine equivalent octopoles. The method was applied to both the untransposed line and the transposed version. The values obtained by calculation are similar to those obtained by measurement, which confirms the phenomenon produced and validates the correctness of the mathematical model.

Observed Paradox
Specialists of a Romanian distribution system operator have observed on a few occasions a phenomenon that they could not explain, which was produced on HVAC overhead lines of 110 kV voltage, operating in no-load conditions.
In such a regime, the line takes capacitive reactive energy from the network to which it is connected. This is equivalent to the delivery of inductive reactive energy in the network, which is in line with the following expectations: currents with a pronounced capacitive character are flowing on the line, determined by the natural capacities of the line. These currents have relatively low rms (root mean square) values and produce reactive energy losses on inductive longitudinal equivalent impedances, which are, however, much lower than the reactive energy absorbed by the natural capacities of the line. So, the reactive energy absorbed at the beginning of the line is capacitive.
The paradox refers to active energy. For many transmission lines operating in no-load conditions, the active energy at the beginning of the line is negative, which means that the line injects active energy into the network to which it is connected, behaving as an active power source. This phenomenon is contrary to expectations, because the line operating in no-load conditions is an active energy receiver. In this regime, the line should consume active energy, mainly due to the thermal effect of the current (Joule-Lenz), and to a very small extent due to the imperfection of the insulation and the corona effect. Measuring the active powers on the phases, it was found that the sign of the active energy comes from the sign of the active powers: in some cases on one phase and in other cases on two phases, the active power has the negative sign, and the sum of the three active phases is also negative.
For example, take the case of 110 kV OHL, located in a network area to which the transformer substations are connected for the supply of electric railway traction. The simplified single-wire electrical diagram of this area is shown in Figure 1.

Observed Paradox
Specialists of a Romanian distribution system operator have observed on a few occasions a phenomenon that they could not explain, which was produced on HVAC overhead lines of 110 kV voltage, operating in no-load conditions.
In such a regime, the line takes capacitive reactive energy from the network to which it is connected. This is equivalent to the delivery of inductive reactive energy in the network, which is in line with the following expectations: currents with a pronounced capacitive character are flowing on the line, determined by the natural capacities of the line. These currents have relatively low rms (root mean square) values and produce reactive energy losses on inductive longitudinal equivalent impedances, which are, however, much lower than the reactive energy absorbed by the natural capacities of the line. So, the reactive energy absorbed at the beginning of the line is capacitive.
The paradox refers to active energy. For many transmission lines operating in no-load conditions, the active energy at the beginning of the line is negative, which means that the line injects active energy into the network to which it is connected, behaving as an active power source. This phenomenon is contrary to expectations, because the line operating in no-load conditions is an active energy receiver. In this regime, the line should consume active energy, mainly due to the thermal effect of the current (Joule-Lenz), and to a very small extent due to the imperfection of the insulation and the corona effect. Measuring the active powers on the phases, it was found that the sign of the active energy comes from the sign of the active powers: in some cases on one phase and in other cases on two phases, the active power has the negative sign, and the sum of the three active phases is also negative.
For example, take the case of 110 kV OHL, located in a network area to which the transformer substations are connected for the supply of electric railway traction. The simplified single-wire electrical diagram of this area is shown in Figure 1.   I 1av = 7.132 A; I 2av = 7.653 A; I 3av = 7.574 A When analyzing the recorded values of the active powers on the three phases, the paradox described above can be observed. On two of the phases (1 and 2), the active powers have the negative sign, so they are delivered from the line to the network to which it is connected, and on the third phase (phase 3), the active power has the positive sign, being consumed (taken) from the network. The average values over the time interval of the active power's simultaneous measurement per phase are as follows: P 1av = −39.996 kW, P 2av = −1.190 kW, P 3av = +34.082 kW. Thus, the average value of the sum of the active powers over the same time period is as follows: P tot av = −2.368 kW. Regarding the reactive powers per phase, there is a variation in the absolute values higher than the active power. The average values over the time interval of the simultaneous measurement of the reactive powers for each phases were as follows: Q 1av = −483.905 kvar, Q 2av = −524.903 kvar, Q 3av = −514.780 kvar. The average value of the sum of the reactive powers over the same time period is as follows: Q tot av = −1523.588 kvar. When analyzing the recorded values of the active powers on the three phases, the paradox described above can be observed. On two of the phases (1 and 2), the active powers have the negative sign, so they are delivered from the line to the network to which it is connected, and on the third phase (phase 3), the active power has the positive sign, being consumed (taken) from the network.  The fact that the sum of the measured values of the active power on the phases is negative has a minor effect on the line operation and on the measurement and protection systems. However, if the section concerned is at the edge of the two operators' power management areas, the results may become important. Such a situation can be accepted as an effect of errors in the measurement systems, namely current transformer (CT) errors. It is known that at very low currents in relation to the rated current, CTs have very high relative errors for both the rms values and the angles of the secondary currents. For the studied OHL, the rms values of the rated currents are 800/1 A, so in the no-load operating conditions, the primary currents are less than 1% of the primary CT rated current.
However, it is difficult to accept that the active powers on the phases are positive, and the measuring systems indicate negative values. The causes of this phenomenon must be searched elsewhere, more precisely in the asymmetry of the line equivalent phase parameters.

Calculation of the Equivalent Phase Parameters for the Considered Overhead Line
The line considered only has one three-phase circuit, but it is installed on towers for doublecircuit lines. The second circuit, which was to be installed with the development of the network in that area, has not yet been built.
As mentioned above, the rated voltage of the line is 110 kV, and its length is 47.4 km. It was built with 86 metal towers, most of them having the geometrical characteristics shown in Figure 3. The three phase conductors (indices 1, 2, 3 corresponding to the phases a, b, c) and the protective conductor (index 4, corresponding to conductor n, grounded-g) are aluminum conductor steelreinforced cable (ACSR)-types and have rated sections of 185/32 mm 2 and 95/55 mm 2 , respectively. The insulators of the line are ball-and-socket-type glass strings, with a length of 1.8 m. The fact that the sum of the measured values of the active power on the phases is negative has a minor effect on the line operation and on the measurement and protection systems. However, if the section concerned is at the edge of the two operators' power management areas, the results may become important. Such a situation can be accepted as an effect of errors in the measurement systems, namely current transformer (CT) errors. It is known that at very low currents in relation to the rated current, CTs have very high relative errors for both the rms values and the angles of the secondary currents. For the studied OHL, the rms values of the rated currents are 800/1 A, so in the no-load operating conditions, the primary currents are less than 1% of the primary CT rated current.
However, it is difficult to accept that the active powers on the phases are positive, and the measuring systems indicate negative values. The causes of this phenomenon must be searched elsewhere, more precisely in the asymmetry of the line equivalent phase parameters.

Calculation of the Equivalent Phase Parameters for the Considered Overhead Line
The line considered only has one three-phase circuit, but it is installed on towers for double-circuit lines. The second circuit, which was to be installed with the development of the network in that area, has not yet been built.
As mentioned above, the rated voltage of the line is 110 kV, and its length is 47.4 km. It was built with 86 metal towers, most of them having the geometrical characteristics shown in Figure 3. The three phase conductors (indices 1, 2, 3 corresponding to the phases a, b, c) and the protective conductor (index 4, corresponding to conductor n, grounded-g) are aluminum conductor steel-reinforced cable (ACSR)-types and have rated sections of 185/32 mm 2 and 95/55 mm 2 , respectively. The insulators of the line are ball-and-socket-type glass strings, with a length of 1.8 m.
The three-phase modeling of the line consists of the determination of the primitive equivalent parameters and the equivalent series and shunt phase parameters. Along these lines, it should be taken into account that each conductor is located in both in its own electric and magnetic fields and in the electric and magnetic fields of the neighboring conductors, with the whole conductor assembly being located near the ground, the potential of which is null. The three-phase modeling of the line consists of the determination of the primitive equivalent parameters and the equivalent series and shunt phase parameters. Along these lines, it should be taken into account that each conductor is located in both in its own electric and magnetic fields and in the electric and magnetic fields of the neighboring conductors, with the whole conductor assembly being located near the ground, the potential of which is null.

Calculation of the Primitive Impedance Matrix
For the calculation of the self and mutual impedances of an electric line, the Carson's relations are universally accepted [1,[26][27][28][29]. For their deduction, the method of conductor images was applied in the conditions of ground modeling using the following hypotheses: the ground is a uniform solid of infinite dimensions, with the outer side a uniform surface that is perfectly flat, and whose electrical resistivity is constant. It is also considered that the end effect introduced by the grounding of the neutral (the protective) conductor is reduced for the frequency values found in the network, which is why it can be neglected. By applying additional mathematical simplifications, Carson's modified equations were obtained, which have proved to be accurate enough and therefore widely used both for OHL and underground line modeling.
In Figure 4a, the unit length segment of a three-phase OHL, having a protective conductor grounded at both ends to which the unitary primitive impedances are associated [1], is presented.

Calculation of the Primitive Impedance Matrix
For the calculation of the self and mutual impedances of an electric line, the Carson's relations are universally accepted [1,[26][27][28][29]. For their deduction, the method of conductor images was applied in the conditions of ground modeling using the following hypotheses: the ground is a uniform solid of infinite dimensions, with the outer side a uniform surface that is perfectly flat, and whose electrical resistivity is constant. It is also considered that the end effect introduced by the grounding of the neutral (the protective) conductor is reduced for the frequency values found in the network, which is why it can be neglected. By applying additional mathematical simplifications, Carson's modified equations were obtained, which have proved to be accurate enough and therefore widely used both for OHL and underground line modeling.
In Figure 4, the unit length segment of a three-phase OHL, having a protective conductor grounded at both ends to which the unitary primitive impedances are associated [1], is presented.  Carson's modified equations for calculating the unitary primitive self and mutual impedances are as follows: Carson's modified equations for calculating the unitary primitive self and mutual impedances are as follows: where z p ii is the self-impedance of conductor i in Ω/mile; z p ij is the mutual impedance between conductors i and j in Ω/mile; r i is the resistance of conductor i in Ω/mile; f is the frequency in Hz; D ij is the distance between conductor i and conductor j in ft; GMR i is the geometric mean radius of conductor i in ft; ρ is the resistivity of the earth in Ω·m; and k 1 , k 2 , k 3 are constants with the following values: According to Kirchhoff's voltage law applied to the circuit in Figure 4, one can write the matrix relation as follows:

Calculation of the Primitive Impedance Matrix
The matrix of the self and mutual phase impedances is obtained from the primitive impedance matrix by applying Kron's reduction method [30]. For this, the four matrices above are partitioned in such a way that Equation (3) can be written as follows: The partial matrices are therefore: Kron's reduction method can apply if the following is true: For a line such as the one above, Equation (6) is valid because the neutral (or protective) conductor is grounded so that the voltages between the two points are null at both ends of the considered segment. Under these conditions Equation (3) is reduced to the following: Energies 2019, 12, 734 9 of 17 In Equation (7) one can identify the phase impedance matrix, which is square, of dimension 3 × 3, and it can be determined using the primitive impedances partial matrices: [z abc ] =    z aa z ab z ac z ba z bb z bc z ca z cb z cc For the studied OHL, the results of the calculation of the primitive impedances matrix, and of the phase impedance matrix, respectively, expressed in SI units, are as follows: The calculation of the equivalent parameters was done under simplifying conditions, according to which the line conductors are considered to be parallel to the ground and the towers and the insulation has the same configuration along the whole length of the line.

Calculation of the Potential Coefficients Matrix
In order to determine the phase admittance matrix, the potential of a conductor belonging to an OHL must first be determined, located in its own electric field and in the electric fields determined by the other conductors.
Applying the image conductor method, one can determine the voltage drop between the conductor and the ground [1]. For each conductor, the electrical loads of all the conductors influence its voltage drop to the ground. This dependence corresponds to self and mutual primitive potential coefficients.
For a three-phase line having three phase conductors and a neutral (or protective) conductor, the equations between the potentials of the four conductors and their electrical loads can be written by using the primitive potential coefficients matrix: The primitive capacities matrix can be associated to the primitive potential coefficients. These capacities correspond to the capacitive couplings formed between each pair of elements of the four conductors and ground assembly ( Figure 5). Partitioning the three matrices, Equation (9) can also be written as follows: In Equation (10), the partial matrices are as follows: For the calculation of the values of the primitive potential coefficients per unit length line, Carson's equations, using the same concept of conductor images, is applied [1]: where is the self-potential coefficient, in mile/μF; is the mutual potential coefficient, in mile/μF; is the distance from conductor i and its image, in ft; is the distance from conductor i and the image of conductor j, in ft; is the radius of conductor i, in ft; is the distance from conductor i to conductor j, in ft; and = 11.17689 is the constant resulting from considering the value of the air relative permittivity = 1.424 • 10 μF/mile.

Calculation of the Potential Coefficients Matrix
The determination of the phase potential coefficients is done by observing that by grounding the protective conductor, its potential becomes null: Equation (10) is reduced to the following: Applying Kron's reduction method [30], the primitive potential coefficient matrix is reduced to the phase potential coefficient matrix, which can be calculated by an expression in which the partial primitive coefficients matrices intervene: Applying Kron's reduction method [30], the primitive potential coefficient matrix is reduced to the phase potential coefficient matrix, which can be calculated by an expression in which the partial primitive coefficients matrices intervene: The phase capacities matrix, calculated for the unit length of the line (in µF/miles), immediately results as the inverse of the phase potential coefficients matrix: According to the commonly accepted practice, most of the time, the shunt conductance of an OHL can be neglected, so that for the phase admittance matrix, expressed in µS/miles, the following relationship results: For the analyzed OHL, the calculation results of the primitive potential coefficients matrix, phase potential coefficient matrix, phase capacitances matrix, and phase capacitive admittances matrix, expressed in SI units, are as follows:  The calculation of the equivalent parameters was done under the same simplification conditions as those mentioned above.

Power Flow Calculation in Phase Coordinates
In the case of the OHLs with relatively small lengths (tens of kilometers), only the values of the electrical amounts at the ends of the line are of concern. In these cases, the calculation of the power flow is accurate enough if the line is modeled by means of the phase total impedances and admittance matrices. These matrices are calculated by multiplying the unitary phase impedance and admittance matrices, with the length of the line as follows: where L is the line length. The four-wire electric line is reduced to the three-phase conductors, and the steady-state operating conditions can be analyzed by performing the power flow calculation in phase coordinates, using an equivalent octopole with concentrated parameters (Figure 6).

= •
where L is the line length. The four-wire electric line is reduced to the three-phase conductors, and the steady-state operating conditions can be analyzed by performing the power flow calculation in phase coordinates, using an equivalent octopole with concentrated parameters (Figure 6).  The equations between phase currents and phase voltages (between phases and ground conductors), depending on the phase total impedances and admittance matrices, are as follows: Equations (21) and (22) allow for the back forward power flow calculation.
In the case of long OHLs or situations where the variation of the values of electrical amounts along the line are monitored, it should be divided into segments of small lengths, which are modeled by appropriate equivalent octopoles. The limit of two successive octopoles correspond to the sections of the line where the values of electrical amounts are monitored.

Power Flow Calculation for the Studied OHL and Comparison with Measured Values
Mathematical modeling in the phase coordinates of the no-load operating conditions of the studied OHL was performed in order to determine if the paradox described above can be achieved. For this purpose, the values of the previously calculated phase parameters were used. In order to determine the variation in electrical amounts along the line, modeling was achieved by dividing it into nine segments of equal length, which highlight the 10 sections. Each of the nine segments was modeled by an equivalent octopole with concentrated parameters. Then, Equations (21) and (22) were successively applied for the back forward power flow calculation (section k = 10-the end of the line to section k = 1-the beginning of the line). For each section the voltages, currents, and active and reactive powers were determined.
In the second part of the numerical study, a hypothetical transposition of the line was considered to evaluate the effect of this operation and to verify the correctness of the mathematical model. For this purpose, two permutations of the three-phase conductors were considered, and they were done in Sections 4 and 7, respectively. The line was thus divided into three parts of equal length (each containing three out of the nine equal segments) without considering the tower's position in the field. Phase conductors' permutations were introduced into the calculus by appropriate permutations in matrices D (distances between the four conductors' matrix) and S (distances between the conductors and their images' matrix). Figure 7 presents the results of the electrical amounts obtained by numerical calculation in the two situations: the non-transposed line (a1-4)), and the transposed line (b1-4), respectively.
The numerical study was performed by applying the Mathcad software tool under the following conditions: • Equivalent parameters of type R, L, and C are considered ideal and perfectly linear; • The voltages applied to the three-phase circuit are perfectly sinusoidal (the harmonic regime caused by the railway traction is neglected). The numerical study was performed by applying the Mathcad software tool under the following conditions: • Equivalent parameters of type R, L, and C are considered ideal and perfectly linear; • The voltages applied to the three-phase circuit are perfectly sinusoidal (the harmonic regime caused by the railway traction is neglected). Only the values calculated for Section 1 can be compared to those measured. Therefore, in the numerical study carried out for both the untransposed line and the transposed version, we aimed to obtain values for the voltages as close as possible to the measured values and to their average value, respectively. From the set of recorded values, the ones from 2 October 2018 13:01:19 '584 were retained, corresponding to a very short period in which the operating regime was quasi-stationary.
For an easier comparison, the measured and calculated values were grouped in Table 1. Considering the calculated values (Figure 6), the comparison between them, and the measured values (Table 1), the following observations can be made:

•
In both situations, the calculated voltages at the end of the line were higher than those at the beginning of the line, which is in line with expectations. This is the consequence of producing the Ferranti effect, typical for AC power lines operating in no-load conditions.  Only the values calculated for Section 1 can be compared to those measured. Therefore, in the numerical study carried out for both the untransposed line and the transposed version, we aimed to obtain values for the voltages as close as possible to the measured values and to their average value, respectively. From the set of recorded values, the ones from 2 October 2018 13:01:19 '584 were retained, corresponding to a very short period in which the operating regime was quasi-stationary.
For an easier comparison, the measured and calculated values were grouped in Table 1. Considering the calculated values (Figure 6), the comparison between them, and the measured values (Table 1), the following observations can be made:

•
In both situations, the calculated voltages at the end of the line were higher than those at the beginning of the line, which is in line with expectations. This is the consequence of producing the Ferranti effect, typical for AC power lines operating in no-load conditions. • The variation along the line of rms currents and reactive powers is in line with expectations.
They are null at the end of the line and have their maximum values at the beginning of the line, showing a basically linear variation. The currents have a pronounced capacitive behavior, which is reflected in the character of the reactive powers and in their much higher values than those of the active powers.

•
The numerical study of the untransposed line confirms the paradox found by measurements to the studied OHL. At the beginning of the line, on two of the phases, the active power is negative (delivered to the network to which it is connected), and on the third one, it is positive (taken from the network). In absolute values, the variation along the line of the active power is also linear. However, the sum of the values at the beginning of the line (Section 1) is positive (consumed), which is in line with expectations and confirms the correctness of the mathematical model.  As a result of the numerical study, it can be stated that there is a high similarity between the measured values and those calculated for the untransposed line. The lower or higher differences are primarily determined by the errors of the measurement systems, but also by the simplifications introduced in the mathematical model. The errors occurred in the measured values due to the presence of the harmonic regime determined by the railway in the network to which the studied OHL is connected, must not be forgotten.

The Paradox Explanation
The paradox of the appearance of negative active powers at the beginning of an untransposed OHL is not actually a paradox. This regime is determined by the asymmetry of the equivalent capacities that are formed between the phase conductors. For the studied line, this asymmetry was accentuated as a result of both the asymmetry of the mutual position of the phase conductors and the asymmetry of the position of the phase conductors with respect to the ground and the protective conductor.
It is known that when three single-phase capacitor banks having unequal capacitances are connected between the phase conductors of a three-phase circuit (in delta connection), they will absorb capacitive reactive powers on the three phases of the network, and in addition, they will determine an unequal active power flow [24]. Depending on the type and magnitude of the asymmetry of the three equivalent capacitances, they will take over the active power from the network during one or two phases and deliver this power back into the network over the other phases (phase). There is a redistribution of active powers between the phases of the network. The operation of unbalanced capacitive compensators designed to load balancing in three-phase networks is based on this phenomenon [31].
The same phenomenon also occurs between the phases of an untransposed OHL operating in no-load conditions, although the equivalent capacitances are uniformly distributed along the line. By making the algebraic sum of the three active powers on the phases, a positive value is always obtained. This value is much lower than the values of the active powers flowing through the phases. This is nothing other than the consumption of active power in no-load operating conditions on all the three phases, corresponding to losses through the Joule-Lenz effect and the corona effect. The transposition of phase conductors of an OHL leads to the permutation of the redistribution of active powers between the phases, so that over the entire length of the line, this redistribution is nullified and a symmetrization of all the electrical amounts at the ends of the line is obtained.

Conclusions
In practice, there are many situations in OHL operation for which solving and understanding is not useful to apply symmetrical coordinates modeling. This is especially true the case of untransposed OHLs, for which is necessary to return to modeling in phase coordinates.
Carson's modified equations are applied to the calculation of the primitive equivalent impedances and the primitive potential coefficients. Kron's transformation method allows for the determination of equivalent phase impedances, phase potential coefficients, and phase capacitances, respectively. Both are practical solutions for the phase coordinate modeling of an OHL, and they are accurate enough and are thus universally accepted. In the paper, these solutions were used to obtain a mathematical confirmation of a so-called paradox found for OHLs operating in no-load conditions. For the case study, a real OHL was considered, for which the mathematical modeling was performed in phase coordinates. The results of the numerical study of the no-load regime led to similar results to those measured under real conditions.
In fact, the asymmetry of equivalent phases capacitances, even if they are uniformly distributed along the line, determines for an untransposed OHL operating in no-load conditions, the same effect as in the case of an asymmetrical capacitor bank used for shunt capacitive compensation and load balancing. Besides the capacitive reactive power consumption from the network to which it is connected, it determines the redistribution of the active power between the network phases.