Permeability (Earth sciences)

Last updated

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.

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.

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.

<span class="mw-page-title-main">Porous medium</span> Material containing fluid-filled voids

In materials science, a porous medium or a porous material is a material containing pores (voids). The skeletal portion of the material is often called the "matrix" or "frame". The pores are typically filled with a fluid. The skeletal material is usually a solid, but structures like foams are often also usefully analyzed using concept of porous media.

In physics and engineering, a constitutive equation or constitutive relation is a relation between two physical quantities that is specific to a material or substance, and approximates the response of that material to external stimuli, usually as applied fields or forces. They are combined with other equations governing physical laws to solve physical problems; for example in fluid mechanics the flow of a fluid in a pipe, in solid state physics the response of a crystal to an electric field, or in structural analysis, the connection between applied stresses or loads to strains or deformations.

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 fractures 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.

In physics and engineering, permeation is the penetration of a permeate through a solid. It is directly related to the concentration gradient of the permeate, a material's intrinsic permeability, and the materials' mass diffusivity. Permeation is modeled by equations such as Fick's laws of diffusion, and can be measured using tools such as a minipermeameter.

In multiphase flow in porous media, the relative permeability of a phase is a dimensionless measure of the effective permeability of that phase. It is the ratio of the effective permeability of that phase to the absolute permeability. It can be viewed as an adaptation of Darcy's law to multiphase flow.

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.

The Kozeny–Carman equation is a relation used in the field of fluid dynamics to calculate the pressure drop of a fluid flowing through a packed bed of solids. It is named after Josef Kozeny and Philip C. Carman. The equation is only valid for creeping flow, i.e. in the slowest limit of laminar flow. The equation was derived by Kozeny (1927) and Carman from a starting point of (a) modelling fluid flow in a packed bed as laminar fluid flow in a collection of curving passages/tubes crossing the packed bed and (b) Poiseuille's law describing laminar fluid flow in straight, circular section pipes.

In fluid dynamics, Luke's variational principle is a Lagrangian variational description of the motion of surface waves on a fluid with a free surface, under the action of gravity. This principle is named after J.C. Luke, who published it in 1967. This variational principle is for incompressible and inviscid potential flows, and is used to derive approximate wave models like the mild-slope equation, or using the averaged Lagrangian approach for wave propagation in inhomogeneous media.

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">Diffusion</span> Transport of dissolved species from the highest to the lowest concentration region

Diffusion is the net movement of anything generally from a region of higher concentration to a region of lower concentration. Diffusion is driven by a gradient in Gibbs free energy or chemical potential. It is possible to diffuse "uphill" from a region of lower concentration to a region of higher concentration, as in spinodal decomposition. Diffusion is a stochastic process due to the inherent randomness of the diffusing entity and can be used to model many real-life stochastic scenarios. Therefore, diffusion and the corresponding mathematical models are used in several fields beyond physics, such as statistics, probability theory, information theory, neural networks, finance, and marketing.

<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.

<span class="mw-page-title-main">Morris Muskat</span> American petroleum engineer

Morris Muskat was an American petroleum engineer. Muskat refined Darcy's equation for single phase flow, and this change made it suitable for the petroleum industry. Based on experimental results worked out by his colleagues, Muskat and Milan W. Meres also generalized Darcy's law to cover multiphase flow of water, oil and gas in the porous medium of a petroleum reservoir. The generalized flow equation provides the analytical foundation for reservoir engineering that exists to this day.

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