Full text loading...
-
Julia Application in Calculating Reflectivity from High Order Anisotropic Interface
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, 78th EAGE Conference and Exhibition 2016, May 2016, Volume 2016, p.1 - 5
Abstract
Elastic anisotropy in the Earth causes many artifacts in an isotropic seismic processing and interpretation. To date, these issues have been ameliorated by using approximations to the full solutions for wave propagation and reflectivity for special material symmetries. This algorithm has been written to extend these capabilities to the general case of reflectivity from the interface between two anisotropic slabs of arbitrary symmetry and orientation. To achieve this, the algorithm solves for polarization, amplitude and slowness of all the wave modes generated by a plane wave incident to the interface. In the first step, the plane-wave velocities and polarizations of all three orthogonal wave modes are calculated for a given incidence angle. Second, the algorithm determines the reflection and transmission angles of all of the possible scattered modes followed by their respective velocities and polarization vectors. With this information, the algorithm solves system of equations incorporating the imposed boundary conditions to arrive at the scattered wave amplitudes. We have tested this algorithm in some of the most computationally difficult models to ensure there is no energy leakage in the system of calculations. This algorithm could be used in seismic processing, migration, and inversion of seismic data in anisotropic media.