3-D TAU-P interpolation
View Patent ↗Methods and apparatus for removing noise from signal traces collected during a seismic gather, particularly removing the aliased energy in the time-slowness (tau-P) domain so that aliasing noise does not lead to high noise levels in the inverse transformed data are provided. Removal of the aliasing noise may provide for quieter seismic traces and may increase the likelihood for successful three-dimensional (3-D) tau-P interpolation and detection of seismic events.
1. A computer-readable medium containing a program for removing aliasing noise from three-dimensional (3-D) time-slowness (tau-P) data, which, when executed by a processor, performs operations comprising:
decomposing the tau-P data into a plurality of frequency bands, wherein a template frequency band in the plurality of frequency bands does not contain aliasing noise and at least a portion of the other frequency bands constitute a subset of frequency bands for processing;
calculating an envelope for the template frequency band and for each of the frequency bands in the subset to generate a plurality of envelope cuboids;
for each frequency band in the subset from the lowest to the highest frequency band in the subset, multiplying an envelope cuboid associated with the present frequency band with an envelope cuboid associated with an adjacent lower frequency band or with the template frequency band to produce an envelope mask for the present frequency band, saving the envelope mask as the envelope cuboid associated with the present frequency band, and multiplying the envelope mask for the present frequency band with the tau-P data for the present frequency band to produce processed tau-P data for the present frequency band such that the aliasing noise is removed; and
summing the tau-P data for the template frequency band and the processed tau-P data for each of the frequency bands in the subset to produce tau-P data without aliasing noise.
2. The computer-readable medium of claim 1 , wherein the operations further comprise:
normalizing each time slice within each of the plurality of envelope cuboids; and
applying an amplitude filter to each normalized time slice within each of the plurality of envelope cuboids such that envelope amplitudes below a threshold value are attenuated or zeroed and envelope amplitudes above the threshold value are set to one before multiplying the envelope cuboids.
3. The computer-readable medium of claim 1 , wherein the operations further comprise removing edge effects from the template envelope cuboid before multiplying the envelope cuboids.
4. The computer-readable medium of claim 3 , wherein the operations further comprise:
calculating a center of amplitude for each time slice in the envelope cuboid associated with the template frequency band;
calculating a rotational moment of inertia with respect to the center of amplitude for each time slice in the template envelope cuboid; and
for each time slice in the template envelope cuboid, attenuating or zeroing the time slice if the moment of inertia is greater than an inertia threshold value before multiplying the envelope cuboids.
5. The computer-readable medium of claim 4 , wherein the operations further comprise:
normalizing each time slice within the template envelope cuboid; and
applying an amplitude filter to each normalized time slice within the template envelope cuboid such that envelope amplitudes below a template threshold value are attenuated or zeroed and envelope amplitudes above the template threshold value are set to one before calculating the center of amplitude.
6. The computer-readable medium of claim 1 , wherein the operations further comprise performing a linear Radon transform on a window of time-offset data to form the 3-D tau-P data.
7. The computer-readable medium of claim 1 , wherein the operations further comprise performing a linear inverse Radon transform on the tau-P data without aliasing noise to produce time-offset data without aliasing noise.