Microseismic Monitoring and Geomechanical Modelling of CO2 - bris
Microseismic Monitoring and Geomechanical Modelling of CO2 - bris
Microseismic Monitoring and Geomechanical Modelling of CO2 - bris
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
CHAPTER 6.<br />
GENERATING ANISOTROPIC SEISMIC MODELS BASED ON GEOMECHANICAL SIMULATION<br />
This approach does have its strengths, in particular the reduction <strong>of</strong> a complex system to 3<br />
empirically determined constants, which facilitates the population <strong>of</strong> any models we may wish to<br />
build. For this reason the 3rd-order approach is becoming increasingly popular in the industrial<br />
sector (Herwanger, 2007). However, in order to fully define the 3rd order tensor for isotropic samples,<br />
velocity measurements are required at non-hydrostatic stresses (Prioul et al., 2004). As a result, there<br />
is little experimental data <strong>of</strong> the sort required to determine the 3rd-order terms, <strong>and</strong> as they are<br />
empirical constructs, they tell us little about the processes that lead to the development <strong>of</strong> non-linear<br />
elasticity. For this reason I prefer an approach that considers the micro-structure <strong>of</strong> the rocks in<br />
question.<br />
6.3 A micro-structural model for nonlinear elasticity<br />
A number <strong>of</strong> micro-structural models exist in the literature (e.g., Zatsepin <strong>and</strong> Crampin, 1997; Shapiro<br />
<strong>and</strong> Kaselow, 2005) where variations in micro-structural parameters are defined as a function <strong>of</strong> stress<br />
<strong>and</strong> related to the overall elastic properties <strong>of</strong> the rock via an effective medium model. In this section<br />
we consider the generalised effective medium approach <strong>of</strong> Schoenberg <strong>and</strong> Sayers (1995). In this<br />
approach the elasticity <strong>of</strong> a rock is evaluated in terms <strong>of</strong> the stiffness <strong>of</strong> its mineral components <strong>and</strong><br />
the presence <strong>of</strong> low volume displacement discontinuities, which serve to increase compliance. This<br />
model is highly generalised, where few assumptions need to be made about the discontinuities. By<br />
assuming that the discontinuities can be considered as rotationally invariant cracks, we can extend<br />
the model by using a number <strong>of</strong> methods available in the literature to describe such discontinuities<br />
(e.g., Sayers <strong>and</strong> Kachanov, 1995; Hudson et al., 1996; Hall, 2000) <strong>and</strong> how they might vary with<br />
pressure (e.g., Tod, 2002). Following the approach <strong>of</strong> Tod (2002), I develop a simple model to describe<br />
the change in elasticity <strong>of</strong> a rock as a function <strong>of</strong> the stress applied to it. This approach is capable<br />
<strong>of</strong> considering anisotropy that develops due to both intrinsic rock properties <strong>and</strong> as a result <strong>of</strong> nonhydrostatic<br />
stresses. It is also capable <strong>of</strong> providing a framework within which we might consider<br />
damage due to coring or thermal effects.<br />
6.3.1 Theoretical background<br />
Schoenberg <strong>and</strong> Sayers (1995) introduce an effective medium approach to describe the compliance, S,<br />
<strong>of</strong> a damaged rock. This approach is defined in terms <strong>of</strong> a matrix material <strong>and</strong> a r<strong>and</strong>om distribution<br />
<strong>of</strong> low volume, poorly bonded discontinuities. When a stress is applied across such a discontinuity,<br />
there will be a difference in displacement between the faces - a displacement discontinuity, [u i ] - that<br />
is proportional to the traction, t i = σ ij n j , on the discontinuity surface s. Hence, for a discontinuity<br />
with normal n, the total displacement discontinuity [u i ] is given by<br />
∫<br />
[u i ]ds ∝ σ ij n j . (6.4)<br />
s<br />
108