ToolFusion Civil
ToolFusion Civil
Calculator D4

Geophysical Methods in Site Investigation: Seismic Refraction vs. Resistivity

Seismic refraction measures how fast sound waves travel through the ground to find layers and their depths; resistivity measures how much the ground resists electricity to detect water, clay, or rock types.

⚠️ Why It Matters

1
Inaccurate bedrock depth estimation
↓
2
Overdesign of foundations or underdesign of piles
↓
3
Excessive settlement or differential movement
↓
4
Structural damage to superstructure
↓
5
Costly remediation or litigation

πŸ“˜ Definition

Seismic refraction is a geophysical method that interprets the travel times of compressional (P-) waves refracted at subsurface acoustic impedance boundaries to infer layered velocity structure and depth to bedrock or competent strata. Electrical resistivity imaging (ERI) measures spatial variations in bulk electrical resistivity by injecting direct current between electrode pairs and measuring resulting potentials, enabling discrimination of lithology, saturation state, and fracture zones based on charge-transfer properties of subsurface materials.

🎨 Concept Diagram

Ground SurfaceBedrock (Vp, ρ)Weathered ZoneSoil/OverburdenResistivity currentSeismic rayElectrodesGeophone

AI-generated illustration for visual understanding

πŸ’‘ Engineering Insight

Never treat seismic refraction or resistivity as standalone truth β€” both are non-unique inverse problems. A 10% error in first-break picking can shift interpreted bedrock depth by 2–5 m in low-velocity overburden; always tie results to at least two independent boreholes with P-wave logging and resistivity probes. Resistivity anomalies require hydrogeologic context: a 100 Ω·m zone may be clay *or* saline groundwater β€” never assign lithology without supporting lab tests.

πŸ“– Detailed Explanation

Geophysical methods like seismic refraction and resistivity provide indirect, large-scale views of subsurface conditions where drilling alone is sparse and costly. Refraction relies on wavefronts bending upward when crossing into faster layers β€” the 'head wave' arrival time vs. distance plot yields velocities and layer depths via intercept-time and slope analysis. Resistivity uses Ohm’s law applied across the ground: injected current flows preferentially through conductive paths (e.g., water-filled fractures), and measured potential differences allow reconstruction of subsurface resistivity distribution.

Both methods face fundamental limitations: refraction fails where velocity inversions occur (e.g., stiff clay over soft sand), and resistivity struggles with equivalence (e.g., thin conductive layer vs. thick moderate-resistivity layer yielding identical ρa). Modern practice uses joint inversion β€” coupling seismic and resistivity datasets under shared geological constraints β€” to reduce non-uniqueness and improve confidence in interpreted interfaces and material properties.

At advanced levels, time-lapse (4D) resistivity monitors moisture migration during dewatering or grouting, while seismic interferometry extracts virtual refraction responses from ambient noise β€” eliminating need for active sources. Integration with LiDAR-derived topography, GIS-based geological modeling, and machine-learning-assisted inversion (e.g., physics-informed neural networks) is now standard in major infrastructure projects, but only when validated against rigorous ground-truthing protocols per ASTM G304 and ISO 22475-1.

πŸ”„ Engineering Workflow

Step 1
Step 1: Define geotechnical objectives (e.g., bedrock depth, fault location, aquifer geometry)
β†’
Step 2
Step 2: Conduct desktop study (geology maps, borehole logs, historical geophysics)
β†’
Step 3
Step 3: Select method(s) and acquisition parameters (array geometry, electrode spacing, shot points, sampling rate)
β†’
Step 4
Step 4: Acquire field data with QA/QC (repeatability checks, noise monitoring, polarity verification)
β†’
Step 5
Step 5: Process and invert data (time-term analysis for refraction; smoothness-constrained least-squares for resistivity)
β†’
Step 6
Step 6: Integrate with borehole control (Vp-log correlation, ρ-log calibration, lithology mapping)
β†’
Step 7
Step 7: Deliver interpreted section with uncertainty bounds and engineering recommendations

πŸ“‹ Decision Guide

Rock/Field Condition Recommended Design Action
Shallow, highly variable overburden (e.g., glacial till with boulders, >5 m thick) Use combined seismic refraction + resistivity; deploy short-offset refraction (≀50 m) with high-density ERI (1–2 m electrode spacing) to resolve lateral heterogeneity
Deep (>30 m), low-velocity weathered zone over competent bedrock with suspected karst features Prioritize high-resolution resistivity (dipole-dipole, 64+ electrodes) with tomographic inversion; supplement with long-offset refraction (β‰₯200 m) to constrain Vp of bedrock
Urban site with high electromagnetic noise and limited access Use passive seismic (microtremor H/V) and resistivity with gradient array to minimize current injection; avoid refraction due to poor signal-to-noise ratio

📊 Key Properties & Parameters

P-wave Velocity (Vp)

300–6,000 m/s (soil: 150–800 m/s; weathered rock: 800–2,500 m/s; intact bedrock: 2,500–6,000 m/s)

The propagation speed of compressional seismic waves through a material, governed by elastic moduli and density.

⚡ Engineering Impact:

Directly correlates with stiffness and bearing capacity β€” used to estimate Young’s modulus and classify rock mass quality per ASTM D7400.

Bulk Resistivity (ρ)

1–10,000 Ω·m (clay: 1–50 Ω·m; saturated sand: 50–500 Ω·m; limestone: 500–5,000 Ω·m; dry granite: 2,000–10,000 Ω·m)

The intrinsic opposition of a volume of earth material to electrical current flow, expressed in ohm-meters.

⚡ Engineering Impact:

Enables identification of groundwater-bearing zones, clay lenses, or saline contamination β€” critical for drainage design and corrosion risk assessment.

Critical Distance (Xc)

10–200 m (depends on layer thickness and velocity contrast)

The minimum offset distance at which the refracted head wave first arrives before the direct wave, marking the transition to layer-controlled arrivals.

⚡ Engineering Impact:

Determines minimum survey length required to resolve target layer β€” undersized arrays yield ambiguous or non-unique velocity models.

Apparent Resistivity (ρa)

Same as bulk resistivity but spatially variable; standard Wenner array ρa ranges: 1–5,000 Ω·m

The resistivity value computed from measured voltage and current assuming a homogeneous half-space, used as input for inversion modeling.

⚡ Engineering Impact:

Misinterpretation of ρa without 2D/3D inversion leads to false stratigraphic boundaries and erroneous soil-structure interaction assumptions.

πŸ“ Key Formulas

Refraction Intercept Time

t_i = (2h / V1) * √((V2Β² βˆ’ V1Β²) / (V2Β²))

Calculates intercept time on travel-time curve to determine depth h to refractor with velocities V1 (overburden) and V2 (refractor)

Variables:
Symbol Name Unit Description
t_i Refraction Intercept Time s Intercept time on the travel-time curve
h Depth to Refractor m Vertical depth from surface to refracting layer
V1 Overburden Velocity m/s Seismic wave velocity in the overlying (upper) layer
V2 Refractor Velocity m/s Seismic wave velocity in the refracting (lower) layer
Typical Ranges:
Shallow clay over chalk
20–60 ms
Gravel over bedrock
10–40 ms
⚠️ t_i < 100 ms ensures reliable layer resolution; >150 ms indicates insufficient offset or excessive noise
πŸ”§ Open Calculator

Wenner Apparent Resistivity

ρ_a = 2Ο€aR

Computes apparent resistivity from measured resistance R and electrode spacing a

Variables:
Symbol Name Unit Description
ρ_a Apparent Resistivity Ω·m Resistivity derived from measured resistance and electrode spacing
a Electrode Spacing m Distance between adjacent electrodes in the Wenner array
R Measured Resistance Ξ© Resistance measured between the inner electrodes
Typical Ranges:
Urban clayey soils
10–100 Ω·m
Carbonate bedrock
500–3,000 Ω·m
⚠️ R < 10 kΞ© required for stable measurement; use current β‰₯100 mA in high-resistivity terrain
πŸ”§ Open Calculator

🏭 Engineering Example

Crossrail Bond Street Station Box, London, UK

London Clay (overlying Chalk bedrock)
Vp (Clay)
650 m/s
ρ (Clay)
35 Ω·m
Vp (Chalk)
2,400 m/s
ρ (Chalk)
220 Ω·m
Bedrock Depth (refraction)
18.2 Β± 0.7 m
Bedrock Depth (resistivity)
17.9 Β± 1.1 m

πŸ—οΈ Applications

  • Foundation design for tall buildings
  • Tunnel alignment and excavation support planning
  • Landfill liner integrity assessment
  • Groundwater resource characterization

πŸ”§ Try It: Interactive Calculator

πŸ”§
Geotechnical Site Investigation Calculator
Multi-tab calculator for Geotechnical Site Investi...
Open β†’

πŸ“‹ Real Project Case

Urban Transit Tunnel Alignment Through Mixed-Soil Stratigraphy

3.2 km cut-and-cover metro extension in Jakarta, Indonesia

Challenge: Variable soil profile (soft clay β†’ weathered volcanic tuff β†’ dense sand) causing differential settle...
Dense Sand (Ο†β€²=36.4Β°, Kβ‚€=0.41)Weathered Volcanic TuffSoft Clay (Cv=0.82 mΒ²/yr)InclinometerSecant PilesJet-grouted secant piles (staged excavation)Differential settlement & excavation instabilitySoil Stratigraphy Survey:SPT + CPT + Seismic RefractionDesign Parameters:Ο†β€² = 36.4Β° | Kβ‚€ = 0.41 | Cv = 0.82 mΒ²/yr
Read full case study β†’
High-Rise Foundation Design on Residual Lateritic Soil Tailings Storage Facility (TSF) Stability Assessment Post-Earthquake Offshore Wind Turbine Monopile Foundation in Glacial Till Historic Masonry Bridge Retrofit on Alluvial Floodplain

🎨 Technical Diagrams

SurfaceBedrock (Vp=2400 m/s)Refraction ray pathShot
ρ = 35 Ω·m (Clay)ρ = 220 Ω·m (Chalk)ρ = 5,000 Ω·m (Limestone)Depth

πŸ“š References

[1]
[2]
[3]
Guidelines for the Use of Geophysical Methods in Engineering Site Investigations β€” International Society for Soil Mechanics and Geotechnical Engineering (ISSMGE)