12.3 Ground water hydrology

Achievable FE Civil

12. Water resources engineering

Our FE Civil course is currently in development and is a work-in-progress.

Ground water hydrology

3 min read

Font

Discuss

Feedback

This chapter covers the following:

Darcy’s law
Unconfined aquifer - Dupuit’s formula
Confined aquifer - Theim equation

Darcy’s law

Darcy’s law is the core relationship used to describe groundwater flow through saturated porous media. It says that the volumetric flow rate is proportional to the hydraulic gradient, with hydraulic conductivity capturing how easily water moves through the material (based on both fluid properties and the permeability of the soil or rock).

Darcy’s law applies when flow is laminar, which is typical in most groundwater systems. You’ll use it as the starting point for analyzing seepage, aquifer behavior, and contaminant transport.

$Q = - K A (\frac{d h}{d x})$

$Q$ = discharge rate (ft³/sec or m³/s)
$K$ = hydraulic conductivity (ft/sec or m/s)
$h$ = hydraulic head (ft or m)
$A$ = cross-sectional area of flow (ft² or m²)

The negative sign indicates that flow goes from higher head to lower head.

Specific discharge:

$q = - K (\frac{d h}{d x})$

Average seepage velocity:

$v = \frac{q}{n} = \frac{- K}{n} (\frac{d h}{d x})$

$n$ = effective porosity

Unconfined aquifer - Dupuit’s formula

Dupuit’s formula is used for steady-state groundwater flow in an unconfined aquifer, where the top of the saturated zone is the water table. The approach relies on the Dupuit assumptions:

Vertical flow components are neglected.
The hydraulic gradient is approximated by the slope of the water table.

These assumptions simplify the flow field, but the result is often accurate enough for shallow, laterally extensive unconfined aquifers. You can use the equation to estimate discharge to a well, drawdown, or the shape of the water table.

$Q = \frac{π K ( h _{1}^{2} - h _{2}^{2} )}{ln ( \frac{r _{2}}{r _{1}} )}$

Where:

$Q$ = flowrate of water drawn from well (cfs)
$K$ = hydraulic conductivity (ft/sec)
$h_{1}$ = water height at perimeter (ft)
$h_{2}$ = water height at radius $r_{2}$ (ft)
$r_{1}$ = radius of well (ft)
$r_{2}$ = distance from well centerline to $h_{2}$ point (ft)
$Q / D_{w}$ = specific capacity
$D_{w}$ = well drawdown (ft)

Confined aquifer - Theim equation

The Thiem equation is a classical steady-state solution for radial flow toward a pumping well in a confined aquifer. It connects the pumping rate to aquifer properties and shows that hydraulic head changes logarithmically with distance from the well.

The equation assumes:

The aquifer is homogeneous and isotropic.
The well fully penetrates the confined aquifer.
Flow is steady (equilibrium conditions).

These conditions make the Thiem equation especially useful for interpreting pumping test data and estimating transmissivity.

$Q = \frac{2 π T ( h _{1} - h _{2} )}{ln ( \frac{r _{2}}{r _{1}} )}$

Where:

$T = K b$ = transmissivity (ft²/sec)
$b$ = thickness of confined aquifer (ft)
$h_{1}, h_{2}$ = head at radius $r_{1}, r_{2}$ (ft)
$r_{1}, r_{2}$ = radii from pumping well (ft)
$H$ = height of piezometric surface prior to pumping (ft)

Transmissivity, $T$

Transmissivity is the product of hydraulic conductivity and the confined aquifer thickness, $b$ :

$T = K b [L^{2} / T]$

Storativity or storage coefficient of an aquifer, $S$

Storativity (storage coefficient) is the volume of water taken into or released from storage per unit surface area per unit change in potentiometric (piezometric) head.

Darcy’s law

Governs groundwater flow in saturated porous media
Key formula: $Q = - K A (\frac{d h}{d x})$
- $Q$ : discharge, $K$ : hydraulic conductivity, $A$ : area, $d h / d x$ : hydraulic gradient
Specific discharge: $q = - K (\frac{d h}{d x})$
Seepage velocity: $v = \frac{q}{n} = \frac{- K}{n} (\frac{d h}{d x})$
- $n$ : effective porosity

Unconfined aquifer - Dupuit’s formula

Used for steady-state flow in unconfined aquifers
Assumes:
- Negligible vertical flow
- Hydraulic gradient ≈ water table slope
Key formula: $Q = \frac{π K ( h _{1}^{2} - h _{2}^{2} )}{l n ( \frac{r _{2}}{r _{1}} )}$
- $h_{1}$ , $h_{2}$ : water heights, $r_{1}$ , $r_{2}$ : radii

Confined aquifer - Thiem equation

Describes steady-state radial flow to a well in confined aquifer
Assumes:
- Homogeneous, isotropic aquifer
- Well fully penetrates aquifer
- Steady flow
Key formula: $Q = \frac{2 π T ( h _{1} - h _{2} )}{l n ( \frac{r _{2}}{r _{1}} )}$
- $T$ : transmissivity, $h_{1}$ , $h_{2}$ : heads, $r_{1}$ , $r_{2}$ : radii

Transmissivity, $T$

$T = K b$ (hydraulic conductivity × aquifer thickness)
Units: $[L^{2} / T]$ (e.g., ft²/sec)

Storativity (storage coefficient), $S$

Volume of water released or stored per unit area per unit head change
Important for aquifer response to pumping

Next | 13.1 Water quality

Ground water hydrology

This chapter covers the following:

Darcy’s law
Unconfined aquifer - Dupuit’s formula
Confined aquifer - Theim equation

Darcy’s law

Darcy’s law applies when flow is laminar, which is typical in most groundwater systems. You’ll use it as the starting point for analyzing seepage, aquifer behavior, and contaminant transport.

$Q = - K A (\frac{d h}{d x})$

$Q$ = discharge rate (ft³/sec or m³/s)
$K$ = hydraulic conductivity (ft/sec or m/s)
$h$ = hydraulic head (ft or m)
$A$ = cross-sectional area of flow (ft² or m²)

The negative sign indicates that flow goes from higher head to lower head.

Specific discharge:

$q = - K (\frac{d h}{d x})$

Average seepage velocity:

$v = \frac{q}{n} = \frac{- K}{n} (\frac{d h}{d x})$

$n$ = effective porosity

Unconfined aquifer - Dupuit’s formula

Dupuit’s formula is used for steady-state groundwater flow in an unconfined aquifer, where the top of the saturated zone is the water table. The approach relies on the Dupuit assumptions:

Vertical flow components are neglected.
The hydraulic gradient is approximated by the slope of the water table.

$Q = \frac{π K ( h _{1}^{2} - h _{2}^{2} )}{ln ( \frac{r _{2}}{r _{1}} )}$

Where:

$Q$ = flowrate of water drawn from well (cfs)
$K$ = hydraulic conductivity (ft/sec)
$h_{1}$ = water height at perimeter (ft)
$h_{2}$ = water height at radius $r_{2}$ (ft)
$r_{1}$ = radius of well (ft)
$r_{2}$ = distance from well centerline to $h_{2}$ point (ft)
$Q / D_{w}$ = specific capacity
$D_{w}$ = well drawdown (ft)

Confined aquifer - Theim equation

The equation assumes:

The aquifer is homogeneous and isotropic.
The well fully penetrates the confined aquifer.
Flow is steady (equilibrium conditions).

These conditions make the Thiem equation especially useful for interpreting pumping test data and estimating transmissivity.

$Q = \frac{2 π T ( h _{1} - h _{2} )}{ln ( \frac{r _{2}}{r _{1}} )}$

Where:

$T = K b$ = transmissivity (ft²/sec)
$b$ = thickness of confined aquifer (ft)
$h_{1}, h_{2}$ = head at radius $r_{1}, r_{2}$ (ft)
$r_{1}, r_{2}$ = radii from pumping well (ft)
$H$ = height of piezometric surface prior to pumping (ft)

Transmissivity, $T$

Transmissivity is the product of hydraulic conductivity and the confined aquifer thickness, $b$ :

$T = K b [L^{2} / T]$

Storativity or storage coefficient of an aquifer, $S$

Storativity (storage coefficient) is the volume of water taken into or released from storage per unit surface area per unit change in potentiometric (piezometric) head.

Achievable FE Civil

12. Water resources engineering

Our FE Civil course is currently in development and is a work-in-progress.