Permeability (Earth sciences)

Permeability in fluid mechanics and the Earth sciences (commonly symbolized as k) is a measure of the ability of a porous material (often, a rock or an unconsolidated material) to allow fluids to pass through it.

Symbol used to represent in situ permeability tests in geotechnical drawings

Permeability

Permeability is a property of porous materials that is an indication of the ability for fluids (gas or liquid) to flow through them. Fluids can more easily flow through a material with high permeability than one with low permeability.<ref>"Reading: Porosity and Permeability | Geology". courses.lumenlearning.com. Retrieved 2022-01-14.</ref> The permeability of a medium is related to the porosity, but also to the shapes of the pores in the medium and their level of connectedness.<ref>Fu, Jinlong; Thomas, Hywel R.; Li, Chenfeng (January 2021). "Tortuosity of porous media: Image analysis and physical simulation" (PDF). Earth-Science Reviews. 212: 103439. Bibcode:2021ESRv..21203439F. doi:10.1016/j.earscirev.2020.103439. S2CID 229386129.</ref> Fluid flows can also be influenced in different lithological settings by brittle deformation of rocks in fault zones; the mechanisms by which this occurs are the subject of fault zone hydrogeology.<ref name="faultzone">Bense, V.F.; Gleeson, T.; Loveless, S.E.; Bour, O.; Scibek, J. (2013). "Fault zone hydrogeology". Earth-Science Reviews. 127: 171–192. Bibcode:2013ESRv..127..171B. doi:10.1016/j.earscirev.2013.09.008.</ref> Permeability is also affected by the pressure inside a material.

Units

The SI unit for permeability is the square metre (m²). A practical unit for permeability is the darcy (d), or more commonly the millidarcy (md) (1 d <math>\approx</math> 10⁻¹² m²). The name honors the French Engineer Henry Darcy who first described the flow of water through sand filters for potable water supply. Permeability values for most materials commonly range typically from a fraction to several thousand millidarcies. The unit of square centimetre (cm²) is also sometimes used (1 cm² = 10⁻⁴ m² <math>\approx</math> 10⁸ d).

Applications

The concept of permeability is of importance in determining the flow characteristics of hydrocarbons in oil and gas reservoirs,<ref>Guerriero V, et al. (2012). "A permeability model for naturally fractured carbonate reservoirs". Marine and Petroleum Geology. 40: 115–134. Bibcode:1990MarPG...7..410M. doi:10.1016/j.marpetgeo.2012.11.002.</ref> and of groundwater in aquifers.<ref>Multiphase fluid flow in porous media From Transport in porous media</ref>

For a rock to be considered as an exploitable hydrocarbon reservoir without stimulation, its permeability must be greater than approximately 100 md (depending on the nature of the hydrocarbon – gas reservoirs with lower permeabilities are still exploitable because of the lower viscosity of gas with respect to oil). Rocks with permeabilities significantly lower than 100 md can form efficient seals (see petroleum geology). Unconsolidated sands may have permeabilities of over 5000 md.

The concept also has many practical applications outside of geology, for example in chemical engineering (e.g., filtration), as well as in Civil Engineering when determining whether the ground conditions of a site are suitable for construction.

Description

Permeability is part of the proportionality constant in Darcy's law which relates discharge (flow rate) and fluid physical properties (e.g. viscosity), to a pressure gradient applied to the porous media:<ref>Controlling Capillary Flow, an application of Darcy's law, at iMechanica</ref>

<math>v = \frac {k}{\eta} \frac{\Delta P}{\Delta x}</math> (for linear flow)

Therefore:

<math>k = v \frac{\eta \Delta x}{\Delta P}</math>

where:

<math>v</math> is the fluid velocity through the porous medium (i.e., the average flow velocity calculated as if the fluid was the only phase present in the porous medium) (m/s)

<math>k</math> is the permeability of a medium (m²)

<math>\eta</math> is the dynamic viscosity of the fluid (Pa·s)

<math>\Delta P</math> is the applied pressure difference (Pa)

<math>\Delta x</math> is the thickness of the bed of the porous medium (m)

In naturally occurring materials, the permeability values range over many orders of magnitude (see table below for an example of this range).

Relation to hydraulic conductivity

The global proportionality constant for the flow of water through a porous medium is called the hydraulic conductivity (K, unit: m/s). Permeability, or intrinsic permeability, (k, unit: m²) is a part of this, and is a specific property characteristic of the solid skeleton and the microstructure of the porous medium itself, independently of the nature and properties of the fluid flowing through the pores of the medium. This allows to take into account the effect of temperature on the viscosity of the fluid flowing though the porous medium and to address other fluids than pure water, e.g., concentrated brines, petroleum, or organic solvents. Given the value of hydraulic conductivity for a studied system, the permeability can be calculated as follows:

<math> k = K \frac {\eta} {\rho g}</math>

where

• <math>k</math> is the permeability, m²
• <math>K</math> is the hydraulic conductivity, m/s
• <math>\eta</math> is the dynamic viscosity of the fluid, Pa·s
• <math>\rho</math> is the density of the fluid, kg/m³
• <math>g</math> is the acceleration due to gravity, m/s².

Anisotropic permeability

Tissue such as brain, liver, muscle, etc can be treated as a heterogeneous porous medium. Describing the flow of biofluids (blood, cerebrospinal fluid, etc.) within such a medium requires a full 3-dimensional anisotropic treatment of the tissue. In this case the scalar hydraulic permeability is replaced with the hydraulic permeability tensor so that Darcy's Law reads<ref>Sowinski, Damian (2021). "Poroelasticity as a Model of Soft Tissue Structure: Hydraulic Permeability Reconstruction for Magnetic Resonance Elastography in Silico". Frontiers in Physics. 8: 637. arXiv:2012.03993. Bibcode:2021FrP.....8..637S. doi:10.3389/fphy.2020.617582. PMC 9635531. PMID 36340954.</ref>

<math>\boldsymbol q = -\frac{1}{\eta}\boldsymbol \kappa \cdot\nabla P </math>

• <math>\boldsymbol q</math> is the Darcy flux, or filtration velocity, which describes the bulk (not microscopic) velocity field of the fluid, <math>[\text{Length}][\text{Time}]^{-1}</math>
• <math>\eta</math> is the dynamic viscosity of the fluid, <math>[\text{Mass}][\text{L}]^{-1}[T]^{-1}</math>
• <math>\boldsymbol \kappa </math> is the hydraulic permeability tensor, <math>[\text{L}]^2</math>
• <math>\nabla </math> is the gradient operator, <math>[\text{L}]^{-1}</math>
• <math>P</math> is the pressure field in the fluid, <math>[\text{M}][\text{L}][\text{T}]^{-1}</math>

Connecting this expression to the isotropic case, <math>\boldsymbol \kappa = k\mathbb 1</math>, where k is the scalar hydraulic permeability, and 1 is the identity tensor.

Determination

Permeability is typically determined in the lab by application of Darcy's law under steady state conditions or, more generally, by application of various solutions to the diffusion equation for unsteady flow conditions.<ref> "CalcTool: Porosity and permeability calculator". www.calctool.org. Retrieved 2008-05-30.</ref>

Permeability needs to be measured, either directly (using Darcy's law), or through estimation using empirically derived formulas. However, for some simple models of porous media, permeability can be calculated (e.g., random close packing of identical spheres).

Permeability model based on conduit flow

Based on the Hagen–Poiseuille equation for viscous flow in a pipe, permeability can be expressed as:

<math>k_{I}=C \cdot d^2</math>

where:

<math>k_{I}</math> is the intrinsic permeability [length²]

<math>C</math> is a dimensionless constant that is related to the configuration of the flow-paths

<math>d</math> is the average, or effective pore diameter [length].

Absolute permeability (aka intrinsic or specific permeability<ref>"Chapter 2: Physical Properties and Principles | Freeze and Cherry Groundwater Book". 2016-09-08. Retrieved 2023-05-02.</ref>)

Absolute permeability denotes the permeability in a porous medium that is 100% saturated with a single-phase fluid. This may also be called the intrinsic permeability or specific permeability. These terms refer to the quality that the permeability value in question is an intensive property of the medium, not a spatial average of a heterogeneous block of material equation 2.28; and that it is a function of the material structure only (and not of the fluid). They explicitly distinguish the value from that of relative permeability.

Permeability to gases

Sometimes permeability to gases can be somewhat different than those for liquids in the same media. One difference is attributable to "slippage" of gas at the interface with the solid<ref>L. J. Klinkenberg, "The Permeability Of Porous Media To Liquids And Gases", Drilling and Production Practice, 41-200, 1941 (abstract).</ref> when the gas mean free path is comparable to the pore size (about 0.01 to 0.1 μm at standard temperature and pressure). See also Knudsen diffusion and constrictivity. For example, measurement of permeability through sandstones and shales yielded values from 9.0×10⁻¹⁹ m² to 2.4×10⁻¹² m² for water and between 1.7×10⁻¹⁷ m² to 2.6×10⁻¹² m² for nitrogen gas.<ref>J. P. Bloomfield and A. T. Williams, "An empirical liquid permeability-gas permeability correlation for use in aquifer properties studies". Quarterly Journal of Engineering Geology & Hydrogeology; November 1995; v. 28; no. Supplement 2; pp. S143–S150. (abstract)</ref> Gas permeability of reservoir rock and source rock is important in petroleum engineering, when considering the optimal extraction of gas from unconventional sources such as shale gas, tight gas, or coalbed methane.

Permeability tensor

To model permeability in anisotropic media, a permeability tensor is needed. Pressure can be applied in three directions, and for each direction, permeability can be measured (via Darcy's law in 3D) in three directions, thus leading to a 3 by 3 tensor. The tensor is realised using a 3 by 3 matrix being both symmetric and positive definite (SPD matrix):

• The tensor is symmetric by the Onsager reciprocal relations
• The tensor is positive definite because the energy being expended (the inner product of fluid flow and negative pressure gradient) is always positive

The permeability tensor is always diagonalizable (being both symmetric and positive definite). The eigenvectors will yield the principal directions of flow where flow is parallel to the pressure gradient, and the eigenvalues represent the principal permeabilities.

Ranges of common intrinsic permeabilities

These values do not depend on the fluid properties; see the table derived from the same source for values of hydraulic conductivity, which are specific to the material through which the fluid is flowing.<ref>Bear, Jacob, 1972. Dynamics of Fluids in Porous Media, Dover. ISBN 0-486-65675-6</ref>

Permeability	Pervious				Semi-pervious				Impervious
Unconsolidated sand and gravel	Well sorted gravel		Well sorted sand or sand and gravel			Very fine sand, silt, loess, loam
Unconsolidated clay and organic					Peat		Layered clay			Unweathered clay
Consolidated rocks	Highly fractured rocks				Oil reservoir rocks			Fresh sandstone		Fresh limestone, dolomite		Fresh granite
k (cm2)	0.001	0.0001	10−5	10−6	10−7	10−8	10−9	10−10	10−11	10−12	10−13	10−14	10−15
k (m2)	10−7	10−8	10−9	10−10	10−11	10−12	10−13	10−14	10−15	10−16	10−17	10−18	10−19
k (millidarcy)	10+8	10+7	10+6	10+5	10,000	1,000	100	10	1	0.1	0.01	0.001	0.0001

See also

• Fault zone hydrogeology
• Hydraulic conductivity
• Hydrogeology
• Permeation
• Petroleum geology
• Relative permeability
• Klinkenberg correction
• Electrical resistivity measurement of concrete

References

• Wang, H. F., 2000. Theory of Linear Poroelasticity with Applications to Geomechanics and Hydrogeology, Princeton University Press. ISBN 0-691-03746-9

External links

• Defining Permeability
• Tailoring porous media to control permeability
• Permeability of Porous Media
• Graphical depiction of different flow rates through materials of differing permeability
• Web-based porosity and permeability calculator given flow characteristics
• Multiphase fluid flow in porous media
• Florida Method of Test For Concrete Resistivity as an Electrical Indicator of its Permeability Archived 2011-06-16 at the Wayback Machine