Telegraph equations are derived from the equations of transmission line theory. They describe the relationships between the currents and voltages on a portion of an electric line as a function of the linear constants of the conductor (resistance, conductance, inductance, capacitance). Their resolution makes it possible to determine the variation of the current and the voltage as a function of time at each point of the line. By adopting a general sinusoidal form, we propose a new exact solution to the telegraphers’ partial differential equations. Different simulations have been carried out considering the parameter of the 12/20 (24) kV Medium Voltage Cable NF C 33,220. The curves of the obtained solution better fit the real voltage curves observed in the electrical networks in operation.
Electrical transmission lines conduct single-phase, two-phase, or three-phase electrical voltage and current in opposite directions to each other [
M. K. Smail [
Other authors have proposed analytical solutions, notably C. R. Paul [
All of the above solutions have given interesting results in their respective fields of application. As for the studies presented in [
Many of the analytical solutions are more mathematical than physical. Their representative curves do not always fit the actual voltage and current curves seen on power system control and supervision systems.
To this end, a new exact solution could improve electricity management by minimizing energy losses on the lines.
This paper proposes a new exact solution of telegraph equations for better electricity management. This solution describes, at a given location, the time-variation of the voltage amplitude in an electric cable. It is suitable for single-phase, two-phase, and three-phase voltage. The resulting solution allows the generation of the electrical wave and the monitoring of its propagation along the line. It allows following the variation of the wave as a function of time, space, and phase. Its curves better fit the shape of the instantaneous voltage in an electrical network in operation than those of the previous solutions. This new solution presents the variation of the signal shape as a function of time and the angle of the phase shift between the voltage and current signals in the network.
The remainder of the paper is structured in three parts: The first part presents the methodical approach adopted to describe the electrical model and the different equations. The second part deals with formulating the proposed exact solution of the telegraph equations and presents and discusses simulation results. The third part presents the conclusion and perspectives for further work.
Electrical transmission lines interconnect a source to an electrical charge by at least two conducting wires composed of a charge wire and a neutral one. The charge wire conducts the electrical signal while the neutral wire is used to form a circuit loop.
A portion of a two-wire line is usually modeled by a quadrupole as shown in
In the quasi-steady state, the equations relating voltage and current will be of the partial derivative type concerning time and space. Thus, by applying Kirchhoff’s laws to
v ( x , t ) − R ⋅ d x ⋅ i ( x , t ) − L ⋅ d x ⋅ ∂ i ( x , t ) ∂ t − v ( x + d x , t ) = 0 (1.1)
i ( x , t ) − G ⋅ d x ⋅ v ( x , t ) − C ⋅ d x ⋅ ∂ v ( x , t ) ∂ t − i ( x + d x , t ) = 0 (1.2)
By introducing the notion of finite differences, Equations (1.1) and (1.2) become (1.3) and (1.4):
− v ( x + d x , t ) − v ( x , t ) d x = R ⋅ i ( x , t ) + L ⋅ ∂ i ( x , t ) ∂ t (1.3)
− i ( x + d x , t ) − i ( x , t ) d x = G ⋅ v ( x , t ) + C ⋅ ∂ v ( x , t ) ∂ t (1.4)
Assuming that dx is infinitesimal, we obtain Equations (1.5) and (1.6):
− ∂ v ( x , t ) ∂ x = R ⋅ i ( x , t ) + L ⋅ ∂ i ( x , t ) ∂ t (1.5)
− ∂ i ( x , t ) ∂ x = G ⋅ v ( x , t ) + C ⋅ ∂ v ( x , t ) ∂ t (1.6)
The second derivative of Equations (1.5) and (1.6) concerning space gives Equations (1.7) and (1.8):
− ∂ 2 v ( x , t ) ∂ x 2 = R ⋅ ∂ i ( x , t ) ∂ x + L ⋅ ∂ 2 i ( x , t ) ∂ x ∂ t (1.7)
− ∂ 2 i ( x , t ) ∂ x 2 = G ⋅ ∂ v ( x , t ) ∂ x + C ⋅ ∂ 2 v ( x , t ) ∂ x ∂ t (1.8)
By combining Equations (1.1), (1.5), and (1.7), respectively (1.2), (1.6), and (1.8), we find Equations (1.9) and (1.10) below. These equations correspond to the one-dimensional wave equations called the telegraph equations of voltage (1.9) and current (1.10).
∂ 2 v ( x , t ) ∂ x 2 = L ⋅ C ⋅ ∂ 2 v ( x , t ) ∂ t 2 + ( R ⋅ C + L ⋅ G ) ⋅ ∂ v ( x , t ) ∂ t + R ⋅ G ⋅ v ( x , t ) (1.9)
∂ 2 i ( x , t ) ∂ x 2 = L ⋅ C ⋅ ∂ 2 i ( x , t ) ∂ t 2 + ( R ⋅ C + L ⋅ G ) ⋅ ∂ i ( x , t ) ∂ t + R ⋅ G ⋅ i ( x , t ) (1.10)
In this paper, we propose an exact solution to Equation (1.9). This latter solution depends on time, position, and phase. It can be used to simulate the variation of single-phase, two-phase, or three-phase voltage on electrical power lines.
Equations (1.9) and (1.10) have similar forms. They well describe the propagation of voltage and current in electric cables.
For Equation (1.9), we propose a solution of the form:
v ( x , φ , t ) = u ( x , φ ) e ( − j ω t + φ ) (1.11)
where t is the time variable, x is the position and φ is the phase. u ( x , φ ) is the wave amplitude and ω is the wave pulsation.
The solution of Equation (1.10) can be deduced from the expression:
v ( x , φ , t ) = Z ( x , φ ) i ( x , φ , t ) (1.12)
where Z ( x , φ ) is the characteristic impedance of the transmission line. Therefore, finding the solution to one of the equations allows us to deduce the other solution.
Let’s pose a = LC; b = RC + LG and c = RG, then Equation (1.9) becomes:
∂ 2 v ( x , t ) ∂ x 2 = a ∂ 2 v ( x , t ) ∂ t 2 + b ∂ v ( x , t ) ∂ t + c v ( x , t ) (1.13)
Let v ( x , φ , t ) given by Equation (1.11) be the general form of the solution of Equation (1.9). By derivating two times v ( x , φ , t ) according to space and time, we obtain:
∂ 2 u ( x , φ ) ∂ x 2 + [ ( a ω 2 − c ) + b j ω ] u ( x , φ ) = 0 (1.14)
Equation (1.14) is a second-order differential equation without a second member whose characteristic equation is given by (1.15)
r 2 = ( c − a ω 2 ) − j b ω (1.15)
Let r = α + j β such that r 2 = α 2 − β 2 + 2 j α β , and | r | 2 = α 2 + β 2 .
We derive the system of equations given by:
{ α 2 − β 2 = c − a ω 2 ( a ) α 2 + β 2 = ( c − a ω 2 ) 2 + ( b ω ) 2 ( b ) 2 α β = − b ω ( c )
whose solution is given by:
α = ± c − a ω 2 + ( c − a ω 2 ) 2 + ( b ω ) 2 2 (d)
β = ± a ω 2 − c + ( c − a ω 2 ) 2 + ( b ω ) 2 2 (e)
The product α ⋅ β = − b w 2 < 0 . Since α denotes the attenuation coefficient we choose α < 0 and β > 0 . We derive a real solution of Equation (1.14) given by
u ( x , φ ) = ( A ⋅ cos ( β x + φ ) + B ⋅ sin ( β x + φ ) ) e ( α ⋅ x − φ ) (1.16)
Then a solution to Equation (1.9) is:
v ( x , φ , t ) = ( A cos ( β x + φ ) + B sin ( β x + φ ) ) e [ α ⋅ x − j ω t ] (1.17)
where A andB are time-dependent parameters to be determined with boundary conditions on position, phase, and time.
The parameters A and B are time-dependent. To determine them let us apply the boundary conditions to x and φ . When x = 0 and φ = 0 , the line is disconnected from the distribution source, therefore the voltage is equal to that of the production source which is v 0 ( t ) .
Similarly, if x = 0 , and t = 0 , the voltage depends on the phase and is noted v 0 ( φ ) .
Thus the expressions of the parameters A and B are determined as follows:
· If x = 0 et φ = 0 then v ( x , φ , t ) = v ( 0 , 0 , t ) = v 0 ( t ) = A e − j ω t therefore A = v 0 ( t ) e j ω t
· If x = 0 et t = 0 then
v 0 ( φ ) = A cos ( φ ) + B sin ( φ ) = A 2 + B 2 cos ( β + φ )
with
cos ( β ) = A A 2 + B 2 et sin ( β ) = B A 2 + B 2
Also d v 0 ( φ ) d φ = − A sin ( φ ) + B cos ( φ ) = 0 therefore
v 0 ( φ ) = ( A + B A ) cos ( φ )
Assuming that β is very small compared to φ , we have:
A + B A = A 2 + B 2 , hence
B = 2 A 2 A 2 − 1 = 2 v 0 2 ( t ) ⋅ e j 2 ω t v 0 2 ( t ) ⋅ e j 2 ω t − 1 .
Given A and B the expression of the voltage is:
v ( x , φ , t ) = [ v 0 ( t ) cos ( β x + φ ) + 2 v 0 2 ( t ) ⋅ e j ω t v 0 2 ( t ) ⋅ e j 2 ω t − 1 sin ( β x + φ ) ] e α ⋅ x (1.18)
To find a real solution to the telegraph equations, we decompose Equation (1.18) into its real and imaginary parts. To do this, we write the expression of B e − j ω t in its real and imaginary forms.
B e − j ω t = v 0 2 ( t ) cos ω t − cos ω t − j [ v 0 2 ( t ) sin ω t + sin ω t ] v 0 2 ( t ) 2 + 1 2 ⋅ v 0 2 ( t ) − cos 2 ω t (1.19)
By replacing this expression of B e − j w t in Equation (1.18) we deduce that:
v ( x , φ , t ) = [ v 0 ( t ) cos ( β x + φ ) + v 0 2 ( t ) cos ω t − cos ω t − j [ v 0 2 ( t ) ⋅ sin ω t + sin ω t ] v 0 2 ( t ) 2 + 1 2 ⋅ v 0 2 ( t ) − cos 2 ω t sin ( β x + φ ) ] e α ⋅ x
(1.20)
When we decompose the voltage expression into its real part ( R e [ v ] ) and imaginary one ( I m [ v ] ) where v ( x , φ , t ) = R e [ v ] + j I m [ v ] , we identify:
R e [ v ] = [ v 0 ( t ) cos ( β x + φ ) + v 0 2 ( t ) cos ω t − cos ω t v 0 2 ( t ) 2 + 1 2 ⋅ v 0 2 ( t ) − cos 2 ω t ] e α x
I m [ v ] = − v 0 2 ( t ) sin ω t + sin ω t v 0 2 ( t ) 2 + 1 2 v 0 2 ( t ) − cos 2 ω t sin ( β x + φ ) e α ⋅ x
From this solution, we run simulations on R e [ v ] and I m [ v ] as a function of cable length, time, and the angle θ of the cosine of the electrical network.
The solution of the telegraph equations depends on the linear constants R, C, L, and G of the power line. In this paper, we consider the Medium Voltage Cables NF C 33,220 standards of 12/20 (24) kV. It is an aluminum cable of a nominal cross-section of 150 mm2 (square millimeter), whose linear constants in AC transmission current of frequency 50 Hz and temperature 90˚C are R = 0.265 Ω/km, C = 0.24 μF/km, L = 0.41 mH/km and G = 1/R Siemens/m. In the literature the phase shifts often considered are: φ = 60˚; 90˚; 120˚.
For simulations, the source voltage v 0 ( t ) considered is
v 0 ( t ) = V max sin ( ω t + θ ) .
where:
- V max is the maximum voltage on a medium voltage transmission power line. In Burkina Faso the effective voltages used for medium voltage transmission are 15,000 V; 20,000 V and 33,000 V with a tolerance of ±5%. For this study, we consider the RMS voltage of 20,000 V i.e. a maximum voltage V max = 20000 2 V ± 1000 2 V .
- θ is the phase shift between voltage and current in the power system. It is determined by the cosine of the network which is cos θ = 0.8 or θ = 0.20 π rad.
- The pulsation is ω = 100 π rad / s with a transmission frequency f = 50 Hz .
Three situations are considered for the simulations:
Situation 1: Single-phase signal in three dimensions (3D) space.
We plot here the variation of the real and imaginary parts of the voltage v ( x , φ , t ) as a function of time and length.
Situation 2: Single-phase, two-phase, and three-phase signals in the plane.
Firstly, we plot v ( x , φ , t ) according to the positionx, when the time is set to t = 3 s .
Secondly, we plot v ( x , φ , t ) as a function of time. For this purpose, we set the wave speed to v = 5 m / s and fix the length x of the cable at 0.05 km.
Situation 3: Single-phase, two-phase, and three-phase signals in three dimensions projected on the plane.
In this third part, we vary v ( x , φ , t ) as a function of time and length simultaneously and then plot a projection in the time plane. This corresponds to the curves visualized by electrical network supervision systems.
These two figures show a sinusoidal character of the voltage according to the length and time.
This paper proposes an exact solution to telegraph equations. For the founded solution, the amplitude of the voltage varies with position and time. By taking the phase into account, the proposed solution allows the simulation of single-phase, two-phase and three-phase voltages. The choice of the initial condition of the voltage is a factor that influences the shape of the solution. The simulations showed that the voltage variation corresponds well to the patterns observed by electrical network supervision systems.
This work could be extended to a multi-branch power line section to describe the voltage behavior in the medium voltage (MV) and low voltage (LV) networks. Also, this exact solution could be explored to study the stability of an electrical network in the case of an energy mix. The a priori knowledge of the voltage level as a function of position and time, provided by the proposed solution, could also be exploited to detect and locate faults on a power line.
The authors declare no conflicts of interest regarding the publication of this paper.
Konane, D., Ouedraogo, W.Y.S.B., Guingane, T.T., Zongo, A., Koalaga, Z. and Zougmoré, F. (2022) An Exact Solution of Telegraph Equations for Voltage Monitoring of Electrical Transmission Line. Energy and Power Engineering, 14, 669-679. https://doi.org/10.4236/epe.2022.1411036