Grounding grid sizing: earth resistance, GPR, step and touch voltage
Sizing a grounding grid crosses soil resistivity, fault current and grid geometry to estimate the earth resistance, the ground potential rise (GPR) and to compare it with the step and touch voltages the human body can tolerate (IEEE 80 / NBR 15751).
When to use
Use it whenever you need to design or verify the ground grid of a substation, machine room, medium-voltage switchgear room or industrial plant. The method sets the area and length of buried conductor, estimates the earth resistance by the Sverak or Laurent methods, computes the GPR (ground potential rise during the fault) and compares the result with the maximum step and touch voltages the human body can withstand for the fault clearing time. It is the first step of every grounding project: it shows whether the grid protects people against shock, sizes the conductor by the adiabatic criterion and anticipates the need for a gravel layer before construction. It also serves to diagnose existing grids that fail to reach the design target resistance.
What grounding grid sizing is
Sizing a grounding grid is not simply burying copper conductors and driving rods: it is ensuring that, at the instant of an earth fault, no person inside or around the installation is exposed to a dangerous voltage. The design crosses three groups of data — the soil resistivity, the fault current and the grid geometry — to answer two distinct questions: what earth resistance the grid achieves, and whether the voltages that arise in the soil during the fault stay below what the human body can tolerate.
The most common conceptual mistake is to treat grounding as merely a “lower the ohms” problem. A grid can have an extremely low resistance and still be dangerous, because safety depends on the surface potential gradients, not only on the total resistance. That is why the IEEE 80 and NBR 15751 method clearly separates the resistance calculation from the safety criterion.
Geometry and the buried conductor length
It all starts with the grid geometry. The grid is a mesh of horizontal conductors of sides Lx by Ly, with spacing D between parallel cables, supplemented by vertical rods. The number of conductors in each direction is:
n_x = ⌊Lx/D⌋ + 1 ; n_y = ⌊Ly/D⌋ + 1
and the total buried length sums the mesh plus the rods:
L = (n_x·Ly + n_y·Lx) + n_rods·L_rod
This length L is decisive: the more buried conductor, the lower the resistance and the more uniform the potential. The area A = Lx·Ly defines the equivalent radius r = √(A/π), which appears in the resistance floor imposed by the soil.
Grid resistance: Sverak and Laurent
The tool computes the resistance by two methods. The Laurent-Niemann one is the simplified estimate:
R = ρ/(4·r) + ρ/L
The first term, ρ/(4r), is the resistance floor — the minimum a plate of area A would reach, regardless of how much cable is buried. The second term decreases with length. Laurent is fast and conservative, but it ignores the burial depth.
The Sverak method, adopted by IEEE 80, corrects this by capturing the depth h:
R = ρ·[ 1/L + 1/√(20·A) · (1 + 1/(1 + h·√(20/A))) ]
It is the expression recommended for design, because a grid buried deeper has lower resistance — an effect Laurent does not see. In both methods, the soil resistivity ρ is the dominant variable: doubling ρ almost doubles the resistance.
Ground potential rise (GPR)
Once the resistance is known, the GPR (Ground Potential Rise) is computed — the potential the whole grid reaches relative to remote earth during the fault:
GPR = I_G · R , with I_G = Df · Sf · I
The grid design current I_G is not necessarily the symmetrical short circuit I. The decrement factor Df (IEEE 80 §15) corrects the asymmetry in the first cycles, and the division factor Sf represents the fraction of the short circuit that actually returns through the grid (part flows via shield wires and neutrals). By default Df = Sf = 1 is adopted, throwing the whole current into the grid — the most conservative assumption.
The GPR can reach thousands of volts. It is not, by itself, the failure criterion — but it is the starting point for comparison with the safety limits.
Conductor section by the adiabatic criterion
The grid conductor must withstand the fault current thermally without melting. By the IEEE 80 adiabatic criterion (all the heat stays in the copper during the short circuit):
A_mm² = I_kA · Kf · √t · 0.5067
where Kf = 7.06 for hard-drawn copper and the factor 0.5067 converts kcmil to mm². A 50 mm² copper floor is also applied over the thermal value, required for mechanical robustness and corrosion resistance, since the conductor will stay buried for 30 to 40 years. The adopted section is the larger of the two.
The safety criterion: step and touch voltages
Here lies the heart of the design. IEEE 80 defines the maximum tolerable voltages of the human body as a function of the exposure time:
E_touch = (1000 + 1.5·Cs·ρs) · k/√t
E_step = (1000 + 6·Cs·ρs) · k/√t
The touch voltage (between the hand on a grounded structure and the feet) is the most critical limit, almost always lower than the step voltage. The coefficient k is 0.116 for a 50 kg body and 0.157 for 70 kg, and t is the fault clearing time — the longer it is, the lower the tolerable voltage.
The factor Cs represents the effect of the surface gravel layer:
Cs = 1 − 0.09·(1 − ρ/ρs)/(2·hs + 0.09)
With no gravel, Cs = 1. With a high-resistivity layer (ρs ≈ 2,000–5,000 Ω·m), Cs drops below 1, which raises the tolerable voltages — often exactly what makes a grid pass the safety criterion.
How to interpret the result
The calculation ends with a screening: it compares the GPR with the tolerable touch voltage. If GPR ≤ E_touch, the grid is classified as ok. If GPR > E_touch, the result is analyse mesh — a sign that, despite the low resistance, the potentials may be dangerous and the design needs refinement. The typical fixes are:
- Add or thicken the gravel layer (raises Cs and the tolerable voltages);
- Densify the grid (a smaller spacing D reduces the local mesh voltage);
- Increase the area or the conductor length (lowers the resistance and the GPR);
- Reduce the fault clearing time with faster protection.
Practical design considerations
- Measure the resistivity in the field (Wenner method, IEEE 81) rather than adopting a table value; the soil is stratified and ρ varies with the season.
- Do not rely on the ohms alone: the final approval is the voltage criterion, not the resistance in isolation.
- Consider the gravel layer from the start: it is usually the cheapest way to pass the touch and step voltages.
- Align standard and method: IEEE 80 and NBR 15751 govern the resistance calculation (Sverak), the GPR, the step and touch voltages and the conductor adiabatic criterion; IEEE 81 guides the resistivity measurement.
Following this chain — geometry, conductor length, Sverak resistance, GPR, adiabatic section and verification of the step and touch voltages — yields a grounding grid that protects people and agrees with the worked examples of IEEE 80.
Formulas and fundamentals
R = ρ·[ 1/L + 1/√(20·A) · (1 + 1/(1 + h·√(20/A))) ] Earth resistance of a grid buried at depth h. ρ is the soil resistivity [Ω·m], L the total length of buried conductor (grid + rods) [m], A the area occupied by the grid [m²] and h the burial depth [m]. This is the Sverak expression adopted in IEEE 80; it captures the depth effect that Laurent ignores.
R = ρ/(4·r) + ρ/L Simplified estimate (Laurent-Niemann). ρ is the resistivity [Ω·m], r the equivalent radius of the area √(A/π) [m] and L the total conductor length [m]. The first term is the floor imposed by the area (equivalent plate); the second decreases with the buried cable length.
L = (n_x·Ly + n_y·Lx) + n_rods·L_rod ; n_x = ⌊Lx/D⌋+1 ; n_y = ⌊Ly/D⌋+1 Total buried length. Lx and Ly are the grid sides [m], D the spacing between parallel conductors [m], n_x and n_y the number of conductors in each direction, and n_rods·L_rod the total of vertical rods [m]. The more buried conductor, the lower the resistance.
GPR = I_G·R ; I_G = Df·Sf·I Potential the grid reaches relative to remote earth during the fault. I is the single-line-to-ground fault current [A], R the grid resistance [Ω], Df the decrement factor (asymmetry, IEEE 80 §15) and Sf the current division factor. By default Df = Sf = 1 (the whole current flows into the grid — conservative).
A_mm² = max( I_kA·Kf·√t · 0.5067 ; 50 ) Minimum grid conductor section. I_kA is the fault current [kA], Kf the material constant (7.06 for hard-drawn copper), t the fault clearing time [s] and 0.5067 converts kcmil to mm². A 50 mm² copper floor is also imposed for mechanical robustness and corrosion over a 30–40 year burial.
E_step = (1000 + 6·Cs·ρs)·k/√t ; E_touch = (1000 + 1.5·Cs·ρs)·k/√t Safety limits for the human body. ρs is the resistivity of the surface layer (gravel) [Ω·m], Cs the surface derating factor (=1 with no layer), t the fault time [s] and k = 0.116 for a 50 kg body or 0.157 for 70 kg. Cs = 1 − 0.09·(1 − ρ/ρs)/(2·hs + 0.09), with hs the gravel thickness [m].
Standards & methods
- IEEE Std 80 — Guide for Safety in AC Substation Grounding (Sverak method, step and touch voltages, conductor adiabatic criterion)
- ABNT NBR 15751 — Substation grounding systems — Requirements
- ABNT NBR 5419 — Protection against lightning (LPS grounding)
- IEEE Std 81 — Guide for Measuring Earth Resistivity, Ground Impedance and Earth Surface Potentials
- IEC 60364-5-54 — Low-voltage electrical installations: earthing arrangements and protective conductors
Typical reference values
| Quantity | Typical range | Note |
|---|---|---|
| Soil resistivity (ρ) | 10 to 1000 Ω·m | Moist clay ≈ 50–100 Ω·m; dry sandy/rocky soil > 500 Ω·m. |
| Gravel resistivity (ρs) | 2,000 to 5,000 Ω·m | A 0.10–0.15 m layer over the soil greatly raises the tolerable voltage (Cs). |
| Target grid resistance | ≤ 1 to 10 Ω | HV substations ≤ 1 Ω; common industrial installations ≤ 5–10 Ω. |
| Burial depth (h) | 0.5 to 0.8 m | Above the frost line and below the seasonal moisture variation layer. |
| Conductor spacing (D) | 3 to 10 m | Smaller spacing evens out the potential and lowers the mesh voltage. |
| Minimum conductor section | ≥ 50 mm² (copper) | Mechanical/corrosion floor; the adiabatic value may require more at high current. |
| Fault clearing time (t) | 0.1 to 1.0 s | Longer time raises the conductor section and lowers the tolerable voltage. |
Worked example
30 × 40 m substation grid in 100 Ω·m soil
Inputs
- Soil resistivity
- ρ = 100 Ω·m
- Ground fault current
- I = 10 kA
- Fault clearing time
- t = 0.5 s
- Grid dimensions
- 30 × 40 m
- Conductor spacing
- D = 5 m
- Vertical rods
- 20 × 3 m
- Method
- Sverak (no gravel) —
Results
- Grid area
- A = 1,200 m²
- Total buried length
- L ≈ 610 m
- Grid resistance (Sverak)
- R ≈ 1.42 Ω
- Conductor section
- 50 mm²
- GPR (potential rise)
- ≈ 14,157 V
- Tolerable touch voltage
- ≈ 189 V
- Tolerable step voltage
- ≈ 262 V
The geometry yields 7 conductors along 40 m and 9 along 30 m, for 550 m of grid plus 60 m of rods — L ≈ 610 m. By Sverak, R = 100·[1/610 + 1/√24000·(1 + 1/(1 + 0.5·√(20/1200)))] ≈ 1.42 Ω, close to the floor ρ/(4r) = 1.28 Ω imposed by the area. With I_G = 10 kA (Df = Sf = 1), the GPR reaches ≈ 14,157 V — far above the tolerable touch voltage of ≈ 189 V computed WITHOUT a gravel layer (Cs = 1, 50 kg body, t = 0.5 s). The screening returns 'analyse mesh': low resistance alone does not guarantee safety. The typical fix is to add a high-resistivity gravel layer (ρs ≈ 3,000 Ω·m, 0.10–0.15 m), which raises Cs and the tolerable voltages, and to densify the grid (smaller D) to reduce the local mesh voltage.
Common mistakes
- Confusing grid resistance with safety: a grid can have low R and still fail the step and touch voltages — the final criterion is GPR × tolerable voltage, not the ohms alone.
- Omitting the gravel layer: with no gravel Cs = 1 and the tolerable voltage collapses; a 0.10–0.15 m layer with high ρs is often what makes the grid pass.
- Using Laurent (no depth) in deep soil: the Sverak method is more accurate because it captures h; Laurent only serves as a quick, conservative estimate.
- Undersizing the conductor by the adiabatic criterion alone and ignoring the 50 mm² mechanical/corrosion floor for long-duration burial.
- Adopting the fault current without the decrement factor Df and the division factor Sf: the current that actually flows through the grid (I_G = Df·Sf·I) can differ greatly from the symmetrical short circuit.
- Computing a single soil resistivity value, ignoring stratification (two-layer model) and the seasonal moisture variation that changes ρ.
Frequently asked questions
What is the difference between the Sverak and Laurent methods?
Laurent-Niemann is a simplified estimate, R = ρ/(4r) + ρ/L, that depends only on the equivalent radius of the area and the cable length — it ignores the burial depth. The Sverak method, adopted in IEEE 80, adds the depth term h and is more accurate for real grids. Use Laurent as a quick, conservative check; use Sverak for the design.
Does low resistance guarantee a safe grid?
No. The resistance sets the ground potential rise (GPR = I·R), but the safety of people depends on the GPR and the soil potential gradients staying below the tolerable step and touch voltages of the body. A grid with low R can fail the voltages if the gravel layer is missing or the mesh is coarse. The final criterion is the comparison against the IEEE 80 limits.
What is the gravel layer for?
Gravel (crushed stone) is a high-resistivity surface layer (ρs ≈ 2,000–5,000 Ω·m) over the soil. It raises the contact resistance between the feet and the ground, increasing the step and touch voltages the body tolerates. It enters the calculation through the factor Cs: with no gravel Cs = 1; with gravel Cs < 1, which raises the tolerable limit and is often what makes the grid pass.
How is the conductor section sized?
By the IEEE 80 adiabatic criterion: A_mm² = I_kA·Kf·√t·0.5067, where Kf depends on the material (7.06 for hard-drawn copper) and t is the fault time. A 50 mm² copper floor is then applied over that thermal value for mechanical robustness and corrosion resistance, since the conductor stays buried for 30–40 years.
What are the Df and Sf factors?
Df is the decrement factor (IEEE 80 §15), which corrects the asymmetry of the fault current in the first cycles; Sf is the current division factor, representing the fraction of the short circuit that actually returns through the grid (the rest flows via shield wires, neutrals, etc.). The grid design current is I_G = Df·Sf·I. By default the tool adopts Df = Sf = 1 (the whole short circuit into the grid), which is conservative.
Which soil resistivity should I use?
The resistivity ρ should come from field measurement (Wenner method, IEEE 81), because it ranges from ~10 Ω·m in moist clay to >1,000 Ω·m in dry sandy or rocky soil. Real soils are stratified; when the variation is large, model it as two layers. Also account for the seasonal moisture variation, which can double ρ in the dry season.
Glossary
- Soil resistivity (ρ)
- Ground property that measures the opposition to current flow, in Ω·m. It is the main input variable and governs the grid resistance directly.
- GPR (Ground Potential Rise)
- The potential the grid reaches relative to remote earth during the fault, GPR = I·R. The higher it is, the greater the risk of dangerous voltages.
- Touch voltage
- Potential difference between the hand (in contact with a grounded structure) and the feet of a person. It is the most critical limit, usually lower than the step voltage.
- Step voltage
- Potential difference between a person's two feet, ~1 m apart, on the soil during the fault. The tolerable limit is given by IEEE 80 as a function of ρs, Cs, t and the body weight.
- Cs factor (gravel)
- Surface-layer derating factor. Cs = 1 with no gravel; with a high-resistivity layer Cs < 1, which raises the voltages the body tolerates.
- Adiabatic criterion
- Assumption that all the fault-current heat stays in the conductor (no dissipation to the soil) during the short circuit, used to size the minimum grid conductor section.
- Equivalent radius (r)
- Radius of the circle with the same area as the grid, r = √(A/π). It appears in the resistance floor ρ/(4r) and in the Laurent method.
- Decrement factor (Df)
- IEEE 80 multiplier that converts the symmetrical short-circuit current into the effective asymmetrical current of the first cycles, between 1 and ~2.