P213 Di050016013 Marquee

Advanced Multi-Well Pumping Test Analysis Tool

WELLS offers detailed aquifer drawdown simulations and variable-rate pumping test analysis using analytical methods and parameter estimation integration.

WELLS is a versatile C-based code for analyzing multi-well, variable-rate pumping tests. It simulates drawdown in confined, unconfined, and leaky aquifers, accounting for wellbore storage, complex pumping rates, and boundary conditions.

Download

The software is open source and available on GitLab.

Julia version of WELLS is available on GitHub.

Analytical Solutions

WELLS computes the drawdown/hydraulic head at one or more predefined locations due to pumping at one or multiple locations in an aquifer. The computation is based on a variety of available analytical solutions.

The analytical solutions implemented in WELLS can be broadly classified in the following categories:

Multi-well Variable-Rate Pumping-Test Analysis 

System Geometry

Wells Schematic Confined

The schematic of system geometry for a confined aquifer (with an overlying and underlying aquitard) pumping test is shown in the adjacent figure.

The pumping well of finite radius rw is partially penetrating the confined aquifer of thickness b between depths l and d below the top impermeable boundary.

The pumping well has wellbore storage coefficient Cw (volume of water released from well storage per unit drawdown in it).

The observation well penetrating between depths z1 and z2 below the top impermeable boundary is located at distance r from the axis of the pumping well.

Analytical solutions for confined aquifers

The following analytical solutions are implemented in WELLS for confined aquifers:

  • Theis [1935] Solution : For fully and partially penetrating well ( d=0, l=b ) of zero radius (rw = Cw =0 ) in an isotropic (Kr = Kz) confined aquifer.
  • Hantush [1964] Solution : For partially penetrating well (0 ≤ d,l ≤ b ) of zero radius (rw = Cw =0 ) in an anisotropic (Kr ≠ Kz) confined aquifer.
  • Papadopulos and Cooper [1967] Solution : For fully penetrating well (d=0, l=b ) of finite radius (rw ≠ 0 ) with a finite wellbore storage capacity (Cw ≠ 0 ) in an isotropic (Kr = Kz) confined aquifer solution :
  • Yang et.al. [2006] Solution : For partially penetrating well (0 ≤ d,l ≤ b ) of finite radius (rw ≠ 0 ) with no wellbore storage capacity (Cw = 0 ) in an anisotropic (Kr ≠ Kz) confined aquifer.
  • Mishra et.al. [2011] Solution : For partially penetrating well (0 ≤ d,l ≤ b ) of finite radius (rw ≠ 0 ) with a finite wellbore storage capacity (Cw ≠ 0 ) in an anisotropic (Kr ≠ Kz) confined aquifer.

References

 

Multi-well Variable Rate Pumping Test Analysis

System Geometry

Wells Schematic Unconfined

The schematic of system geometry for unconfined aquifer pumping test is shown in adjacent figure.

The pumping well of finite radius rw is partially penetrating the confined aquifer of initial saturated thickness b between depths l and d below the initial water table.

The pumping well has wellbore storage coefficient Cw (volume of water released from well storage per unit drawdown in it).

The observation well penetrating between depths z1 and z2 above the bottom impermeable boundary is located at distance r from the axis of the pumping well.

Analytical solutions for unconfined aquifer

The following analytical solutions are implemented in WELLS for unconfined aquifers:

  • Transformed Theis [1935] Solution for unconfined aquifer : For fully penetrating well of zero radius
  • Tartakovsky and Neuman [2007] solution : For partially penetrating well of zero radius in unconfined aquifer having infinitely thick unsaturated zone (L → ∞).
  • Mishra and Neuman [2010] solution : For partially penetrating well of zero radius in unconfined aquifer having finitely thick unsaturated zone.
  • Mishra and Neuman [2011] solution : For partially penetrating well of finite radius in unconfined aquifer having finitely thick unsaturated zone.

References

System Geometry

Wells Schematic Leaky Confined

The schematic of system geometry for leaky confined aquifer pumping test is shown in adjacent figure.

The pumping well of finite radius rw is partially penetrating the leaky confined aquifer of thickness b between depths l and d below the top aquitard.

The pumping well has wellbore storage coefficient Cw (volume of water released from well storage per unit drawdown in it).

The observation well penetrating between depths z1 and z2 above the bottom aquitard and is located at distance r from the axis of the pumping well.

Analytical solutions for leaky confined aquifer

The following analytical solutions are implemented in WELLS for leaky confined aquifers:

  • Hantush [1960] Solution : For fully and partially penetrating well ( d=0, l=b ) of zero radius (rw = Cw =0 ) in an isotropic (Kr = Kz) confined aquifer.
  • Mishra and Vessilinov [2011] Solution : For partially penetrating well (0 ≤ d,l ≤ b ) of finite radius (rw ≠ 0 ) with a finite wellbore storage capacity (Cw ≠ 0 ) in an anisotropic (Kr ≠ Kz) confined aquifer.

References

System Geometry

Wells Schematic Leaky Unconfined

The schematic of system geometry for leaky unconfined aquifer pumping test is shown in adjacent figure.

The pumping well of finite radius rw is partially penetrating the leaky unconfined aquifer having initial saturated thickness b between depths l and d below the initial water table.

The pumping well has wellbore storage coefficient Cw (volume of water released from well storage per unit drawdown in it).

The observation well penetrating between depths z1 and z2 above the top boundary of underlying aquitard is located at distance r from the axis of the pumping well.

Analytical solutions for leaky unconfined aquifer

The following analytical solutions are implemented in WELLS for leaky unconfined aquifers:

  • Mishra and Vessilinov [2011] Solution : For partially penetrating well of finite radius with wellbore storage

References

  • Mishra P.K. and V.V. Vessilinov, 2011. Semi-analytical solution for Radial flow towards large diameter well in a leaky confined aquifer, Journal of Hydrology, In review.

Applications

The code was originaly developed at the University of Mining & Geology, Sofia, Bulgaria in 1992 by Velimir V. Vesselinov (Monty).

Contact Us

For inquires please fill out this form