Hydrogeology: Porosity and Permeability, Aquifer Properties, Darcy’s Law, Groundwater Flow, Well Hydraulics and Exploration

This chapter answers the Part B1 heading “Hydrogeology”: “Porosity and permeability; Groundwater storage and aquifer properties; Darcy’s law; Hydraulic conductivity and groundwater flow dynamics; Groundwater exploration, well hydraulics.” It is the second subject of Part B1 in the Geology and Geophysics (GG) paper. It teaches pore space and conductivity, aquifer types and storage coefficients, Darcy’s law with its limits, flow in layered and regional systems, steady and transient flow to a well, and the geological, remote-sensing and geophysical search for groundwater. The numericals are layered conductivity, specific yield and recharge, seepage velocity and travel time, Thiem and Cooper-Jacob analysis and Archie’s law.

1. Porosity and permeability

Porosity n = V_v/V is the fraction of a rock or soil that is void. Primary porosity forms with the rock (spaces between grains, vesicles in basalt); secondary porosity comes later (joints, fractures, solution cavities, weathering). Porosity depends on sorting, packing and grain shape but not on grain size for uniform spheres: cubic packing gives about 47.6 per cent and rhombohedral packing about 26 per cent. Poor sorting lowers n because fines fill the gaps between large grains. Effective porosity n_e counts only interconnected pores that can transmit water, so n_e ≤ n. After drainage, n = S_y + S_r: specific yield (water that drains under gravity) plus specific retention (water held by molecular attraction and capillarity).

Permeability is the ease of transmission. Intrinsic permeability k (m², or darcy, 1 darcy ≈ 10⁻¹² m²) belongs to the medium alone; hydraulic conductivity K = kρg/μ (m/s) also depends on the fluid. Typical orders: sand about 10⁻⁵ to 10⁻³ m/s, silt 10⁻⁹ to 10⁻⁶ m/s, clay below 10⁻⁹ m/s. Fractured, karstic and vesicular rocks owe high K to secondary porosity, while sound crystalline rock has very low K.

Layering. For flow parallel to layers of thickness b_i and conductivity K_i, K_h = ΣK_i b_i/Σb_i (arithmetic, controlled by the best layer); for flow across them, K_v = Σb_i/Σ(b_i/K_i) (harmonic, controlled by the worst layer). Example: layers (2 m, 10 m/d), (3 m, 1 m/d), (5 m, 0.1 m/d) give K_h = 23.5/10 = 2.35 m/d and K_v = 10/(0.2 + 3 + 50) = 0.19 m/d. K_h always exceeds K_v in layered ground.

⚠️ High porosity is not high permeability
Clay can have a higher total porosity than sand, yet its pores are tiny, water clings to the grains and K is smaller by orders of magnitude. A productive aquifer needs a high K and a large specific yield.

2. Groundwater occurrence, storage and aquifer properties

Below the surface the unsaturated (vadose) zone holds soil water, intermediate water and the capillary fringe; below the water table (pressure equal to atmospheric) the saturated zone holds groundwater. An aquifer yields usable water; an aquitard is a poorly permeable layer that leaks; an aquiclude (clay) transmits negligible water. An unconfined aquifer has a free water table. A confined aquifer lies between aquitards, its water is under pressure and rises in a well to the piezometric surface; where that surface is above ground the well flows (artesian). A perched water table sits on a local aquitard above the main one.

  • Specific yield S_y: volume of water an unconfined aquifer yields per unit area per unit fall of the water table, by draining pores; typically 0.01 to 0.3.
  • Specific storage S_s: volume released per unit volume of aquifer per unit fall of head by the elastic expansion of water and compaction of the skeleton, S_s = ρg(α + nβ). Storativity S = S_s b is dimensionless: roughly 10⁻⁵ to 10⁻³ for confined aquifers, while for unconfined ones S ≈ S_y.
  • Transmissivity T = K b (m²/d) is the flow per unit width under unit gradient through the full saturated thickness b, so Q = T i w across a width w.

Water released is V = S·A·Δh. For an unconfined aquifer of 50 km² with S_y = 0.12, a 2 m fall of the water table releases 50 × 10⁶ × 0.12 × 2 = 12 million m³. A confined aquifer of the same area with S = 2 × 10⁻⁴ releases only 0.2 million m³ for a 20 m head fall: it gives water by pressure relief, not by draining pores, and stays fully saturated. In the water-table fluctuation method the recharge of a season is R = S_y Δh per unit area (plus the draft, if the rise observed is net of pumping).

⚠️ Pumping a confined aquifer does not dewater it
Head falls but the aquifer stays saturated, so T is unchanged unless head drops below the aquifer top. In an unconfined aquifer pumping dewaters the upper part, the saturated thickness b falls and so does T.

3. Darcy’s law and its limits

Henry Darcy (1856), working on the sand filters of Dijon, found that discharge is proportional to the area and the head loss and inversely proportional to the length: q = Q/A = −K dh/dl = K i, and Q = K i A. The hydraulic head is h = z + p/(ρg) (elevation plus pressure head; velocity head is negligible). Water moves from high to low total head, not necessarily from high to low pressure. The gradient i = Δh/L is dimensionless.

The Darcy flux q is a volume per unit area per time, not a speed of water. Water moves only through the effective pores, so the seepage (average linear) velocity is v_s = q/n_e, larger than q, and the travel time over a distance L is L/v_s. Example: heads 105.2 m and 104.7 m in wells 250 m apart give i = 0.002; with K = 20 m/d and n_e = 0.25, q = 0.04 m/d, v_s = 0.16 m/d, and 400 m takes 2500 days (about 6.8 years). Three wells are needed to fix the gradient and its direction: flow is down-gradient, perpendicular to the head contours.

  • Limits. Darcy’s law holds for laminar flow, up to a Reynolds number (based on grain size) of about 1 to 10. It fails at high velocities in coarse gravel or open fractures, especially near a pumped well, and applies to saturated flow, since K in the unsaturated zone depends on moisture.
  • Laboratory K. Constant-head permeameter (sands): K = VL/(A h t). Falling-head permeameter (silts, clays): K = (aL/At) ln(h₁/h₂). Small samples miss fractures and layers, so field pumping tests give the representative K.
  • Flow nets and Dupuit. In a flow net, flow lines cross equipotentials at right angles and q = K H N_f/N_d per unit width (N_f flow channels, N_d head drops). For unconfined flow between two water levels h₁ and h₂ a distance L apart, Dupuit gives q = K(h₁² − h₂²)/(2L).
⚠️ Pressure is not head
Always compare total head, never pressure or depth to water alone. A deeper piezometer has a larger pressure head, yet the flow still goes from the larger total head to the smaller. Also, the Darcy flux is not the speed at which a tracer moves; the tracer moves at v_s.

4. Hydraulic conductivity and groundwater flow dynamics

Refraction. At a boundary between media of conductivity K₁ and K₂, flow lines bend such that tan θ₁/tan θ₂ = K₁/K₂ (θ measured from the normal to the boundary): they turn towards the normal in the less conductive layer. Flow concentrates in the high-K layers, which is why lenses of sand and gravel, buried channels and fracture zones carry most of the water.

Governing equations. For steady flow in a homogeneous isotropic medium, head obeys Laplace’s equation, ∇²h = 0. For transient flow in a confined aquifer, ∂²h/∂x² + ∂²h/∂y² = (S/T) ∂h/∂t. The ratio T/S is the hydraulic diffusivity: it is large in confined aquifers (S is small), so a pressure change at one point is felt far away quickly, whereas in an unconfined aquifer the response is slow because S ≈ S_y is large.

Regional flow. Recharge (rain, influent streams, canals, irrigation return) enters at topographic highs and discharges at lows as springs, effluent streams, wetlands and evapotranspiration. The water table is a subdued copy of the topography, and flow systems nest as local, intermediate and regional (Tóth). In a recharge area head falls with depth (downward flow); in a discharge area head rises with depth (upward flow, artesian conditions possible). Natural velocities are small, from centimetres to a few metres per day, so residence times run from years to millennia.

🧠 Reading a piezometer nest
Total head = screen elevation + water level above it. Compute it at each depth, then subtract. If the deeper piezometer shows the lower total head, flow is downward (recharge setting); if higher, flow is upward (discharge setting). Never decide from pressure head alone.

5. Well hydraulics

Pumping lowers the head around a well to form a cone of depression; drawdown s = h₀ − h. For steady radial flow to a fully penetrating well (Thiem, equilibrium) in a confined aquifer, with heads h₁ and h₂ at radii r₁ and r₂, Q = 2πT(h₂ − h₁)/ln(r₂/r₁), that is s₁ − s₂ = [Q/(2πT)] ln(r₂/r₁). For an unconfined aquifer, Q = πK(h₂² − h₁²)/ln(r₂/r₁). Example: Q = 1800 m³/d, T = 600 m²/d, observation wells at 5 m and 50 m give s₁ − s₂ = (1800/(2π × 600)) × ln 10 = 1.10 m.

Theis (non-equilibrium) gives the drawdown at radius r and time t: s = Q W(u)/(4πT), u = r²S/(4Tt). It assumes a confined, homogeneous, isotropic aquifer of uniform thickness and infinite extent, a fully penetrating well of negligible diameter, an initially flat piezometric surface, constant Q and instantaneous release from storage. For u < 0.05 (later times, or small r), the Cooper-Jacob form is a straight line on a semi-log plot: s = [2.303Q/(4πT)] log₁₀(2.25Tt/(r²S)). The drawdown per log cycle Δs gives T = 2.303Q/(4πΔs) ≈ 0.183Q/Δs, and the time t₀ at which the line cuts s = 0 gives S = 2.25Tt₀/r². Example: Q = 1200 m³/d and Δs = 0.5 m give T = 440 m²/d.

Boundaries are handled by image wells. An impermeable boundary (fault, bedrock wall) acts like an image well pumping the same rate, so drawdown rises faster and the semi-log slope doubles. A recharge boundary (perennial river) acts like an image injection well, so drawdown flattens and stabilises. Well performance: the specific capacity is Q/s_w, and the total drawdown in the pumped well is s_w = BQ + CQ², where BQ is the linear aquifer loss and CQ² the well loss from turbulent flow through screen and casing. Drawdowns of neighbouring wells add (superposition). A recovery test, with residual drawdown against t/t′, gives T independently. In an unconfined aquifer the Theis equation is usable only if drawdown is small against the saturated thickness b, after correcting s′ = s − s²/(2b).

⚠️ Which log, which constant
The Thiem difference uses the natural logarithm, ln(r₂/r₁); the Cooper-Jacob slope is per log₁₀ cycle, hence the 2.303. Cooper-Jacob is valid only for small u, so early data lie above the straight line and must not be fitted. Thiem also needs steady drawdown, which is reached only when a recharge source balances pumping.

6. Groundwater exploration

Exploration proceeds from cheap to costly: a desk study of geology, geomorphology, rainfall and existing wells; remote sensing and GIS (lineaments, drainage density, geomorphic units and land use combined by weighted overlay into potential zones); field hydrogeology (well inventory, water-level contours, springs, quality); geophysics; test drilling with logging; and finally a pumping test. In alluvium and valley fills the targets are thick sand and gravel and buried channels. In karst they are solution-widened joints and conduits. In basalt they are vesicular and jointed zones, weathered tops and inter-flow beds. In crystalline hard rock (granite-gneiss) the water is held in the weathered mantle and the fractures below it, so favourable sites are thick weathering in topographic lows, lineament intersections and fracture zones, while fresh massive rock is dry and a dyke may act as a barrier.

Geophysical methods for groundwater
MethodWhat it sensesGroundwater use
Electrical resistivity (VES)Resistivity layering with depthDepth to aquifer and bedrock, thickness of weathering, fresh-saline contrast
Seismic refractionP-wave velocity layeringWater table and bedrock in unconsolidated cover
Electromagnetic and VLFLateral conductivity changesFracture and shear zones, clay lenses
Gravity and magneticDensity and susceptibility contrastsBuried valleys and basins; dykes and faults

In a vertical electrical sounding the current electrodes are spread symmetrically about a fixed centre. Apparent resistivity is ρ_a = 2πa ΔV/I for the Wenner array (spacing a) and ρ_a = π(L² − l²)/(2l) · ΔV/I for the Schlumberger array (L = AB/2, l = MN/2). Archie’s law links resistivity to the pores of a clean saturated formation: F = ρ₀/ρ_w = a φ^(−m), with a ≈ 1 and m ≈ 2 for clean sands (use the values a question states). For φ = 0.25 and water resistivity 20 Ωm, F = 16 and ρ₀ = 320 Ωm. Resistivity falls with higher porosity, saline water and clay content. Borehole logs complete the picture: resistivity and spontaneous-potential logs identify beds, and the gamma log marks clay.

⚠️ Low resistivity is not a good aquifer
A conductive layer may be a saline-water aquifer or a clay, which holds water but yields almost none. Equivalent models also fit the same sounding curve. Resolve the ambiguity with borehole logs, lithology and water-quality data, never from resistivity alone.

Key takeaways

  • Porosity n = V_v/V, with n = S_y + S_r and n_e ≤ n; K = kρg/μ. Layered ground has K_h = ΣK_i b_i/Σb_i and K_v = Σb_i/Σ(b_i/K_i). High porosity does not mean high permeability.
  • Water released is V = S A Δh: unconfined S ≈ S_y (0.01–0.3) drains pores, confined S (about 10⁻⁵–10⁻³) comes from elastic storage. T = Kb, and recharge by water-table fluctuation is S_y Δh plus draft.
  • Darcy: q = Ki and Q = KiA; flow follows total head h = z + p/ρg. Seepage velocity v_s = q/n_e. The law holds for laminar flow (Re below about 1–10).
  • Flow refracts as tan θ₁/tan θ₂ = K₁/K₂; transient flow follows ∂²h/∂x² + ∂²h/∂y² = (S/T)∂h/∂t. Recharge areas have downward and discharge areas upward head gradients.
  • Thiem (confined): s₁ − s₂ = [Q/(2πT)] ln(r₂/r₁). Theis: s = QW(u)/(4πT), u = r²S/(4Tt). Cooper-Jacob (u < 0.05): T = 2.303Q/(4πΔs), S = 2.25Tt₀/r². An impermeable boundary doubles the semi-log slope; a recharge boundary flattens it.
  • Explore by geology, geomorphology, lineaments and GIS, then geophysics (VES: ρ_a = 2πaΔV/I for Wenner) and logs. Archie F = ρ₀/ρ_w = aφ^(−m); low resistivity may be clay or saline water, not a good aquifer.

Practice questions (18)

Attempt each one before opening the answer. Every explanation names the tempting wrong option as well as the right one, because that is where marks are lost.

  1. A bed of clay and a bed of well-sorted coarse sand are compared. Which statement is correct?

    1. The clay can have the higher total porosity yet a far lower conductivity, because its pores are tiny and hold water tightly
    2. The clay has the higher hydraulic conductivity, because its larger total porosity gives more pore space for water to flow through
    3. The sand has the higher specific retention, because coarse grains hold more water against gravity than fine grains do
    4. The two beds have equal hydraulic conductivity, because conductivity depends only on porosity and not on pore size
    Show answer

    Answer: A — The clay can have the higher total porosity yet a far lower conductivity, because its pores are tiny and hold water tightly

    Total porosity of clay (often 40 per cent or more) can exceed that of sand, but the pores are tiny and water is bound to the grains, so K is smaller by many orders of magnitude. Conductivity depends on pore size and connection, not on porosity alone, so the claim of higher K for clay is wrong. Fine grains hold more water against gravity, so clay, not sand, has the larger specific retention.
  2. A saturated sand has a porosity of 0.38. After free drainage under gravity it still holds water equal to 9 per cent of its total volume. Its specific yield is

    1. 0.09, the volume held by molecular attraction and capillarity
    2. 0.47, the porosity plus the water that stays in the sand
    3. 0.29, the porosity minus the specific retention
    4. 0.38, equal to the total porosity of the sand
    Show answer

    Answer: C — 0.29, the porosity minus the specific retention

    Porosity splits as n = S_y + S_r. The water that stays after drainage (0.09) is the specific retention, so S_y = 0.38 − 0.09 = 0.29. The value 0.09 is the retention itself, 0.38 would mean that no water is retained at all, and adding the two quantities breaks the relation n = S_y + S_r.
  3. When the head in a confined aquifer falls, the water released from storage comes mainly from

    1. gravity drainage of the pores, as in an unconfined aquifer, so that the aquifer progressively desaturates
    2. evaporation of water from the capillary fringe lying above the piezometric surface of the aquifer
    3. downward drainage of the overlying unsaturated zone into the top of the confined aquifer
    4. expansion of the water and compaction of the skeleton as pore pressure falls, the aquifer staying saturated
    Show answer

    Answer: D — expansion of the water and compaction of the skeleton as pore pressure falls, the aquifer staying saturated

    In a confined aquifer the pores stay full: a head fall lowers the pore pressure, so the water expands slightly and the skeleton compacts, which defines the small storativity S = S_s b. Gravity drainage of pores, with specific yield, belongs to unconfined aquifers. Evaporation from a capillary fringe and drainage from above do not describe release from confined storage.
  4. Two piezometers at the same site have their screens at elevations of 90 m and 70 m above a datum. Water stands at 98 m above the datum in the shallow piezometer and at 95 m above the datum in the deep one. The vertical flow between the two screens is

    1. downward, because total head is 98 m at the shallow screen and 95 m at the deep screen, so head falls with depth
    2. upward, because the pressure head is 25 m at the deep screen against 8 m at the shallow one, and water moves towards lower pressure
    3. zero, because both piezometers sit in the same vertical column and so share the same total head
    4. upward, because the deeper screen lies at the lower elevation, and groundwater always moves from low elevation to high elevation
    Show answer

    Answer: A — downward, because total head is 98 m at the shallow screen and 95 m at the deep screen, so head falls with depth

    Total head is the water level above the datum: 98 m at the shallow screen and 95 m at the deep one. Flow goes from high to low total head, so it is downward, a recharge setting. The pressure heads (8 m and 25 m) are a distractor: pressure is higher at depth, yet flow does not follow pressure. Elevation alone does not set the direction either, and different total heads mean there is flow.
  5. Which statement about the Darcy flux q and the seepage velocity v_s in a saturated sand is correct?

    1. v_s = q·n_e, so the pore water moves more slowly than the Darcy flux because the grains obstruct it
    2. v_s = q/n_e, so pore water moves faster than q suggests, as flow uses only part of the cross-section
    3. v_s equals q, because the Darcy flux is the average speed at which the water moves through the pores
    4. q = v_s/n_e, so the Darcy flux is larger than the actual speed of the water in the pores
    Show answer

    Answer: B — v_s = q/n_e, so pore water moves faster than q suggests, as flow uses only part of the cross-section

    The Darcy flux is the discharge divided by the whole cross-sectional area, but water can pass only through the effective pores, so the true average speed is v_s = q/n_e, which exceeds q because n_e < 1. Travel times and tracer speeds therefore use v_s. The flux is not the pore speed, and the relations that multiply by n_e invert the geometry of the flow.
  6. Darcy’s linear law is least likely to hold in which of the following situations?

    1. Slow laminar flow through fine sand under a small natural hydraulic gradient of about 0.001
    2. Fast flow through coarse gravel or an open fracture near a pumped well, at Reynolds number above 10
    3. Slow seepage of water through a thick silt layer under a small difference of head across it
    4. Regional laminar flow in a sandstone aquifer at a speed of only a few centimetres per day
    Show answer

    Answer: B — Fast flow through coarse gravel or an open fracture near a pumped well, at Reynolds number above 10

    Darcy’s law assumes laminar, viscous-dominated flow, valid up to a Reynolds number of roughly 1 to 10. Near a pumped well in gravel or an open fracture, velocities are high enough for inertial and turbulent losses, so the relation between flux and gradient becomes non-linear. Slow flow in fine sand, silt or a regional sandstone aquifer is the standard domain of the law.
  7. In a pumping test on a confined aquifer, drawdown plotted against log time first follows one straight line and later a second straight line with twice the slope. This is best explained by

    1. a perennial river acting as a recharge boundary, equivalent to an image injection well
    2. a single impermeable boundary, equivalent to an image well pumping at the same rate
    3. a leaky overlying aquitard that begins to supply water to the aquifer
    4. the aquifer becoming unconfined at late times, so that the storativity rises towards the specific yield
    Show answer

    Answer: B — a single impermeable boundary, equivalent to an image well pumping at the same rate

    An impermeable boundary reflects the cone of depression, which is equivalent to a second well of the same sign pumping at the same rate, so after its effect arrives the drawdown rate doubles and the semi-log slope doubles. A recharge boundary or a leaky aquitard supplies water and flattens the curve, and a delayed-yield effect flattens the middle part rather than doubling the final slope.
  8. A fully penetrating well pumps 1500 m³/d from a confined aquifer of transmissivity 300 m²/d. Assume steady radial flow. The difference in drawdown between observation wells at 20 m and 200 m from the well is closest to

    1. 0.80 m, from using log₁₀(200/20) = 1 in place of the natural logarithm
    2. 0.92 m, from using 4π in place of 2π in the denominator
    3. 1.83 m, from [Q/(2πT)] ln(200/20)
    4. 11.5 m, from Q/T multiplied by ln(200/20), omitting 2π
    Show answer

    Answer: C — 1.83 m, from [Q/(2πT)] ln(200/20)

    Thiem’s equation gives s₁ − s₂ = [Q/(2πT)] ln(r₂/r₁) = [1500/(2π × 300)] × ln 10 = 0.7958 × 2.3026 = 1.83 m. The value 0.80 uses a base-10 logarithm, which is wrong in this form of the equation, 0.92 uses the Theis denominator 4π, and 11.5 forgets the 2π.
  9. A vertical electrical sounding in an alluvial plain shows a layer of low resistivity between 25 m and 60 m depth. Which interpretation is most defensible before drilling?

    1. It is a fresh-water sand aquifer, because water always lowers resistivity, so fresh-water sand has the lowest resistivity of any layer
    2. It is dry sand, because air-filled pores conduct electric current better than water-filled pores
    3. It is unweathered crystalline bedrock, because such rock has lower resistivity than saturated sediments
    4. It may be clay or sand saturated with saline water, so logs, lithology and water-quality data are needed to tell them apart
    Show answer

    Answer: D — It may be clay or sand saturated with saline water, so logs, lithology and water-quality data are needed to tell them apart

    Both clay and saline-water-saturated sand have low resistivity, and only one of them is a useful fresh-water aquifer, so the sounding cannot decide by itself. Fresh-water sand has moderate to high resistivity, not the lowest. Dry sand is highly resistive, not conductive, and fresh unweathered crystalline rock is far more resistive than saturated sediments.
  10. In a granite-gneiss terrain, which site is the best first target for a bore well?

    1. A hill crest where fresh massive granite is exposed and no joints or weathering are visible
    2. A flat area with fresh massive gneiss at the surface and no lineament mappable on the imagery
    3. A steep scarp with thin soil over massive rock, where runoff is high and infiltration is poor
    4. A topographic low at a lineament intersection, over thick weathering on fractured rock
    Show answer

    Answer: D — A topographic low at a lineament intersection, over thick weathering on fractured rock

    In crystalline rock the water is held in the weathered mantle and in fractures below it, so the best sites are lows with thick weathering and intersecting lineaments, which concentrate recharge and fracturing. Fresh massive rock, a hill crest or a steep scarp with thin soil offers little storage, little infiltration and few open fractures.
  11. Which of the following statements about aquifer storage properties are correct?

    1. For the same head decline an unconfined aquifer releases far more water per unit area than a confined one, because S_y is much larger than confined storativity
    2. Specific retention is the volume of water that drains freely from a unit volume of the formation under gravity
    3. Transmissivity T = Kb has the dimensions of length squared per unit time
    4. Storativity of a confined aquifer is expressed in units of inverse metres
    Show answer

    Answer: A — For the same head decline an unconfined aquifer releases far more water per unit area than a confined one, because S_y is much larger than confined storativity; C — Transmissivity T = Kb has the dimensions of length squared per unit time

    S_y (0.01 to 0.3) is orders of magnitude larger than confined S (about 10⁻⁵ to 10⁻³), and T = K × b has units m²/d. Specific retention is the water held against gravity, which is the opposite of what drains freely, and storativity S = S_s b is dimensionless; the quantity with units of inverse length is the specific storage S_s.
  12. Which of the following statements about groundwater flow are correct?

    1. Groundwater always flows from a point of higher pressure to a point of lower pressure, whatever their elevations
    2. For flow parallel to layers the effective conductivity is the thickness-weighted arithmetic mean, and for flow across the layers it is the harmonic mean
    3. At a boundary between two media, flow lines refract so that tan θ₁/tan θ₂ = K₁/K₂, with the angles measured from the normal
    4. The seepage velocity is always smaller than the Darcy flux, because the pore water is slowed by the grains
    Show answer

    Answer: B — For flow parallel to layers the effective conductivity is the thickness-weighted arithmetic mean, and for flow across the layers it is the harmonic mean; C — At a boundary between two media, flow lines refract so that tan θ₁/tan θ₂ = K₁/K₂, with the angles measured from the normal

    The thickness-weighted arithmetic mean (parallel flow) and harmonic mean (flow across layers) are the standard layered-medium results, and the tangent law of refraction is tan θ₁/tan θ₂ = K₁/K₂. Water moves down total head, so a deeper point may have the higher pressure and still receive no flow from the shallower one. Seepage velocity v_s = q/n_e is larger than the Darcy flux, not smaller.
  13. Which of the following statements about well hydraulics are correct?

    1. The Theis solution assumes a confined, homogeneous, isotropic aquifer of infinite extent and a fully penetrating well of negligible diameter
    2. The Cooper-Jacob straight-line method is valid when u is large, at the earliest times after pumping begins
    3. Well loss is the part of the drawdown proportional to Q, whereas aquifer loss is proportional to Q²
    4. A recharge boundary such as a perennial river acts like an image injection well and makes the drawdown flatten
    Show answer

    Answer: A — The Theis solution assumes a confined, homogeneous, isotropic aquifer of infinite extent and a fully penetrating well of negligible diameter; D — A recharge boundary such as a perennial river acts like an image injection well and makes the drawdown flatten

    The listed Theis assumptions are the standard ones, and a recharge boundary behaves as an image injection well so drawdown flattens. The Cooper-Jacob line is valid for small u (later times), not large u. In s_w = BQ + CQ², the linear term BQ is the aquifer loss and the quadratic term CQ² is the well loss, the reverse of the statement given.
  14. Two observation wells 300 m apart in a sand aquifer have heads of 52.4 m and 51.2 m. The hydraulic conductivity is 15 m/d and the effective porosity is 0.20. The time taken by groundwater to travel 600 m along the flow line at the seepage velocity is ____ days.

    Numerical answer — type the value.

    Show answer

    Answer: 2000

    Gradient i = (52.4 − 51.2)/300 = 0.004. Darcy flux q = K i = 15 × 0.004 = 0.06 m/d. Seepage velocity v_s = q/n_e = 0.06/0.20 = 0.30 m/d. Time = 600/0.30 = 2000 days. Using q instead of v_s would give 10 000 days.
  15. An unconfined aquifer of area 25 km² has specific yield 0.08. Between the pre-monsoon and post-monsoon surveys the water table rose by 2.4 m. This rise is net of a seasonal pumping draft of 1.3 million m³ taken from the aquifer in the same period; other inflows and outflows are negligible. By the water-table fluctuation method, gross recharge (net storage gain plus draft) is ____ million m³, to one decimal place.

    Numerical answer — type the value.

    Show answer

    Answer: 6.1

    Net storage gain = S_y × Δh × A = 0.08 × 2.4 × 25 × 10⁶ = 4.8 × 10⁶ m³ = 4.8 million m³. Because that rise is after the draft was withdrawn, gross recharge = 4.8 + 1.3 = 6.1 million m³. Reporting 4.8 would understate recharge by the amount pumped out.
  16. Three horizontal layers lie one above another: 4 m thick with K = 12 m/d, 2 m thick with K = 0.3 m/d and 4 m thick with K = 6 m/d. The effective vertical hydraulic conductivity of the sequence, to two decimal places, is ____ m/d.

    Numerical answer — type the value.

    Show answer

    Answer: 1.30

    For flow across layers K_v = Σb_i/Σ(b_i/K_i). Total thickness = 10 m. Σ(b_i/K_i) = 4/12 + 2/0.3 + 4/6 = 0.3333 + 6.6667 + 0.6667 = 7.6667 d. K_v = 10/7.6667 = 1.30 m/d. The arithmetic mean, 7.26 m/d, is the horizontal conductivity and is wrong for vertical flow, where the thin 0.3 m/d layer controls.
  17. Cooper-Jacob analysis of a pumping test at Q = 2400 m³/d gives a drawdown of Δs = 0.8 m per log₁₀ cycle of time. The straight line cuts zero drawdown at t₀ = 0.002 d in an observation well at r = 30 m. Use T = 2.303Q/(4πΔs) and S = 2.25Tt₀/r². The storativity S, expressed as a multiple of 10⁻³ and to two decimal places, is ____ × 10⁻³.

    Numerical answer — type the value.

    Show answer

    Answer: 2.75

    T = 2.303 × 2400/(4π × 0.8) = 5527.2/10.053 = 549.8 m²/d. S = 2.25 × 549.8 × 0.002/30² = 2.474/900 = 2.749 × 10⁻³, i.e. 2.75 × 10⁻³. The slope gives T and only then does the intercept give S; using r² = 30 instead of 900 would give a value 30 times too large.
  18. A clean saturated sand has porosity 0.20 and pore-water resistivity 15 Ωm. Use Archie’s law ρ₀ = ρ_w a φ^(−m) with a = 1 and m = 2. The bulk resistivity of the sand is ____ Ωm.

    Numerical answer — type the value.

    Show answer

    Answer: 375

    Formation factor F = a φ^(−m) = 0.20^(−2) = 25. Bulk resistivity ρ₀ = F × ρ_w = 25 × 15 = 375 Ωm. Using φ^m instead of φ^(−m) would give 0.6 Ωm, which is far below the pore-water resistivity and therefore impossible for a rock that also contains insulating grains.