HUBER VELOCITY SPACE MULTIPLE ELIMINATION

Transforming seismic data into another domain and then muting is a common multiple suppression technique. The common method is the parabolic radon transform Foster and Mosher (1992). This method has the advantage of having a analytic inverse (and is therefore faster), but involves approximating moveouts by parabolas. Lumley et al. (1994) used the more expensive hyperbolic transform and went a step further by forming it as an inversion problem,  

$$\bf d\approx {\bf H} {\bf V} \bf m ,$$ (5)

where:

$\bf d$

is the CMP gather,

$\bf H$

is a half derivative operator Prucha (1999),

$\bf V$

is a velocity transform operator.

They then used an L₂ conjugate gradient algorithm to estimate a ($\tau,v$) model and then muted out the multiples. Unfortunately both methods suffer from multiple and primary energy overlapping.