Effective permeability ($k_{\mathrm{eff}}$)

Effective permeability is a bulk permeability of a porous medium that reproduces the total Darcy flux through a heterogeneous sample under specified boundary conditions. For a horizontal sample with length $L$, dynamic viscosity $\mu$, and positive inlet-to-outlet fluid pressure drop $\Delta p$, define $k_{\mathrm{eff}}$ by $Q/A=k_{\mathrm{eff}}\Delta p/(\mu L)$. The description is useful when individual layers or passages are smaller than the scale of the calculation. It depends on direction, geometry and boundary conditions; it is not generally the arithmetic mean of local intrinsic permeabilities.

Layered porous permeability illustrates the mechanism: when every path must cross a resistant layer, successive hydraulic resistances add; when more permeable paths span the whole sample alongside less permeable ones, their discharges add. This is the same organizing idea as series and parallel resistors. For fixed fluid properties the corresponding hydraulic conductivities average in the same way.

Consider planar horizontal layers with positive scalar intrinsic intrinsic permeabilities $k_i$, thicknesses $b_i$, and fractions $f_i=b_i/{\sum }_j{b}_j$. Assume each layer is uniform and isotropic, fully saturated by the same Newtonian fluid of constant mass density and dynamic viscosity, with perfect hydraulic contact and steady incompressible flow. Use boundaries that impose pure flow normal to the layers or a common gradient of hydraulic head along them; additional interface resistance and finite-edge effects are excluded.

For vertical flow normal to the layers, every layer carries the same normal Darcy flux by mass conservation. The hydraulic head drops add. Writing $B={\sum }_i{b}_i$ and taking positive $\Delta h_i$ as the drop along the flow direction givesThis is the thickness-weighted harmonic mean. It is hydraulic head, rather than fluid pressure alone, that accounts correctly for gravity in the vertical calculation. A thin layer of very small permeability of a porous medium can control the whole discharge because all flow must pass through it.

For horizontal flow parallel to the layers, each layer experiences the same gradient of hydraulic head. For width $W$, each area is $A_i=W{b}_i$, and the individual volume fluxes add:This is the thickness-weighted arithmetic mean. A highly permeable layer can carry most of the horizontal flow because it supplies a continuous route from inlet to outlet. Woessner and Poeter derive both cases for hydraulic conductivity; identical fluid properties allow the common factor $\rho g/\mu$ to cancel.^[ref]

For two equally thick layers, take $k_1=k_0$ and $k_2=9k_0$, where $k_0=10^{-12}{\mathrm{m}}^2$. The values 1 and 9 are therefore permeabilities in units of $k_0$, rather than dimensionless intrinsic intrinsic permeabilities. ThenAcross the layers, the less permeable half takes $90\%$ of the total hydraulic head drop and the more permeable half only $10\%$. Along the layers, the more permeable half carries $90\%$ of the total volume flux. At the same magnitude of imposed gradient of hydraulic head, horizontal Darcy flux is $5/1.8=25/9\simeq 2.78$ times the vertical Darcy flux.

The arrangement thus creates physical anisotropy from locally isotropic ingredients. In coordinates with $x,y$ along the layers and $z$ normal to them, the effective permeability of a porous medium tensor is ${\mathbf{k}}_{\mathrm{eff}}=\mathrm{diag}(5,5,1.8)k_0$. A driving gradient of hydraulic head with both horizontal and vertical components generally produces a discharge in a different direction. This distinction matters when deciding whether a geological bed conducts aquifer flow along itself or acts as a barrier to flow across it.

Regenerate the diagram by running python ../layered-porous-permeability.py from _media; it writes its PNG into the current directory.

At small areal number density $\varphi$, identical aligned cracks contribute independently the fixed-throughflow dissipation reduction of a conductive crack. Equating the mean dissipation density to $\mu U^2/k_{xx}^{\ast }$ gives the displayed first-order enhancement. Writing $k/[1-\pi \varphi a^2\alpha /(1+\alpha )]$ is an algebraic noninteraction estimate, not a claim of accuracy at higher density. The enhancement saturates as crack conductance becomes large compared with transport through the surrounding porous matrix.

In a slender layer, two parallel sublayers of fractions $x/L$ and $1-x/L$ have local effective permeability $K(x)=k_{1}x/L+k_2(1-x/L)$. Integrating their successive hydraulic resistance gives the logarithmic mean displayed above. A leading pressure gradient uniform across the thickness does not make parcel paths horizontal; a small transverse velocity preserves mass conservation as the cross-sectional conductance changes.