Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information

# Rock density and porosity

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

Understanding the density and porosity of reservoir rocks is a key factor in estimating their hydrocarbon potential. Density and porosity are related.

## Density

Density is defined as the mass per volume of a substance.

....................(1)

typically with units of g/cm3 or kg/m3. Other units that might be encountered are lbm/gallon or lbm/ft3 (see Table 1).

For simple, completely homogeneous (single-phase) material, this definition of density is straightforward. However, earth materials involved in petroleum engineering are mixtures of several phases, both solids (minerals) and fluids. Rocks, in particular, are porous, and porosity is intimately related to density.

## Porosity

Porosity (Φ) is defined as the nonsolid or pore-volume fraction.

....................(2)

Porosity is a volume ratio and thus dimensionless, and is usually reported as a fraction or percent. To avoid confusion, particularly when variable or changing porosities are involved, it is often reported in porosity units (1 PU = 1%).

Several volume definitions are required to describe porosity:

....................(3)

From these we can define the various kinds of porosity encountered:

....................(4)

## Relationship of density and porosity

Fig. 1 shows the appearance of these types of porosity in a sandstone.

Similarly, the definitions of the standard densities associated with rocks then follows:

....................(5)

where Ms, Md, Msat, Mb, and Mfl are the mass of the solid, dry rock, saturated rock, buoyant rock, and fluid, respectively.

The density of a composite such as rocks (or drilling muds) can be calculated from the densities and volume fraction of each component. For a two-component system,

....................(6)

where ρmix is the density of the mixture; ρA is the density of Component A; ρB is the density of B; A and B are the volume fractions of A and B respectively (and so B = 1− A).

Expanding this into a general system with n components,

....................(7)

For example, exploiting Eqs. 4, 5, and 6 for a rock made up of two minerals, m1 and m2, and two fluids, f 1 and f 2, we find

....................(8) and ....................(9)

Eq. 8 is a fundamental relation used throughout the earth sciences to calculate rock density. Given a porosity and specific fluid, density can be easily calculated if the mineral or grain density is known. Grain densities for common rock-forming minerals are shown in Table 2. The result of applying Eq. 9 is shown in Fig. 1.

Note in Table 2 that there are several densities reported for the same mineral group, such as feldspar or clay. The density will change systematically as composition varies. For example, in the plagioclase series, the density increases as sodium (albite, ρ = 2.61 g/cm3) is replaced by calcium (anorthite, ρ = 2.75 g/cm3). The most problematic minerals are clays, particularly expanding clays (montmorillonite or smectite) capable of containing large and variable amounts of water. In this case, densities can vary 40% or more. This is a particular problem, because clays are among the most common minerals in sedimentary rocks.

Reservoir rocks often contain significant amounts of semisolid organic material such as bitumen. These will typically have light densities similar in magnitude to those of coals.

Pore fluid densities are covered in detail in Pore fluid properties.

## In-situ density and porosity

In general, density increases and porosity decreases monotonically with depth. This is expected, because differential pressures usually increase with depth. As pressure increases, grains will shift and rotate to reach a more dense packing. More force will be imposed on the grain contacts. Crushing and fracturing is a common result. In addition, diagenetic processes such as cementation work to fill the pore space. Material may be dissolved at point contacts or along styolites and then transported to fill pores. Some of the textures resulting from these processes were seen in the photomicrographs of Rock types. In Fig. 2[1][2][3][4][5][6][7][8], generalized densities as a function of depth for shales are plotted. The shapes and overall behaviors for these curves are similar, even though they come from a wide variety of locations with different geologic histories. These kinds of curves are often fit with exponential functions in depth to define the local compaction trend.

Differential or effective pressures do not always increase with increasing depth. Abnormally high pore fluid pressures ("overpressure") can occur because of:

• Rapid compaction
• Low permeability
• Mineral dewatering
• Migration of high-pressure fluids

The high pore pressure results in an abnormally low differential of effective pressure. This can retard or even reverse the normal compaction trends. Such a situation is seen in Fig. 3[9]. Porosities for both shales and sands show the expected porosity loss with increasing depth in the shallow portions. However, at about 3500 m, pore pressure rises and porosity actually increases with depth. This demonstrates why local calibration is needed. It also indicates the pressure dependence of rock properties.

## Measurement techniques

### Laboratory

Numerous methods can be used in the laboratory to determine porosity and density. The most common are by saturation weight and Boyle’s law. For rocks without sensitive minerals such as smectites, the porosity and dry, grain, and saturated densities can be derived from the saturated mass, dry mass, and volume (or buoyant weight). These measurements allow calculation of saturated, dry, and grain density as well as porosity and mineral and pore volume by employing Eqs. 3 through 5.

The Boyle’s law technique measures the relative changes in gas pressures inside a chamber with and without a rock specimen. The internal (connected) pore volume is calculated from these variations in pressure, from which porosities and densities are extracted.

### Logging

Several logging techniques are available to measure density or porosity.[10][11] These indirect techniques can have substantial errors depending on borehole conditions, but they do provide a measure of the in-situ properties. Gamma ray logs bombard the formation with radiation from an active source. Radiation is scattered back to the logging tool, depending on the electron density of the material. Formation density is extracted from the amplitude of these back-scattered gamma rays. The neutron log estimates porosity by particle interaction with hydrogen atoms. Neutrons lose energy when colliding with hydrogen atoms, thus giving a measure of the hydrogen content. Because most of the hydrogen in rocks resides in the pore space (water or oil), this is then related to the liquid-filled porosity. Note that the neutron log will include bound water within clays as porosity. In addition, when relatively hydrogen-poor gas is the pore fluid, the neutron log will underestimate porosity. In a similar fashion, the nuclear magnetic resonance (NMR) log will resolve the hydrogen content. This tool, however, has the ability to differentiate between free bulk water and bound water. Sonic logs are also used for porosity measurements, particularly when anomalous minerals (such as siderite) or borehole conditions render other tools less accurate. The technique involves inverting velocity to porosity using one of the relationships provided in Elastic wave velocities. Gravimetry has also been used downhole to measure variations in density. Although this tool is insensitive to fine-scale changes, it permits density measurement far out into the formation.

### Seismic

On a coarse scale, densities can sometimes be extracted from seismic data. This method requires separating the density component of impedance. This normally requires an analysis of the seismic data as a function of offset or reflection angle. This technique will probably see more use as seismic data improves and is further incorporated into reservoir description.

## Nomenclature

 Vpor = total pore volume, m3 or cm3 Vp-con = connected pore volume, m3 or cm3 Vp-iso = isolated pore volume, m3 or cm3 Φ = porosity Φfx = fracture porosity Φp-e = effective porosity Φp-iso = isolated, ineffective porosity ρ = density, kg/m3 or g/cm3 ρb = bulk density, kg/m3 or g/cm3 ρB = brine density, kg/m3 or g/cm3 ρd = dry density, kg/m3 or g/cm3 ρfl = fluid density, kg/m3 or g/cm3 ρg = grain or mineral density, kg/m3 or g/cm3 ρG = gas density, kg/m3 or g/cm3 ρO = oil density, kg/m3 or g/cm3 ρsat = saturated density, kg/m3 or g/cm3 ρW = water density, kg/m3 or g/cm3 M = molecular weight, g/mole Af 1, Af 2 = fraction fluid component 1, 2, etc. Am1, Am2 = fraction mineral component 1, 2, etc. A1, A2 = fraction component 1, 2, etc.

## References

1. Castagna, J.P., Batzle, M.L., and Kan, T.K. 1993. Rock physics—the link between rock properties and AVO response. In Offset-Dependent Reflectivity—Theory and Practice of AVO Analysis, ed. P. Castagna and M.M. Backus, No. 8, 124–157. Tulsa, Oklahoma: Investigations in Geophysics series, Society of Exploration Geophysicists.
2. Rieke III, H.H. and Chilingarian, G.V. 1974. Compaction of Argillaceous Sediments. Amsterdam, The Netherlands: Elsevier Scientific Publishing Company.
3. McCulloh, T.H. 1967. Mass properties of sedimentary rocks and gravimetric effects of petroleum and natural-gas reservoirs. USGS Professional Paper 528-A, Department of the Interior, United States Geological Survey, Washington, DC http://pubs.usgs.gov/pp/0528a/report.pdf.
4. Storer, D. 1959. Compaction of the argiuaceous sediments in the Padano Basin. In The Gasiferous Deposits of Western Europe, Vol. 2, 519–536. Roma, Italy: Accademia Nazionale dei Lincei.
5. Dickinson, G. 1953. Geological Aspects of Abnormal Reservoir Pressures in Gulf Coast Louisiana. AAPG Bull. 37 (2): 410-432.
6. Eaton, B.A. 1969. Fracture gradient Prediction and Its Application in Oilfield Operations. J Pet Technol 21 (10): 1353–1360. SPE-2163-PA. http://dx.doi.org/10.2118/2163-PA.
7. Dallmus, K.F. 1958. Mechanics of basin evolution and its relation to the habitat of oil in the basin. In Habitat of Oil, L.G. Weeks, No. 36, 2071–2174. Tulsa, Oklahoma: AAPG Memoir, AAPG.
8. Skeels, C. 2001. Hydrol. Processes 15 (25 October 2001): 3073.
9. Stuart, C.A. 1970. Geopressures. Presented at the Louisiana State University 1970 Symposium on Abnormal Subsurface Pressure, Baton Rouge, Louisiana, USA.
10. Schlumberger. 1985. Log Interpretation Charts, 1985 edition. Sugar Land, Texas: Schlumberger.
11. Hurst, A., Griffiths, C.M., and Worthington, P.F. 1992. Geologic Applications of Wireline Logs II, No. 65. Bath, UK: Geological Society Publishing House.