We propose a set of local operators to deal with internal discontinuities that do not coincide with collocated Cartesian grids in 2D transversely isotropic media. We use globally optimally accurate operators in order to obtain high accuracy. We derive modified operators by extrapolating wavefields from nearby grid points to the discontinuous point, introducing boundary conditions; we then distribute those conditions to the nearby grid points. Numerical examples suggest that the operators improve the coherency of the wavefront. We would like to optimise the local operators to control all the error propagation during the modelling.
|Original language||English (US)|
|Title of host publication||78th EAGE Conference and Exhibition 2016 - Workshops|
|Publisher||EAGE Publications BV|
|State||Published - Mar 13 2017|
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledgements: We thank GPX consortium, composed of French National Research Agency (ANR) and industrial sponsors, CGG, Schlumberger and Total, who supported O.O. during his master research. C.C. was supported by Laboratoire de Sismologie at Institut de Physique du Globe de Paris.