Seismic wavefield decomposition using recursive radon transform

A non-transitory computer-readable memory having computer-executable instructions stored thereon that, when executed by a computer processor, perform steps comprising:

• receiving a frequency-domain seismic wavefield associated to a plurality of grid points,
• wherein the frequency-domain seismic wavefield is organized into a plurality of slices,
• wherein each of the plurality of slices is organized into a sequence of overlapping windows,
• wherein each of the sequence of overlapping windows comprises grid points among the plurality of grid points,
• wherein each n-th overlapping window in the sequence of overlapping windows comprises incremental grid points among the grid points and common grid points among the grid points, and
• wherein n is an integer greater than one;
• for each of the plurality of slices in turn:
• for each of the grid points in a first overlapping window in the sequence of overlapping windows:
• determining a value of a first intermediate function using, at least in part, the frequency-domain seismic wavefield,
• for each n-th overlapping window in turn:
• for each of the incremental grid points:
• determining a value of the first intermediate function using, at least in part, the frequency-domain seismic wavefield,
• for each of the common grid points:
• determining a value of the first intermediate function by selecting the previously-determined value of the first intermediate function;
• determining a Radon-transformed seismic wavefield using, at least in part, the value of the first intermediate function for the plurality of grid points.

A system comprising:

• a seismic processing system configured to:
• receive a frequency-domain seismic wavefield associated to a plurality of grid points,
• wherein the frequency-domain seismic wavefield is organized into a plurality of slices,
• wherein each of the plurality of slices is organized into a sequence of overlapping windows,
• wherein each of the sequence of overlapping windows comprises grid points among the plurality of grid points,
• wherein each n-th overlapping window in the sequence of overlapping windows comprises incremental grid points among the grid points and common grid points among the grid points, and
• wherein n is an integer greater than one,
• for each of the plurality of slices in turn:
• for each of the grid points in a first overlapping window in the sequence of overlapping windows:
• determine a value of a first intermediate function using, at least in part, the frequency-domain seismic wavefield;
• for each n-th overlapping window in turn:
• for each of the incremental grid points:
• determine a value of the first intermediate function using, at least in part, the frequency-domain seismic wavefield;
• for each of the common grid points:
• determine a value of the first intermediate function by selecting the previously-determined value of the first intermediate function,
• determine a Radon-transformed seismic wavefield using, at least in part, the value of the first intermediate function for the plurality of grid points, and
• determine a migrated seismic image using, at least in part, the frequency-domain seismic wavefield and the Radon-transformed seismic wavefield;
• and
• a seismic interpretation workstation configured to:
• determine a location of a hydrocarbon reservoir within a subterranean region of interest using, at least in part, the migrated seismic image,
• wherein the subterranean region of interest is represented by the plurality of grid points.