Permeability (materials science)

Last updated

Permeability in fluid mechanics, materials science and 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.

Contents

Symbol used to represent in situ permeability tests in geotechnical drawings In situ permeameter test.svg
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. [1] 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. [2] 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. [3] Permeability is also affected by the pressure inside a material.

Units

The SI unit for permeability is the square metre (m2). A practical unit for permeability is the darcy (d), or more commonly the millidarcy (md) (1 d 10−12 m2). 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 (cm2) is also sometimes used (1 cm2 = 10−4 m2 108 d).

Applications

The concept of permeability is of importance in determining the flow characteristics of hydrocarbons in oil and gas reservoirs, [4] and of groundwater in aquifers. [5]

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: [6]

(for linear flow)

Therefore:

where:

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)
is the permeability of a medium (m2)
is the dynamic viscosity of the fluid (Pa·s)
is the applied pressure difference (Pa)
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: m2) 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:

where

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 [7]

Connecting this expression to the isotropic case, , 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. [8]

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:

where:

is the intrinsic permeability [length2]
is a dimensionless constant that is related to the configuration of the flow-paths
is the average, or effective pore diameter [length].

Absolute permeability (aka intrinsic or specific permeability [9] )

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 [ clarification needed ][ further explanation needed ]; 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 [10] 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−19 m2 to 2.4×10−12 m2 for water and between 1.7×10−17 m2 to 2.6×10−12 m2 for nitrogen gas. [11] 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 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. [12]

PermeabilityPerviousSemi-perviousImpervious
Unconsolidated sand and gravelWell sorted gravel Well sorted sand or sand and gravelVery fine sand, silt, loess, loam
Unconsolidated clay and organic Peat Layered clay Unweathered clay
Consolidated rocksHighly fractured rocks Oil reservoir rocksFresh sandstone Fresh limestone, dolomite Fresh granite
k (cm2)0.0010.000110−510−610−710−810−910−1010−1110−1210−1310−1410−15
k (m2)10−710−810−910−1010−1110−1210−1310−1410−1510−1610−1710−1810−19
k (millidarcy)10+810+710+610+510,0001,0001001010.10.010.0010.0001

See also

Footnotes

  1. "Reading: Porosity and Permeability | Geology". courses.lumenlearning.com. Retrieved 2022-01-14.
  2. 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.
  3. 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.
  4. 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.
  5. Multiphase fluid flow in porous media From Transport in porous media
  6. Controlling Capillary Flow, an application of Darcy's law, at iMechanica
  7. 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.
  8. "CalcTool: Porosity and permeability calculator". www.calctool.org. Retrieved 2008-05-30.
  9. "Chapter 2: Physical Properties and Principles | Freeze and Cherry Groundwater Book". 2016-09-08. Retrieved 2023-05-02.
  10. L. J. Klinkenberg, "The Permeability Of Porous Media To Liquids And Gases", Drilling and Production Practice, 41-200, 1941 (abstract).
  11. 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)
  12. Bear, Jacob, 1972. Dynamics of Fluids in Porous Media, Dover. ISBN   0-486-65675-6

Related Research Articles

In fluid mechanics, the Rayleigh number (Ra, after Lord Rayleigh) for a fluid is a dimensionless number associated with buoyancy-driven flow, also known as free (or natural) convection. It characterises the fluid's flow regime: a value in a certain lower range denotes laminar flow; a value in a higher range, turbulent flow. Below a certain critical value, there is no fluid motion and heat transfer is by conduction rather than convection. For most engineering purposes, the Rayleigh number is large, somewhere around 106 to 108.

A Newtonian fluid is a fluid in which the viscous stresses arising from its flow are at every point linearly correlated to the local strain rate — the rate of change of its deformation over time. Stresses are proportional to the rate of change of the fluid's velocity vector.

Darcy's law is an equation that describes the flow of a fluid through a porous medium. The law was formulated by Henry Darcy based on results of experiments on the flow of water through beds of sand, forming the basis of hydrogeology, a branch of earth sciences. It is analogous to Ohm's law in electrostatics, linearly relating the volume flow rate of the fluid to the hydraulic head difference via the hydraulic conductivity. In fact, the Darcy's law is a special case of the Stokes equation for the momentum flux, in turn deriving from the momentum Navier-Stokes equation.

In physics, Washburn's equation describes capillary flow in a bundle of parallel cylindrical tubes; it is extended with some issues also to imbibition into porous materials. The equation is named after Edward Wight Washburn; also known as Lucas–Washburn equation, considering that Richard Lucas wrote a similar paper three years earlier, or the Bell-Cameron-Lucas-Washburn equation, considering J.M. Bell and F.K. Cameron's discovery of the form of the equation in 1906.

The darcy and millidarcy are units of permeability, named after Henry Darcy. They are not SI units, but they are widely used in petroleum engineering and geology. The unit has also been used in biophysics and biomechanics, where the flow of fluids such as blood through capillary beds and cerebrospinal fluid through the brain interstitial space is being examined. A darcy has dimensional units of length2.

In science and engineering, hydraulic conductivity, is a property of porous materials, soils and rocks, that describes the ease with which a fluid can move through the pore space, or fracture network. It depends on the intrinsic permeability of the material, the degree of saturation, and on the density and viscosity of the fluid. Saturated hydraulic conductivity, Ksat, describes water movement through saturated media. By definition, hydraulic conductivity is the ratio of volume flux to hydraulic gradient yielding a quantitative measure of a saturated soil's ability to transmit water when subjected to a hydraulic gradient.

Used in hydrogeology, the groundwater flow equation is the mathematical relationship which is used to describe the flow of groundwater through an aquifer. The transient flow of groundwater is described by a form of the diffusion equation, similar to that used in heat transfer to describe the flow of heat in a solid. The steady-state flow of groundwater is described by a form of the Laplace equation, which is a form of potential flow and has analogs in numerous fields.

The Richards equation represents the movement of water in unsaturated soils, and is attributed to Lorenzo A. Richards who published the equation in 1931. It is a quasilinear partial differential equation; its analytical solution is often limited to specific initial and boundary conditions. Proof of the existence and uniqueness of solution was given only in 1983 by Alt and Luckhaus. The equation is based on Darcy-Buckingham law representing flow in porous media under variably saturated conditions, which is stated as

The void ratio of a mixture of solids and fluids, or of a porous composite material such as concrete, is the ratio of the volume of the voids filled by the fluids to the volume of all the solids. It is a dimensionless quantity in materials science and in soil science, and is closely related to the porosity, the ratio of the volume of voids to the total volume, as follows:

The air permeability specific surface of a powder material is a single-parameter measurement of the fineness of the powder. The specific surface is derived from the resistance to flow of air through a porous bed of the powder. The SI units are m2·kg−1 or m2·m−3.

<span class="mw-page-title-main">HydroGeoSphere</span>

HydroGeoSphere (HGS) is a 3D control-volume finite element groundwater model, and is based on a rigorous conceptualization of the hydrologic system consisting of surface and subsurface flow regimes. The model is designed to take into account all key components of the hydrologic cycle. For each time step, the model solves surface and subsurface flow, solute and energy transport equations simultaneously, and provides a complete water and solute balance.

Volume viscosity is a material property relevant for characterizing fluid flow. Common symbols are or . It has dimensions, and the corresponding SI unit is the pascal-second (Pa·s).

<span class="mw-page-title-main">Viscosity</span> Resistance of a fluid to shear deformation

The viscosity of a fluid is a measure of its resistance to deformation at a given rate. For liquids, it corresponds to the informal concept of "thickness": for example, syrup has a higher viscosity than water. Viscosity is defined scientifically as a force multiplied by a time divided by an area. Thus its SI units are newton-seconds per square meter, or pascal-seconds.

Porosity or void fraction is a measure of the void spaces in a material, and is a fraction of the volume of voids over the total volume, between 0 and 1, or as a percentage between 0% and 100%. Strictly speaking, some tests measure the "accessible void", the total amount of void space accessible from the surface.

<span class="mw-page-title-main">Radiation stress</span> Term in physical oceanography

In fluid dynamics, the radiation stress is the depth-integrated – and thereafter phase-averaged – excess momentum flux caused by the presence of the surface gravity waves, which is exerted on the mean flow. The radiation stresses behave as a second-order tensor.

<span class="mw-page-title-main">Saffman–Taylor instability</span>

The Saffman–Taylor instability, also known as viscous fingering, is the formation of patterns in a morphologically unstable interface between two fluids in a porous medium, described mathematically by Philip Saffman and G. I. Taylor in a paper of 1958. This situation is most often encountered during drainage processes through media such as soils. It occurs when a less viscous fluid is injected, displacing a more viscous fluid; in the inverse situation, with the more viscous displacing the other, the interface is stable and no instability is seen. Essentially the same effect occurs driven by gravity if the interface is horizontal and separates two fluids of different densities, the heavier one being above the other: this is known as the Rayleigh-Taylor instability. In the rectangular configuration the system evolves until a single finger forms, whilst in the radial configuration the pattern grows forming fingers by successive tip-splitting.

The viscous stress tensor is a tensor used in continuum mechanics to model the part of the stress at a point within some material that can be attributed to the strain rate, the rate at which it is deforming around that point.

In fluid mechanics, fluid flow through porous media is the manner in which fluids behave when flowing through a porous medium, for example sponge or wood, or when filtering water using sand or another porous material. As commonly observed, some fluid flows through the media while some mass of the fluid is stored in the pores present in the media.

The shear viscosity of a fluid is a material property that describes the friction between internal neighboring fluid surfaces flowing with different fluid velocities. This friction is the effect of (linear) momentum exchange caused by molecules with sufficient energy to move between these fluid sheets due to fluctuations in their motion. The viscosity is not a material constant, but a material property that depends on temperature, pressure, fluid mixture composition, local velocity variations. This functional relationship is described by a mathematical viscosity model called a constitutive equation which is usually far more complex than the defining equation of shear viscosity. One such complicating feature is the relation between the viscosity model for a pure fluid and the model for a fluid mixture which is called mixing rules. When scientists and engineers use new arguments or theories to develop a new viscosity model, instead of improving the reigning model, it may lead to the first model in a new class of models. This article will display one or two representative models for different classes of viscosity models, and these classes are:

The porous medium equation, also called the nonlinear heat equation, is a nonlinear partial differential equation taking the form:

References