Fluid flow through porous media

{{Short description|Manner in which fluids behave when flowing through a porous medium}}

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.

Classical flow mechanics in porous media assumes that the medium is homogenous, isotropic, and has an intergranular pore structure. It also assumes that the fluid is a Newtonian fluid, that the reservoir is isothermal, that the well is vertical, etc. Traditional flow issues in porous media often involve single-phase steady state flow, multi-well interference, oil-water two-phase flow, natural gas flow, elastic energy driven flow, oil-gas two-phase flow, and gas-water two-phase flow.{{Citation |title=Classical flow theory in porous media |date=2020-02-28 |url=https://www.worldscientific.com/doi/abs/10.1142/9789811219535_0003 |work=Fluid Flow in Porous Media |pages=69–197 |publisher=WORLD SCIENTIFIC |doi=10.1142/9789811219535_0003 |isbn=978-981-12-1952-8 |s2cid=242538024 |access-date=2022-07-01}}

The physicochemical flow process will involve various physical property changes and chemical reactions in contrast to the basic Newtonian fluid in the classical flow theory of porous system. Viscosity, surface tension, phase state, concentration, temperature, and other physical characteristics are examples of these properties. Non-Newtonian fluid flow, mass transfer through diffusion, and multiphase and multicomponent fluid flow are the primary flow issues.{{Citation |title=Physicochemical flow in porous media |date=2020-02-28 |url=https://www.worldscientific.com/doi/abs/10.1142/9789811219535_0005 |work=Fluid Flow in Porous Media |pages=249–285 |publisher=WORLD SCIENTIFIC |doi=10.1142/9789811219535_0005 |isbn=978-981-12-1952-8 |s2cid=241660780 |access-date=2022-07-01}}

Governing laws

The movement of a fluid through porous media is described by the combination of Darcy's law with the principle of conservation of mass in order to express the capillary force or fluid velocity as a function of various other parameters including the effective pore radius, liquid viscosity or permeability.[https://hal.archives-ouvertes.fr/hal-02310722/document#page=6 Tailoring porous media for controllable capillary flow] Journal of colloid and interface science (2019), 539, 379-387 However, the use of Darcy's law alone does not produce accurate results for heterogeneous media like shale, and tight sandstones, where there is a huge proportion of nanopores. This necessitates the use of a flow model that considers the weighted proportion of various flow regimes like Darcy flow, transition flow, slip flow, and free molecular flow.{{cite journal | last1=Saurabh | first1=Suman | title=Flow regime evolution and stress-dependent permeability in nanoporous rocks | date=2022-02-15 | url=https://www.sciencedirect.com/science/article/abs/pii/S0016236121022857| journal=Fuel| volume=310 | page=122413 |publisher=Elsevier|doi=10.1016/j.fuel.2021.122413 | bibcode=2022Fuel..31022413S | s2cid=243474407 }}

=Darcy's law=

align="right" width="20%" class="wikitable"
SymbolDescription
QVolumetric flow rate [m3/s]
kPermeability of porous medium [m2]. The permeability is a function of material type, and also varies with stress, temperature, etc.
\muFluid viscosity [Pa.s]
ACross-sectional area of Porous medium [m2]
(p_b-p_a)Pressure drop across medium [Pa]
LLength of sample [m]

The basic law governing the flow of fluids through porous media is Darcy's Law, which was formulated by the French civil engineer Henry Darcy in 1856 on the basis of his experiments on vertical water filtration through sand beds.{{cite journal | doi=10.1007/BF01036523 | volume=1 | title=Flow in porous media I: A theoretical derivation of Darcy's law | journal=Transport in Porous Media | pages=3–25|year = 1986|last1 = Whitaker|first1 = Stephen| issue=1 | bibcode=1986TPMed...1....3W | s2cid=121904058 }}

According to Darcy's law, the fluid's viscosity, effective fluid permeability, and fluid pressure gradient determine the flow rate at any given location in the reservoir.{{cite journal |last1=Johns |first1=Russell T. |title=Oil Recovery |journal=Encyclopedia of Energy |date=1 January 2004 |pages=701–713 |doi=10.1016/b0-12-176480-x/00255-2 |url=https://doi.org/10.1016/B0-12-176480-X/00255-2 |publisher=Elsevier |isbn=9780121764807 |language=en}}

{{Equation box 1

|indent=::

|equation=:Q=\frac{-k A}{\mu} \frac{(p_b - p_a)}{L}

|cellpadding

|border

|border colour = #0073CF

|background colour=#F5FFFA

}}

For transient processes in which the flux varies from point to-point, the following differential form of Darcy’s law is used.

{{Equation box 1

|indent=::

|equation=:Q= \frac {-k A}{\mu} \left( \frac{d p}{d x} \right)

|cellpadding

|border

|border colour = #0073CF

|background colour=#F5FFFA

}}

Darcy's law is valid for situation where the porous material is already saturated with the fluid. For the calculation of capillary imbibition speed of a liquid to an initially dry medium, Washburn's or Bosanquet's equations are used.

=Mass conservation=

Mass conservation of fluid across the porous medium involves the basic principle that mass flux in minus mass flux out equals the increase in amount stored by a medium.{{cite book|last1=Bear|first1=Jacob|title=Dynamics of fluids in Porous Media|url=https://books.google.com/books?id=fBMeVSZ_3u8C&dq=Mass+conservation+in+porous+media&pg=PP1|isbn=9780486131801|date=2013-02-26|publisher=Courier Corporation }} This means that total mass of the fluid is always conserved. In mathematical form, considering a time period from t to \Delta t , length of porous medium from x to \Delta x and m being the mass stored by the medium, we have

: [A\rho (x)q(x) - A\rho (x + \Delta x)q(x + \Delta x )]\Delta t = m(t + \Delta t) - m(t) .

Furthermore, we have that m = \rho V_p, where V_p is the pore volume of the medium between x and x+\Delta x and \rho is the density. So m = \rho V_p = \rho \phi V = \rho \phi A\Delta x, where \phi is the porosity. Dividing both sides by A\Delta x, while \Delta x \rightarrow 0, we have for 1 dimensional linear flow in a porous medium the relation

: \frac{-d(\rho q)}{d x} = \frac{d(\rho \phi)}{d t} ~~~~ (i)

In three dimensions, the equation can be written as

: \frac{d(\rho q)}{d x} + \frac{d(\rho q)}{d y} + \frac{d(\rho q)}{d z} = \frac{-d(\rho \phi)}{d t}

The mathematical operation on the left-hand side of this equation is known as the divergence of \rho q, and represents the rate at which fluid diverges from a given region, per unit volume.

Diffusion Equation

align="right" class="wikitable"
Material TypeCompressibility (m2N−1 or Pa−1)https://eng.ucmerced.edu/people/jfisher/.../EnveEss110_20081110.pdf{{Dead link|date=December 2019 |bot=InternetArchiveBot |fix-attempted=yes }}
Clay10−6 - 10−8
Sand10−7 - 10−9
Gravel10−8 - 10−10
Jointed rock10−8 - 10−10
Sound Rock10−9 - 10−11
Water (beta)4.4 x 10−10

Using product rule(and chain rule) on right hand side of the above mass conservation equation (i),

: \frac{d(\rho \phi)}{d t} = \rho \frac{d \phi}{d t} + \phi \frac{d \rho}{d t} = \rho \frac{d \phi}{d P} \frac{d P}{d t} + \phi \frac{d \rho}{d P} \frac{d P}{d t} = \rho \phi \left [\frac{1}{\phi} \frac{d \phi}{d P} + \frac{1}{\rho} \frac{d \rho}{d P} \right] \quad \frac{d P}{d t} = \rho \phi [c_\phi + c_f]\frac{d P}{d t} ~~~~ (ii)

Here, c_f = compressibility of the fluid and c_\phi = compressibility of porous medium.{{cite journal|last1=Ohirhian|first1=Peter|title=Steady State Compressible Fluid Flow in Porous Media}} Now considering the left hand side of the mass conservation equation, which is given by Darcy's Law as

: \frac{-d(\rho q)}{d x} = \frac{-d}{dx} \left [\frac{- \rho k}{\mu} \frac{d P}{d x} \right] \quad = \frac{k}{\mu} \left [\rho \frac{d^2 P}{d x^2} + \frac{d \rho}{d P} \frac{d P}{d x} \frac{d P}{d x} \right] \quad = \frac{\rho k}{\mu} \left[\frac{d^2 P}{d x^2} + \left (\frac{1}{\rho} \frac{d \rho}{d P}\right)\left(\frac{d P}{d x} \right)^2 \right ] \quad = \frac{\rho k}{\mu} \left [\frac{d^2 P}{d x^2} + c_f \left(\frac{d P}{d x}\right)^2 \right] \quad ~~~~ (iii)

Equating the results obtained in (ii) & (iii), we get

: \frac{d^2 P}{d x^2} + c_f\left(\frac{d P}{d x}\right)^2 = \frac{\phi \mu (c_f + c_\phi)}{k} \frac{d P}{d t}

The second term on the left side is usually negligible, and we obtain the diffusion equation in 1 dimension as

: \frac{d P}{d t} = \frac{k}{\phi \mu c_t} \frac{d^2 P}{d x^2}

where c_t = c_f + c_\phi .{{cite web|last1=Zimmerman|first1=Dr. R.W.|title=Flow in porous media|url=http://imperial.ac.uk}}

References

{{Reflist}}

Further reading

  • {{cite book|last=Acheson|first=D. J.|title=Elementary Fluid Dynamics|publisher=Clarendon Press|year=1990|isbn=0-19-859679-0}}
  • {{cite book|last=Lamb|first=Horace|author-link=Horace Lamb|title=Hydrodynamics|edition=6th|publisher=Cambridge University Press|year=1994|isbn=0-521-45868-4}} Originally published in 1879, the 6th extended edition appeared first in 1932.
  • {{cite book|last=Milne-Thompson|first=L. M.|title=Theoretical Hydrodynamics|edition=5th|publisher=Macmillan|year=1968}} Originally published in 1938.
  • {{cite book|last=Shinbrot|first=M.|title=Lectures on Fluid Mechanics|publisher=Gordon and Breach|year=1973|isbn=0-677-01710-3}}
  • {{cite book|last=Chanson|first=H.|author-link=Hubert Chanson|title=Applied Hydrodynamics: An Introduction to Ideal and Real Fluid Flows|publisher=CRC Press, Taylor & Francis Group, Leiden, The Netherlands, 478 pages|year=2009|isbn=978-0-415-49271-3}}
  • {{cite book|last=Clancy|first=L. J.|title=Aerodynamics|publisher=Pitman Publishing Limited|location=London|year=1975|isbn=0-273-01120-0}}
  • {{citation | last1=Nazarenko | first1=Sergey | year=2014 | title=Fluid Dynamics via Examples and Solutions | publisher=CRC Press (Taylor & Francis group) | isbn=978-1-43-988882-7 }}