5.4 Methods of coupling flow and geomechanics

In conventional fluid flow formulations, the pore volume variation only depends on the pore volume compressibility coefficient. The rock compressibility is assumed to be constant and the reservoir permeability is unaffected by pore pressure changes. However, the injection of CO2, especially into highly compacted, faulted and fractured formations, causes a strain localisation on fracture and fault planes and results in a change in permeability or transmissibility. To account for geomechanical effects due to stress changes in and around the injection formation, the fluid flow problem must be solved with a geomechanical model that can predict the evolution of stress dependent parameters, such as porosity, rock compressibility, and permeability. The coupling can be carried out by integrating the mechanical concepts in reservoir simulation. The geomechanical equilibrium equation and the fluid mass balance equation should be solved iteratively. In the case of highly compacted, faulted and fractured reservoirs, the coupling may also lead to a modification of the transmissibility matrix due to fracture and fault permeability enhancement resulting from rock deformation (Longuemare et al., 2002).

The fully coupled and partially coupled approaches can be used to solve the stress dependent CO2 geological storage problem.

The fully coupled approach simultaneously solves the whole set of equations in one simulator. The fully coupled method offers internal consistency for the simultaneous resolution of both flow and stress equations, but the hydraulic or geomechanical mechanisms are often simplified by comparison with conventional uncoupled geomechanical and reservoir approaches. TOUGH-FRAC, a simulator for non-isothermal multiphase flow in porous media with geomechanical coupling, is an example of such a code which models plume dispersion and impact of stresses due to CO2 interactions.