# A Practical Algorithm for Calculating the Impulse Earthing Resistances of Vertical Earthing Electrodes

Yongzheng Zhang1*, Xiaoqing Zhang2

1Global Energy Interconnection Research Institute Co., Ltd., Beijing, China.

2School of Electrical Engineering, Beijing Jiaotong University, Beijing, China.

*Corresponding author:Yongzheng Zhang, Global Energy Interconnection Research Institute Co., Ltd., Beijing, China. E-mail: zhangyz2@126.com

Citation: Zhang Y, Zhang X. A Practical Algorithm for Calculating the Impulse Earthing Resistances of Vertical Earthing Electrodes. J Electron Adv Electr Eng. 2020;1(2):15-20.

https://dx.doi.org/10.47890/JEAEE/2020/YongzhengZhang/11120008

Received Date: October 25, 2020; Accepted Date: December 10, 2020; Published Date: December 15, 2020

## Abstract

An algorithm is proposed in this paper for calculating the impulse earthing resistances of vertical earthing electrodes. The proposed algorithm employs the average potential method to derive the formula of the low current earthing resistance. Unlike the previous algorithm, the soil ionization effect under high impulse current is taken into account by introducing a nonlinear characteristic to represent the relationship between the electric field and current density in the ionization zone around the earthing electrode. On the basis of the nonlinear characteristic, the effective radius is evaluated for the equivalent earthing electrode. Then, the impulse earthing resistance can be calculated by substituting the effective radius into the formula of the low current earthing resistance. A comparison is also made between calculated and measured results to confirm the validity of the proposed algorithm.

Keywords: Earthing Resistance; Vertical Earthing Electrode; Average Potential Integral; Soil Ionization; Current Density;

## Introduction

The impulse earthing resistance of earthing electrodes is a very important factor which has to be taken into consideration in the lightning protection design of civil buildings and electric substations. An appropriate choice of protection measures against lightning overvoltage depends to a large degree on the knowledge of the values of the impulse earthing resistances of earthing electrodes. As the simplest form of earthing electrodes, vertical erathing electrodes are widely used in practical earthing systems. The impulse resistances of vertical earthing electrodes have been investigated for many years. Although a few algorithms were presented by different authors [1−3], their modeling for soil ionization effect under high impulse current is still a problem. These previous algorithms usually utilized a linear characteristic to represent the relationship between the electric field and current density in the ionization zone around the earthing electrode. In fact, this relationship has pronounced nonlinearity for typical kinds of soils in terms of the experimental investigation [4]. The nonlinear characteristic should therefore betaken into account in order to perform a more accurate calculation of the impulse earthing resistance.

The aim of this paper is to propose an efficient algorithm for calcu-lating the impulse resistances of vertical earthing electrodes. For the sake of engineering application, the inhomogeneity of soil is neglected. In the algorithm, the vertical earthing electrode is divided into a series of segments and the average potential method is employed to calculate the resistance matrix. The formula for the low current earthing resistance is derived from the resistance matrix. The soil ionization effect under high impulse current is further considered by using a nonlinear characteristic to represent the relationship between the electric field and current density in the ionization zone. With a simplified treatment made for the ionization zone, the effective radius of the equivalent earthing electrode is evaluated from the nonlinear characteristic. The impulse earthing resistance is then obtained by substituting the effective radius into the formula of the low current earthing resistance. Validity of the proposed algorithm has been verified by comparing calculated and measured results.

## Low Current Earthing resistance

Consider a vertical earthing electrode in a homogeneous soil with the resistivity ρ, as shown in Figure 1. When a low current I diffuses from it into the soil, the soil ionization cannot be caused. According to the electromagnetic field theory, the presence of the interface between

Figure 1: Vertical earthing

air and soil can be taken into account by installing the image electrode[5], as illustrated in Figure 2. The image electrode depicted by the dotted line is installed above the earth surface; meanwhile the actual and image electrodes are symmetrical about the earth surface. In consideration of

Figure 2: Division of actual and image electrodes into N segments.

the diffused current distribution along the electrode length, the actual electrode and its image are divided into N segments, respectively (see Figure 2）. At an arbitrary point on the surface of jth actual segment, as shown in Figure 3, the potential generated by the currents of jth actual segment and its image segment (j′ th segment) is calculated as [5]

where τj and τj are the linear current densities of jth segment and its image, respectively. As the segment is rather short, the diffused current distribution on the segment is approximately considered to be uniform. Letting Ij and Ij denote the diffused currents of jth segment and its image, respectively,

Figure 3: Sketch for calculating self resistance of jth segment.

τj and τj become τj= Ij/(z2-z1) and τj= Ij/(z2-z1). Owing to the image symmetry, we have Ij= Ij. Thus, (2.1) is rewritten as

By taking the integral average of φj(zs) over the length of jth segment, the average potential is obtained as

where

The self resistance of jth segment is given as

In a similar manner, the average potential generated by the currents of kth segment and its image on the surface of jth segment, as shown in Figure 4, is found as

where Ik is the current of kth segment. The mutual resistance between jth and kth segments is given as

Figure 4: Sketch for calculating mutual resistance between jth and kth segments.

By taking the self and mutual resistances as the matrix elements, the resistance matrix of the vertical earthing electrode is formed as

Using R gives the relationship between voltage and current

where Uj (j=1, 2, ..., N) is the voltage of jth segment. As the conducing surfaces are equipotential, we have U1=U2= ··· =UN=U [5]. Calculation of the inverse matrix R−1 gives the following expression

where Gjk (j, k=1, 2, ..., N) is the element of inverse matrix R−1. In terms of (9), the total current is found as

As a result, the earthing resistance can be obtained as

As shown in Figure 1, a vertical earthing electrode of r0=0.015m is considered here. The earthing electrode is made of plain carbon steel (con-ductivity σ=0.5×107 S•m-1). The soil is kind of clay. Its resistivity is ρ=42Ω·m and water content is about 20%. The earthing resistances calculated from (11) are given in Figure 5. The corresponding measured values are given simultaneously in Figure 5, which was obtained by the fall of potential method [6]. As indicated in Figure 5, the calculated values are close to the measured ones.

Figure 5: Calculated and measured values of low current earthing resistance.

## Impulse Earthing Resistance

For a high impulse current with crest value I, representative of lightning, the impulse resistance Rie of the earthing electrode is defined as a ratio of the crest value Um of the impulse potential on the earthing electrode to I, i.e. Rie=Um/I [7]. The impulse current can produces great current density and high electric field intensity near the earthing electrode. When the electric field intensity on the surface of the earthing electrode exceeds the critical value Ec of soil ionization gradient, the breakdown will occur. This process can be illustrated by Figure 6 [8].

Figure 6: Impulse breakdown of soil around a vertical earthing electrode, 1—arc zone, 2—streamer zone, 3—semiconductive zone, 4—constant resistivity zone.

As the current increases, streamers are developed and in turn arcs are generated. Within the streamer and arc zones, the resistivity decreases from its original value to a limit approaching conductor [8]. In addition, there is a semiconductive zone between the streamer zone and the constant resistivity zone. For the purpose of simplifying calculation, this process can be described by an equivalent model shown in Figure 7. In the equivalent model, the semiconductive zone is neglected since it is small.

Figure 7: Ionization zone, (a) non-uniform along electrode length, (b) uniform along electrode length.

The streamer and arc zones are equivalently modeled as an ionization zone [8,9]. As the vertical earthing electrodes used in the actual cases usually have a shorter length (less than 5m), the non-uniformity of the ionization zone along the electrode length, as shown in Figure 7(a), is not pronounced and may be neglected. Therefore, the ionization zone is approximately considered to be uniform along the electrode length, as shown in Figure 7(b) [8−10]. The border of theionizationzoneis delimitedbythecriticalelectric field value Ec. A nonlinear characteris-tic is introduced to represent the relationship between the electric field andcurrent density in theionizationzone. The nonlinear characteristic is given as [4]

where a and b are constants and J is the current density

Substituting (12) into (13), the electric field intensity is rewritten as

On the border of the ionization zone, E should satisfy the boundary condition.

The values of a, b and Ec can be found in [4] for typical kinds of soils. Thus, the boundary radius of the ionization zone can be derived from (14) and (15)

where

The soil ionization is basically equivalent to an increase in the dimension of the earthing electrode, which can be taken into account by the effective radius re of the equivalent earthing electrode in constant resistivity zone, asshowninFigure 8. The voltage between the earthing electrode and the border of the ionization zone, as shown in Figure 8(a), is expressed as

where

On the other hand, the voltage between the equivalent earthing electrode and the border of the ionization zone, as shown in Figure 8(b), is expressed as

Figure 8: Sketch for evaluation of the effective radius, 1—ionization zone, 2—constant resistivity zone, (a) ionization zone around the actual earthing electrode, (b) constant resistivity zone around the equivalent earthing electrode.

Because Ui and Ue must be equal, the effective radius of the equivalent eathing electrode can be determined from (18) and (20)

where λ=exp(χ) and

As the integration for (19) is difficult to be evaluated analytically, a numerical solution is given to obtain W

where Δr is the radial step and N=(ri-r0)/Δr.

By replacing r0 with re in (11), the impulse earthing resistance can be obtained. To check the validity of the proposed algorithm, two numerical examples are given here. The data for the first example are: r0=0.0127m, l=3.05m, ρ=87.2Ω·m, Ec=127kV/m, a=3094.6 and b=0.51 [4]. At different current crest values,the calculated impulse earthing resistances are shown in Figure 9.

Figure 9: Calculated and measured impulse earthing resistances.

The data for the second example are: r0=0.025m, l=1m, ρ=43.5Ω·m, Ec=350kV/m, a=219.05 and b=0.82 [4]. The calculated crest potentials on the earthing electrode (Um=RieI) are shown in Figure 10. Furthermore, the corresponding measured results [11,12] are given in Figures 9 and 10, respectively, for comparison. It can be seen from Figures 9 and 10 that a better agreement appears between calculated and measured results.

Figure 10: Calculated and measured crest potentials.

## Conclusions

An algorithm for calculating the impulse earthing resistances of vertical earthing electrodes has been proposed in this paper. It uses the average potential method to derive the formula of the low current earthing resistances. In the ionization zone, a nonlinear characteristic is intro-duced in representingtherelationship between electric field and current density and has the capability of giving an appropriate consideration for the soil ionization effect under high impulse current. On the basis of the nonlinear characteristic, a simplified procedure has been developed to evaluate the effective radius of the equivalent earthing electrode. The impulse earthing resistance can then be obtained by substituting the effective radius into the formula of the low current earthing resistance. The calculated results are compared whith the measured ones and an agreement appears between them, which shows the applicability of the proposed algorithm in practical lightning protection design.

## References

1. Salgado LA, Guardado JL, Torres J and Hernandez EO. Transient analysis of grounding systems under lightning strikes considering soil ionization. IEEE Industry Applications Society Annual Meeting (IAS) .2010;1-7. doi: 10.1109/IAS.2010.5614431
2. Cooray V. Lightning protection. IET power and energy series. 2010.
3. Geri A, Behavior of grounding systems exited by high impulse currents: the model and its validation. IEEE Transactions on Power Delivery. 1999;14(3):1008–1017. doi: 10.1109/61.772347
4. Annenkov AZ. The Volt–Ampere characteristic of soil under lightning surge currents. Electrichestvo, No. 2, pp.57–60, 1988.
5. Chow YL, Elsherbiny MM and Salama MMA. Efficient computation of rodbed grounding resistance in a homogeneous earth by Galerkin’s moment method. IEE Proceedings–Generation. Transmission and  Distribution. 1995;142(6):653–660. doi: 10.1049/ip-gtd:19952103
6. Chinese national standard GB /T 17949.1-2000. Guide for measuring resistivity, grounding impedance and earth surface potentials of a ground system. China Planning Press. Beijing. 2001.
7. Sekioka S, Sonoda T and Ametani A. Experimental study of current-dependent grounding resistance of rod electrode.  IEEE Transactions on Power Delivery. 2005;20(2):1569–1576. doi: 10.1109/TPWRD.2004.838660
8. Hileman AR. Insulation coordination for power systems. New York, Marcel Dekker Inc. 1999
9. Annenkov AZ. A method for calculating the impulse earthing resisitance of a grounding rod in sandy soil. Electrichestvo. 1996;14–21.
10. Mohamad NNor. Simulation of earthing systems under high impulse. 8th International Conference on Properties and applications of Dielectric Materials. 2006;347–352. doi: 10.1109/ICPADM.2006.284187
11. Nabkova EI. Investigation on the impulse earthing characteristics based on experimental data. Electrichestvo. 1984;46–49.
12. Bellaschi PL, Armington RE and Snowden AE. Impulse and 60 cycle characteristics of driven grounds–Ⅱ. Transactions of American Institute of Electrical Engineers. 1942;61(6):349–363. doi: 10.1109/T-AIEE.1942.5058539