18-Env-A3 Geotechnical and Hydrogeological Engineering · December 2014
Question 6 of 6: Unconfined Aquifer — Well Drawdown & Tracer Travel Time
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
National Exams — December 2014 — 04-Env-A3 / Geotechnical & Hydrogeological Engineering. 3 hours duration; open book exam, any non-communicating calculator permitted. The first five questions as they appear in the answer book are marked (20 marks each, 100 marks total); all six are solved below for completeness.
Reference texts. Braja M. Das, Principles of Geotechnical Engineering (9th ed.) — weight–volume relations, compaction, seepage/flow nets, effective stress and shear strength chapters; Craig & Knappett, Craig's Soil Mechanics (8th ed.) — cross-reference for flow nets, method of fragments and shear strength; Freeze & Cherry, Groundwater (1979) — Darcy's law and the Dupuit–Thiem equation for unconfined radial flow to a well.
Question 6: Unconfined Aquifer — Well Drawdown & Tracer Travel Time (20 marks)
Find. (a) drawdown at $r_1$ and $r_2$; (b) travel time for a conservative tracer from $r_1$ to $r_2$ (or $r_2$ to $r_1$, the direction of flow) at this pumping rate.
Approach. Apply the Dupuit–Forchheimer steady, unconfined, radial form of Darcy's law between the far-field radius of influence and each observation well to get the saturated thickness (hence drawdown) at $r_1$ and $r_2$, then integrate the porosity-corrected (actual) seepage velocity along the radial flow path between the two wells to get the travel time.
Figure 3. Unconfined aquifer with the pumping well and the two observation wells; the true drawdown is millimetre-scale given the very high $K$, so the curve is drawn exaggerated for visibility.
Governing equation. Steady radial (Dupuit) unconfined flow between the far-field boundary and radius $r$ gives
$$Q=\frac{\pi K\left(h_0^2-h(r)^2\right)}{\ln(R/r)}\ \Rightarrow\ h(r)=\sqrt{h_0^2-\frac{Q}{\pi K}\ln\!\left(\frac{R}{r}\right)}.$$
Part (a) — drawdown at $r_1$ and $r_2$. With $K=4320\ \text{m/day}$, $h_0=7\ \text{m}$, $R=4000\ \text{m}$, $Q=100\ \text{m}^3/\text{day}$:
$$h(10)=6.99685\ \text{m}\ \Rightarrow\ s_1=h_0-h(10)=\boxed{3.15\ \text{mm}},$$
$$h(50)=6.99769\ \text{m}\ \Rightarrow\ s_2=h_0-h(50)=\boxed{2.31\ \text{mm}}.$$
Part (b) — tracer travel time. The actual (seepage) velocity at radius $r$ is the Darcy (superficial) velocity divided by porosity, $v(r)=\dfrac{Q}{2\pi r\,h(r)\,n}$, so the travel time between $r_1$ and $r_2$ is
$$t=\int_{r_1}^{r_2}\frac{2\pi n\,r\,h(r)}{Q}\,dr.$$
Evaluating this integral numerically (with $h(r)$ from Step 1) gives
$$t=\boxed{184.7\ \text{days}}\ (\approx6.1\ \text{months}).$$
Check: $K=5\ \text{cm/s}$ is at the very high end of natural aquifer permeability (clean gravel/karst range), which is why a normal municipal-scale pumping rate of 100 m³/day produces only millimetre-scale drawdowns at these radii — the arithmetic is correct, but a real well test at this Q would likely be too small a stress to resolve reliably against seasonal water-table noise, and a much higher Q would normally be used to develop a supply well in an aquifer this permeable.