1. Introduction
With the development of urbanization, transmission and distribution corridors become increasingly tense. Therefore, the construction of urban transmission lines is of great significance to meet the needs of urban economic development [
1]. In the construction of urban transmission lines, due to the scarcity of urban land resources, the transmission corridors of overhead lines are often tense, and the construction of overhead lines affects the urban landscape. In the construction of urban transmission lines, transmission cables with a high utilization rate of transmission corridors and high environmental friendliness are often used. In this case, the safe operation of power cables becomes critical. In terms of distribution cables, compared with single-core cables, three-core symmetrical power cables are favored in medium-low voltage (below 35 kV) distribution projects due to their lower cost [
2]. The current is one of the most important electrical state quantities, and the current can reflect the running state of the three-core cable, so the phase current measurement of cables is of great significance [
3].
Current monitoring has been investigated by different researchers. With the development of magnetic measurement technology, the use of magnetic sensors to measure current has become a method recognized by the power industry. Reference [
4] designed a circular magnetic field sensor array with a 45° difference for the arc current of a single-core cable. In reference [
5], a current sensor composed of three Hall sensors with a difference of 120° is constructed, and a method of single-conductor current measurement is proposed. Reference [
6] designed a Hall sensor with a phase difference of 90° to measure the current of a single conductor. Reference [
7] proposed a method to measure the magnetic field of an 11 kV single-core cable with a linear giant reluctance array. Compared with the ferromagnetic current transformer, this method improved the accuracy. However, the measurement object of the above literature is the current of single conductor, and it is not suitable for multi-conductor systems such as three-core cables.
In general, phase A, phase B, and phase C are distributed symmetrically with a 120° difference in the three-core power cable [
8]. Due to such geometry, traditional current measurement methods based on an electromagnetic induction principle are not suitable, since the magnetic field cancels out around the surface of the cable.
For the three-core cable, reference [
9] designed a giant magnetoresistance sensor to monitor the charging status, but the method of current measurement is not given in the paper. Reference [
10] proposed an analytical calculation method for the surface magnetic field of three-core cables. However, this method requires three magnetic sensors to be installed at specific positions, which limits its application.
In this paper, the mathematical model of the phase current and magnetic field around three-core cables is established based on the Ampere theorem. By arranging multiple sensors around the cable to measure the magnetic field, the phase current in the three-core cable can be effectively measured via an inversion model. The proposed methodology is validated through a test.
2. Materials and Methods
The structure of the general three-core cable is shown in
Figure 1, and the parameters are shown in
Table 1.
The surface magnetic flux density of the three-core cable rotates with the phase current. And the surface magnetic field is the vector superposition of the three-phase current magnetic field. Since the distance from a particular point on the surface of the cable to the three cores is different, the superimposed magnetic field is not perfectly zero [
11].
In this paper, four magnetic sensors are used for phase current measurement. The sensors are numbers S
1–S
4 and placed on the surface of the three-core cable. As shown in
Figure 2, point O is the center of the three-core cable while points A, B, and C are the centers of phase A, B, and C cores, respectively. The distance between point O and the center of each core is
,
, and
. The radius of the cable is
R. The angle between OA, OB, OC, and the
x axis forward is
α,
β,
γ, respectively, with 120° difference. The coordinates of points A, B, and C can be obtained as (
cos
α,
sin
α), (
cos
β,
sin
β), and (
cos
γ,
sin
γ), respectively.
The magnetic flux density
B1–
B4 along the tangential direction generated at S
1–S
4 can be obtained, as shown in Equations (1)–(4) [
12]. In the equations,
μ represents the permeability of the vacuum. The three-phase current
,
,
, and rotation angle
α are unknown. The cable radius
R, the distance between the cable center and the three-phase core center
,
,
, and the magnetic flux density measured using the magnetic sensor are known. The number of unknowns in the equations is smaller than the number of equations, which makes the equations overdetermined. And the phase current of the cable can be obtained by solving the equations. The surface magnetic field of a three-core cable is a vector superposition of each current magnetic field. The equations of tangential magnetic flux density can be derived from the Ampere-loop theorem [
12]. Taking Equation (1) as an example, the three parts in Equation (1) are the magnetic flux density of each current at S
1, and the total magnetic flux density can be obtained by superimposing them.
Because of the sinusoidal term of rotation angle α, the equations are nonlinear. It is impossible to obtain a specific expression for the phase current based on magnetic flux density. The problem can be solved numerically to obtain the current. The nonlinear least square method is used to approximate the optimal solution. The phase currents for the three-core cable can be obtained using this method.
In this paper, Newton–Raphson algorithm is used for the nonlinear least square method. The algorithm is the result of Taylor expansion and the approximation of nonlinear equations. For a given function
f(
x) (least square between actual magnetic flux density and estimated magnetic flux density),
x0 is the approximate root of the function and
f′(
x) is the first derivative of the function. The function is considered to be locally linearized near
x0. The specific equation is shown below. The final result is obtained through continuous iteration, where
n is the iteration number.
In order to show the magnetic field distribution, the finite element method is used for the magnetic field analysis. The parameters in
Table 1 are used for the simulation. Furthermore, the frequency in the simulation is 50 Hz. The current amplitude is 707.1 A. The three–phase currents differ by 120° in the phase. The mesh diagram with triangles is shown in
Figure 3. The governing equations to analyze the magnetic field are shown in Equations (8)–(11). The result is shown in
Figure 4. Along the surface red route indicated in
Figure 4, the magnetic flux density around the cable is extracted and shown in
Figure 5.
where
H represents the magnetic intensity,
J represents the conducted current density,
D represents the electric displacement,
E represents the electric field intensity,
A represents the magnetic potential, and
represents the conductivity.
As can be seen from
Figure 4, the magnetic field can be observed outside the cable, since the distances from the cable outer surface to three cores are different. The magnetic flux density is not evenly distributed at the measuring positions of the magnetic sensors. It can be seen from
Figure 5 that the magnetic flux density of the surface rotates as the location changes, with a maximum value occurring every 120°. Eight different points are selected in
Figure 5. The magnetic flux density results derived using Equations (1)–(4) and the finite element method are shown in
Table 2. The comparison shows that the error is about 2%, which confirms the accuracy of Equations (1)–(4). The main reason for the error is that there is a skin effect and a proximity effect, which makes the cable current distribution uneven.
As for the influence of steel wire armor on magnetic field distribution, the magnetic field simulation after removing the steel wire armor is performed in this paper, as shown in
Figure 6.
The magnetic flux density of the three-core cable without steel wire armor is shown in
Figure 6. Since the three-core cable has a balanced current, the steel wire armor has little influence on the magnetic field distribution of the three-core cable, and it has little influence on the measurement results of the magnetic sensor. So, the magnetic field distribution in
Figure 6 is almost identical to
Figure 4.
3. Validation of Sensor’s Performance via Single-Core Cable
In the follow-up experiment, the magnetic sensors will be used for the current measurement of the single-core cable. This aim of this section is to test the performance of the magnetic sensor and then verify whether the related experiments of the three-phase cables can be carried out. In the experiment, the cable is long enough to be considered as an infinite straight wire. The magnetic sensor used is the commercial TMR2905 model [
13]. It is produced by MultiDimension Technology Co., Ltd., which located in Jiangsu, China.
For a single-phase wire, the relationship between its magnetic flux density and current is shown in Equation (12) [
14].
where
μ represents the vacuum permeability,
I represents the current energized in the wire, and
x represents the distance between the measuring point and the wire. The sensitivity of the magnetic sensor used is 2 V/mT, and the corresponding magnetic field strength can be calculated by the voltage measured from the magnetic sensor. The magnetic flux density can then be used to derive the current based on Equation (12). The deduced current is compared with the actual current and the relative error can be calculated. Based on these results,
Figure 7 is plotted for better visualization.
It can be seen from
Figure 7 that the error between the actual current results and the calculated results is less than 10%, indicating that this method is effective for the measurement of cable current, and when the current is greater than or less than 10 A, the larger the current, the smaller the experimental error. The errors are mainly instrumental errors and environmental errors. There are errors in the calibration and sensitivity of the instrument. In addition, the ambient temperature and air pressure will also affect the experimental results. This experiment lays a foundation for the phase current measurement of the following three-phase cable.