Prestack first-break traveltime tomography using the double-square-root eikonal equation

Traveltime tomography with shot-based eikonal equation fixes shot positions then relies on inversion to resolve any contradicting information between independent shots and achieve a possible cost-function minimum. On the other hand, the double-square-root (DSR) eikonal equation that describes the whole survey, while providing the same first-arrival travel-times, allows not only the receivers but also the shots to change position and therefore leads to faster convergence in tomo-graphic inversion. The DSR eikonal equation can be solved by a version of the fast-marching method (FMM) with special treatment for its singularity at horizontally traveling waves. For inversion, we use an upwind finite-difference scheme and the adjoint-state method to avoid explicit calculation of Fréchet derivatives. The proposed method generalizes to the 3D case straightforwardly.