Knowledge / Geology / Geophysical exploration

Geophysical exploration

Geology (B.Sc) · Applied geology · free study notes for exams, placements and interviews

In short

Geophysical exploration measures contrasts in physical properties of rocks - density, magnetic susceptibility, electrical resistivity, elastic wave velocity, radioactivity - from the surface, air or boreholes, and inverts those measurements into a picture of the subsurface. It is the only non-destructive way to "see" rock below cover before drilling, so it decides where expensive boreholes go in mineral, groundwater, petroleum and civil projects. In India it is practised by the Geological Survey of India (airborne and ground surveys), NGRI Hyderabad, CGWB, ONGC, Oil India, AMD, MECL, NMDC, Coal India/CMPDI and private consultants. It carries heavy weight in GATE Geology & Geophysics (GG), UPSC Geologist/Geophysicist (Combined Geo-Scientist) and JAM/university Applied Geology papers, and almost every campus technical round on exploration starts with "which method would you choose and why".

Explained

Fundamentals: what geophysics actually measures

Seismic methods

Gravity and magnetic methods

Electrical and electromagnetic methods

Radiometric, well-logging and other methods

Indian practice, standards and regulation

Key points

Formulas and examples

ρa = K · ΔV / I — apparent resistivity (Ω·m); ΔV in volts, I in amperes; K in metres
K(Wenner) = 2πa — a = electrode spacing (m)
K(Schlumberger) = π (L² − l²) / (2l) — L = AB/2, l = MN/2 (m)
Snell's law: sin i / V₁ = sin r / V₂; critical angle ic = sin⁻¹(V₁/V₂)
Refraction, two layers (shot-point intercept): z₁ = (ti · V₁V₂) / (2√(V₂² − V₁²)) — ti = intercept time (s), z₁ = layer thickness (m)
Refraction, crossover distance: xc = 2 z₁ √((V₂+V₁)/(V₂−V₁))
Reflection (flat layer, CMP): t² = t₀² + x²/V²; t₀ = 2z/V — x = offset (m), t = travel time (s)
Reflection coefficient R = (ρ₂V₂ − ρ₁V₁)/(ρ₂V₂ + ρ₁V₁); Z = ρV (kg·m⁻²·s⁻¹)
Newton: g = G M / r²; G = 6.674 × 10⁻¹¹ N·m²/kg²; 1 mGal = 10⁻⁵ m/s²
Free-air correction = +0.3086 h (mGal, h in m); Bouguer slab = −0.04193 ρ h (mGal, ρ in g/cm³, h in m)
Bouguer anomaly = g(obs) − g(theoretical) + 0.3086h − 0.04193ρh + terrain correction
Buried sphere, maximum half-width: depth z ≈ 1.305 x½ (x½ = half-width of the anomaly profile)
Total magnetic field F = √(H² + Z²); tan I = Z/H — I = inclination, H and Z in nT
Magnetic induction in a body: J = κ H — κ = magnetic susceptibility (dimensionless, SI)
Skin depth δ = 503 √(ρ / f) metres — ρ in Ω·m, f in Hz
Archie: Sw = [ (a · Rw) / (φm · Rt) ]1/n — Sw = water saturation, φ = porosity, Rw = formation water resistivity, Rt = true resistivity, a ≈ 1, m ≈ 2, n ≈ 2
Chargeability (time domain) M = (1/(t₂−t₁)) ∫ V(t) dt / Vp (mV/V)
Seismic velocity in a rod/solid: Vp = √((K + 4µ/3)/ρ) ; Vs = √(µ/ρ)

Example 1 - Seismic refraction depth to bedrock (dam site). A 12-geophone refraction spread over weathered alluvium gives direct-wave velocity V₁ = 600 m/s, refractor velocity V₂ = 2,400 m/s, and an intercept time ti = 0.030 s.

Example 2 - Apparent resistivity, Wenner and Schlumberger. (a) Wenner array, a = 10 m, I = 0.5 A, ΔV = 0.040 V.

(b) Schlumberger, AB/2 = 100 m, MN/2 = 5 m, I = 1.2 A, ΔV = 0.018 V.

Example 3 - Bouguer anomaly reduction. Observed g = 978,520.40 mGal at a station of elevation 420 m; theoretical (latitude) value = 978,540.00 mGal; assumed crustal density ρ = 2.67 g/cm³; terrain correction = +1.20 mGal.

Example 4 - Reflection coefficient and NMO. Shale (ρ = 2.40 g/cm³, V = 2,500 m/s) overlies limestone (ρ = 2.65 g/cm³, V = 4,000 m/s).

Example 5 - Skin depth and Archie saturation. (a) An EM system at f = 1,000 Hz over ground of ρ = 50 Ω·m: δ = 503 √(50/1,000) = 503 × 0.2236 = 112 m. Dropping to f = 10 Hz: δ = 503 × √5 = 503 × 2.236 = 1,125 m - depth needs low frequency. (b) A sandstone log gives φ = 0.25, Rt = 20 Ω·m, Rw = 0.5 Ω·m, a = 1, m = n = 2.

MethodProperty sensedTypical depth reachMain targetsCost/speed
Reflection seismicAcoustic impedance100 m to >10 kmOil, gas, coal, basin structureVery high cost, slow
Refraction seismicP velocity5-100 m (engineering)Bedrock depth, rippability, Vs30Moderate
GravityDensity10 m to crustalChromite, iron, basins, cavitiesLow cost, slow stations
MagneticsSusceptibility, remanenceSurface to basementMagnetite, kimberlite, dykes, basementCheapest areal
Resistivity/ERTResistivity1-300 mGroundwater, salinity, foundationsLow cost
IPChargeability50-500 mDisseminated sulphidesModerate, slow
EM/TEMConductivity10 m to >500 mMassive sulphides, aquifers, saline zonesModerate
SPNatural potential<100 mSulphides, graphite, dam seepageVery cheap
RadiometricsK, U, ThCentimetresUranium, lithology, placersCheap, airborne
GPRDielectric permittivity0.5-20 mUtilities, rebar, shallow layers, cavitiesModerate
ArrayGeometric factor KStrengthWeakness
Wenner2πaStrong signal, good vertical resolutionNeeds all four electrodes moved; slow
Schlumbergerπ(L² − l²)/2lDeep sounding, fewer moves, standard for VESWeak signal at large L; MN must be increased
Dipole-dipoleπn(n+1)(n+2)aExcellent lateral resolution, used for IPLow signal, high noise sensitivity

Where it is used

Interview questions

Common mistakes

Related topics

Hydrogeology · Economic geology and ore deposits · Petroleum geology · Engineering geology · Natural hazards · Remote sensing and GIS in geology

Open in the Knowledge reader →Find jobs →