In order to perform large scale numerical simulation of wave propagation in 3D heterogeneous multiscale viscoelastic media, Finite Difference technique and its parallel implementation based on domain decomposition is used. A couple of typical statements of borehole geophysics are dealt with—sonic log and cross well measurements. Both of them are essentially multiscales, which claims to take into account heterogeneities of very different sizes in order to provide reliable results of simulations. Locally refined spatial grids help us to avoid the use of redundantly tiny grid cells in a target area, but cause some troubles with uniform load of Processor Units involved in computations. We present results of scalability tests together with results of numerical simulations for both statements performed for some realistic models.