Theis Recovery Solution for Nonleaky Confined Aquifers
- Data requirements
- Solution options
- Estimated parameters
- Curve matching tips
Related Solution Methods
- Theis (1935) for recovery tests
- Cooper and Jacob (1946)
- Dougherty and Babu (1984)
- Theis (1935) for step tests
Charles Vernon Theis (1900-1987) was the first groundwater hydrologist to develop a rigorous mathematical model of transient flow of water to a pumping well by recognizing the physical analogy between heat flow in solids and groundwater flow in porous media.
The Theis (1935) solution is useful for determining the transmissivity of nonleaky confined aquifers from recovery tests. Analysis involves matching a straight line to residual drawdown data collected after the termination of a pumping test. The solution assumes a line source for the pumped well and therefore neglects wellbore storage.
If your observation data contain both pumping and recovery measurements, you may use the Theis recovery method to analyze just the recovery data.
- aquifer has infinite areal extent
- aquifer is homogeneous, isotropic and of uniform thickness
- control well is fully penetrating
- flow to control well is horizontal
- aquifer is nonleaky confined
- flow is unsteady
- water is released instantaneously from storage with decline of hydraulic head
- diameter of pumping well is very small so that storage in the well can be neglected
- values of are small (i.e., is small and is large)
C.V. Theis, a groundwater hydrologist with the U.S. Geological Survey, derived the following approximate linear equation to predict residual drawdown in a homogeneous, isotropic and nonleaky confined aquifer assuming a fully penetrating line sink that discharged at a constant rate prior to recovery:
- is pumping rate [L³/T]
- is residual drawdown [L]
- is storativity during pumping [dimensionless]
- is storativity during recovery [dimensionless]
- is elapsed time since start of pumping [T]
- is elapsed time since pumping stopped [T]
- is transmissivity [L²/T]
- is the slope of the fitted line (change in residual drawdown per log cycle equivalent time)
is found from the intersection of the line with the axis of the plot. In the absence of boundary effects, should be close to unity. A value of suggests recharge during the test; whereas may indicate a no-flow boundary.
For uninterrupted variable-rate pumping, one may use the principle of superposition in time to compute residual drawdown for constant-rate steps (Birsoy and Summers (1980)):
- is adjusted time [T]
- is pumping rate in the constant-rate step [L³/T]
- is pumping rate in the last constant-rate step prior to recovery [L³/T]
- is residual drawdown [L]
- is elapsed time when the constant-rate step starts [T]
- pumping and observation well locations
- pumping rate(s)
- observation well measurements (time and displacement)
- variable pumping rates
- (ratio of storativity during pumping to storativity during recovery)
Curve Matching Tips
- Choose Match>Visual to perform visual curve matching using the procedure for straight-line solutions.
- Use parameter tweaking to perform visual curve matching and sensitivity analysis.
- Perform visual curve matching prior to automatic estimation to obtain reasonable starting values for the aquifer properties.
Theis, C.V., 1935. The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage, Am. Geophys. Union Trans., vol. 16, pp. 519-524.