Microseismic Monitoring and Geomechanical Modelling of CO2 - bris
Microseismic Monitoring and Geomechanical Modelling of CO2 - bris
Microseismic Monitoring and Geomechanical Modelling of CO2 - bris
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 6.<br />
GENERATING ANISOTROPIC SEISMIC MODELS BASED ON GEOMECHANICAL SIMULATION<br />
Define S r (from EBSD or high pressure behaviour)<br />
⇓<br />
For each velocity measurement at specified pressure:<br />
Define C obs from velocities<br />
⇓<br />
Loop over B N /B T values<br />
Initialise crack density tensor<br />
α kk = 0<br />
⇓<br />
Loop over number <strong>of</strong> iterations<br />
Define C model = f(S r , α kk , B N /B T )<br />
(equations 6.22 – 6.24)<br />
⇓<br />
Define model misfit:<br />
δb = C obs − C model<br />
⇓<br />
Define Jacobean: J = δC<br />
δα<br />
⇓<br />
Invert δb = Jδm<br />
⇓<br />
Update model: α kk = α kk + δm<br />
End iteration loop<br />
⇓<br />
Back calculate velocities using final value <strong>of</strong> C model<br />
⇓<br />
Compute misfit between modelled <strong>and</strong> observed velocities<br />
δV = (V obs −V model ) 2<br />
V obs<br />
End B N /B T grid-search<br />
⇓<br />
Select the value <strong>of</strong> B N /B T that minimises δV<br />
⇓<br />
Select the values <strong>of</strong> α kk computed with this B N /B T value<br />
⇓<br />
Move on to next pressure measurement<br />
Figure 6.3: Workflow for inverting velocity measurements for B N /B T <strong>and</strong> the crack density tensor.<br />
Modified from Hall et al. (2008).<br />
116