REVERSE TIME MIGRATION NOISE REMOVAL USING A DEPTH-DEPENDENT WAVENUMBER FILTER
Examples of methods and systems are disclosed. The methods may include obtaining a seismic data regarding a subsurface region of interest, wherein the seismic data comprises a plurality of time-space waveforms. The methods may also include obtaining a seismic velocity model. The methods may further include determining a migrated seismic image based on the plurality of time-space waveforms and the seismic velocity model. The method may still further include generating a filtered seismic image by applying a depth-dependent attenuation operator to the migrated seismic image and determining a drilling target in the subsurface region based on the filtered seismic image.
1 . A method, comprising:
obtaining a seismic data regarding a subsurface region of interest, wherein the seismic data comprises a plurality of time-space waveforms;
obtaining a seismic velocity model;
determining a migrated seismic image based on the plurality of time-space waveforms and the seismic velocity model;
generating a filtered seismic image by applying a depth-dependent attenuation operator to the migrated seismic image; and
determining a drilling target in the subsurface region based on the filtered seismic image.
2 . The method of claim 1 , further comprising planning, using a wellbore planning system, a planned wellbore trajectory to intersect the drilling target.
3 . The method of claim 2 , further comprising drilling, using a drilling system, a wellbore guided by the planned wellbore trajectory.
4 . The method of claim 1 , wherein the depth-dependent attenuation operator comprises a function that decreases with depth.
5 . The method of claim 1 , wherein generating a filtered seismic image further comprises:
generating a transformed seismic image by transforming the migrated seismic image into components of a vertical wavenumber;
multiplying the transformed seismic image by a depth-dependent weight; and
performing an inverse transform.
6 . The method of claim 1 , wherein determining the migrated seismic image comprises performing reverse time migration.
7 . The method of claim 1 , wherein the depth-dependent attenuation operator comprises a power function of a vertical wavenumber.
8 . The method of claim 5 , wherein generating a transformed seismic image comprises performing short-window Fourier Transforms.
9 . A non-transitory computer-readable medium comprising computer-executable instructions stored thereon that, when executed on a processor, cause the processor to perform:
obtaining seismic data regarding a subsurface region of interest, wherein the seismic data comprises a plurality of time-space waveforms;
obtaining a seismic velocity model;
determining a migrated seismic image based on the plurality of time-space waveforms and the seismic velocity model;
generating a filtered seismic image by applying a depth-dependent attenuation operator to the migrated seismic image; and
determining a drilling target in the subsurface region based on the filtered seismic image.
10 . The non-transitory computer-readable medium of claim 9 , further comprising computer-executable instructions that cause the processor to perform:
generating a transformed seismic image by transforming the migrated seismic image into components of a vertical wavenumber;
multiplying the transformed seismic image by a depth-dependent weight; and
performing an inverse transform.
11 . The non-transitory computer-readable medium of claim 9 , further comprising computer-executable instructions that cause the processor to perform:
determining the migrated seismic image by performing reverse time migration.
12 . The non-transitory computer-readable medium of claim 9 , further comprising computer-executable instructions that cause the processor to perform:
generating a transformed seismic image by performing short-window Fourier Transforms.
13 . A system, comprising:
a seismic acquisition system configured to record seismic data regarding a subsurface region of interest, wherein the seismic data comprises a plurality of time-space waveforms; and
a seismic processor configured to receive the seismic data and to:
obtain a seismic velocity model,
determine a migrated seismic image based on the plurality of time-space waveforms and the seismic velocity model,
generate a filtered seismic image by applying a depth-dependent attenuation operator to the migrated seismic image, and
determine a drilling target in the subsurface region based on the filtered seismic image.
14 . The system of claim 13 , further comprising a wellbore planning system configured to plan a planned wellbore trajectory based, at least in part, on the filtered seismic image.
15 . The system of claim 14 , further comprising a drilling system configured to drill a portion of a wellbore guided by the planned wellbore trajectory.
16 . The system of claim 13 , wherein the depth-dependent attenuation operator comprises a function that decreases with depth.
17 . The system of claim 13 , wherein the seismic processor is further configured to:
generate a transformed seismic image by transforming the migrated seismic image into components of a vertical wavenumber;
multiply the transformed seismic image by a depth-dependent weight; and
perform an inverse transform.
18 . The system of claim 13 , wherein the seismic processor is further configured to determine the migrated seismic image by performing reverse time migration.
19 . The system of claim 13 , wherein the depth-dependent attenuation operator comprises a power function of a vertical wavenumber.
20 . The system of claim 13 , wherein the seismic processor is further configured to generate a transformed seismic image by performing short-window Fourier Transforms.