Development of an Inversion Method for Low Velocity Medium


  • A. Afnimar Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology
  • Andri Dian Nugraha Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology
  • Ahmad Syahputra Faculty of Mining and Petroleum Engineering, Bandung Institute of Technology



Eikonal-solver method, Fermat principle, knot-point parameterization, low velocity structure, pseudo-bending method


The main problem with the inversion of a low velocity medium is the application of an appropriate ray tracing method after choosing a suitable model parameterization. Block parameterization is not suitable, because it is not capable of representing the velocity model well. A large amount of blocks with a small grid size are needed to express the model well, but in that case, a ray coverage problem will be encountered. A knot-point parameterization model is better suited than a block model, because it can express the velocity model well, while the number of variables is much smaller. Ray calculation using the pseudo-bending method is not appropriate for the velocity model because of an instability problem at high velocity gradients. The crucial problem of this method involves the initial ray-path that is optimized in order to obtain the "true" ray, but does not satisfy the Fermat principle. These problems can be solved by applying the eikonal-solver method, because this can handle high-velocity gradients and does not need an initial ray path. Using a suitable model parameterization and appropriate ray tracing method, the inversion can obtain good results that fit the desired output. Applying a block model and the pseudo-bending method will not produce the desired output.


Zelt, C.A. & R.B. Smith, Seismic Traveltime Inversion for 2-D Custal Velocity Structure, Geophys. J. Int, 108, pp. 16-34, 1992.

Afnimar, Koketsu, K. & Nakagawa, K., Joint Inversion of Refraction and Gravity Data for Three-Dimensional Topography of A Sediment Basement-Interface, Geophys. J. Int., 151, pp. 243-254, 2002.

Nugraha, A.D. & Mori, J., Three-Dimensional Velocity Structure in The Bungo Channel and Shikoku Area, Japan, and Its Relationship to Low-Frequency Earthquakes, Geophys. Res. Lett, 33, L24307. doi:10.1029/ 2006GL028479, 2006.

Zhao, D., Hasegawa, A. & Horiuchi, S., Tomographic Imaging of Pand S Wave Velocity Structure Beneath Northeastern Japan, J. Geophys. Res., 97, pp. 19909-19928, 1992.

Miyazawa, M. & Kato, M., On Interpolation Functions in Traveltime Tomography, Geophys. J. Int., 158, pp.169-178, 2004.

Um, J. & Thurber, C.H., A Fast Algorithm for Two-Point Seismic Ray Tracing, Bull. Seism. Soc. Am., 77, pp. 972-986, 1987.

Afnimar & Koketsu, K., Finite-Difference Traveltime Calculation for Head Waves Traveling Along an Irregular Interface, Geophys. J. Int, 143, pp. 729-734, 2000.

Vidale, J., Finite-Difference Calculation of Traveltimes, Bull. Seism. Soc. Am., 78, pp. 2062-2076, 1988.

Vidale, J., Finite-Difference Calculation of Traveltimes in Three Dimensions, Geophysics, 55, pp. 521-526, 1990.

Podvin, P. & Lecomte, I., Finite-Difference Computation of Traveltimes in Very Contrasted Velocity Models: A Massively Parallel Approach and Its Associated Tools, Geophys. J. Int, 105, pp. 271-284, 1991.

Paige, C.C. & Saunders, M.A., LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares, ACM Trans. Math. Soft., pp. 43-71, 195-209, 1982.

Grandis, H., Pengantar Pemodelan Inversi Geofisika, Himpunan Ahli Geofisika Indonesia (HAGI), 2009.