Skip to main content

Seepage in anisotropic soil

Short answer

Transform horizontal dimensions by √(kz/kx) to make the section Laplace‑governed; an equivalent isotropic permeability k'=√(kx kz). Draw a square flow net on the transformed section; seepage q = √(kx kz) h (Nf/Nd).

Derive the relationship for the seepage discharge through anisotropic soil.

In anisotropic soil the permeability differs in the two directions: kxk_x (horizontal, usually larger) and kzk_z (vertical). The flow net cannot be drawn directly with square fields, so the section is transformed.

Governing equation

Continuity for steady flow, ∂vx∂x+∂vz∂z=0\dfrac{\partial v_x}{\partial x} + \dfrac{\partial v_z}{\partial z} = 0, with Darcy's law vx=kx∂h∂xv_x = k_x\dfrac{\partial h}{\partial x} and vz=kz∂h∂zv_z = k_z\dfrac{\partial h}{\partial z} (sign neglected), gives

kx∂2h∂x2+kz∂2h∂z2=0k_x\frac{\partial^2 h}{\partial x^2} + k_z\frac{\partial^2 h}{\partial z^2} = 0

Transformation

Put X=xkzkxX = x\sqrt{\dfrac{k_z}{k_x}}. Then ∂∂x=kzkx∂∂X\dfrac{\partial}{\partial x} = \sqrt{\dfrac{k_z}{k_x}}\dfrac{\partial}{\partial X} and the equation becomes

kz∂2h∂X2+kz∂2h∂z2=0  ⇒  ∂2h∂X2+∂2h∂z2=0k_z\frac{\partial^2 h}{\partial X^2} + k_z\frac{\partial^2 h}{\partial z^2} = 0 \;\Rightarrow\; \frac{\partial^2 h}{\partial X^2} + \frac{\partial^2 h}{\partial z^2} = 0

which is the Laplace equation for an isotropic medium. So the section is redrawn with all horizontal dimensions multiplied by kz/kx\sqrt{k_z/k_x} (vertical dimensions unchanged), and the flow net with square fields is drawn on this transformed section.

Discharge

The flow across a vertical plane in the real section, through a vertical height dzdz:

Δq=vx dz=kx∂h∂xdz=kxkzkx∂h∂Xdz=kxkz ∂h∂X dz\Delta q = v_x\,dz = k_x\frac{\partial h}{\partial x}dz = k_x\sqrt{\frac{k_z}{k_x}}\frac{\partial h}{\partial X}dz = \sqrt{k_xk_z}\,\frac{\partial h}{\partial X}\,dz

This is Darcy's law in the transformed section with the equivalent permeability

k′=kx kzk' = \sqrt{k_x\,k_z}

For a flow net of square fields with total head loss hh:

q=kxkz  h NfNd\boxed{q = \sqrt{k_x k_z}\;h\,\frac{N_f}{N_d}}

The same expression follows by considering the flow across a horizontal plane, so the result is independent of the direction.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Soil Mechanics chapter 7 questions →