Next: SIMULTANEOUS ESTIMATION OF WAVELET
Up: DECONVOLUTION WITH A KNOWN
Previous: Building a more complex
Although I have not yet explored the practical implications of this,
the operator
, previously defined, can be described
as a particular case of a more general operator
:
| ![\begin{displaymath}
{\bf r} = {\cal D}_L^{\alpha}\{{\bf c}\} \;\;\; \longrightar...
...r
{\displaystyle \sum_{j=i-n}^{i+n}}
\mid c_j^{\alpha} \mid },\end{displaymath}](img29.gif) |
(10) |
![\begin{displaymath}
\mbox{\hspace{-4.1cm} where \hspace{1.0cm}}
n = {L + 1 \ove...
...d \hspace{0.6cm}
$\alpha$ \hspace{0.6cm} is any odd integer}. \end{displaymath}](img30.gif)
The operator
can be defined as
| ![\begin{displaymath}
\delta^L = \lim_{\alpha \rightarrow \infty} \: {\cal D}_L^{\alpha}.\end{displaymath}](img31.gif) |
(11) |
Figure 7 shows that using this operator
with finite values of
is equivalent to the introduction of
the damping factor in equation (7). Also we can see that
the results obtained with small values of
are
similar to those obtained with the predictive deconvolution.
An Lp norm deconvolution can be formulated as an iterative,
reweighted least-squares inversion with weighting factors given
by
, where
are the residuals of the previous iteration (Darche, 1989).
Since the value of the exponent
will directly affect the relative weighting of the residuals,
we can consider the ICP decon as somehow equivalent to
an
norm deconvolution.
knowt3
Figure 7 (a) Original reflectivity series. The other traces
correspond to the results of the ICP algorithm using the
generalized operator of equation (10) with
different values of
: (b)
,(c)
, (d)
,and (e)
.
Next: SIMULTANEOUS ESTIMATION OF WAVELET
Up: DECONVOLUTION WITH A KNOWN
Previous: Building a more complex
Stanford Exploration Project
1/13/1998