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.
C in In Chapter 6 (<strong>and</strong> Angus et al., 2009) we have found that the aspect ratio rarely varies between<br />
CHAPTER 7.<br />
FORWARD MODELLING OF SEISMIC PROPERTIES<br />
7.2.2 Initialisation<br />
In order to ensure that the relative velocities computed using our rock physics model match the relative<br />
stiffnesses <strong>of</strong> the materials used in the geomechanical modelling, I use an initialisation procedure such<br />
that the initial stiffness is set to be equal to that used to compute the geomechanical deformation,<br />
C mech =C in . These values can be increased by a specified percentage to represent the fact that<br />
dynamic stiffness (the stiffness used to calculate seismic velocities, which are low strain <strong>and</strong> high<br />
strain rate) is generally observed empirically to be larger than static stiffness (used to calculate<br />
geomechanical deformation, which has high strain at a low strain rate). By doing so I preserve<br />
the relative stiffness differences between reservoir <strong>and</strong> non-pay units, although absolute values are<br />
increased for the dynamic stiffness. I will refer to the stiffness used to compute the initial seismic<br />
velocities as C in , while recognising that it may be equal to, or tied to, the geomechanical stiffness<br />
C mech . Alternatively, where the seismic velocities <strong>of</strong> layers are known, these can be used to define<br />
sedimentary rocks, so I specify that the initial average aspect ratio a 0 is fixed. This means that there<br />
are two parameters that can be varied to ensure that the dynamic stiffness tensor, C, at the initial<br />
ρ d<br />
ρ fl<br />
Φ<br />
K fl<br />
P fl<br />
C mech<br />
C in<br />
σ ij<br />
β w<br />
ξi<br />
0<br />
a 0<br />
ξ f<br />
a f<br />
θ f<br />
ϕ f<br />
ω<br />
M g<br />
Results from MORE-ELFEN<br />
Density <strong>of</strong> the dry rock<br />
Effective density <strong>of</strong> the multiphase pore fluid<br />
Porosity<br />
Effective bulk modulus <strong>of</strong> the multiphase pore fluid<br />
Pore fluid pressure<br />
Elastic stiffness used to compute geomechanical deformation<br />
Elastic stiffness used to compute the initial seismic velocities<br />
Stress tensor<br />
User-defined inputs<br />
Biot-Willis parameter<br />
Initial crack density tensor at zero stress<br />
Initial crack aspect ratio<br />
Number density <strong>of</strong> any user-defined fracture sets<br />
Aspect ratio <strong>of</strong> any user-defined fracture sets<br />
Azimuth <strong>of</strong> normals to user-defined fracture sets<br />
Inclination <strong>of</strong> normals to user-defined fracture sets<br />
Dominant frequency <strong>of</strong> incident seismic energy, used to compute squirtflow<br />
effects<br />
Characteristic grain size, used to compute squirt-flow effects<br />
Table 7.1: List <strong>of</strong> SeisModel c⃝ input parameters<br />
136