Proceedings of the 14th International Middle East Power Systems Conference (MEPCON’10), Cairo University, Egypt, December 19-21,

2010, Paper ID 229.

Earth Surface Potentials and GPR for Grids Having Different

Conductors' Distributions Buried in different Layer Soil Structures
Prof.Dr.Osama Elsayed Gouda Eng Hanan Mohamed Ahmed Ibrahim
Department of Electrical Power and Machines National Energy Control Center
Faculty of Engineering, Cairo University Egyptian Electricity Transmission Company
[email protected] [email protected]

Abstract - This paper presents simplified approaches for the

convenient design of grounding grid configurations buried in Lik
uniform or non-uniform soil. The parallel grid conductors with S ik = × 100% . (1)
different spaces are distributed using these approaches according
L
to soil parameters which grid is buried in, hence, the equalization
of current density and then decreasing of GPR are achieved. The
image method is suggested to calculate the GPR and the earth
S ik = b1e − b 2 i + b3 . (2)
surface potentials with accurate results for the grids buried in where Sik is the grid percentage in each segment. The
uniform or non-uniform soil structure. Charts of voltage profiles coefficients b1, b2 and b3 are dependent on number of
have been designed in 2D and 3D to make the selection of suitable segments of the conductor [5].
earthing system easier.
B. OPTIMAL COMPRESSION RATIO (OCR)
Index Terms – Touch voltage, earth surface potential,
The optimal design of grounding systems for substations
unequally spaces, conductor distribution.
is to suitably arrange the conductors of grounding systems to
equalize the leakage current distribution and the potential
I. INTRODUCTION
gradient of ground surface and at the same time to keep the
Since the optimal design of the grounding system maximum touch voltages within the safe limit; this would
conditional with the safety criteria for the substations and ensure making all grounding conductors sufficiently utilized
persons, thus, the safety of personnel in and around electric so certified as a safe and economic design method.
power installations has been a prime concern. The safety of The conductor compression ratio “C” is related to the grid
persons is affected by the ground potential rise of grounded dimensions and the maximum conductor separation dmax which
structures during electric power faults and then subjecting him occurs toward the center of the grid and is given by [6]:
or herself to voltages. So, in order to limit these voltages to L(1 − C )
safe values, it is necessary to design a right grounding system d max = N even. (3)
1 + C − 2C ( N / 2+1)
that effectively connects all metallic structures of the
L(1 − C )
substations to ensure the safety criteria for both human and d max = N odd. (4)
substations taking into account the economical purposes. 2(1 − C ( N −1) / 2 )
However, the unequally spaced grid can make current where N is the desired number of conductors perpendicular to
density which is emanating from grid conductors more a given grid axis of length L.
uniform, and therefore remarkably decrease potential When the grounding conductors are arranged according
gradients on the earth’s surface [1-3]. Then, the maximum to an exponent regularity, the conductors span decreases
touch voltage for this grid will reach its lowest values and thus gradually from the center to the side of the grounding grid.
the safety conditions for human above the substation earth The ith conductor span from the center is [6, 7]
surface will be ensured.
d i = d max C i . (5)

II. METHOLOTOGY (UNEQUALLY SPACES GRIDS DESIGN

For the optimal design of grounding system, the optimal
PROCEDURES)
compression ratio could be determined, which is defined as
A. UNEQUALLY SPACED GROUNDING GRIDS PROCEDURE the compression ratio, where the touch voltage reaches its
To get the best voltage gradient distribution, the minimum if the grounding grid is designed under this
unequally grid conductor spacing must be increase compression ratio [6]. An empirical expression is obtained to
progressively from the periphery to the center. To get the best calculate OCR as follows [7, 8]:
distribution for this conductor spacing, if a conductor spacing
is divided by k segments, the length of the ith segment, Lik can OCR=a0+a1exp (0.0001h) +a2exp (bh). (6)
be expressed in percent of the total length of the conductor L where h is the grid burial depth in m and
[4, 5]: b=-0.3503- 9.6311exp (-0.03666L). (7)

544
 a0=a01+a02k+a03k2. (8) cases is given in figure 1. Current emanating from the grid
a1=a11+a12k+a13k2. (9) into the soil is I=10KA. In this study the soil ionization effects
a2=a21+a22k. (10) are ignored. By applying image technique [8, 9], the potential
K is number of segments of the conductor in any differences in volts and as a percentage of GPR for all grid
direction. The relationship between a0, a1... and L is cleared in modules according to number of parallel conductors are
[8]. calculated, and saved in table I for the grid module of (A=105
x 105 m2, N=8 x 8) as an example.

III. EFFECT OF CONDUCTORS' DISTRIBUTIONS ON VOLTAGE

PROFILE IN CASE OF GRIDS BURIED IN UNIFORM SOIL

A. Using Unequally Spaced Grounding Grids Procedure

To illustrate the effect of equally and unequally grid

conductor distribution on the earth surface potentials above
the grid, a square grid of 105 x 105 m2 with 9, 13 and 17
equally and unequally spaces parallel conductors in each Fig. 1 An example for grid module of N=64 with equally, unequally spaced
having denser conductors at the edges and the middle respectively.
direction is taken as an example. An example of the studied

TABLE I
EARTH SURFACE POTENTIALS IN VOLTS AND AS PERCENTAGES OF GPR ALONG CENTER CONDUCTOR
FOR EQUALLY AND UNEQUALLY SPACED GRID OF (A=105 x 105 M2, N=8 x 8)

Equally spaces grid conductors

x, (m) -52.5 -39.38 -26.25 -13.13 0 13.1 26.3 39.4 52.5
V, (volt) 4003 4732 5076 5249 5297 5249 5076 4732 4003
V, ( % GPR) 0.76 0.89 0.96 0.99 1 0.99 0.96 0.89 0.76
Unequally spaces having denser conductors at the grid edges
x, (m) -52.4 -46.8 -35.1 -19.1 0.0 19.1 35.1 46.8 52.4
V, (volt) 4080 4431 4629 4762 4799 4762 4629 4431 4080
V, ( % GPR) 0.85 0.92 0.96 0.99 1 0.99 0.96 0.92 0.85
Unequally spaces having denser conductors at the grid middle
x, (m) -52.4 -33.3 -17.3 -5.6 0.0 5.6 17.3 33.3 52.4
V, (volt) 4156 5265 5785 6138 6197 6138 5785 5265 4156
V, ( % GPR) 0.66 0.85 0.93 0.99 1 0.99 0.93 0.85 0.66

From the obtained results, it's concluded that as Noteworthy from the table, and to confirm the previous
increasing of equally spacing parallel grid conductors conclusion, that the increasing of unequally grid conductors
included in the same area, the GPR for the grid decreases and having denser conductors at the grid middle increases the
the earth surface potentials for grid increases, and then the GPR for this grid.
potential difference between each two conductors will Colored graphical representations for earth surface
decrease and the voltage profile will be more uniform. potentials above the mentioned grid in all conductors'
Regarding to grid conductors' distribution with unequally distribution cases for every one meter in both x and y
spaces, the percentages of the potential differences between direction and saved in 3D contours plots as shown in figures 2
center and corner conductors for this grid of (A=105 x 105 to 4.
m2) whatever number of conductors are around 26 %, 37% The potential difference between center and corner
and 15% in cases of equally spaces, unequally spaces having conductors is an important factor to determine the most
denser conductors at the middle and edges respectively. So, efficient design of grounding system, because of this
and as a conclusion, the conductor distribution with unequally difference representing the maximum touch voltage for the
spaces having denser conductors at the edges provides the grid. Figures 5 to 7 illustrate the earth surface potentials and
most efficient design because it satisfies the lowest GPR and touch voltage along the center conductor of the mentioned
the lowest potential differences between center and corner grid modules in different cases of conductors' distribution.
conductors comparing with the results that occurred in other From the figures, and for all grid modules whatever number of
cases of conductors' distribution included with the same area, conductors, it is cleared that the maximum touch voltage
thus the voltage distribution will be more uniform. along the center conductor reaches its lowest value in case of
unequally spaced grid having denser conductors at the edges.

545
Fig. 2 Voltage profile along all grid conductors in case of equally spaced grid. Fig. 5 Earth surface potentials and touch voltages over the center conductor
Number of meshes =64 for mentioned grid having N=64

Fig. 6 Earth surface potentials and touch voltages over the center conductor
Fig. 3 Voltage profile along all grid conductors in case of unequally spaced
for mentioned grid having N=144
grid having denser conductors at the middle. Number of meshes =64

Fig. 4 Voltage profile along all grid conductors in case of unequally spaced Fig. 7 Earth surface potentials and touch voltages over the center conductor
grid having denser conductors at the edges. Number of meshes = 64 for mentioned grid having N=256.

546
B. By Using Optimal Compression Ratio (OCR)

To illustrate the effect of unequally conductors'

distribution by using OCR on the earth surface potentials and
touch voltage profiles for the grid buried in uniform soil and
to satisfy the optimal grounding grid design, the grid module
of (A=105 x 105 m2, N=8 x 8) as an example, buried in a
uniform soil of 100 ohm.m resistivity, the earth surface
potentials and touch voltages are calculated along the center
conductor for this grid under different cases of compression
ratios by applying equations (3) to (5) as presented in the
figure 8.
By the analysis of calculated data represented in the
figures, it is stated that the maximum touch voltage along the
center conductor reaches its lowest value when the
compression ratio for this grid reaches 0.538, it means and as Fig. 9 Maximum touch voltage over center conductor for mentioned grid as a
a conclusion, at the optimal compression ratio obtained by function of conductor compression ratio.
applying equations (6) to (10), the best arrangement of
grounding conductor span is achieved, so the minimum touch
voltage and thus, the optimal design for grounding grid are
satisfied. This result is influenced if compared with the
different cases of conductors' distribution mentioned in the
previous section for the same grid module of (A=105 x 105
m2, N=8 x 8).

Fig. 10 Earth surface potentials over all grid conductors for grid of (A=105 x
105 m2, N=8 x 8) in case of C=0.

Fig. 8 Earth surface potentials and touch voltages over the center conductor
for mentioned grid having different cases of compression ratios.

Figure 9 illustrate the relationship between the maximum

touch voltages over center conductor for the mentioned
unequally spaced grid as a function of conductor compression
ratio. Remarkable, this relationship usually has an obvious ''U''
shape.

Colored graphical representations for earth surface

potentials above the mentioned grid in different cases of
compression ratio for every one meter in both x and y
direction are saved in 3D contours plots as shown in figures Fig. 11 Earth surface potentials over all grid conductors for grid of (A=105 x
105 m2, N=8 x 8) in case of C=0.2.
10 to 14. The difference between figures 10 and 11 on one
side and figures 12 to 14 on the other side is referred to the
wide range compression ratio differences as shown in figures.

547
 IV. EFFECT OF CONDUCTORS' DISTRIBUTIONS ON VOLTAGE
PROFILE IN CASE OF GRIDS BURIED IN DOUBLE LAYER SOIL

To illustrate the effect of changing the first layer depth of

double layer soil structure with variation of the upper and
lower soil resistivities on the voltage profile, 2D graphically
representation for earth surface potentials is plotted as in
figure 15 along the center conductor for a 105 x 105 m2 grid
area with 9 equally spaces conductors in each direction buried
in non-uniform soil having specifications that given in table
II. The current emanating for all cases is 10,000 Amp.

TABLE II
SPECIFICATIONS OF MENTIONED NON UNIFORM SOIL

ρ1 100 100 100 100 100

ρ2 100 40 40 250 250
Fig. 12 Earth surface potentials over all grid conductors for grid of (A=105 x H infinity 10 20 10 20
105 m2, N=8 x 8) in case of C=0.6.
From the figure, it is cleared that the ground potential rise
(GPR) for the grid buried in the high resistivity upper layer
soil will be lower than that resulted above the grid if it is
buried in the upper layer soil of low resistivity. Also, the GPR
increases as the depth of the high resistivity upper layer soil
increases and the opposite results in case of increasing the
depth of the upper layer of low resistivity.

Fig. 13 Earth surface potentials over all grid conductors for grid of (A=105 x
105 m2, N=8 x 8) in case of C=1.

Fig. 15 Earth surface potentials along center conductor for equally spaced grid
(A= 105 x 105 m2, N=8 x 8) with various soil structures.

To illustrate the effect of soil non-uniformity on the

grounding grid system having different conductors'
distribution by using OCR, table III shows the earth surface
potentials and maximum touch voltages along center
conductor for the mentioned grid of (A=105 x 105 m2, N=8 x
8) buried in upper layer of double layer soil having the same
parameters structure mentioned above.
Figure 16 clears that the compression ratio which
Fig. 14 Earth surface potentials over all grid conductors for grid of (A=105 x determine the arrangement of grid conductors satisfies the
105 m2, N=8 x 8) in case of C=0.54. lowest touch voltage for this grid according to the soil
parameters in which grid is buried in. From the figure, it's
concluded that in cases of (ρ1<ρ2), the compression ratio that

548
satisfies the lowest maximum touch voltage will takes values 0.4, it means that the grid conductors will be denser at the
less than the optimal compression ratio (OCR) determined by edges and having widely conductor spans at the middle. The
equation (6) for this grid. The results will be in the contrary of spaces between conductors will increase progressively from
this in cases of (ρ1>ρ2). the periphery to the center and the middle spans will decrease
For the mentioned grid buried in double layer soil gradually in cases of (ρ1>ρ2) where the maximum touch
structures, in cases of (ρ1<ρ2) and since the maximum touch voltage reaches its lowest values at the compression ratio of
voltage reaches its lowest values at the compression ratio of 0.6.

TABLE III
EARTH SURFACE POTENTIALS AND MAXIMUM TOUCH VOLTAGES IN VOLTS ALONG CENTER CONDUCTOR FOR MENTIONED GRID
HAVING DIFFERENT CASES OF CONDUCTORS' DISTRIBUTION BURIED IN DIFFERENT CASES OF SOIL STRUCURE BY USING OCR

uniform soil ρ2=40, H=10m ρ2=40, H=20m ρ2=250, H=10m ρ2=250, H=20m
max max max max max
Type of unequally grid GPR touch GPR touch GPR touch GPR touch GPR touch
voltage voltage voltage voltage voltage
equally spaced grid, means C=1 5297 1346 2925 779 3435 982 9956 2302 8582 1851
unequally spaced grid with C=0.8 4933 982 2708 549 3147 681 9404 1778 8146 1435
unequally spaced grid with C=0.6 4603 582 2597 432 2978 445 8851 1186 7718 962
unequally spaced grid with OCR=0.54 4544 479 2584 483 2955 506 8688 989 7593 802
unequally spaced grid with C=0.4 4489 563 2655 680 2987 704 8346 518 7331 409
unequally spaced grid with C=0.2 4670 988 2898 1053 3233 1125 8284 777 7360 776
unequally spaced grid with C=0 5140 1608 3395 1624 3730 1721 8682 1453 7783 1417

• The conductors' distribution with unequally spaces

having denser conductors at the edges provides the
most efficient design.
• As a conclusion, for grids buried in double layer soil,
the compression ratio that satisfies the lowest touch
voltage will differs according to soil parameters.
REFERENCES
[1] A. Puttarach N.Chakpitak T. Kasirawat and C. Pongsriwat, Department of
Electrical Engineering, ChiangMai University, “Ground Layer Depth and
the Effect of GPR”, IEEE PES Power Africa 2007 Conference and
Exposition, July 2007.
[2] Yves Rajotte, Jean De Steve, Jaues Fortin, Richard Lehoux, Georges
Simard, “Earth Potential Rise Influence near HV Substations in Rural
Areas”, CIRED (Turin) 2005. June 2005.
[3] BHARAT SANCHAR NIGAM LIMITED, A Govt. of India Enterprise,
''Measurement Earth Potential Rise (EPR) Zone'', No. EPR/A-002, CGM
and Technical and Development Circle, Issue No. 01, No of Pages. 06,
Fig. 16 Maximum touch voltage over center conductor for mentioned grid 18-03-2005.
buried in various soil structures as a function of compression ratio. [4] B. Phithakwong, N. Kraisnachinda, S. Banjongjit, C. Chompoo-Inwai, M.
Kando, “New Techniques the Computer-Aided Design for Substation
Grounding ”, Power Engineering Society Winter Meeting, 2000. IEEE,
V. CONCLUSION Vol. 3, pp. 2011-2015, Jan. 2000.
[5] L. Huang, X. Chen, and H. Yan, “Study of Unequally Spaced Grounding
The following points are briefly analyzing the main Grids,” IEEE Trans. Power Del., Vol. 10, No. 2, pp. 716–722, Apr. 1995.
conclusions of this paper: [6] J. He, R. Zeng, Y. Gao, Y. Tu, W. Sun, J. Zou, and Z. Guan, “Optimal
Design of Grounding System Considering the Influence of Seasonal
• As increasing of equally spacing parallel grid
Frozen Soil Layer”, IEEE Trans. Power Del., vol. 20, no. 1, Jan. 2005.
conductors included in the same area, the GPR for [7] W.Sun, JinglingHe, Yanqing Gao, R.Zeng, W.Wu and Qi Su: “Optimal
the grid decreases and the minimum earth surface Design Analysis of Grounding Grids for Substations Built in Non-uniform
potentials for grid increases, and then the potential Soil”, IEEE, Trans. Power apparatus and system, 2000, pp.1455-1460.
difference between each two conductors will [8] A. Thabet: “Grounding Systems of Electric Substations in Non-uniform
Earth Structure with New Analysis”, M. Sc. Thesis, High Institute of
decrease, so, the voltage profile will be more Energy, Aswan, 2002.
uniform. [9] E.A. Ibrahim: “Grounding of the Nuclear Power Plants Using Grounding
Grids”, Ph. D. thesis, Faculty of Engineering, Cairo University, 1993.

549