1) Constant Velocity Scans

2) Radon plus KL to eliminate Radon artifacts

3) Spiking decon followed by Maximum Kurtosis Phase Correction

1) Constant Velocity Scans

A script used to compute Constant Velocity Scans for velocity analysis.

• See the script cvs.m
• Download the script cvs.m

You will also need the following script:

• See the script plot_scans.m
• Download the script plot_scans.m

Constant Velocity Scans, NMO velocity increases from left to right. [scans.jpg]

2) Radon plus KL to eliminate Radon artifacts

We use Radon demultiple and the KL transform to eliminate multiples from a NMO corrected gather.
Parabolic Radon Transform: Hampson, D., 1986, Inverse velocity stacking for multiple elimination, Journal of the CSEG, vol 22, no 1., 44-55.
KL transform: Jones, I.F. and Levy, S. 1987, Signal-to-noise ratio enhancement in multichannel seismic data via the Karhunen-Loeve transform, Geophysical Prospecting 35,12-32.

• See the script radon_kl.m
• Download the script radon_kl.m

Panels from left to right: Original gather after NMO correction, Radon panel, estimated primaries (see residual energy from the multiples), estimated primaries after KL filtering, artifacts removed by the KL filter.
[radon_kl.jpg]

3) Spiking decon followed by Maximum Kurtosis Phase Correction

Spiking deconvolution followed by phase rotation and band-pass filtering were applied to a small data set from the Western Canadian Basin.

• See the script wcb_decon.m
• Download the script wcb_decon.m
• Download the data used for this example wcb_small_stack.su

The script should reproduce the following figures:

After spiking deconvolution we calculate the variation of the Kurtosis with phase rotation. We choose the phase correction that maximizes the kurtosis (Longbottom, J., Walden, A.T. and White, R.E., 1988, Principles and application of maximum kurtosis phase estimation, Geophysical Prospecting, 36, 115-138)

Spectra before and after spiking deconvolution and band-pass filtering.

From left to right: Original data, data after spiking deconvolution, data after spiking deconvolution followed by maximum kurtosis phase rotation and band-pass filtering