Geomechanical reservoir models

Geomechanical models simulate changes in rock properties with time and pressure. In performance-based reservoir analysis, a geomechanical model may be coupled with a flowing material balance, typecurve, analytical, numerical, or unconventional reservoir module analysis. In addition, geomechanical models are useful for, and may be necessary additions to, analyses of overpressured reservoirs. The primary properties of interest are porosity, formation compressibility (pore compressibility), and permeability. The impact of these properties on the fluid flow model are discussed in the following sections.

Porosity and formation compressibility

Pore volume and pore volume compressibility are related through the following equation:

In normally pressured reservoirs, the pore volume change with pressure is considered minimal, and thus the pore volume (formation) compressibility retains a very small, constant value. However, in overpressured reservoirs, the natural compaction process is incomplete, as a large portion of the overburden remains supported by high internal pore pressure. As this pressure is released through fluid production, the pore space may reduce significantly. Thus, under overpressured conditions, formation compressibility is relatively large, and may have significant variation with pressure. This becomes important to the material balance equation, as pore compressibility is a significant energy term in overpressured reservoirs.

In normally pressured gas reservoirs, the energy of the formation is usually negligible, compared to the energy of the gas.  In oil reservoirs, the formation compressibility may be significant at any pressure. The graph below compares the energies (compressibilties) of formation and fluid (gas) over a large pressure range. In areas where they are near the same order of magnitude, the formation compressibility cannot be ignored.

For overpressured reservoirs, there are very few published correlations. Most appear to be highly specific to area and rock type, and thus are not useful for universal application.

Hall

The most well-known and widely used correlation for formation compressibility was developed by Hall, and is a function of porosity. The Hall correlation is based on laboratory data and is considered reasonable for normally pressured sandstones. It tends to under-predict formation compressibility under high pressure conditions. This correlation is used when calculating initial formation compressibility on the Basic Reservoir Properties view.

Dobrynin

The Dobrynin (1962) correlation can predict a relative change in formation compressibility over a range of net overburden pressures. It is designed for overpressured reservoirs. The correlation assumes a semi-logarithmic relationship between formation compressibility and net overburden pressure between two limiting pressures, pmin and pmax. However, it does not predict the initial formation compressibility.


Laurent et al. (1993) propose an empirical correlation for Biot Coefficient (α):

The Dobrynin correlation is considered a reasonable predictor of variation in formation compressibility, provided that the following conditions are met:

  • Net overburden pressure must lie within the range of pmin and pmax (usually 150 to 30,000 psi)
  • Overburden pressure must be higher than the product of the initial pressure and Biot’s number (initial net overburden pressure must be greater than zero)

The Dobrynin correlation may be used to generate a continuous compressibility profile using a known formation compressibility value:

  • Initial formation compressibility (usually between 10 and 80 microcips for overpressured reservoirs)
  • Overburden pressure gradient (Pgrd, psi/ft)
  • Reservoir depth (ft)

Formation compressibility is then solved for any value of pressure.

Pressure-dependent permeability

In the standard pressure transient equations, permeability is usually considered to be constant. However, there are several situations in which this assumption may not be a valid:

  • Compaction in overpressured reservoirs
  • Very low permeability reservoirs (in general)
  • Unconsolidated and/or fractured formations

Variable permeability is an important (and sometimes critical) consideration, with regards to the fluid flow model. As a result of pore volume reduction during fluid withdrawal, the available flow area is reduced, and thus the permeability decreases with pressure. As with formation compressibility, the variable permeability effect is greatest in overpressured conditions. The process of pore volume collapse in overpressured reservoirs is referred to as compaction.

One way to account for a variable permeability over time is to modify the definition of pseudo-pressure and pseudo-time. Another way to account for variable permeability is through the use of formation compressibility correlations.

Yilmaz and Nur

Yilmaz and Nur present a generic correlation, suitable for the pressure variable permeability observed in extremely low permeability reservoirs. The correlation introduces a "permeability modulus" (γ) which has a form identical to compressibility. The magnitude of the permeability modulus is a strong function of the formation compressibility.

Dobrynin

Dobrynin presents a more detailed correlation, designed for overpressured reservoirs, which relates the fractional change in permeability due to formation compressibility, overburden pressure, and a dimensionless correlation factor (γ). Between a certain minimum and maximum pressure range (pmin and pmax), this correlation produces a semi-log straight line.

where:

γ is 0.33 for uniform pore distribution or very high pore compressibility, higher for poorly sorted sandstones with low pore compressibility.

The equation below relates the permeability at any pore pressure to that at conditions of zero net overburden pressure:

It is more suitable to have an equation that relates permeability at any pressure to permeability at initial conditions, thus we rewrite the above as follows.

where:

Biot’s Number default: 0.85 to 1

Maximum Formation Compressibility default: 30 to 80 microcips

Pmax = 30,000 psi

Pmin = 150 psi

Geomechanical effects for oil analytical models

This content applies to water analytical models as well.

Background

Geomechanical effects are useful for (and may be necessary additions to) the analysis of over-pressured reservoirs. The major geomechanical effects are:

  • Pressure-dependent permeability effect (Yilmaz and Nur correlation)
  • Pressure-dependent formation compressibility effect (Dobrynin correlation)

Permeability changes are incorporated into a special formulation of pseudo-pressure for oil, but they also enter into the definition of pseudo-time. In addition, formation compressibility changes are incorporated into the definition of pseudo-time.

Implementation notes: 2014v2 release

For the limited scope of the 2014v2 release, only the pseudo-pressure aspect of pressure-dependent permeability has been implemented. To be consistent with standard analytical model formulations, the pseudo-pressure for oil assumes a constant formation volume factor and constant viscosity.

Pseudo-time (which is affected by both changing permeability and by compressibility) has not been implemented.

Consequences:

  • By using regular time instead of pseudo-time for oil, forecasts tend to be more conservative. The degree of conservativeness is directly related to the combined effect of permeability, and formation compressibility variations with pressure.
  • By using constant oil properties, slight changes in pressure, volume, temperature (PVT) are not be tracked (this is a relatively minor issue).

Recommended Usage:

  • For oil reservoirs that stay above the bubble point
  • For oil reservoirs that are at or below the bubble point, but with total compressibility that accounts for gas coming out of solution.