15.01.2015 Views

Microseismic Monitoring and Geomechanical Modelling of CO2 - bris

Microseismic Monitoring and Geomechanical Modelling of CO2 - bris

Microseismic Monitoring and Geomechanical Modelling of CO2 - bris

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!