Abstract: A seismic velocity analysis method includes tying velocity parameter values such as residual velocity values to geological horizons (reflectors) within a seismic exploration volume. Common image gathers (CIGs) such a common reflection point (CRP) gathers or angle-domain common image gathers (ACIGs) are generated for a set of CIG grid points. Computed best-fit residual velocity values are then snapped to a neighboring horizon or vertically interpolated to the horizon, to generate residual velocity values along the horizon. The residual velocity values for points along the horizon are then selectively employed in updating the velocity model for the volume of interest.

Abstract: Migration velocity analysis is performed using Angle-Domain Common Image Gathers (ACIGs). When the correct velocity model is employed for migration, all ACIG events corresponding to a subsurface location are aligned along a horizontal line. Residual moveout can be performed on each ACIG with a suite of trial residual velocity values, according to an angle-domain residual moveout equation. A best-fit residual velocity value that leads to horizontally-aligned events upon moveout can be selected by generating a distribution of semblance (amplitude summed over a given depth) over residual velocity. Best-fit residual velocity values corresponding to selected subsurface points can be employed to update the initial velocity model using a vertical update, normal ray update, or tomographic update method.

Abstract: Seismic traveltimes are computed using a conditional high-order method: the traveltime computation operator is second- or higher-order if enough suitable upwind traveltimes are available, and first-order otherwise. Typically, a first-order operator is employed only around singularities, e.g. corners and cusps; a high-order operator is employed for the vast majority of grid points in the target volume. Selectively switching to first-order makes the method relatively simple and computationally efficient, as compared to a pure high-order method. At the same time, most traveltimes are computed to high-order accuracy. In the preferred embodiment, the traveltime front is selectively advanced at its minimum traveltime grid point, using a finite-difference approximation to the eikonal equation. A narrow band propagation zone is used to advance the finite-difference stencil. Tentative traveltimes for the narrow band adjacent to the traveltime front are computed using the eikonal equation and arranged on a heap.