In this section, the CGG method is tested on a velocity-stack inversion, which is useful not only for velocity analysis but also for various data processing. The applications of the velocity-stack inversion include separating multiples from the signal in the velocity domain Kostov and Nichols (1995); Lumley et al. (1995), multiple-attenuation techniques using parabolic Radon transforms Herrmann et al. (2000); Kabir and Marfurt (1999), missing-trace interpolation in the CMP domain Ji (1994), and so on. In these applications, the velocity-stack panels obtained by inversion are usually required to be as spiky and sparse as possible. Then the hyperbolic events represented by the isolated peaks in the velocity-stack panel are more easily distinguished from the rest of the noise.
The conventional velocity stack is performed by summing or estimating semblance Taner and Koehler (1969) along the various hyperbolas in a CMP gather, resulting in a velocity-stack panel. Ideally a hyperbola in a CMP gather should be mapped onto a point in a velocity-stack panel. Summation along a hyperbola, or hyperbolic Radon transform (HRT), does not give such resolution. To obtain a velocity-stack panel with better resolution, Thorson and Claerbout 1985 formulated it as an inverse problem in which the velocity domain is the unknown space. If we find an operator that transforms a point in a model space (velocity-stack panel), , into a hyperbola in data space (CMP gather), ,
(10) |
(11) |
The least-squares solution has some attributes that may be undesirable. If the model space is overdetermined and has bursty noise in data, the least-squares solutions usually will be spread over all the possible solutions. Other methods may be more useful if we desire a parsimonious representation. To obtain a more robust solution, Nichols 1994 used the IRLS method for L^{1}-norm minimization, and Guitton and Symes 2003 used the L-BFGS method for Huber-norm minimization. Another possibility is the CGG method proposed in the preceding section. In the next subsections the results of the CGG method for the velocity-stack inversion are compared with the results of conventional LS and L^{1}-norm IRLS.