Darcy's law


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.

Background

Darcy's law was first determined experimentally by Darcy, but has since been derived from the Navier–Stokes equations via homogenization methods. It is analogous to Fourier's law in the field of heat conduction, Ohm's law in the field of electrical networks, and Fick's law in diffusion theory.
One application of Darcy's law is in the analysis of water flow through an aquifer; Darcy's law along with the equation of conservation of mass simplifies to the groundwater flow equation, one of the basic relationships of hydrogeology.
Morris Muskat first refined Darcy's equation for a single phase flow by including viscosity in the single phase equation of Darcy. This change made it suitable for researchers in the petroleum industry. Based on experimental results by his colleagues Wyckoff and Botset, Muskat and Meres also generalized Darcy's law to cover a multiphase flow of water, oil and gas in the porous medium of a petroleum reservoir. The generalized multiphase flow equations by Muskat and others provide the analytical foundation for reservoir engineering that exists to this day.

Description

Darcy's law, as refined by Morris Muskat, in the absence of gravitational forces and in a homogeneously permeable medium, is given by a simple proportionality relationship between the instantaneous flow rate through a porous medium, the permeability of the medium, the dynamic viscosity of the fluid, and the pressure drop over a given distance, in the form
This equation, for single phase flow, is the defining equation for absolute permeability.
With reference to the diagram to the right, the flux, or discharge per unit area, is defined in units, the permeability in units, the cross-sectional area in units, the total pressure drop in units, the dynamic viscosity in units, and is the length of the sample in units. A number of these parameters are used in alternative definitions below. A negative sign is used in the definition of the flux following the standard physics convention that fluids flow from regions of high pressure to regions of low pressure. Note that the elevation head must be taken into account if the inlet and outlet are at different elevations. If the change in pressure is negative, then the flow will be in the positive direction. There have been several proposals for a constitutive equation for absolute permeability, and the most famous one is probably the Kozeny equation.
The integral form of the Darcy law is given by:
where is the total discharge. By considering the relation for static fluid pressure :
one can deduce the representation
where ν is the kinematic viscosity.
The corresponding hydraulic conductivity is therefore:
This quantity, often referred to as the Darcy flux or Darcy velocity, is not the velocity at which the fluid is traveling through the pores. The flow velocity is related to the flux by the porosity and takes the form
Darcy's law is a simple mathematical statement which neatly summarizes several familiar properties that groundwater flowing in aquifers exhibits, including:
A graphical illustration of the use of the steady-state groundwater flow equation is in the construction of flownets, to quantify the amount of groundwater flowing under a dam.
Darcy's law is only valid for slow, viscous flow; however, most groundwater flow cases fall in this category. Typically any flow with a Reynolds number less than one is clearly laminar, and it would be valid to apply Darcy's law. Experimental tests have shown that flow regimes with Reynolds numbers up to 10 may still be Darcian, as in the case of groundwater flow. The Reynolds number for porous media flow is typically expressed as
where is the kinematic viscosity of water, is the specific discharge, is a representative grain diameter for the porous media.

Derivation

For stationary, creeping, incompressible flow, i.e., the Navier–Stokes equation simplifies to the Stokes equation, which by neglecting the bulk term is:
where is the viscosity, is the velocity in the direction, is the gravity component in the direction and is the pressure. Assuming the viscous resisting force is linear with the velocity we may write:
where is the porosity, and is the second order permeability tensor. This gives the velocity in the direction,
which gives Darcy's law for the volumetric flux density in the direction,
In isotropic porous media the off-diagonal elements in the permeability tensor are zero, for and the diagonal elements are identical,, and the common form is obtained
The above equation is a governing equation for single phase fluid flow in a porous medium.

Use in petroleum engineering

Another derivation of Darcy's law is used extensively in petroleum engineering to determine the flow through permeable media — the most simple of which is for a one-dimensional, homogeneous rock formation with a single fluid phase and constant fluid viscosity.
Almost all oil reservoirs have a water zone below the oil leg, and some have also a gas cap above the oil leg. When the reservoir pressure drops due to oil production, water flows into the oil zone from below, and gas flows into the oil zone from above, and we get a simultaneous flow and immiscible mixing of all fluid phases in the oil zone. The operator of the oil field may also inject water in order to improve oil production. The petroleum industry is therefore using a generalized Darcy equation for multiphase flow that was developed by Muskat et alios. Because Darcy's name is so widespread and strongly associated with flow in porous media, the multiphase equation is denoted Darcy's law for multiphase flow or generalized Darcy equation or simply Darcy's equation or simply flow equation if the context says that the text is discussing the multiphase equation of Muskat et alios. Multiphase flow in oil and gas reservoirs is a comprehensive topic, and one of many articles about this topic is Darcy's law for multiphase flow.

Additional forms

Quadratic law

For flows in porous media with Reynolds numbers greater than about 1 to 10, inertial effects can also become significant. Sometimes an inertial term is added to the Darcy's equation, known as Forchheimer term. This term is able to account for the non-linear behavior of the pressure difference vs flow data.
where the additional term is known as inertial permeability.
The flow in the middle of a sandstone reservoir is so slow that Forchheimer's equation is usually not needed, but the gas flow into a gas production well may be high enough to justify use of Forchheimer's equation. In this case the inflow performance calculations for the well, not the grid cell of the 3D model, is based on the Forchheimer equation. The effect of this is that an additional rate-dependent skin appears in the inflow performance formula.
Some carbonate reservoirs have many fractures, and Darcy's equation for multiphase flow is generalized in order to govern both flow in fractures and flow in the matrix. The irregular surface of the fracture walls and high flow rate in the fractures, may justify use of Forchheimer's equation.

Correction for gases in fine media (Knudsen diffusion or Klinkenberg effect)

For gas flow in small characteristic dimensions, the particle-wall interactions become more frequent, giving rise to additional wall friction. For a flow in this region, where both viscous and Knudsen friction are present, a new formulation needs to be used. Knudsen presented a semi-empirical model for flow in transition regime based on his experiments on small capillaries. For a porous medium, the Knudsen equation can be given as
where is the molar flux, is the gas constant, is the temperature, is the effective Knudsen diffusivity of the porous media. The model can also be derived from the first-principle-based binary friction model. The differential equation of transition flow in porous media based on BFM is given as
This equation is valid for capillaries as well as porous media. The terminology of the Knudsen effect and Knudsen diffusivity is more common in mechanical and chemical engineering. In geological and petrochemical engineering, this effect is known as the Klinkenberg effect. Using the definition of molar flux, the above equation can be rewritten as
This equation can be rearranged into the following equation
Comparing this equation with conventional Darcy's law, a new formulation can be given as
where
This is equivalent to the effective permeability formulation proposed by Klinkenberg:
where is known as the Klinkenberg parameter, which depends on the gas and the porous medium structure. This is quite evident if we compare the above formulations. The Klinkenberg parameter is dependent on permeability, Knudsen diffusivity and viscosity.

Darcy's law for short time scales

For very short time scales, a time derivative of flux may be added to Darcy's law, which results in valid solutions at very small times,
where is a very small time constant which causes this equation to reduce to the normal form of Darcy's law at "normal" times. The main reason for doing this is that the regular groundwater flow equation leads to singularities at constant head boundaries at very small times. This form is more mathematically rigorous, but leads to a hyperbolic groundwater flow equation, which is more difficult to solve and is only useful at very small times, typically out of the realm of practical use.

Brinkman form of Darcy's law

Another extension to the traditional form of Darcy's law is the Brinkman term, which is used to account for transitional flow between boundaries,
where is an effective viscosity term. This correction term accounts for flow through medium where the grains of the media are porous themselves, but is difficult to use, and is typically neglected. For example, if a porous extracellular matrix degrades to form large pores throughout the matrix, the viscous term applies in the large pores, while Darcy's law applies in the remaining intact region. This scenario was considered in a theoretical and modelling study. In the proposed model, the Brinkman equation is connected to a set of reaction-diffusion-convection equations.

Validity of Darcy's law

Darcy's law is valid for laminar flow through sediments. In fine-grained sediments, the dimensions of are small and thus flow is laminar. Coarse-grained sediments also behave similarly but in very coarse-grained sediments the flow may be turbulent. Hence Darcy's law is not always valid in such sediments.
For flow through commercial circular pipes, the flow is laminar when Reynolds number is less than 2000 and turbulent when it is more than 4000, but in some sediments it has been found that flow is laminar when the value of Reynolds number is less than 1.