Next: Synthetic data tests
Up: Guitton & Symes: The
Previous: Introduction
The velocity domain representation of seismic data is an alternative to the standard
CMP presentation. Transformation of CMP data into the velocity domain (producing a velocity
model or panel of the data) exhibits clearly the moveout inherent in the data and therefore,
forms a convenient basis for velocity analysis as a linear inverse problem.
The velocity transform A from the model space (velocity domain)
into the data space (CMP gathers) stretches the velocities back
in the offset plane (superposition of hyperbolas) whereas the adjoint operation (A')
squeezes the data (summation over hyperbolas):
![\begin{displaymath}
\bold{A}=\bold{HS},
\end{displaymath}](img7.gif)
with
![\begin{displaymath}
\bold{S}m(t,x) = \displaystyle \left.
\sum_{s}\frac{t}{\tau}w(s,x,\tau)m(\tau,s)
\right\vert _{\tau=\sqrt{t^2-s^2x^2}}\ \end{displaymath}](img8.gif)
![\begin{displaymath}
\bold{A'}=\bold{S'H'},
\end{displaymath}](img9.gif)
with
![\begin{displaymath}
\bold{S'}d(\tau,s) = \displaystyle \left.
\sum_{x}w(s,x,\tau)d(t,x)
\right\vert _{t=\sqrt{\tau^2+s^2x^2}}\ \end{displaymath}](img10.gif)
where
is a weighting function,
is a filter that we define later.
is related to the velocity stack as defined by Taner and Koehler (1969).
The problem is: given a CMP gather can we find a velocity panel which synthesizes it ![$\it{via}$](img14.gif)
? In equations, given data
, we want to solve for model
:
![\begin{displaymath}
\bold{Am} \approx \bold{d}.
\end{displaymath}](img17.gif)
A simple way to solve this problem is to find a model
that minimizes the mean
square misfit
![\begin{displaymath}
(\bold{Am}-\bold{d})'(\bold{Am}-\bold{d}).
\end{displaymath}](img18.gif)
This optimization problem is equivalent to the linear system
(``normal equations'')
![\begin{displaymath}
\bold{A'Am} = \bold{A'd}.
\end{displaymath}](img19.gif)
This system is easy to solve if
, i.e if
is close
to unitary: then
.In general,
is far from an unitary operator for many reasons.
However, the choice of a weighting function compensates to some extent for
geometrical spreading and other effects Claerbout and Black (1997):
![\begin{displaymath}
w(s,x,\tau)=\frac{1}{\sqrt{(\tau^2+s^2x^2)^{1/2}}}\frac{\tau}{\sqrt{\tau^2+s^2x^2}}\sqrt{xs}.
\end{displaymath}](img22.gif)
The summation in the velocity space boosts low frequencies. Claerbout and Black (1997) suggest that a good choice
of filter
is a half derivative operator (
).
These choices for
and
bring
closer to being
an unitary operator.
Since the data is noisy, the modeling operator is not unitary and the numbers of equations
and unknowns may be large, an iterative data-fitting approach seems reasonable:
![\begin{displaymath}
\rm{min}_\bold{m}(\it{E}(\bold{Am - d}))
\end{displaymath}](img24.gif)
where
is the model,
, the data we want to fit,
the modeling operator,
and E a misfit measurement function we have to choose.
We have already presented one possibility, namely that E is the l2 norm (least squares inversion).
A convenient iterative method in this case is to solve the normal equation using conjugate gradient iteration.
We refer to this approach as ``CG'' or ``l2''. An alternative approach is to take for E the Huber function
introduced in the first section. With this Huber misfit measure, the velocity transform inverse problem is no
longer equivalent to a linear system. We choose to solve it using a general-purpose nonlinear
optimizer, as mentioned above, rather than one of several special-purpose methods invented for this
type of problem. We refer to this approach as''Huber'' or the ``Huber solver''.
The next two parts of this paper compare the performance of the CG algorithm to Huber in
the velocity analysis application.
Next: Synthetic data tests
Up: Guitton & Symes: The
Previous: Introduction
Stanford Exploration Project
4/20/1999