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.
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.
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.[1] 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.[2]
| Symbol | Description | |
|---|---|---|
Q | Volumetric flow rate [m<sup>3</sup>/s] | |
k | Permeability of porous medium [m<sup>2</sup>]. The permeability is a function of material type, and also varies with stress, temperature, etc. | |
\mu | Fluid viscosity [Pa.s] | |
A | Cross-sectional area of Porous medium [m<sup>2</sup>] | |
(pb-pa) | Pressure drop across medium [Pa] | |
L | Length of sample [m] |
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.[4]
For transient processes in which the flux varies from point to-point, the following differential form of Darcy’s law is used.
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 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.[5] This means that total mass of the fluid is always conserved. In mathematical form, considering a time period from
t
\Deltat
x
\Deltax
m
[A\rho(x)q(x)-A\rho(x+\Deltax)q(x+\Deltax)]\Deltat=m(t+\Deltat)-m(t).
Furthermore, we have that
m=\rhoVp
Vp
x
x+\Deltax
\rho
m=\rhoVp=\rho\phiV=\rho\phiA\Deltax,
\phi
A\Deltax
\Deltax → 0
| -d(\rhoq) | |
| dx |
=
| d(\rho\phi) | |
| dt |
~~~~(i)
In three dimensions, the equation can be written as
| d(\rhoq) | |
| dx |
+
| d(\rhoq) | |
| dy |
+
| d(\rhoq) | |
| dz |
=
| -d(\rho\phi) | |
| dt |
The mathematical operation on the left-hand side of this equation is known as the divergence of
\rhoq
| Material Type | Compressibility (m2N−1 or Pa−1)[6] | |
|---|---|---|
| Clay | 10−6 - 10−8 | |
| Sand | 10−7 - 10−9 | |
| Gravel | 10−8 - 10−10 | |
| Jointed rock | 10−8 - 10−10 | |
| Sound Rock | 10−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),
| d(\rho\phi) | |
| dt |
=\rho
| d\phi | |
| dt |
+\phi
| d\rho | |
| dt |
=\rho
| d\phi | |
| dP |
| dP | |
| dt |
+\phi
| d\rho | |
| dP |
| dP | |
| dt |
=\rho\phi\left[
| 1 | |
| \phi |
| d\phi | |
| dP |
+
| 1 | |
| \rho |
| d\rho | |
| dP |
\right]
| dP | |
| dt |
=\rho\phi[c\phi+
| c | ||||
|
~~~~(ii)
Here,
cf
c\phi
| -d(\rhoq) | |
| dx |
=
| -d | |
| dx |
\left[
| -\rhok | |
| \mu |
| dP | |
| dx |
\right] =
| k | |
| \mu |
\left[\rho
| d2P | |
| dx2 |
+
| d\rho | |
| dP |
| dP | |
| dx |
| dP | |
| dx |
\right] =
| \rhok | \left[ | |
| \mu |
| d2P | |
| dx2 |
+\left(
| 1 | |
| \rho |
| d\rho | \right)\left( | |
| dP |
| dP | |
| dx |
\right)2\right] =
| \rhok | |
| \mu |
\left[
| d2P | |
| dx2 |
+cf\left(
| dP | |
| dx |
\right)2\right] ~~~~(iii)
Equating the results obtained in
(ii)
(iii)
| d2P | |
| dx2 |
+
| c | ||||
|
\right)2=
| \phi\mu(cf+c\phi) | |
| k |
| dP | |
| dt |
The second term on the left side is usually negligible, and we obtain the diffusion equation in 1 dimension as
| dP | |
| dt |
=
| k | |
| \phi\muct |
| d2P | |
| dx2 |
where
ct=cf+c\phi