Fast 3D “wave-consistent” ray tomography based on a polynomial representation of a model

In the paper, we describe an original 3D travel-time tomography approach. It is based on the new realization of the bending method which to some extent takes into account the band-limited nature of real seismic signals propagation. As a result, two-point ray tracing provides more reliable ray trajectories and travel times in complex media. Another original feature of the proposed tomography is that the model is represented using the Chebyshev polynomials. Such parameterization allows analytical calculation of travel times and their derivatives with respect to model parameters and significantly reduces the number of parameters to be recovered during inversion compared to more common grid tomography. In certain situations, the proposed approach provides significant computational advantages. Numerical examples prove its efficiency.