The efficient computation of finite-frequency traveltime and amplitude sensitivity kernels for velocity and attenuation perturbations in global seismic tomography poses problems both of numerical precision and of validity of the paraxial approximation used. We investigate these aspects, using a local model parameterization in the form of a tetrahedral grid with linear interpolation in between g...

We adopt a recent work in Chung, Qian, Uhlmann and Zhao (Inverse Problems, 23(2007) 309-329) to develop a phase space method for reconstructing pressure wave speed and shear wave speed of an elastic medium from travel time measurements. The method is based on the so-called Stefanov-Uhlmann identity which links two Riemannian metrics with their travel time information. We design a numerical algo...

A higher-order parabolic equation is used to compute the traveltime (phase) and the amplitude in constant density acoustic media. This approach is in the frequency domain, thereby avoiding the high frequency approximation inherent in the Eikonal equation. Intrinsic attenuation can be naturally incorporated into the calculation. The error at large angles of propagation caused by the expansion of...

seismic tomography is an imaging technique which creates maps of subsurface elastic properties such as p/s wave velocity, density and attenuation, based on observed seismograms and use of sophisticated inversion algorithms. amongst different acquisition geometries, seismic cross-hole tomography has a special position in geophysical surveys with many applications in hydrocarbons, coal and other ...

Pseudo-depth monotonicity condition is an important assumption of the inverse scattering internal multiple attenuation algorithm. Analysis reveals that this condition is equivalent to a vertical-time monotonicity condition which is different than the total traveltime monotonicity suggested in recent literature/discussions. For certain complex media, the monotonicity condition can be too restric...

The Stolt migration extension to a variable velocity case describes the velocity heterogeneity with a constant parameter, which is related to the stretch transformation of the time axis. We exploit a connection between modified dispersion relations and nonhyperbolic traveltime approximations to derive an explicit expression for the stretch parameter. This analytical expression allows one two ac...

The point source traveltime eld has an upwind singularity at the source point. Consequently, all formally high-order nite-diierence eikonal solvers exhibit rst-order convergence and relatively large errors. Adaptive upwind nite-diierence methods based on high-order Weighted Essentially NonOscillatory (WENO) Runge-Kutta diierence schemes for the paraxial eikonal equation overcome this diiculty. ...

We present a method of ray tracing that is based on a system of differential equations equivalent to the eikonal equation, but formulated in the ray coordinate system. We use a first-order discretization scheme that is interpreted very simply in terms of the Huygens’ principle. The method has proved to be a robust alternative to conventional ray tracing, while being faster and having a better a...