Parallel numerical simulation of 3D acoustic logging

Finite-difference method for simulation of sonic waves propagating in a vicinity of a borehole filled with fluid and surrounded by a 3D heterogeneous elastic medium is developed. The method utilizes an explicit second order of approximation FD scheme on staggered grids that approximates the elastodynamic system of equations in cylindrical coordinates. A computational domain is surrounded by special Perfectly Matched Layer for cylindrical coordinate system designed to attenuate waves reflected from outer boundaries. Parallelization is based on domain decomposition approach and implemented with the help of MPI library. Results of numerical experiments are presented.