Modeling data error during deconvolution |

The decon filter , parameterized by , we take as noncausal. The constraint is no longer a spike at zero lag, but a filter whose log spectrum vanishes at zero lag, , so we are now constraining the mean of the log spectrum. This is a fundamental change which we confess to being somewhat mysterious.

The single regression for
including noise
now becomes two.

The notation means the data fitting is done under a hyperbolic penalty function. The regularization need not be . To save clutter I leave it as until the last step when I remind how it can easily be made hyperbolic.

Under the constraint of a causal filter with
,
traditional auto regression for
with its regularization looks more like

Comparing equations 1-2 with 3-4 you see we are not simply rehashing traditional methodology but seem to be off in a wholly new direction! We are here because solved our non-minimum phase problem, and seeing sea swell in our estimated shot wavelets told us we need to replace by .

Antoine noticed the quasi-Newton method of data fitting requires gradients but not knowledge of how to update residuals so the only thing we really need to think about is getting the gradient. The gradient wrt is the same as before (Claerbout et al. (2011)) except that replaces . The gradient wrt is the new element here.

Let
,
, and
be time functions (data, noise, and filter).
Let
be the residual.
Let
.
Expressing our two regressions in the time domain we minimize

(5) |

A scaling factor is required between the terms. We expect to learn it by experimentation.

Now we go after the gradient, the derivative of the penalty function wrt each component of noise
.
Let the derivative of the penalty function
wrt its argument
be
called the softclip and be denoted
.
Let
denote the FT of
.
Let
be the time reverse of
while in Fourier space
is the conjugate of
.

(6) | |||

(7) | |||

(8) | |||

(9) | |||

(10) | |||

(11) | |||

(12) |

For simplicity I set out with a quadratic penalty function on the noise, but it is easy to make it hyperbolic. Simply use softclip on . Change to .

Now having the gradient we should be ready to code.

Modeling data error during deconvolution |

2012-10-29