compCmp
Figure 6 Portions of CMP gathers after radon demultiple. Left shows decimated CMP gather. Right side shows reinterpolated CMP gather, with interpolated traces thrown out after radon demultiple. Note that in the left panel aliased multiples at far offset are not well suppressed (A), and that some interesting energy in the mostly primary region around 3 seconds is attenuated (B).

 

compStack
Figure 7 Subsalt portion of stacked section. The left panel shows decimated data after demultiple and stack, the right panel shows same data after interpolation and stack. Multiples (A) are arguably better attenuated in the interpolated data. The decimated data shows attenuated primary events (B), and layer-like artifacts (C).

 

 

Short Note
Decon and interpolation with nonstationary filters

Sean Crawley, Robert Clapp, and Jon Claerbout

sean@sep.stanford.edu

ABSTRACT

Building on the notions of time-variable filtering and the helix coordinate system, we develop software for filters that are smoothly variable in multiple dimensions. Multiscale prediction-error filters (PEFs) can estimate dips from recorded data and use the dip information to fill in unrecorded shot or receiver gathers. The data are typically divided into patches with approximately constant dips. Instead, we estimate a set of smoothly varying filters, up to one PEF for every data sample. They are more memory-intensive to estimate, but the smoothly varying filters do give more accurate interpolation results than discrete patches. Finally, we offer an improved method of controlling the smoothness of the filters. We design filters like directional derivatives that we call ``flag filters''. They destroy dips in easily adjusted directions. We use them in residual space to encourage dips in the specified directions. We develop the notion of ``radial-flag filters'', i.e., flag filters oriented in the radial direction (lines of constant x/t in (t,x) space). Break a common-midpoint gather into pie shaped regions bounded by various values of x/t. Such a pie-shaped region tends to have constant dip spectrum throughout the region so it is a natural region for smoothing estimates of dip spectra or of gathering statistics (via 2-D PEFs). In this paper we use smoothly variable PEFs to interpolate missing traces, though obviously they have many other uses.