Home / Science / Earth sciences / Hydrogeology / Well field drawdownTool #129

Well field drawdown calculator: Theis superposition for several wells, boundaries by image wells, and capture zones

Add up the drawdown from several pumping or injection wells that start at different times, with an optional river or no-flow boundary. Draw a contour map at any day, follow one point over time, and find the capture zone of a well in regional flow.

Keep it one tap away for everyday use.

Bookmark Well field drawdown?

Press Ctrl + D now. Your browser opens its own bookmark box, where you choose the name and folder.

  1. Wells
  2. Superpose
  3. Contour map
  4. Capture zone
Nothing is uploadedWorks offline once open

Aquifer and wells

Uses the Theis solution for each well and adds the drawdowns: a confined aquifer, even and endless apart from the one boundary you choose. A boundary is handled by a mirror-image well. Drawdowns larger than about a quarter of the aquifer's thickness break the assumptions.

Drawdown

Capture zone of one well in regional flow

The steady capture zone of a fully penetrating well in an even aquifer with uniform flow. All water inside the line ends up in the well; water outside flows past. The example values are made up.

Download

What it does

When several wells pump from the same aquifer, their cones of depression overlap and the drawdowns add up. A town well field, a dewatering scheme on a building site or a farm with a few boreholes all work this way. This tool adds up the Theis drawdown from each well. Each well can start on a different day, and an injection well is simply a negative rate.

It can also include one straight boundary:

  • A no-flow boundary, such as a rock wall or a sealed fault, makes the drawdown larger.
  • A constant-head boundary, such as a river or a lake that feeds the aquifer, makes it smaller and lets it level off.

The boundary is handled the classic way, with a mirror-image well on the other side of the line.

The second card draws the capture zone of one well in a regional flow: the area whose water ends up in the well. That matters for protecting a drinking-water well, or for designing a pump-and-treat clean-up.

How to use it

  1. Enter the aquifer’s transmissivity T (m²/day), storativity S and the well radius.
  2. Paste or open the wells: name, x (east) and y (north) in metres, the pumping rate in m³/day and the day each starts.
  3. Optional: add a boundary along a north–south or east–west line.
  4. Choose the day for the map and an observation point to follow over time.
  5. For the capture zone, enter the well’s rate, the aquifer thickness, K and the regional gradient.

Worked example

All numbers are made up. In an aquifer with T = 500 m²/day and S = 0.0001, PW1 pumps 800 m³/day from day 0. PW2 starts at 500 m³/day on day 10, and PW3 at 300 m³/day on day 20. On day 30 the drawdown at the point (200, 100) is 2.17 m. Add a river along x = −300 m as a constant-head boundary, and it falls to 0.749 m, because the river feeds the aquifer.

For the capture zone, a well pumping 800 m³/day from a 20 m thick aquifer with K = 25 m/day and a gradient of 0.002 gives a Darcy flux of 0.05 m/day. Its stagnation point is 127 m downstream and the zone is 800 m wide far upstream.

Good to know

  • Theis assumes a confined aquifer of even thickness. For an unconfined aquifer the answers are fair while the drawdown stays small compared with the saturated thickness.
  • Drawdown right at a pumping well also depends on well losses, which this tool leaves out.
  • Real capture zones are shaped by recharge, layering and changing flow, so treat these as first estimates.

FAQ

Can I model recovery when a pump stops? Yes. Add the same well again with a negative rate starting on the day it stops.

How many wells? Up to 50. The map is computed on a 70 × 70 grid, so it stays quick.

Where the numbers come from

  • Theis solution, superposition in space and time, and the method of image wells for straight boundaries: standard hydrogeology, written by fibub from the equations.
  • Capture zone of one well in uniform flow (stagnation point Q/(2πBq), width Q/(Bq)): the classic analytical solution, written from the stream function.
  • Example data are made up by fibub for illustration.

Standard scientific values, such as constants, atomic weights and melting points, are facts found in every textbook; we name where ours were checked. If you are an author or publisher and want a value, a credit or a table changed or removed, tell us and we will act quickly.

Last updated 2026-11-08.

Change log for Well field drawdown

  1. Add the well field drawdown and capture zone tool (#129) Details on GitHub

Every change to fibub

Related tools

Piper diagram maker #125 Paste or open major-ion analyses in mg/L, meq/L or mmol/L. Get a Piper diagram with groups, Stiff diagrams, the ion balance error for each sample, the water type (such as Ca-HCO₃ or Na-Cl), the diamond field and total hardness. Download as PNG, PDF or Word. Pump test analyser #126 Paste or open time–drawdown readings from a pumping test. Get transmissivity T and storativity S from a Theis curve fit and the Cooper-Jacob straight line (with a check that it applies), T from recovery data, and T and S from several observation wells. Log plots, PNG, PDF and Word downloads. Slug test analyser #127 Paste or open water level readings from a slug or bail test. Choose the fit window, get the basic time lag T0 and hydraulic conductivity K by Hvorslev, or Bouwer-Rice with your own ln(Re/R), with a filter pack correction, warnings and a semi-log plot. Darcy flow calculator #128 Groundwater flow from Darcy's law: flux, discharge, seepage velocity and travel time. Find the hydraulic gradient and flow direction from three or more wells, convert K to intrinsic permeability, and estimate K from a sieve curve (d10, Cu, Cc and Hazen). 1D Schrödinger solver #923 Choose a textbook potential, type a formula or draw your own, and get the lowest energy levels and wavefunctions of the one-dimensional Schrödinger equation by the matrix (finite-difference) method, with the exact levels alongside for the box, harmonic oscillator, finite well, Morse, linear and Pöschl–Teller potentials. Adduct and neutral loss calculator #984 Type a formula or a neutral monoisotopic mass and get the m/z of 25 common electrospray adducts in positive and negative mode, the ions left after common neutral losses (water, ammonia, CO2, sugars, phosphate and your own), and match a list of observed m/z values within a ppm tolerance. With no molecule, it turns observed m/z values into candidate neutral masses.