Vertical Electrical Sounding Lab
Schlumberger and Wenner soundings to a layered earth: depths, true resistivities, aquifers and a geoelectric section.
Turns the field readings of a resistivity sounding into apparent resistivity, joins the curve segments measured with different potential electrode spacings, and inverts the curve to a model of two to six horizontal layers by damped least squares with an exact forward solution of the layered earth. Each layer gets its thickness, depth and true resistivity with an uncertainty, an interpretation in terms of rock, clay and water, and the Dar Zarrouk parameters. Soundings saved along a line are joined into a geoelectric section.
Overview
Vertical electrical sounding is the most used geophysical method in Iraq for groundwater, for the depth of gypsum and salt, and for foundation and dam sites. The field work is simple, but the interpretation is usually done with old curve-matching charts or desktop programs that most students and field geologists do not have. This lab does the full interpretation on a phone, without a connection, from the moment the last reading is taken.
Readings can be entered as apparent resistivity, as resistance ΔV/I, or as voltage and current, typed or pasted from a spreadsheet. The number of layers can be chosen or found automatically, any thickness or resistivity known from a borehole can be fixed, and the model can be edited by hand to see how the curve responds.
What you enter · What you get
What you enter
- Array (Schlumberger or Wenner) and form of the readings
- AB/2 (or a), MN/2 and the reading of each point
- Number of layers, or automatic
- Thickness or resistivity known from boreholes (optional)
- Groundwater resistivity or conductivity (optional)
- Station name, distance along the line and elevation
What you get
- Apparent resistivity and joined curve segments
- Layer thickness, depth and true resistivity with uncertainty
- Curve type, fit error and depth of investigation
- Interpretation of each layer and likely aquifers, with porosity and water salinity estimates
- Dar Zarrouk parameters and anisotropy
- Sounding curve, geoelectric column and geoelectric section
How it works in the field
Method and equations
Apparent resistivity
For a symmetric array with current electrodes at ±L and potential electrodes at ±l the geometric factor is K = π(L² − l²)/(2l); Schlumberger uses L = AB/2 and l = MN/2, Wenner uses L = 1.5a and l = 0.5a so that K = 2πa.
Layered earth
The resistivity transform of the layers is built upward from the half-space by Pekeris’ recurrence, and the potential is its Hankel transform, evaluated with the 61-point digital linear filter of Guptasarma and Singh (1997). The exact finite MN of each reading is used, so segments are modelled as measured.
Inversion
The logarithms of the resistivities and thicknesses are adjusted by the Levenberg–Marquardt method to minimise the misfit of log ρa, from several starting models; the parameter uncertainty comes from the covariance matrix at the solution.
Interpretation
Dar Zarrouk longitudinal conductance S = Σh/ρ and transverse resistance T = Σhρ (Maillet 1947); Archie’s law F = ρ/ρw = φ⁻ᵐ for clean water-bearing sand; water salinity TDS ≈ 0.64 EC.
Accuracy and error limits
- The forward solution agrees with the exact image series of a two-layer earth to better than 0.0001 %. The uncertainty of the result is that of the data and of the method: with a fit of 2–5 % the resistivity of thick layers is usually within ±10–20 % and their depth within ±10–25 %, as shown for each layer.
- What is not covered: equivalence (a thin conductive layer is resolved only by its conductance h/ρ and a thin resistive layer only by its transverse resistance hρ, so different models fit equally well), suppression of a thin layer whose resistivity lies between its neighbours, dipping or laterally changing layers, and anything deeper than about a third of the largest AB/2. The lithology is a suggestion from resistivity ranges; clay and saline water give the same low values.
- To reduce the error: use six to ten readings per decade of AB/2, overlap two readings when MN is changed, keep AB/2 at least five times MN/2, extend AB/2 to three times the target depth, fix layers known from nearby boreholes, compare the result with the logs, and repeat soundings across lateral changes rather than trusting a single one.
Who uses it
Standards and references
- Koefoed, O. (1979). Geosounding Principles 1: Resistivity Sounding Measurements. Elsevier, Amsterdam.
- Guptasarma, D. & Singh, B. (1997). New digital linear filters for Hankel J0 and J1 transforms. Geophysical Prospecting, 45(5), 745–762. Filter coefficients from the libdlf library (CC BY 4.0).
- Pekeris, C.L. (1940). Direct method of interpretation in resistivity prospecting. Geophysics, 5(1), 31–42.
- Marquardt, D.W. (1963). An algorithm for least-squares estimation of nonlinear parameters. Journal of the Society for Industrial and Applied Mathematics, 11(2), 431–441.
- Maillet, R. (1947). The fundamental equations of electrical prospecting. Geophysics, 12(4), 529–556.
- Archie, G.E. (1942). The electrical resistivity log as an aid in determining some reservoir characteristics. Transactions of the AIME, 146, 54–62.
- Telford, W.M., Geldart, L.P. & Sheriff, R.E. (1990). Applied Geophysics, 2nd ed. Cambridge University Press.
- Zohdy, A.A.R., Eaton, G.P. & Mabey, D.R. (1974). Application of surface geophysics to ground-water investigations. USGS Techniques of Water-Resources Investigations, Book 2, Chapter D1.
At a glance
- Category
- Geophysics
- Runs on
- Phone, tablet and desktop browser
- Export
- Sounding curve, geoelectric column and section as PNG, JPG and PDF; readings, model and stations as CSV; project backup as JSON; full PDF report
- Access
- Free
Run this survey on your next site
All 64 tools, report export and the AI assistant come with one subscription.
