Open AccessTwo‐Dimensional Seismic Refraction Tomography
Author(s) -
White D. J.
Publication year - 1989
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.1989.tb00498.x
Subject(s) - slowness , inversion (geology) , tomography , jacobian matrix and determinant , seismic tomography , geometry , ray tracing (physics) , refraction , mathematical analysis , inverse problem , perturbation (astronomy) , seismic wave , geology , geodesy , mathematics , optics , physics , geophysics , seismology , tectonics , quantum mechanics
Summary An iterative tomographic inversion scheme is presented for determination of 2‐D velocity structure from seismic refraction first‐arrival traveltimes. the method is suited to refraction profiles where source/receiver spacings are denser than for conventional profiles. the inversion method is based on an iterative solution of the linearized problem, and allows for determination of continuous velocity variations as welt as geometry of subhorizontal interfaces. In each iteration, two‐point ray tracing is performed using a shooting method to construct the linear system. the velocity field is defined using triangular cells within which the velocity gradient is constant, allowing analytic calculation of ray paths. Two different inversion techniques are considered, based on distinct linearized formulations of the forward problem. Inversion using a linearized traveltime‐velocity Jacobian produced better results than a slowness formulation more akin to common series expansion‐techniques. Resolution examples reveal horizontal smearing due to ray geometry, drop‐off in resolution with depth, as well as the effect of source‐receiver geometry and velocity structure on resolution. Inversion examples indicate that a global norm produces solutions closer to the true model than solutions calculated using a smallest perturbation approach, when a good starting model is available. Streak effects caused by inhomogeneous ray coverage and the removal of these effects are demonstrated.