Conjugate guided gradient (CGG) method for robust inversion and its application to velocity-stack inversion |

The CGG method described in the previous section can be used to any inversion problem
whose required properties are the robustness to spiky noise and the parsimony of the model.
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 applications.
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) and Hampson (1986)
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) ,

(11) |

where denotes the common-midpoint (CMP) gather and denotes velocity-stack panel. The slowness-time pair are the coordinate axes of the velocity-stack and the offset-time pair are the coordinate axes of the CMP gather. A straightforward definition for the forward operator is the adjoint of the operator defined by Equation (12). Through the suitable definition of the inner product, turns out to be simply the process of reverse NMO and stacking (Thorson and Claerbout, 1985):

(13) |

(14) |

The least-squares solution has some attributes that may be undesirable. If the model space is overdetermined and has busty 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 of the solution. To obtain a more robust solution, Nichols (1994) and Trad et al. (2003) used the IRLS method for -norm minimization, and Guiiton and Symes (2003) used a quasi-Newton method called limited-memory BFGS Goldfarb (1970); Nocedal (1980); Broyden (1969); Fletcher (1970); Shanno (1970) 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 method and -norm IRLS method.

Conjugate guided gradient (CGG) method for robust inversion and its application to velocity-stack inversion |

2011-06-26