πŸ’§ IGW-NET Β· Quick Tutorial 20 of 31

Tutorial 20: Theis Well Solution

Apply the Theis analytical solution for confined aquifer pumping tests. Fit drawdown data to estimate transmissivity and storativity.

IGW-NET Tutorial 20 Prereq: MAGNET4WATER account 4 sections

This tutorial covers

  1. Part 1: Steady-State Net Drawdown
  2. Part 1b: Steady-State Superimposed Heads
  3. Part 2: Transient Drawdown
  4. What's Next

1Part 1: Steady-State Net Drawdown

Step 1 β€” Enter Synthetic Mode

Navigate to IGW-NET and log in. Use Go to Synthetic 'Go to Synthetic Case Area' (under Utilities) to create a synthetic model domain.

Synthetic mode β€” blank rectangular domain ready for the Theis well experiment
Synthetic mode β€” blank rectangular domain ready for the Theis well experiment

Step 2 β€” Add River Boundaries

Click SaveShape ZoneRect to add rectangular zones along the left and right edges of the domain to represent river features. The head difference between them drives regional groundwater flow:

Left edge: prescribed head = 0 m
Right edge: prescribed head = βˆ’2 m

Step 3 β€” Add a Theis Well

Click DrawWell to add a well near the center of the domain. In the Well Input Options interface:

1. Check the box next to Theis 'Theis' β€” this tells IGW-NET to compute the analytical Theis solution for this well
2. Use the default pumping rate: βˆ’2500 mΒ³/day

The Theis checkbox is the key feature β€” it enables side-by-side comparison of the numerical solution with the exact analytical answer.

Step 4 β€” Set Initial Head and Aquifer Properties

Set the initial head to 0 m (Conceptual Model Tools β†’ DomainAttr β†’ Simulation Settings β†’ Initial & Boundary Condition for Head). Use the default aquifer properties:

Aquifer thickness: 50 m
Hydraulic conductivity: 75 ft/day
Specific storage: 0.00001 1/m
Mode: Steady-state (default)

Step 5 β€” Run the Net Drawdown Simulation

Click Simulate 'SIMULATE'. Two prompts appear:

1. "Theis wells found in your model..." β†’ Click 'OK' to confirm
2. Superposition vs. net drawdown prompt β†’ Click 'Cancel' to compute net drawdown only

Head contours and velocity vectors show the classic cone of depression around the well β€” drawdown relative to the initial head of zero.

Steady-state net drawdown β€” the classic cone of depression around the pumping well
Steady-state net drawdown β€” the classic cone of depression around the pumping well

Step 6 β€” View 3D and Cross-Section Results

Click Display Charts 'Display Charts' (Analysis Tools β†’ Analysis β†’ Display Charts) to see the results in 3D and cross-section view. In the Cross-section Plot, uncheck 'Bot Elevation' and click 'ReDraw' to focus on the drawdown shape near the well.

3D view of the cone of depression β€” the Theis drawdown rendered as a surface
3D view of the cone of depression β€” the Theis drawdown rendered as a surface
Cross-section through the cone of depression β€” drawdown profile at the well
Cross-section through the cone of depression β€” drawdown profile at the well

2Part 1b: Steady-State Superimposed Heads

Step 7 β€” Set Background Head for Superposition

Change the aquifer top elevation to βˆ’10 m (DomainAttr β†’ Aquifer Properties β†’ Top Elevation). This establishes a non-zero background head distribution β€” the regional flow field driven by the two river boundaries.

Step 8 β€” Run the Superimposed Simulation

Click Simulate 'SIMULATE' again. When the superposition prompt appears, click 'OK' this time. IGW-NET now computes: background heads (from regional flow between the two rivers) minus Theis drawdown = the actual head distribution during pumping.

Superimposed heads β€” the cone of depression overlaid on the regional flow field
Superimposed heads β€” the cone of depression overlaid on the regional flow field
Superimposed head distribution in map view β€” showing how pumping distorts the regional flow pattern
Superimposed head distribution in map view β€” showing how pumping distorts the regional flow pattern

What Superposition Reveals

The regional flow field is disrupted: Without pumping, water flows uniformly from the left river (h = 0 m) to the right river (h = βˆ’2 m). With pumping, the cone of depression bends flow lines toward the well. Some water that would have reached the right river is now captured by the well. The superimposed view shows this interaction β€” something net drawdown alone cannot reveal.

Capture zone implications: The distortion of regional flow defines the well's capture zone β€” the area that contributes water to the well. In this simple case, you can see how the well "steals" water from the regional flow. In Tutorial 3 (Particle Tracking) and Tutorial 18 (Probabilistic Capture Zones), you explored this concept in detail.

3Part 2: Transient Drawdown

Why Go Transient?

Steady-state shows the end state: Given infinite time, the drawdown reaches a fixed shape β€” the cone of depression stabilizes when inflow from boundaries balances the pumping rate. But how long does it take to reach steady state? How does drawdown evolve in the first hours, days, and weeks?

Transient shows the process: The Theis solution is inherently transient β€” drawdown grows logarithmically with time. Early-time behavior is dominated by storage release (water squeezed from the aquifer matrix). Late-time behavior approaches steady state as boundary effects arrive. The transition between these regimes is exactly what pumping tests measure β€” and what the Theis solution predicts.

Monitoring wells see the evolution: A monitoring well near the pumping well records drawdown over time β€” the classic pumping test hydrograph. Matching this curve to the Theis solution yields T and S β€” the most common method of aquifer characterization in practice.

Step 9 β€” Enable Transient Flow

Continue from the Part 1 model. Go to DomainAttr β†’ Simulation Settings. Check Transient 'Modeling Transient Flow'. Also check 'Overwrite with Steady State Solution at t=0'. Click 'Save'. The model will now simulate drawdown evolution over time.

Simulation Settings β€” enabling transient flow with steady-state initial condition
Simulation Settings β€” enabling transient flow with steady-state initial condition

Step 10 β€” Add a Monitoring Well

Click DrawWell to add a second well near the pumping well. Check 'Monitoring Well' and set the pumping rate to zero. This well records drawdown over time β€” it's your virtual piezometer for observing the transient response.

Step 11 β€” Run and View Transient Results

Click Simulate 'SIMULATE'. Click 'Cancel' on the superposition prompt to compute net drawdown. The simulation now runs through time β€” watch the cone of depression grow at each time step. Use Display Charts 'Display Charts' to view transient results in 3D and cross-section views.

Transient drawdown β€” the cone of depression evolving over time
Transient drawdown β€” the cone of depression evolving over time
Transient results β€” drawdown distribution at a later time step, approaching steady state
Transient results β€” drawdown distribution at a later time step, approaching steady state

Key Concepts

Verification vs. validation: This tutorial demonstrates verification β€” does the numerical solver reproduce a known analytical answer? This is different from validation (Tutorial 8: Calibration), which asks whether the model reproduces field observations. Verification proves the code works correctly. Validation proves the model represents the real system. Both are necessary for a credible model.

The Theis checkbox is powerful: In traditional workflows, comparing numerical results to analytical solutions requires external scripts, spreadsheet calculations, or separate software. IGW-NET's built-in Theis feature eliminates this β€” one checkbox enables the analytical benchmark alongside the numerical solution. Verification becomes effortless.

From verification to practice: Once you've confirmed the numerical engine reproduces the Theis solution, you can confidently move to complex problems where no analytical solution exists β€” heterogeneous aquifers, irregular boundaries, multiple wells, contaminant transport. The verification gives you the foundation to trust the numerical results.

Superposition principle: The linearity of the confined-aquifer flow equation means drawdowns from multiple wells can be added together, and drawdown can be superimposed on any background flow field. This powerful principle breaks down for unconfined aquifers (nonlinear) and when wells interact with boundaries β€” but it remains a cornerstone of groundwater hydraulics and wellfield design.

4What's Next

Continue to advanced mesh and analysis capabilities:

Tutorial 21: Unstructured Grid (Create) β€” flexible mesh for complex geometries
Tutorial 22: Unstructured Grid (Results) β€” visualize results on unstructured grids
Tutorial 23: MODFLOW Analysis Tool β€” advanced post-processing of MODFLOW output