In this section I will propose a simple parameterization of angle-dependent velocity appropriate for relatively weak anisotropy with a vertical axis of symmetry. Anisotropy will be assumed to derive from layered media that may be isotropic on a fine scale but which will appear to be anisotropic on the scale of larger seismic wavelengths. The anisotropic parameters will be described as a sum of smooth basis functions, with no more spatial variations than necessary to explain the data. The form will allow easy linearization of slowness with respect to the model coefficients.

The parameterization of velocities should have enough degrees of freedom to describe all plausible models; the particular numerical formulation of this parameterization is less important. Our methods should be able to request the velocity at any point and in any direction from an essentially continuous model.

The model should also be differentiable, to allow the efficient optimization of raypaths. If discontinuities are included in the model, then it should be possible to smooth these values numerically to calculate gradients on different scales. I will also require that the resolution of the model should be adjustable at any time, even during optimization of the model. Adjustable resolution will allow the optimization to converge first on the smoothest, most reliable components of the velocity model. As the accuracy of estimated raypaths improves, then more detail will be allowed in velocity models.

No convenient explicit equation exists to describe group velocity as a function of angle. Instead, I will use an approximation with enough degrees of freedom to explain the data well and still adequately span the same range of functions allowed by the exact theory. Tomography will have a limited ability to estimate arbitrary changes in velocity with angle. An approximate equation will have a form that is easy to optimize and yet describe the most important variations in velocity with angle. The errors introduced by analytic approximation are intended to be much smaller than errors introduced by inaccurate traveltimes. We might fit data first with more approximate curves and introduce refinements only when sensitivity improves.

The following equation parameterizes group velocities as a function of group angle , measured from the vertical axis of symmetry:

(1) |

Slowness is integrated as a function
of distance to give traveltimes.
The horizontal velocity *V*_{x} is well defined by a physical
experiment and is measured accurately from surface
or crosswell experiments; the velocity *V*_{z}
along the vertical axis of symmetry is much less well determined.
The first parameter in the square root
is determined second best by surface experiments
Alkhalifah and Tsvankin (1994); Tsvankin and Thomsen (1994)
because it specifies the difference between a normal moveout
(NMO) velocity and *V*_{x} according to
equation (21). Finally, the parameter expresses the third and least well determined part of anisotropy,
giving the fractional change of vertical velocity
according to (19).

The magnitude *V*_{x} will be allowed to change most arbitrarily,
in many dimensions.
Because this anisotropy is assumed to be a layered phenomenon,
will only be allowed to change vertically (or perpendicular
to layering). Because
is so poorly determined by surface data, I will
assume that it has a single global value which can
be chosen to fit well ties when available. Or one may attempt
to predict from the other two parameters by
observing the correlation in values that are produced by
equivalent media calculations from well logs Backus (1962).

We can express the continuous velocity functions as a scaled sum of smooth basis functions. For example, if is an arbitrary Cartesian coordinate, then we can express the slowness as a linear function of discrete parameters , using smooth basis functions :

(2) |

(3) |

Let us designate this discrete set of velocity parameters
as a single vector .
The different elements
of are understood to have different scales,
and will assume different variances during optimization.
The approximate magnitude of *V*_{x} in (2),
and thereby of , is easily
anticipated from a quick glance at traveltimes over certain short
distances.
Theory Backus (1962) can easily
anticipate reasonable magnitudes for and for equivalent layered media.

Now we can write the continuous group velocity explicitly as a function of these discrete parameters

(4) |

(5) |

11/12/1997