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.
7.3. RESULTS FROM SIMPLE GEOMECHANICAL MODELS<br />
7.2.5 Fluid substitution<br />
SeisModel c⃝ also computes the changes in seismic velocity induced by changes in fluid saturation.<br />
SeisModel c⃝<br />
provides two options for doing this - using the low frequency anisotropic Gassmann<br />
equation (Brown <strong>and</strong> Korringa, 1975) given in equation 3.9, or the frequency dependent Chapman<br />
(2003) model which includes squirt flow <strong>and</strong> global flow effects. Because the focus <strong>of</strong> my work is on<br />
the effects <strong>of</strong> stress <strong>and</strong> not fluid substitution, I will use the low frequency Gassmann approach. The<br />
fluid bulk modulus <strong>and</strong> porosity are provided by the MORE-ELFEN results, but the mineral stiffness<br />
is also required to compute the saturated stiffness. In order to estimate this, SeisModel extrapolates<br />
from the background stiffness tensor computed in the initialisation step to the limit <strong>of</strong> zero porosity,<br />
using<br />
C m =<br />
Cr<br />
1 − Φ . (7.6)<br />
The fluid saturated density is computed using the porosity <strong>and</strong> the densities <strong>of</strong> the fluid <strong>and</strong> rock, all<br />
<strong>of</strong> which are provided by MORE-ELFEN, using<br />
ρ = ρ d + Φρ fl . (7.7)<br />
SeisModel c⃝ also provides the option to skip any <strong>of</strong> these steps if we wish to ignore the contribution<br />
from fractures, stress dependence or fluid substitution.<br />
7.2.6 SeisModel c⃝ Output<br />
SeisModel c⃝ generates an output file at each timestep. The header <strong>of</strong> this file contains information<br />
about the grid geometry making up the seismic model, <strong>and</strong> the elevations <strong>of</strong> any surfaces (such as<br />
the top <strong>of</strong> the reservoir) which can be used by ray-tracing algorithms. For each node, the E, N <strong>and</strong> Z<br />
coordinates <strong>of</strong> the node are provided, along with the final stiffness tensor computed by the SeisModel c⃝<br />
workflow in 6×6 Voigt notation. These files can then be used to compute seismic observables such as<br />
travel-time to the top <strong>of</strong> the reservoir, reflection coefficients, <strong>and</strong> shear wave splitting at each timestep.<br />
I will now show the seismic results for the geomechanical models constructed in Chapter 5.<br />
7.3 Results from simple geomechanical models<br />
7.3.1 Overburden travel time-shifts<br />
The normal incidence travel time for a reflection from the top <strong>of</strong> the reservoir can quickly <strong>and</strong> easily<br />
be computed by ray-tracing through the elastic model developed, <strong>and</strong> it is also easily observed on 4-D<br />
seismic surveys. This change to this observable has been the most commonly deployed technique used<br />
to link geomechanical deformation with geophysical observations (e.g., Hatchell <strong>and</strong> Bourne, 2005). I<br />
use the Christ<strong>of</strong>fel equation (Chapter 3) to compute the changes in velocity <strong>of</strong> a vertically propagating<br />
P-wave through the centre <strong>of</strong> the 9 simple reservoir models developed in Chapter 5. The user-defined<br />
parameters that I use are listed in Table 7.2. Having computed the velocities, I compute the change<br />
139