Nonlinear Switching Transient Field Simulation of Cable Joint without Residual Charge

Mingyan Wu¹, Jun Xiong¹, Lei Liao¹, Lu Zhu¹, Ruxin Zhang¹, Zheng Wu¹, Gang Du¹, Xueyou Huang¹, Haiming Li², Jian Zhang², Sizhuo Liao², Binxian Lu^2*

¹Guagnzhou Power Supply CSG, Guangzhou, China.
²Department of Electrical and Electronic Engineering, North China Electric Power University, Beijing, China.

DOI: 10.4236/epe.2020.124B005                

Abstract

The analysis of the impulse voltage on the internal electric field of the cable joint plays a key role in studying the breakdown of the joint. Based on the finite element method, a three-dimensional electromagnetic field simulation model of the cable joint is established in this paper. Simulation results show that the voltage at the head of the cable joint reaches about twice the impulse voltage. The increase of the conductivity of semi-conductive material also leads to the increase of electric field intensity. Then, several points and curves at different positions are selected for further analysis in this paper. Among them, the electric field distortion at the edge of the high voltage shield is the most serious and the electric field in the air gap is the least.

Keywords

Cable Joint, Switching Transient Field, 3D Transient Electromagnetic Simulation Model, Finite Element Method

1. Introduction

With the increasing use of power cables in medium voltage distribution networks, the number of cable joints in use is increasing. The quality of the cable joint manufacturing and installation varies, making the cable joint a weak link in medium voltage cables. Once the insulation performance of the cable joint is deteriorated, it will cause partial discharge and even cause serious accidents such as explosions, which will affect the safe operation of the power grid. During the switching operation, the overvoltage formed during the transient process propagates through the cable joints in the form of electromagnetic waves. Due to the existence of refraction and reflection phenomena, the components of each traveling wave are superimposed, which makes the transient overvoltage rises extremely fast, with large amplitude and high frequency, posing a great threat to the insulation of cable joints. When the conductivity of the semi-conductive material in the cable joint changes, it will affect the electric field distribution of the cable joint and seriously threaten the operational reliability of the cable joint [1].

In order to assure reliable operation of cable joint, a lot of research has been conducted on the space charge characteristics, interface breakdown and flashover characteristics, and optimization of the accessory structure, so as to minimize the possibility of problems in the production, installation and use of cable accessories [2] [3] [4] [5] [6]. These studies are based on steady-state fields and cannot reflect transient processes during switching operations. Therefore, in order to more accurately study the overvoltage and transient electromagnetic fields in cable joints, a full-wave simulation analysis of the three-dimensional electromagnetic field model is needed.

Based on the finite element method, a three-dimensional electromagnetic field simulation model of the cable joint is established in this paper. Under the shock wave form of double exponential wave, the electric field distribution and overvoltage in the cable joint are calculated, and the influence of the conductivity of the semi-conductive material on the electric field distribution and overvoltage is analyzed.

2. Simulation Model Description

2.1. Structure and Material of Model

The structure of the 10 kV AC cable connector is shown in Figure 1. The model used in this paper is obtained by rotating the axisymmetric model around the axis of symmetry. The model consists of 9 parts: conductive copper, inner semi-conductive layer, main cable insulation layer, outer semi-conductive layer, stress cone, high-voltage shield, reinforced insulation layer and insulation shield. The material parameters of each part are shown in Table 1.

Figure 1. Structure of the 10 kV cable joint.

2.2. Control Equation

To calculate the transient field, Maxwell's equations need to be solved. In the region without charge density and current density, the existence of the potential φ can be ignored, and only the magnetic vector potential A is used to represent the electric field strength and magnetic field strength. The control equation is shown in Equations (1) as follow.

$\nabla \times \left(\frac{\text{1}}{{\mu }_{\text{r}}}\nabla \times A\right)+{\mu }_0\sigma \frac{\partial A}{\partial t}+{\mu }_0{\epsilon }_0\frac{\partial }{\partial t}\left({\epsilon }_{\text{r}}\frac{\partial A}{\partial t}\right)=0$ (1)

where t is time, μ_r is relative permeability, μ₀ is the permeability of vacuum, σ is electricity conductivity, ε₀ is the permittivity of vacuum, ε_r is relative permittivity.

The head and end of the cable joint are set as lumped port I and lumped port II, respectively, and the relationship between the voltage U_port and current I_port of the port is expressed as

$Z_{\text{port}}=\frac{U_{\text{port}}}{I_{\text{port}}}$ (2)

where Z_port is the wave impedance of lumped port. The wave impedance of the lumped port at the head is set to 5 Ω. Although the internal impedance of the power supply should actually change with frequency, it is approximated in this paper with a small resistance. The impedance of the lumped port at the end is set to 5 GΩ to simulate an open circuit condition in the laboratory. And the curve of the overvoltage at the end of the cable joint is extracted from this port.

A wave excitation is applied at the head of the joint, and the waveform is shown in Figure 2. The equation of the impulse wave is expressed as

Table 1. Material parameter.

Figure 2. The waveform of the impulse wave.

$U\left(t\right)=10\sqrt{\frac{2}{3}}\left[\exp (-4\cdot {10}^{6}t)-\exp (-4.76\cdot {10}^{8}t)\right]$ (3)

Since the influence of residual charge is not considered, the initial condition of the magnetic vector potential A is set as

$A\left(r,0\right)=0$ (3)

${\frac{\partial A}{\partial t}|}_{t=0}=0$ (4)

3. Results and Discussion

According to the methods and settings above, the overvoltage of lumped ports and the electric field distribution under different conductivity of semi-conductive materials are simulated.

3.1. Simulation Analysis of Overvoltage under Different Conductivity

The voltage of lumped port I and impulse voltage under different conductivity is compared in Figure 3. It can be seen from Figure 3(a) that the voltage of lumped port I is also in the form of a double exponential function, reaching twice the applied voltage. Figure 3(b) is enlarged by the voltage peak time of Figure 3(a). As can be seen from Figure 3(b), when the conductivity changes from 5 S/m to 20 S/m, the voltage gradually decreases with a small change.

And the voltage of lumped port II under different conductivity is compared in Figure 4. Obviously, the voltage oscillation amplitude of lumped port II is larger, and the voltage decreases with the increase of conductivity.

3.2. Analysis of Electric Field Distribution under Different Conductivity

Figure 5 shows the distribution of the electric field under different conductivity when the voltage reaches its peak. As the conductivity increases, the maximum electric field intensity increases, but the growth rate slows down.

For a more detailed analysis of the distribution of the electric field, four points at different locations, namely A, B, C and D, and four curves, l₁, l₂, l₃ and l₄, are

(a) global figure (b) zoomed curve

Figure 3. Voltage of the Lumped Port I vs. Time.

(a) global figure (b) zoomed curve

Figure 4. Voltage of the Lumped Port II vs. Time.

Figure 5. The distribution of the electric field under different conductivity.

selected, as shown in Figure 6.

Figure 7 shows the electric field intensity at point A, B, C and D under different conductivity over time. It can be seen that the conductivity has a greater effect on point A and D, and a smaller effect on point B and C. As the conductivity increases, the electric field intensity at point A increases. But the conductivity has the opposite effect on the electric field intensity at points B, C and D. Among the four points, the electric field distortion is the most serious at point B.

Figure 8 shows the distribution of electric field intensity along the curve l₁, l₂, l₃ and l₄ when the voltage reaches its peak under different conductivity. As the conductivity increases, the growth rate of electric field intensity gradually saturates. Among the four curves, the electric field intensity on curves l₁ and l₂ is larger, while the electric field intensity on curves l₄ is the smallest.

Figure 6. The distribution of the points and curves.

Figure 7. The electric field intensity at point A, B, C and D under different conductivity.

Figure 8. The electric field intensity along the curve l₁, l₂, l₃ and l₄ under different conductivity.

4. Conclusion

Due to the multiple refraction and reflection inside the cable joint, the voltage at the head of the cable joint reaches about twice the impulse voltage. Moreover, the voltage oscillation amplitude at the end of the cable joint is larger than that at the head of the cable joint. As the conductivity increases, the voltage at the head and end decreases slightly. The increase of conductivity also leads to the increase of electric field intensity. After that, different points and curves at different positions are selected for further analysis in this paper. The electric field distortion at the edge of the high voltage shield is the most serious and the electric field in the air gap is the least.

Acknowledgements

This work was supported by the Science and Technology Project of China South Power Gird Co., Ltd., (GZHKJXM20180140).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Mazzanti, G., Montanari, G.C., and Palmieri, F. (2003) Quantities Extracted from Space-charge Measurements as Markers for Insulation Aging. IEEE Transactions on Dielectrics and Electrical Insulation, 10, 198-203. http://dx.doi.org/10.1109/TDEI.2003.1194100
[2]	Tohmine, T., Fujitomi, T., Miyake, H., Tanaka, Y., Ida, Y. and Inoue, Y. (2017) Measurement of Space Charge Accumulated in Multi-layered Samples Composed of Different Insulators Used in the Joints of DC Transmission Cables. 2017 International Symposium on Electrical Insulating Materials (ISEIM), Toyohashi, 2017, 299-302. http://dx.doi.org/10.23919/ISEIM.2017.8088746
[3]	Yaroslavskiy, V., Walker, M., Katz, C. and Keefe, R. (2008) Comparative Laboratory Evaluation of Premolded Joints for Medium Voltage Cables. IEEE Transactions on Power Delivery, 23, 516-522. http://dx.doi.org/10.1109/TPWRD.2007.915909
[4]	Aoran, X. and Yueping, M. (2019) Research Status and Trend of Insulation State Detection Technology for Power Cable Joints. Science and Technology Innovation and Application, 31, 28-31.
[5]	Guo, Q., Xiong, J., Dai, N., Liao, S. and Lu, B. (2018) Influence of the Material Parameters on the Electric Field Intensity Distribution of 10kV Cable Elbow Connector. 2018 2nd IEEE Conference on Energy Internet and Energy System Integration (EI2), Beijing, 2018, 1-5. http://dx.doi.org/10.1109/EI2.2018.8582320
[6]	Heinrich, R., Bonisch, S., Pommerenke, D., Jobava, R. and Kalkner, W. (2000) Broadband Measurement of the Conductivity and the Permittivity of Semiconducting Materials in High Voltage XLPE Cables. 2000 Eighth International Conference on Dielectric Materials, Measurements and Applications (IEE Conf. Publ. No. 473), Edinburgh, 2000, 212-217. http://dx.doi.org/10.1049/cp:20000507