Non-uniqueness with refraction inversion - a syncline model study

• Published in: Geophysical Prospecting Vol. 58; no. 2; pp. 203 - 218
• Main Author: Palmer, Derecke
• Format: Journal Article
• Language: English
• Published: Oxford, UK Blackwell Publishing Ltd 01.03.2010; Blackwell
• Subjects: Applied geophysics; Earth sciences; Earth, ocean, space; Exact sciences and technology; Internal geophysics; Narrows; New York; United States; waves; algorithms; inverse problem; seismic refraction; synclines; accuracy; velocity; North America; amplitude; parametrization; depth; tomography; two-dimensional models; thalwegs; anticlines; programs; travel time
• ISSN: 0016-8025; 1365-2478
• DOI: 10.1111/j.1365-2478.2009.00818.x

ABSTRACT Non‐uniqueness occurs with the 1D parametrization of refraction traveltime graphs in the vertical dimension and with the 2D lateral resolution of individual layers in the horizontal dimension. The most common source of non‐uniqueness is the inversion algorithm used to generate the starting model. This study applies 1D, 1.5D and 2D inversion algorithms to traveltime data for a syncline (2D) model, in order to generate starting models for wave path eikonal traveltime tomography. The 1D tau‐p algorithm produced a tomogram with an anticline rather than a syncline and an artefact with a high seismic velocity. The 2D generalized reciprocal method generated tomograms that accurately reproduced the syncline, together with narrow regions at the thalweg with seismic velocities that are less than and greater than the true seismic velocities as well as the true values. It is concluded that 2D inversion algorithms, which explicitly identify forward and reverse traveltime data, are required to generate useful starting models in the near‐surface where irregular refractors are common. The most likely tomogram can be selected as either the simplest model or with a priori information, such as head wave amplitudes. The determination of vertical velocity functions within individual layers is also subject to non‐uniqueness. Depths computed with vertical velocity gradients, which are the default with many tomography programs, are generally 50% greater than those computed with constant velocities for the same traveltime data. The average vertical velocity provides a more accurate measure of depth estimates, where it can be derived. Non‐uniqueness is a fundamental reality with the inversion of all near‐surface seismic refraction data. Unless specific measures are taken to explicitly address non‐uniqueness, then the production of a single refraction tomogram, which fits the traveltime data to sufficient accuracy, does not necessarily demonstrate that the result is either ‘correct’ or the most probable.