Earthing Fault Simulator (IEEE 80 & IS 3043)

Industrial-grade grounding and earthing fault analyzer. Models substation grid performance per IEEE Std 80-2013 (using Sverak/Schwarz formulas), low/medium voltage installations per IS 3043-2018 (using specific Pipe, Plate, and Strip electrode equations), and high voltage AC stations per IEC 61936-1.

1. Soil & Surface Parameters
2. Grid & Rods Geometry
3. Fault & System Data

Engineering Physics of Substation Earthing

Deep-dive into the thermodynamics, transient electromagnetics, and physiological constraints of industrial ground grids according to IEEE 80-2013 and IEC 60364-4-41 standards.

IEEE 80 / IEC 61936

1. The Invisible Mountain (GPR)

When a massive fault current ($I_G$) flows into the earth, it doesn't just disappear. The earth has resistance ($R_g$). According to Ohm's Law ($V = I_G \times R_g$), pushing thousands of amps into the ground creates a massive pressure spike called the Ground Potential Rise (GPR).

Imagine your substation is sitting on top of an "electrical volcano." At the moment of a fault, the soil directly under the station might jump to 5,000V or 10,000V relative to "remote earth". The danger isn't being at 10,000V; the danger is bridging a voltage difference.

MAX GPR (kV) Remote Earth (0V) Remote Earth (0V)
STEP VOLTAGE TOUCH VOLTAGE
PHYSIOLOGY

2. The Physics of Shock

IEEE 80 defines safe energy limits based on body weight ($50kg$ or $70kg$). The heart enters Ventricular Fibrillation at roughly 50mA-100mA. The formula $I_B = k / \sqrt{t_s}$ dictates that the faster a fault clears, the higher the transient current a human can survive.

Step Voltage ($E_s$): The potential difference between two feet (1m apart). Current flows leg-to-leg, bypassing the heart. Painful, but rarely fatal.
Touch Voltage ($E_t$): The potential difference between a grounded hand and the feet. Current flows directly through the chest. This is the deadly handshake and the primary focus of all compliance testing.

MITIGATION

3. The Crushed Rock Shield ($C_s$)

Substations are covered in crushed rock (gravel) for a critical safety reason. Wet native soil is a good conductor ($\approx 100 \Omega \cdot m$). Clean, wet gravel is a poor conductor ($\approx 3000 \Omega \cdot m$).

By standing on this high-resistance layer, you add a massive "resistor array" in series with your feet. The IEEE 80 Derating Factor ($C_s$) mathematically calculates how this shallow surface layer limits the fault current penetrating the body.

NATIVE SOIL (100 Ω-m) CRUSHED ROCK (3000 Ω-m) SERIES RESISTANCE
Frequency (Hz) Impedance (Z) 50Hz (Grid Fault) 100kHz+ (Lightning)
TRANSIENTS

4. Lightning vs. Phase Faults

Standard engineering designs for 50/60Hz faults. At low frequencies, the earth grid acts primarily as a pure resistor ($Z \approx R_g$). However, a lightning strike is a high-frequency impulse (equivalent to 100kHz or more).

At massive frequencies, the inductance ($L$) of your copper grid becomes the dominating barrier ($X_L = 2\pi f L$). A grid that measures a perfect $0.5\Omega$ at 50Hz might suddenly present $50\Omega$ of impedance to a lightning strike. The "effective area" of the grid shrinks dramatically because the impulse event is over before the fault wave can physically travel to the edges of your copper mesh.

Industrial Earthing Guidelines & FAQs

Advanced engineering Q&A, structural codes, and transient physics for high-voltage and industrial ground design.

1. How is soil resistivity ($\rho$) measured using the Wenner 4-Point method?

The Wenner 4-Point method is the international standard (per IEEE Std 81) for determining soil resistivity profiles. Four auxiliary electrodes are driven into the ground in a straight line at equal spacing $a$. A current $I$ is injected between the two outer electrodes, and the resulting potential difference $V$ is measured between the two inner electrodes.

The apparent resistivity $\rho$ is calculated using the following equation:

$$\rho = \frac{4\pi a R}{1 + \frac{2a}{\sqrt{a^2 + 4b^2}} - \frac{a}{\sqrt{a^2 + b^2}}}$$

Where $R = V/I$ is the measured resistance ($\Omega$), $a$ is the electrode spacing (m), and $b$ is the burial depth of the electrodes (m). If $b \ll a$ (burial depth is small compared to spacing), the formula simplifies to:

$$\rho = 2\pi a R$$
C1 P1 P2 C2 a a a I V

2. What is the physiological basis of the 50 kg vs. 70 kg body weight thresholds?

IEEE Std 80 utilizes Charles Dalziel's empirical formulas for the threshold of ventricular fibrillation. The human heart can tolerate a higher shock current if the duration is short, but the threshold drops exponentially as duration increases. The standard divides populations by average body weight:

50 kg Body Weight: Assumes a lower threshold for fibrillation (conservative, typical default):

$$I_B = \frac{0.116}{\sqrt{t_s}}$$

70 kg Body Weight: Used where the design environment warrants a higher weight limit assumption:

$$I_B = \frac{0.157}{\sqrt{t_s}}$$

During a fault, if the current through the body $I_{\text{body}} = E_{\text{touch}} / (1000 + 1.5\rho_s)$ exceeds $I_B$, fibrillation is statistically probable. Thus, designs are sized to keep body currents strictly below this curve.

Fault Duration (seconds) Tolerable Current I_b (mA) 70 kg Threshold (Ib = 0.157/√t) 50 kg Threshold (Ib = 0.116/√t) 0 2.0s

3. How does the surface layer rock/gravel reduce shock currents, and what is $C_s$?

Crushed rock possesses a high resistivity (typically $2500 - 3000\,\Omega\cdot m$ even when wet). When an operator stands on gravel, the rock layer acts as an insulator, introducing a high contact resistance in series with the feet. Each foot is modeled as a circular plate of radius $0.08\,\text{m}$, giving a ground resistance of $3\rho_s$ per foot.

Because the gravel thickness $h_s$ is finite and the native soil below has lower resistivity $\rho$, the electric field penetrates into the soil, reducing the effective resistance. To account for this, the surface layer derating factor $C_s$ is calculated using the reflection factor $K$:

$$K = \frac{\rho - \rho_s}{\rho + \rho_s}, \quad C_s = 1 - \frac{0.09\left(1 - \frac{\rho}{\rho_s}\right)}{2h_s + 0.09}$$

This derates the apparent resistivity of the rock. The total resistance of the body path for touch voltage is $R_B + 1.5 C_s \rho_s$, where $1.5$ represents two feet in parallel ($3/2 = 1.5$).

Crushed Gravel Layer (Ï s, hs) Native Soil (Ï ) Rf = 3 Cs Ï s Rf = 3 Cs Ï s

4. Why is the Decrement Factor ($D_f$) required, and how does high X/R affect it?

Ground faults are asymmetrical transients. At the instant of contact, the current contains an AC symmetrical component and a decaying DC offset component. The rate of decay of the DC component is governed by the subtransmission/generator system's reactive-to-resistive impedance ratio, or System X/R ratio.

The time constant of decay $T_a$ is defined as:

$$T_a = \frac{X}{2\pi f R}$$

To ensure that the ground conductors do not melt under the thermal stress of this asymmetrical peak and that human safety margins are maintained during the initial cycles of the fault, the symmetrical fault current is multiplied by the Decrement Factor $D_f$:

$$D_f = \sqrt{1 + \frac{T_a}{t_f}\left(1 - e^{-2t_f/T_a}\right)}$$

If the system X/R ratio is high (e.g., 30 near generating stations) and the fault duration is short (e.g. 0.1s), $D_f$ can exceed $1.3$, significantly increasing the effective design current $I_G = I_{\text{sym}} \times D_f$.

Decaying DC Offset Envelope Asymmetrical Fault Current

5. How does high frequency (e.g., lightning strikes) affect earthing impedance?

For power-frequency grounding faults (50/60 Hz), an earthing system acts primarily as a static resistance $R_g$. However, a lightning strike is a transient surge with high-frequency harmonics equivalent to 100 kHz or higher. At these frequencies, the copper grid's self-inductance ($L$) becomes dominant.

The transient impedance $Z_g(f)$ of the grid scales according to:

$$Z_g(f) = \sqrt{R_g^2 + (2\pi f L)^2}$$

Due to propagation delays along conductors, current cannot reach the extremities of a large grid during the transient wavefront. This limits the grid's "effective active radius" to:

$$r_{\text{eff}} \approx \sqrt{\frac{\rho \cdot t_r}{\pi \mu}}$$

Where $t_r$ is the wavefront rise time (e.g. 1.2 $\mu$s). Thus, increasing the grid's physical dimensions beyond $r_{\text{eff}}$ does not reduce surge impedance; instead, short, direct vertical ground rods or low-inductance radial strips must be installed close to surge arresters.

Log Frequency (Hz) Impedance |Z| (Ohms) Resistive Zone (Static Rg) Inductive Zone (2Ï€ f L) 50Hz/60Hz 100kHz (Lightning)

6. What is the design procedure for high-resistivity or rocky soil conditions?

Solid rock has an extremely high soil resistivity ($\rho > 2000\,\Omega\cdot m$), which makes achieving touch/step safety limits difficult. A standard shallow mesh grid is insufficient. The following engineering measures are applied:

  • Deep Boreholes: Deep vertical boreholes (up to 30m) are drilled to tap into moist, lower-resistivity water tables.
  • Ground Enhancement Materials (GEM): Native rocky backfill around rods is replaced with conductive material such as Bentonite clay, Marconite, or chemical gel backfills, which artificially expand the electrical boundary of the conductor.
  • Counterpoise/Radial Wires: Horizontal radial strips are extended outward from the substation perimeter into areas with softer or deeper soil layers.
Dry Rock Layer (Ï 1 = 3000 Ohm-m) Moist Saturated Soil Layer (Ï 2 = 150 Ohm-m) Bentonite / GEM Backfill Copper Ground Rod

7. How does mesh spacing ($D$) impact touch vs. step voltages?

Changing the spacing of grid conductors has different effects on the grid's overall resistance and local safety voltages:

Grid Resistance ($R_g$): The overall grounding resistance is primarily determined by the total area of the grid ($A$) and the total length of the conductors ($L_T$). Decreasing the mesh spacing from 10m to 5m increases the total length of the conductor, which slightly reduces $R_g$ (typically by less than 5%), but this is an inefficient way to lower the overall ground resistance.

Mesh/Touch Voltage ($E_m$): In contrast, touch voltage is a local safety parameters. Shrinking the mesh spacing creates a denser conductor layout. This reduces the electric fields in the ground between conductors, smoothing out the surface potential profile and significantly lowering the actual touch voltage $E_m$ at the mesh center.

Surface Potential Rise Profile Worst Touch Voltage (Em) at center Mesh Spacing D

8. Why are the highest step and touch potentials located at the corners?

Under fault conditions, current leaves the grid and enters the surrounding earth. Conductor elements at the perimeters and corners of the grid have fewer neighboring conductors to share the current discharge. This leads to current crowding at the outer perimeter, resulting in high current densities in the soil at the grid corners.

This localized current crowding creates steep voltage gradients at the edges of the grid. Consequently, the step potential ($E_s$) is highest just outside the perimeter grid loop, and the mesh touch potential ($E_m$) is highest in the corner mesh compartments. To mitigate this risk, ground grids are designed with closely spaced outer conductors and vertical ground rods driven at the grid corners.

Severe Corner Current Crowding Shared Interior Mesh

9. What is Equipotential Bonding, and why is it more critical than low resistance?

A common misconception is that achieving an extremely low grounding resistance (e.g., $< 1\,\Omega$) guarantees safety. If a substation experiences a massive 20 kA fault, even a ground resistance of $0.5\,\Omega$ results in a Ground Potential Rise (GPR) of:

$$GPR = I_G \times R_g = 20000 \times 0.5 = 10000 \text{ V}$$

The substation grid and all connected metal enclosures jump to 10,000 V relative to remote earth. If an operator touches a metal enclosure while standing on the ground, they are exposed to a shock hazard unless Equipotential Bonding is established.

By bonding all metallic structures, fences, operating handles, and concrete reinforcing bars to the ground grid, everything rises to the same electrical potential during a fault. Because there is no potential difference ($\Delta V \approx 0$) across the operator's body, no current flows, and safety is maintained regardless of the GPR magnitude.

Flexible Copper Braid Jumper Main Ground Grid Loop

10. How do IS 3043:2018 and IEEE Std 80-2013 standards compare?

The two standards are applied to different earthing scenarios and have different sizing philosophies:

Design Parameter IS 3043-2018 (Indian Code) IEEE Std 80-2013 (Global Standard)
Primary Scope Low/Medium voltage industrial, commercial, and residential installations. High voltage and extra high voltage transmission/distribution AC substations.
Sizing Criteria Empirical formulas for specific electrode shapes (pipe, plate, strip) to meet a target resistance (e.g., $< 1\,\Omega$ or $< 5\,\Omega$). Safety-based limits where grid geometry is designed to keep touch/step voltages below physiological thresholds.
Electrode Formulas Logs equations based on single electrode dimension profiles ($L$, $d$, depth $h_d$). Multi-mesh grid formulas using Sverak and Schwarz equations to calculate overall ground impedance.
Transients Correction Assumes static power frequencies (50 Hz) for earth resistance sizing. Includes Decrement Factor ($D_f$) for asymmetrical DC offsets and high-frequency checks.
IS 3043 (Pipe/Plate) IEEE 80 (Mesh Grid)