Sie sind hier: ICP » R. Hilfer » Publikationen

VI.C Macroscopic Description

VI.C.1 Macroscopic Equations of Motion

The accepted large scale equations of motion for two phase flow involve a generalization of Darcy’s law to relative permeabilities including offdiagonal viscous coupling terms [349, 270, 322, 321, 350, 351, 352]. The importance of viscous coupling terms has been recognized relatively late [353, 354, 355, 356, 357]. The equations which are generally believed to describe multiphase flow on the reservoir scale as well as on the laboratory scale may be written as [322, 350]

\begin{array}[]{rcl}\overline{\phi}\frac{\displaystyle\partial\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle\partial\overline{t}}&=&\overline{\mbox{\boldmath$\nabla$}}\cdot\overline{\bf v}_{{\scriptscriptstyle{\mathbb{W}}}}\\
\overline{\phi}\frac{\displaystyle\partial\overline{S}_{{\scriptscriptstyle{\mathbb{O}}}}}{\displaystyle\partial\overline{t}}&=&\overline{\mbox{\boldmath$\nabla$}}\cdot\overline{\bf v}_{{\scriptscriptstyle{\mathbb{O}}}}\end{array} (6.38)
\begin{array}[]{rcl}\overline{\bf v}_{{\scriptscriptstyle{\mathbb{W}}}}&=&-\left[{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{W}}}}}\:\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}}\:(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{W}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})+{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{O}}}}}\:\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\:(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{O}}}}-\rho _{{\scriptscriptstyle{\mathbb{O}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})\right]\\
\overline{\bf v}_{{\scriptscriptstyle{\mathbb{O}}}}&=&-\left[{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}}\:\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}}\:(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{W}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})+{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{O}}}}}\:\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\:(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{O}}}}-\rho _{{\scriptscriptstyle{\mathbb{O}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})\right]\end{array} (6.39)
\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}+\overline{S}_{{\scriptscriptstyle{\mathbb{O}}}}=1 (6.40)
\overline{P}_{{\scriptscriptstyle{\mathbb{O}}}}-\overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}=\overline{P}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}) (6.41)

where \overline{A} denotes the macroscopic volume averaged equivalent of the pore scale quantity A. In the equations above {\bf K} stands for the absolute (single phase flow) permeability tensor, {\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{W}}}}} is the relative permeability tensor for water, {\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{O}}}}} the oil relative permeability tensor, and {\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{O}}}}},{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}} denote the possibly anisotropic coupling terms. The relative permeabilities are matrix valued functions of saturation. The saturations are denoted as \overline{S}_{{\scriptscriptstyle{\mathbb{W}}}},\overline{S}_{{\scriptscriptstyle{\mathbb{O}}}} and they depend on the macroscopic space and time variables (\overline{\bf x},\overline{t}). The capillary pressure curve \overline{P}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}) and the relative permeability tensors {\bf K}^{r}_{{ij}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}),i,j=\mathbb{W},\mathbb{O} must be known either from solving the pore scale equations of motion, or from experiment. {\bf K}^{r}_{{ij}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}) and \overline{P}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}) are conventionally assumed to be independent of \overline{\bf v} and \overline{P} and this convention is followed here, although it is conceivable that this is not generally correct [354].

Eliminating \overline{\bf v} and choosing \overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}(\overline{\bf x},\overline{t}) and \overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}(\overline{\bf x},\overline{t}) as the principal unknowns one arrives at the large scale two-phase flow equations

\displaystyle\overline{\phi} \displaystyle\frac{\displaystyle\partial\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle\partial\overline{t}}=\overline{\mbox{\boldmath$\nabla$}}\cdot\left\{{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{W}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}}(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{W}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})\;\right. (6.42)
\displaystyle+ \displaystyle\left.{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{O}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\:\left[(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{W}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})+\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})+(\rho _{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{O}}}})g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z}\right]\right\}
\displaystyle\overline{\phi} \displaystyle\frac{\displaystyle\partial(1-\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})}{\displaystyle\partial\overline{t}}=\overline{\mbox{\boldmath$\nabla$}}\cdot\left\{{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}}\:(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{W}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})\;\right. (6.43)
\displaystyle+ \displaystyle\left.{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{O}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\frac{\displaystyle{\bf K}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\left[(\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{W}}}}g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z})+\overline{\mbox{\boldmath$\nabla$}}\overline{P}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})+(\rho _{{\scriptscriptstyle{\mathbb{W}}}}-\rho _{{\scriptscriptstyle{\mathbb{O}}}})g\:\overline{\mbox{\boldmath$\nabla$}}\overline{z}\right]\right\}

for these two unknowns. Equations (6.42) and (6.43) are coupled nonlinear partial differential equations for the large scale pressure and saturation field of the water phase.

These equations must be complemented with large scale boundary conditions. For core experiments these are typically given by a surface source on one side of the core, a surface sink on the opposite face, and impermeable walls on the other faces. For a reservoir the boundary conditions depend upon the drive configuration and the geological modeling of the reservoir environment, so that Dirichlet as well as von Neumann problems arise in practice [339, 351, 352].

VI.C.2 Macroscopic Dimensional Analysis

The large scale equations of motion can be cast in dimensionless form using the definitions

\overline{\bf x}=\overline{l}\;\widehat{\overline{\bf x}} (6.44)
\overline{\mbox{\boldmath$\nabla$}}=\frac{\widehat{\overline{\mbox{\boldmath$\nabla$}}}}{\overline{l}} (6.45)
\overline{\bf v}=\overline{u}\;\widehat{\overline{\bf v}} (6.46)
\overline{t}=\frac{\overline{l}\;\widehat{\overline{t}}}{\overline{u}} (6.47)
\overline{P}=\overline{P}_{b}\;\widehat{\overline{P}} (6.48)

where as before \widehat{\overline{A}} denotes the dimensionless equivalent of the macroscopic quantity \overline{A}. The length \overline{l} is now a macroscopic length, and \overline{u} a macroscopic (seepage or Darcy) velocity. The pressure \overline{P}_{b} denotes the “breakthrough” pressure from the capillary pressure curve \overline{P}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}). It is defined as

\overline{P}_{b}=\overline{P}_{c}(\overline{S}_{b}) (6.49)

where \overline{S}_{b} is the breakthrough saturation defined as the solution of the equation

\frac{d^{2}\,\overline{P}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})}{d\overline{S}^{2}_{{\scriptscriptstyle{\mathbb{W}}}}}=0. (6.50)

Thus the dimensionless pressure is defined in terms of the inflection point (\overline{P}_{b},\overline{S}_{b}) on the capillary pressure curve, and it gives a measure of the macroscopic capillary pressure. Note that \overline{P}_{b} is process dependent, i.e. it will in general differ between imbibition and drainage. This dependence reflects the influence of microscopic wetting properties [348] and flow mechanisms on the macroscale [358].

The definition (6.48) differs from the traditional analysis [49, 329, 330, 331]. In the traditional analysis the normalized pressure field is defined as

\overline{P}=\frac{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{u}\,\overline{l}}{k}\,\widehat{\overline{P}} (6.51)

which immediately gives rise to three problems. Firstly the permeablity is a tensor, and thus a certain nonuniqueness results in anisotropic situations [339]. Secondly equation (6.51) neglects the importance of microscopic wetting and saturation history dependence. The main problem however is that Eq. (6.51) is not based on macroscopic capillary pressures but on Darcy’s law which describes macroscopic viscous pressure effects. On the other hand the normalization (6.48) is free from these problems and it includes macroscopic capillarity in the same way as the microscopic normalization (6.15) includes microscopic capillarity.

With the normalizations introduced above the dimensionless form of the macroscopic two-phase flow equations (6.42),(6.43) becomes

\displaystyle\overline{\phi}\,\frac{\displaystyle\partial\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle\partial\widehat{\overline{t}}} \displaystyle= \displaystyle\widehat{\overline{\mbox{\boldmath$\nabla$}}}\cdot\left\{{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{W}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\left(\overline{\bf Ca}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\,\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{{\scriptscriptstyle{\mathbb{W}}}}-\overline{\bf Gr}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\,\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{z}}\right)\,\right. (6.52)
\displaystyle+ \displaystyle{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{O}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\,\frac{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\:\left.\left[\left(\overline{\bf Ca}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{{\scriptscriptstyle{\mathbb{W}}}}-\overline{\bf Gr}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\mbox{\boldmath$\nabla$}}\widehat{\overline{z}}\right)\right.\right.
\displaystyle+ \displaystyle\left.\left.\overline{\bf Ca}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})+\left(1-\frac{\rho _{{\scriptscriptstyle{\mathbb{O}}}}}{\rho _{{\scriptscriptstyle{\mathbb{W}}}}}\right)\;\overline{\bf Gr}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{z}}\right]\right\}
\displaystyle\overline{\phi}\,\frac{\displaystyle\partial(1-\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})}{\displaystyle\partial\widehat{\overline{t}}} \displaystyle= \displaystyle\widehat{\overline{\mbox{\boldmath$\nabla$}}}\cdot\left\{{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\left(\overline{\bf Ca}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{{\scriptscriptstyle{\mathbb{W}}}}-\overline{\bf Gr}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{z}}\right)\,\right. (6.53)
\displaystyle+ \displaystyle{\bf K}^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{O}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\,\frac{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\:\left.\left[\left(\overline{\bf Ca}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{{\scriptscriptstyle{\mathbb{W}}}}-\overline{\bf Gr}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\overline{\mbox{\boldmath$\nabla$}}\widehat{\overline{z}}\right)\right.\right.
\displaystyle+ \displaystyle\left.\left.\overline{\bf Ca}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})+\left(1-\frac{\rho _{{\scriptscriptstyle{\mathbb{O}}}}}{\rho _{{\scriptscriptstyle{\mathbb{W}}}}}\right)\;\overline{\bf Gr}^{{-1}}_{{\scriptscriptstyle{\mathbb{W}}}}\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{z}}\right]\right\}.

In these equations the dimensionless tensor

\overline{\bf Ca}_{{\scriptscriptstyle{\mathbb{W}}}}=\frac{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\;\overline{u}\;\overline{l}}{\overline{P}_{b}}\;{\bf K}^{{-1}}=\frac{\mbox{macroscopic viscous pressure drop}}{\mbox{macroscopic capillary pressure}} (6.54)

plays the role of a macroscopic or large scale capillary number. Similarly the tensor

\overline{\bf Gr}_{{\scriptscriptstyle{\mathbb{W}}}}=\frac{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\;\overline{u}}{\rho _{{\scriptscriptstyle{\mathbb{W}}}}\; g}\;{\bf K}^{{-1}}=\frac{\mbox{macroscopic viscous pressure drop}}{\mbox{macroscopic gravitational pressure}} (6.55)

corresponds to the macroscopic gravity number.

If the traditional normalization (6.51) is used instead of the normalization (6.48), and isotropy is assumed then the same dimensionless equations are obtained with

\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}=1 (6.56)

where \overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}} is the macroscopic capillary number. Thus the traditional normalization is equivalent to the assumption that the macroscopic viscous pressure drop always equals the macroscopic capillary pressure. While this assumption is not generally valid, it sometimes is a reasonable approximation as will be illustrated below. First, however, the consequences of the traditional assumption (6.56) for the measurement of relative permeabilities will be discussed.

VI.C.3 Measurement of Relative Permeabilities

For simplicity only the isotropic case will be considered from now on, i.e. let {\bf K}=k\,{\bf 1} where {\bf 1} is the identity matrix. The tensors \overline{\bf Ca}_{{\scriptscriptstyle{\mathbb{W}}}} and \overline{\bf Gr}_{{\scriptscriptstyle{\mathbb{W}}}} then become \overline{\bf Ca}_{{\scriptscriptstyle{\mathbb{W}}}}=\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}\,{\bf 1} and \overline{\bf Gr}_{{\scriptscriptstyle{\mathbb{W}}}}=\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}\,{\bf 1} where \overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}} and \overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}} are the macroscopic capillary and gravity numbers.

The unsteady state or displacement method of measuring relative permeabilities consists of monitoring the production history and pressure drop across the sample during a laboratory displacement process [2, 359, 337]. The relative permeability is obtained as the solution of an inverse problem. The inverse problem consists in matching the measured production history and pressure drop to the solutions of the multiphase flow equations (6.52) and (6.53) using the Buckley-Leverett approximation.

In the present formulation the Buckley-Leverett approximation comprises several independent assumptions. Firstly it is assumed that gravity effects are absent, which amounts to the assumption

\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}\ll\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}. (6.57)

Secondly the viscous coupling terms are neglected, i.e.

k^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{O}}}}}\,\frac{\mu _{{\scriptscriptstyle{\mathbb{W}}}}}{\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\ll\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}\;\mbox{and}\; k^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}}\ll\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}. (6.58)

Finally the resulting equations

\overline{\phi}\,\frac{\displaystyle\partial\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle\partial\widehat{\overline{t}}}=\widehat{\overline{\mbox{\boldmath$\nabla$}}}\cdot\left\{ k^{r}_{{{\scriptscriptstyle{\mathbb{W}}}{\scriptscriptstyle{\mathbb{W}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\,\frac{\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{{\scriptscriptstyle{\mathbb{W}}}}}{\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}}\right\} (6.59)
\overline{\phi}\,\frac{\displaystyle\partial(1-\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})}{\displaystyle\partial\widehat{\overline{t}}}=\widehat{\overline{\mbox{\boldmath$\nabla$}}}\cdot\left\{ k^{r}_{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{O}}}}}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})\,\frac{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{O}}}}}\:\left[\frac{\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{{\scriptscriptstyle{\mathbb{W}}}}}{\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}}+\frac{\widehat{\overline{\mbox{\boldmath$\nabla$}}}\widehat{\overline{P}}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}})}{\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}}\right]\right\}. (6.60)

are further simplified by assuming that the term involving \widehat{\overline{P}}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}) in equation (6.60) may be neglected [332].

Combining (6.57) with the traditional normalization (6.56) yields the consistency condition

\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}\gg 1 (6.61)

for the application of Buckley-Leverett theory in the determination of relative permeabilities. It is now clear from the definition of the macroscopic gravity number, see Eq. (6.55), that the consistent use of Buckley-Leverett theory for the unsteady state measurement of relative permeabilities depends strongly on the flow regime. This is valid whether or not the capillary pressure term \widehat{\overline{P}}_{c}(\overline{S}_{{\scriptscriptstyle{\mathbb{W}}}}) in (6.60) is neglected. In addition to these consistency problems the Buckley-Leverett theory is also plagued with stability problems [360].

VI.C.4 Pore Scale to Large Scale Comparison

The comparison between the macroscopic and the microscopic dimensional analysis is carried out by relating the microscopic and macroscopic velocities and length scales. The macroscopic velocity is taken to be a Darcy velocity defined as (see discussion following equation (5.79))

\overline{u}=\overline{\phi}\, u (6.62)

where \overline{\phi} is the bulk porosity and u denotes the average microscopic flow velocity introduced in the microscopic analysis (Eq. (6.12)). The length scales l and \overline{l} are identical (\overline{l}=l).

Using these relations between microscopic and macroscopic length and time scales together with the assumption of isotropy yields

\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}=\frac{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\phi}\, u\, l}{k\,\overline{P}_{b}}=\frac{u\, l}{\overline{\nu^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}}=\frac{\sigma _{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}}\;\overline{\phi}\; l}{k\;\overline{P}_{b}}\;{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}} (6.63)

as the relationship between microscopic and macroscopic capillary numbers. Similarly one obtains

\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}=\frac{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\phi}\, u}{\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g\, k}=\frac{u}{\overline{u^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}}=\frac{\overline{\phi}\, l^{2}}{k}\,{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}} (6.64)

for the gravity numbers. Taking the quotient gives

\overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}}=\frac{\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}}{\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}}=\frac{\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g\, l}{\overline{P}_{b}}=\frac{l}{\overline{l^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}}=\frac{\sigma _{{\scriptscriptstyle{\mathbb{O}}}}{\scriptscriptstyle{\mathbb{W}}}}{l\,\overline{P}_{b}}\,{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}} (6.65)

for the macroscopic gravillary number. Note that the ratio \sigma _{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}}/(l\,\overline{P}_{b}) is the ratio of the microscopic to the macroscopic capillary pressures. The characteristic numbers

\displaystyle\overline{\nu^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} \displaystyle= \displaystyle\frac{k\,\overline{P}_{b}}{\overline{\phi}\,\mu _{{\scriptscriptstyle{\mathbb{W}}}}} (6.66)
\displaystyle\overline{u^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} \displaystyle= \displaystyle\frac{\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g\, k}{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\phi}} (6.67)
\displaystyle\overline{l^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} \displaystyle= \displaystyle\frac{\overline{P}_{b}}{\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g} (6.68)

are the macroscopic counterparts of the microscopic numbers defined in equations (6.22), (6.30) and (6.34).

An interesting way of rewriting these relationships arises from interpreting the permeability as an effective microscopic cross sectional area of flow, combined with the Leverett J-function. More precisely, let

\Lambda=\sqrt{\frac{k}{\overline{\phi}}} (6.69)

denote a microscopic length which is characteristic for the pore space transport properties. Then equations (6.63), (6.64) and (6.65) may be rewritten as

\displaystyle\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}} \displaystyle= \displaystyle\frac{\overline{l}}{\Lambda}\;\frac{{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}}{J(\overline{S}_{b})\,\cos\theta} (6.70)
\displaystyle\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}} \displaystyle= \displaystyle\frac{\overline{l}^{2}}{\Lambda^{2}}\;{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}} (6.71)
\displaystyle\overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}} \displaystyle= \displaystyle\frac{\Lambda}{\overline{l}}\;\frac{{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}}}{J(\overline{S}_{b})\,\cos\theta} (6.72)

where J(\overline{S}_{b})=(\overline{P}_{b}\sqrt{k/\overline{\phi}})/(\sigma _{{{\scriptscriptstyle{\mathbb{O}}}{\scriptscriptstyle{\mathbb{W}}}}}\cos\theta) is the value of the Leverett-J-function [28, 2] at the saturation corresponding to breakthrough, and \theta is the wetting angle.

The capillary number scales as (\overline{l}/\Lambda) while the gravity number scales as (\overline{l}/\Lambda)^{2}. Inserting (6.71) and (6.72) into (6.57) this implies that the Buckley-Leverett approximation (6.57) becomes invalid whenever \overline{l}<\Lambda{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}}/(J(\overline{S}_{b})\cos\theta).

VI.C.5 Macroscopic Estimates

The present section gives order of magnitude estimates for the relative importance of capillary, viscous and gravity effects at different scales in representative categories of porous media. These estimates illustrate the usefulness of the macroscopic dimensionless ratios for the problem of upscaling.

Three types of porous media are considered: high permeability unconsolidated sand, intermediate permeability sandstone and low permeability limestone. Representative values for \overline{\phi},k and \overline{P}_{b} are shown in Table VI.

Table VI: Representative values for porosity, permeability and breakthrough capillary pressure in unconsolidated sand, sandstone and low permeability limestone.
Quantity Sand Sandstone Limestone
\overline{\phi} 0.36 0.22 0.20  
k 10000 mD 400 mD 3 mD  
\overline{P}_{b} 2000 Pa 10^{4} Pa 10^{5} Pa  
Table VII: Definition and representative values for macroscopic dimensionless numbers in different porous media on laboratory (l_{{\rm lab}}\approx 0.1{\rm m}) and reservoir scale (l_{{\rm res}}\approx 100{\rm m}) under uniform flow conditions (u\approx 3\times 10^{{-6}}\;\mbox{m}\,\mbox{s}^{{-1}}).
Quantity Definition Unconsolidated Sand Sandstone Limestone
Laboratory Reservoir Laboratory Reservoir Laboratory Reservoir
\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}} \frac{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\phi}\, u}{\displaystyle\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g\, k} 0.01 0.13 18.6 
\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}} \frac{\displaystyle\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\phi}\, u\, l}{\displaystyle k\,\overline{P}_{b}} 0.005 4.9 0.015 15.0 0.19 187.5 
\overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}} \frac{\displaystyle\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g\, l}{\displaystyle\overline{P}_{b}} 0.5 492 0.1 115 0.01 10 
\Lambda \sqrt{\frac{\displaystyle k}{\displaystyle\overline{\phi}}} 5.2\,\mu{\rm m} 1.3\,\mu{\rm m} 0.1\,\mu{\rm m} 
{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}/\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}} \frac{\displaystyle\Lambda}{\displaystyle l}\, J(\overline{S}_{b})\cos\theta 1.5\cdot 10^{{-4}} 1.5\cdot 10^{{-7}} 4.8\cdot 10^{{-6}} 4.8\cdot 10^{{-9}} 2.8\cdot 10^{{-7}} 2.8\cdot 10^{{-10}} 
\frac{\displaystyle{\rm Ca}_{c}\,\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}}{\displaystyle{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}} \frac{\displaystyle\Lambda}{\displaystyle l}\, J(\overline{S}_{b})\cos\theta 0.67 --- 0.63 --- 0.71 --- 

To estimate the dimensionless numbers the same microscopic velocity u\approx 3\times 10^{{-6}}\;\mbox{m}\,\mbox{s}^{{-1}} as for the microscopic estimates will be used. The length scale l, however, differs between a laboratory displacement and a reservoir process. l_{{\rm lab}}\approx 0.1{\rm m} and l_{{\rm res}}\approx 100{\rm m} are used as representative values. Combining these values with those in Table IV and VI yields the results shown in Table VII.

The first row in Table VII can be used to check the consistency of the Buckley-Leverett approximation with the traditional normalization. The consistency condition (Eq. (6.61)) is violated for unconsolidated sand and sandstones. Such a conclusion, of course, assumes that the values given in Table VI are representative for these media.

The fifth row in Table VII gives the ratio between macroscopic and microscopic capillary numbers which according to Eq. (6.63) is length scale dependent. The last row in Table VII compares this ratio to the typical critical capillary number {\rm Ca}_{c} reported for laboratory desaturation curves in the different porous media. Using the {\rm Ca}_{c}\approx 10^{{-4}} for sand, {\rm Ca}_{c}\approx 3\cdot 10^{{-6}} for sandstone, and {\rm Ca}_{c}\approx 2\cdot 10^{{-7}} for limestone [28] as before one finds that the corresponding critical macroscopic capillary number is close to 1. This indicates that the macroscopic capillary number is indeed an appropriate measure of the relative strength of viscous and capillary forces.

As a consequence one expects differences between residual oil saturation S_{{{\scriptscriptstyle{\mathbb{O}}}r}} in laboratory and reservoir floods. Given a laboratory measured capillary desaturation curve S_{{{\scriptscriptstyle{\mathbb{O}}}r}}({\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}) as a function of the microscopic capillary number {\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}} the analysis predicts that the residual oil saturation in a reservoir flood can be estimated from the laboratory curve as S_{{{\scriptscriptstyle{\mathbb{O}}}r}}({\rm Ca}_{c}\cdot\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}) [47, 48]. For \overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}>1 the S_{{{\scriptscriptstyle{\mathbb{O}}}r}} value based on macroscopic capillary numbers will in general be lower than the value S_{{{\scriptscriptstyle{\mathbb{O}}}r}}({\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}) expected from using microscopic capillary numbers. Such differences have been frequently observed, and Morrow [361] has recently raised the question why field recoveries are sometimes significantly higher than those observed in the laboratory. The revised macroscopic analysis of [47, 48] suggests a possible answer to this question.

The values of the dimensionless numbers in Table VII allow an assessment of the relative importance of the different forces for a displacement. To illustrate this consider the values \overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}=0.01, \overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}=0.005 and \overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}}=0.5 for unconsolidated sand on the laboratory scale. A moments reflection shows that this implies V\ll G\approx C where V stands for macroscopic viscous forces, C for macroscopic capillary forces, and G for gravity forces. The notation A\ll B indicates that A/B<10^{{-2}} while A<B means 10^{{-2}}<A/B<0.5 and A\approx B stands for 0.5<A/B<2. Repeating this for all cases in Table VII yields the results shown in Table VIII. Table VIII contains also the results from the microscopic dimensional analysis, as well as the results one would obtain from a traditional macroscopic dimensional analysis which assumes \overline{\rm Ca}=1 (see Eq. (6.56)).

Table VIII: Relative importance of viscous (V), gravity (G) and capillary (C) forces in unconsolidated sand, sandstone and limestone. The notation A\ll B (with A,B\in\relax\{ V,G,C\relax\}) indicates that A/B<10^{{-2}} while A<B means 10^{{-2}}<A/B<0.5 and A\approx B stands for 0.5<A/B<2.
Sand Sandstone Limestone
pore scale V\ll G\ll C 
traditional analysis V=C\ll G V=C<G G<V=C 
large scale [47] laboratory scale V\ll G\approx C V<G<C G<V<C 
[48] field scale C<V\ll G C<V<G C<G<V 

Obviously, the relative importance of the different forces may change depending on the type of medium, the characteristic fluid velocities and the length scale. Perhaps this explains part of the general difficulty of scaling up from the laboratory to the reservoir scale for immiscible displacement.

VI.C.6 Applications

The characteristic macroscopic velocities, length scales and kinematic viscosities defined respectively in equations (6.66), (6.67) and (6.68) are intrinsic physical characteristics of the porous media and the fluid displacement processes. These characteristics can be useful in applications such as estimating the width of a gravitational segregation front, the energy input required to mobilize residual oil or gravitational relaxation times.

The macroscopic gravillary number \overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}} defines an intrinsic length scale \overline{l^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} (see Eq. (6.65)). Because \overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}} gives the ratio of the gravity to the capillary forces the length \overline{l^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} directly gives the width of a gravitational segregation front when the fluids are at rest and in gravitational equilibrium, i.e. when viscous forces are negligible or absent. Using the same estimates for \overline{\phi},k and \overline{P}_{b} as those used for Table VII one obtains a characteristic front width of 20\,cm for unconsolidated sand, 1\,m for sandstone, and roughly 10\,m for a low permeability limestone.

Similarly, the macroscopic capillary number defines an intrinsic specific action (or energy input) \overline{\nu^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} via Eq. (6.63) which is the energy input required to mobilize residual oil if gravity forces may be considered negligible or absent. Represenative estimates are given in Table IX below.

The gravitational relaxation time is the time needed to return to gravitational equilibrium after its disturbance. This may be defined from the balance of gravitational forces versus the combined effect of viscous and capillary forces. Analogous to equation (6.35) for the microscopic case the dimensionless ratio becomes

\displaystyle\frac{\overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}}}{\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}} \displaystyle= \displaystyle\frac{(\mbox{macr. gravitational pressure})^{2}}{(\mbox{macr. capillary pressure})\times(\mbox{macr. viscous pressure drop})} (6.73)
\displaystyle= \displaystyle\frac{\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}}}{\overline{\rm Gr}^{2}_{{\scriptscriptstyle{\mathbb{W}}}}}=\frac{\rho^{2}_{{\scriptscriptstyle{\mathbb{W}}}}\, g^{2}\, k\, l}{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\phi}\,\overline{P}_{b}\, u}=\frac{t}{\overline{t^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}}

which defines the gravitational relaxation time \overline{t^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} as

\overline{t^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}=\frac{\overline{l^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}}{\overline{u^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}}=\frac{\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\phi}\,\overline{P}_{b}}{\rho^{2}_{{\scriptscriptstyle{\mathbb{W}}}}\, g^{2}\, k}. (6.74)

Estimated values are given in Table IX. They correspond to gravitational relaxation times of roughly 10\, minutes for unconsolidated sand, 13 hours for a sandstone and 736 days for a low permeability limestone.

Another interesting intrinsic number arises from comparing the strength of macroscopic capillary forces versus the combined effect of viscous and gravity forces

\displaystyle(\overline{\rm Gl}_{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{\rm Ca}_{{\scriptscriptstyle{\mathbb{W}}}})^{{-1}} \displaystyle= \displaystyle\frac{(\mbox{macr. capillary pressure})^{2}}{(\mbox{macr. grav. pressure})\times(\mbox{macr. viscous pressure drop})} (6.75)
\displaystyle= \displaystyle\frac{\overline{\rm Gr}_{{\scriptscriptstyle{\mathbb{W}}}}}{\overline{\rm Ca}^{2}_{{\scriptscriptstyle{\mathbb{W}}}}}=\frac{k\,\overline{P}^{2}_{b}}{\overline{\phi}\,\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g\, u\, l^{2}}=\frac{\overline{Q^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}}{Q}

where Q denotes the volumetric flow rate. Thus \overline{Q^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} defined as

\overline{Q^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}=\overline{l^{*}}^{2}_{{\scriptscriptstyle{\mathbb{W}}}}\,\overline{u^{*}}_{{\scriptscriptstyle{\mathbb{W}}}}=\frac{k\,\overline{P}^{2}_{b}}{\overline{\phi}\,\mu _{{\scriptscriptstyle{\mathbb{W}}}}\,\rho _{{\scriptscriptstyle{\mathbb{W}}}}\, g} (6.76)

is an intrinsic system specific characteristic flow rate. The estimates for \overline{\nu^{*}}_{{\scriptscriptstyle{\mathbb{W}}}},\overline{u^{*}}_{{\scriptscriptstyle{\mathbb{W}}}},\overline{l^{*}}_{{\scriptscriptstyle{\mathbb{W}}}},\overline{t^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} and \overline{Q^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} are summarized in Table IX.

Table IX: Characteristic macroscopic energies, velocities, length scales, time scales and volumetric flow rates for oil-water flow under reservoir conditions in unconsolidated sand, sandstone and low permeability limestone.
Quantity Sand Sandstone Limestone
\overline{\nu^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} 6.1\cdot 10^{{-5}}\,{\rm m}^{2}{\rm s}^{{-1}} 2.0\cdot 10^{{-5}}\,{\rm m}^{2}{\rm s}^{{-1}} 1.6\cdot 10^{{-6}}\,{\rm m}^{2}{\rm s}^{{-1}} 
\overline{u^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} 2.99\cdot 10^{{-4}}\,{\rm ms}^{{-1}} 2.17\cdot 10^{{-5}}\,{\rm ms}^{{-1}} 1.61\cdot 10^{{-7}}\,{\rm ms}^{{-1}} 
\overline{l^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} 0.2\,{\rm m} 1.02\,{\rm m} 10.2\,{\rm m} 
\overline{t^{*}}_{{\scriptscriptstyle{\mathbb{W}}}} 669\,{\rm s} 4.7\cdot 10^{4} s 6.36\cdot 10^{7} s  
\overline{Q^{*}} 1.22\cdot 10^{{-6}}\,{\rm m}^{3}{\rm s}^{{-1}} 2.04\cdot 10^{{-5}}\,{\rm m}^{3}{\rm s}^{{-1}} 1.63\cdot 10^{{-5}}\,{\rm m}^{3}{\rm s}^{{-1}} 

In summary, the dimensional analysis of the upscaling problem for two phase immiscible displacement suggests to normalize the macroscopic pressure field in a way which differs from the traditional normalization. This gives rise to a macroscopic capillary number \overline{\rm Ca} which differs from the traditional microscopic capillary number {\rm Ca} in that it depends on length scale and the breakthrough capillary pressure \overline{P}_{b}. The traditional normalization corresponds to the tacit assumption that viscous and capillary forces are of equal magnitude. With the new macroscopic capillary number \overline{\rm Ca} the breakpoint {\rm Ca}_{c} in capillary desaturation curves seems to occur at \overline{\rm Ca}\approx 1 for all types of porous media. Representative estimates of \overline{\rm Ca} for unconsolidated sand, sandstones and limestones suggest that the residual oil saturation after a field flood will in general differ from that after a laboratory flood performed under the same conditions. Order of magnitude estimates of gravitational relaxation times and segregation front widths for different media are consistent with experiment.