Basic measurements of the core before and after the 1090 days long-term test [3]. **Figure 1** displays SEM images of the core material before and after the flow-through test.

## Abstract

This chapter presents the constitutive equations necessary to interpret laboratory and field data when both solid and pore volume evolve through time due to chemical and mechanical processes. The equations for the porosity evolution that are developed are generic, but the examples presented are acquired from chalk core studies. The processes at play when porosity is subject to change due to volumetric compaction and fluid-rock interactions when porous chalks are continuously flooded are presented here. As the overall solid mass is a conserved quantity, the void space is not. Constitutive equations are therefore required to estimate the time-evolution of the porosity. Laboratory triaxial tests were performed on high-porosity outcrop chalks from Obourg, Liegè, and Mons (Belgium). These tests are being compacted and continuously flooded with MgCl2 brine at elevated temperature and at high stresses. As calcite is replaced by magnesite, the overall mass and solid density change, thereby changing the volume of the solid. At the same time, the bulk volume is changing. Taking both effects into consideration, the pore volume evolution can be determined. We find that the porosity changes in nonintuitive ways as the relative importance of bulk compaction and chemical interaction may vary over time.

### Keywords

- dynamic porosity
- chalk
- dissolution
- precipitation
- deformation
- compaction

## 1. Introduction

Understanding how chemically reactive and mechanically deformable rock-fluid systems consisting of solids and voids evolve is very important to several fields in the Earth sciences. Examples include the lithification of sedimentary strata [1] and long-term creep behavior of crustal rocks [2]. In addition, a range of industrial processes are affected by chemo-mechanical interactions, including, e.g., pharmaceutical and food processing industries, and geotechnical engineering concerning roadwork construction, mass transportation, and slope stability.

The pore volume fraction, the pore size distribution, and the mineral surfaces are key parameters to ensure safe disposal of radioactive waste and captured CO_{2}, and to understand how ores’ deposit evolves with time. In petroleum sciences, chemo-mechanical processes are important to accurately predict the porosity since it is inside the pores where hydrocarbons are stored, and it is through the pores, the hydrocarbons are being produced by miscible and immiscible fluid migration across reactive mineral surfaces that again are subject to change. Both pore volume and production rate are crucial to determine the recoverable hydrocarbon potential.

Reactive pore fluids in nonequilibrium with their host rocks lead to dissolution and precipitation transforming the mineral assembly into another, see for example, [3, 4, 5]. Dissolution and precipitation lead to changes to grain texture and morphology, and the mineral surface’s affinity to oil and water. These are factors, together with porosity, that dictate not only the flow property of the porous rock but also the mechanical parameters that control the stiffness, strength, and the rate at which compaction by grain reorganization and pore collapse occurs. The general processes that are described here is applied to understand how the porosity of chalks develops (Figure 1a, b display unaltered and altered chalks) dynamically in a controlled triaxial cell experiments (Figure 2a), with control of temperature, pore pressure, side stress, and overburden stress of cylindrical samples (Figure 2b).

This chapter deals with some of the constitutive relations that are used to describe the evolution of porous bodies. We incorporate a discussion of how rock-fluid chemistry may impact the grain volume, and review the ways in which total volume reduction may be facilitated in compressive systems. The discussion summarized in the development of a porosity evolution equation in which all effects are included. The usage of the porosity evolution equation is exemplified with references to already published experimental results.

## 2. Constitutive equations for porosity evolution

The basic equations that are used to quantify the porosity evolution through time are presented. The analysis is based on the work presented in Nermoen, et al. [3]. The overall bulk volume of a bi-phase material equals the sum of the solid volume and pore volume

Any changes in solid volume and pore volume lead to changes in the bulk

The pore volume, and hence the porosity, itself is not a conserved quantity. In that case, the bulk volume (size of the object of study) and the solid volume evolution have to be used. Since the volumes are additive by nature, the changes in pore volume can be calculated

At any given time through dynamic porosity evolution, the porosity is given by

When both the bulk volume and the pore volume change dynamically from known measurements before the experiment starts (

Using Eq. (3) enables the determination of the porosity from known quantities

This equation is useful when determining pore volume evolution when considering mechanical and chemical processes that occur at reactive rock-fluid systems exposed to elevated stresses. To simplify the porosity evolution equation further, the volumetric strain and the initial porosity before chemo-mechanical processes occur are introduced

The minus sign in the volumetric strain here are in line with the definition in geotechnical engineering that inward deformation is positive, often different from other fields of sciences. Dividing by the initial bulk volume and employing the definitions Eq. (6) become

Eq. 8 is used to analyze how the pore volume fraction changes as the overall volume and the solid volume changes through time. Typically, it is easier to quantify the changes in the solid volume and total volume because of conservation of mass, but this does not generally apply. In other cases, when the pore volume and solid volume are known, the porosity can be calculated from

This equation could be used when the volumes of injected and produced fluid volumes are monitored and solid volume change can be back-calculated from ion chromatography (IC) of produced fluids. If, however, the bulk volume (e.g., 4D seismic) and the pore volume were obtained from monitoring the injected and produced fluid *volumes*, the porosity is as follows:

## 3. Volumetric strain by imposed stress

In compressive hydrostatic systems, the porous rocks deform by reducing its bulk volume. This may affect the porosity through, for example, Eq. (8). In closed systems, in which the mass and density of the minerals are conserved, the bulk volume reduction equals the pore volume reduction. This is facilitated by grains moving relative to each other, and/or by pressure solution (dissolution of stressed grain contacts and precipitation in unstressed parts of the mineral framework [6]). In open systems subjected reactive flow, both mass and density of the core material may change because of mineral reactions. To evaluate the relative importance of how evolution mechanisms of the solid volume and pore volume dictate the porosity in real systems, a rigorous definition of stresses and strains are required. The aim is to pave the way for quantitative analyses of how stresses impact strains, and how strains impact the porosity chemo-mechanical compaction.

### 3.1. The stress tensor in porous materials

The stress tensor describes the stresses (force per unit area) in a solid porous body. For cylindrical core plug, it is convenient to express the stress tensor

Shear and normal components are abbreviated

In reservoir systems, however, all stress components may apply, and as such, the off-diagonal elements of the stress tensor are nonzero. However, in these cases, the stress tensor can be rotated such that the principal stress notation can be obtained. It is customary procedure to arrange the first, second, and third principal stress directions as

### 3.2. Effective stress

In porous rocks, it is the effective stresses introduced by [7] that drive deformation. The external load applied onto a material that consists of solids and voids is balanced by the interparticle contacts in force networks (material framework) and a fraction

Simultaneously, seepage forces arising from differences in fluid pressure expose a net force onto the solid framework (see Figure 3). In partially consolidated systems, in which the cross area is given by the sum of the consolidated area (solid–solid area) and the area of the fluid-to-solid contact area (

Here, the largest and smallest horizontal stress is abbreviated with an index

In core scale experiments, the directions perpendicular to the

### 3.3. Defining strain

The most commonly used definition of strain, applicable to small finite deformations, is presented here. It is acknowledged that other definitions of strain also exist in the scientific literature. The strain at any time is given by the ratio of elongation divided by the initial length

Stresses may deform Earth materials so that two initially orthogonal directions change by an angle

In three dimensions (cylindrical coordinates), the pairs of shear and normal strains are organized in the strain tensor

Similar to the stress tensor, the shear strains balance each other (

To estimate the porosity evolution, bulk volumetric strain has to be used. The volumetric strain equals the change in volume divided by the initial volume, which is the first strain invariant, given by

If the length and diameter of cylindrical cores are being measured continuously, then the volumetric strain can be estimated. Typically, for small strains, the second and third order terms can be omitted, hence,

## 4. Partitioning time-independent and time-dependent deformation

The volumetric strain can be split into an immediate strain, occurring when the effective stress is being changed, and time-dependent deformation. The two cases are presented briefly in the following sections, even though this is a large area of research. For the time-independent case, Hooke’s law is described before nonlinear models are presented, followed by a short note on plasticity and other failure mechanisms before time-dependent models are described.

### 4.1. Elastic strain: linear elasticity

Hooke’s law is the simplest relation to describe the relation between the stress-strain tensors. It assumes that the deformation is immediate, linear, and reversible. In continuous media, for small stress and strain increments in the linear limit, the

By adding up the three equations expressed in the matrix form earlier

we may use this equation to define the bulk modulus in hydrostatic tests. When omitting higher order terms in the volumetric strain (Eq. (18)), the left hand side of Eq. (20) equals the volumetric strain. For hydrostatic tests, in which the stresses in all spatial directions equal,

In Eq. (21), the bulk modulus (

### 4.2. The effective stress changes that drive deformation

Within the elastic domain, any change in the effective stress drive deformation in the sample, from here on abbreviated with the

In Eq. (23), the underlying assumption is that Young’s modulus and Poisson ratio remain fixed. Furthermore, when pore pressure is included, the effective stress changes due to both axial and radial stress and pore pressure

It is assumed that the Biot stress coefficient remains fixed during loading. Using these definitions into Eq. (32) enables us to fully describe the relation between the stress, pore pressure, and strain in Hooke’s law

### 4.3. Plasticity and irreversible deformation

For a highly porous chalk, a nonzero component of the observed strain is always irreversible when the load is released as exemplified by [8] and [9] where irreversible plasticity is seen also within the ‘elastic’ phase of the

## 5. Time-dependent pore volume reduction processes and compaction

To understand how the mechanical and chemical processes affect the porosity during pore collapse, it is important to take a closer look at how the observed bulk strain can be partitioned. In this section, we consider the simplest possible partition (below) in which the overall strain is partitioned additively into a solid volume and a pore volume component. The relative importance of these mechanisms may be found by the analysis of quantitative measurements of the bulk volume change and the change in the solid volume due to the dissolution/precipitation as the mineral mass and density change over time, while grains reorganize, crush, and solid contacts evolve. In this case, the observed volumetric change can be partitioned additively via

This does not imply that imply that cross terms do not exist in which: (a) the rate of pore volume reduction is sensitive to the reduction in solid volume and (b) how the solid volume rate may depend on how grains reorganize to change the flow pattern and potentially expose new fresh mineral surfaces to the reactive brine. It is likely to assume that based on the accelerated strain presented in [10] minute changes to the solid volume increase the rate of pore collapse (also seen in [4]).

Given the simple partitioning above, a model can be developed to describe the observed creep curve with a few physical parameters (see Eqs. (20)–(23) in [10]). In this model, overall volumetric strain is additively partitioned into a pore and solid volume component in which the pore volume equals,

The solid volume rate is assumed to be constant (^{3}/day determined from ion chromatography data). For Mons chalk at 130°C and 92°C at 1 PV/day of 0.219 MgCl_{2} brine, the solid volume changes approximately by 0.01 and 0.005 cm^{3}/day, respectively. The porosity reduction rate can be proportional to porosity

where

This model takes the initial porosity (

The mathematical models aimed to match observed creep data have a long history, and several, more or less physically based models were reported. Generally, these models do not consider the underlying solid and pore volume contribution, but may still satisfactorily match the observed strains. Three models that have been used are:

Power law with cut-off:

De Waal [11]:

Griggs [12]:

The model parameters

### 5.1. Pore collapse and grain reorganization: the constant solid volume case

The movement of grains relative to each other at high-mean effective pressures causes pore volumes to collapse. It has been experimentally verified that for chalks, the rate of compaction may sometimes accelerate when the fluid composition change, a process termed water weakening. Water weakening has been used to understand reservoir processes [13, 14, 15] and to interpret core experiments as exemplified in Figure 5, where an additional strain of

To understand the immediate additional deformation (i.e., 1–2 days corresponding to a flow of 1-2 pore volumes as seen here), a grain-level approach is required. The grain-grain friction controls and cement bonds binding neighboring grains together control the relative movement of grains. Friction between grains is given by the frictional coefficient times the normal force,

### 5.2. Pressure solution and other grain-reorganization mechanisms

Pressure and temperature are the state variables that control the Gibbs chemical potential [21]. During diagenesis and burial, the chemical stability of mineral phases is altered as the temperature, hydrostatic and lithostatic pressure increases. Pressure solution of stressed grain contacts, and precipitation in unstressed parts of the rock framework, have been used as one of the primary rock-forming mechanisms during diagenesis. Pressure solution can occur in closed systems, in which the overall mass and density remain fixed. For high Biot coefficients, the local stress at particle contacts may become significant (see Figure 7 in [10]), and thus, a stress-dependent production of Ca-ions is observed where more Ca-production for high stress than low stress (see Figure 6 acquired from [10]).

It has been a long-standing discussion how to mathematically describe the relevant thermodynamic pressure for accurate determination of the chemical potential from the stress tensor. Several candidates coexist. The stress tensor in reservoir systems (and core scale experiments) depends on the weight of the overburden (lithostatic weight), side stress (tectonic forces), pore pressure, and the Biot coefficient. The simplest of determining the thermodynamic pressure is using the pore pressure. This way of thinking may seem reasonable at first glance since it is at the interface between the solid and the fluid where the chemical reactions occur. Simultaneously, at the rock-fluid interface, the stresses through the solid framework could also play a role in determining chemical solubility. In that case, the continuum mechanics provide a range of choices for calculating the thermodynamic pressure: (1) the average compressive stress (i.e., the first invariant of the solid framework stress tensor), (2) the principal stresses, thereby leading to different solubility in the different spatial directions. In sedimentary systems, this would often lead to enhanced solubility in the vertical direction as the first principal stress direction is vertical. This may explain the formation of the horizontal stylolites that are sometimes found in calcitic, carbonate, and limestone rocks [6]. (3) The relevant thermodynamic pressure could be related to the stress gradients that have been observed throughout porous materials, termed force chains. At grain-grain contacts, through which the externally imposed loads are being carried, the stresses can be significantly higher than the average. In these regions, the solid-solid stress is given by

Even though pressure solution is a process of chemical nature, it does not necessarily change the solid volume since the mass can be conserved (closed system, i.e., no larger scale mass flow) and the same mineral phase is precipitated as the one dissolved (i.e., same density). In that sense, pressure solution contributes to pore volume reduction rather than the solid volume in the strain partitioning presented here. Hence, pressure solution may fall under mechanical compaction even though the underlying mechanisms of pressure solution are chemically driven.

## 6. Time-dependent solid volume evolution mechanisms

In open, nonequilibrium systems with rock-fluid interactions, the solid volume is subject to change. It has been shown in a range of experiments how additional strain is accumulated during compaction at constant stress conditions when reactive brines are injected [4, 23]. The solid volume varies when solid mass (

The change in solid volume may be evaluated by

Here, the solid volume change is given by the difference between the ratio of the mass and density at a given time and the values before chemo-mechanical processes initiated. The evolution of the solid mass over time is given by the difference between the chemical mass flux in and out of the system, and density changes as new minerals precipitate.

### 6.1. Mass transfer in open systems

When fluids continuously flow and react with the rock, the mass (and hence the solid volume) changes. The chemical flux can be monitored by evaluating the effluent concentration through difference between the ion concentrations in and out of the volume element (Figure 6a, b). This volume element may, in some cases, be between an injector and a producer in an oil field, or a core scale experiment in the laboratory [3]. The concentration of ions can be measured using ion chromatography, and over a time interval

In Eq. (33), the factor

### 6.2. Method to quantify the solid volume evolution

In Eq. (32), the change in solid volume depends on the change in both mass and in density as the minerals dissolve and precipitate. The overall mineral density, as

As the concentrations of different minerals vary, the changes to

It is not always the case that a detailed kinetic chemical model exist tuned to take into account how different mineral mixtures react with fluids in each different case. If the overall density before (

Then, the density at any given time may be estimated using

## 7. Predicting dynamic porosity evolution: an illustrative example

An example of a dynamic porosity development analysis is presented here based on material published in 2015 [3]. The results of an experiment performed over 1090 days where a Liegè (Belgium) chalk sample was exposed to hydrostatic stress (11.1 MPa, approximately 5 MPa above yield) and continuous flow of 0.219 M MgCl_{2} (33 and 99 cm^{3}/day, pore pressure 0.7 MPa and 130°C). Basic sample measurements were performed of dry/saturated mass, pore volume, solid and bulk volume and hence porosity before and after test, hence the mineral density estimated and confirmed using He-pycnometry (Table 1). The bulk core volume was reduced more than 15% and mineral mass is reduced by more than 18% while the density is increased from 2.7 to 2.9 g/cm^{3}. The stresses are sufficient to induce pore collapse for these chalks, and the 0.219 M MgCl_{2}-brine (with equal ion strength as seawater) induced dissolution of the calcium carbonate and precipitation of denser Mg-bearing carbonates (e.g., magnesite and dolomite).

Before test | After test (1090 days) | Change | |
---|---|---|---|

Dry mass (on scale) | 125.57 g | 102.64 g | −22.93 g |

Wet weight (saturated) | 158.56 g | 126.34 g | −32.22 g |

Pore volume | 32.99 cm^{3} | 23.71 cm^{3} | −9.28 cm^{3} |

Solid volume | 46.85 cm^{3} | 35.53 cm^{3} | −11.32 cm^{3} |

Bulk volume | 79.84 cm^{3} | 59.23 cm^{3} | −20.61 cm^{3} |

Mineral density (saturation and pycnometer) | 2.68 and 2.70 g/cm^{3} | 2.89 and 2.90 g/cm^{3} | 0.21 and 0.20 g/cm^{3} |

Porosity (saturation and pycnometer) | 41.3 and 41.7% | 40.0 and 40.1% | −1.3 and −1.6% |

During the test, the axial deformation and the ion concentration of the effluent fluids were measured. In Figure 7a, the Mg and Ca ion concentrations are measured through time. These measurements can be used to find the production rate in g/day using Eq. (33), seen in Figure 7b. A trebling of the inlet flow rate leads to more than a doubling in the calcite dissolution. The mass evolution is used to estimate the dynamic change in density using Eq. (36), and is then combined to estimate the solid volume as function of time as seen by the dotted line in Figure 7c plotted together with the bulk volume, estimated from the axial strain. Thus, the pore volume can be estimated (dashed line in Figure 7c). As can be seen, the pore volume is reduced when bulk compaction dominate the overall process until 200 days. A typical observation from primary creep experiments is that the overall creep rate decreases with time. After 200 days, when the compaction rate has reduced, the flow rate was increased thereby increasing the rate at which dissolution/precipitation occurs, see Eq. (33) where the flow-rate dependency is explicitly shown. At this point of time, the overall porosity dynamics change. In the initial compaction-dominated regime, the overall porosity reduced to a value of as low as 33%, and afterwards it starts increasing (solid line in Figure 7d). At approximately 400 days, the flow rate is then reduced again and the rate of change in porosity is changing accordingly.

From 900 days and onwards, the Ca that was initially found within the core had been produced, after the solid volume was interpreted to be constant and the bulk compaction is facilitated by pore volume reduction, and hence the porosity is decreased to 40.1%.

In the experiment presented here, both pore and solid volume are subject to change. Since only the bulk volume or the solid volume could be determined from axial strain IC data, respectively, the pore volume was determined. As is exemplified in the presented experiment, the porosity evolution dynamics display a complex behavior because of the reduction in pore volume and solid volume. Their rate depends upon stress, the way in which deformation is accumulated and the rate of dissolution/precipitation.

## 8. Summary

Porosity is an important parameter for understanding the diagenetic processes and petrophysical reservoir systems. Its importance to the mechanical stiffness and strength of porous rocks, and to the resource potential, and rate of hydrocarbons produced from reservoirs is evident. The porosity is a dynamic parameter from the strain and chemical reactions from injection of fluids out of equilibrium with the host rock (e.g., seawater brines at elevated temperature in chalks) that induce additional deformation over time. The adsorption of surface-active ions leads to alterations in the forces binding grains together, leading to instantaneous additional deformation.

To understand quantitatively how porosity changes dynamically through time, there are series of processes that needs to be incorporated. This chapter presents some of the ways in which the bulk strain can be partitioned into elastic/plastic components, time-dependent, and time-independent components, and solid volume and pore volume processes. For chalks, the dynamic porosity evolution depends on the relative importance of the different processes at play, that again are functions of the stress, strain, temperature, flow rate, and fluid chemistry.

The methods presented here do not cover all possibilities for porosity evolution determination depending upon measurements that are available. When the bulk volume strain and chemical composition of the effluent fluids are known, the following porosity evolution model applies