Method apparatus and system for migration noise attenuation and image enhancement

Patent metadata

Field	Value
Publication number	US-9575194-B2
Application number	US-201414278039-A
Country	US
Kind code	B2
Filing date	May 15, 2014
Priority date	May 1, 2013
Publication date	Feb 21, 2017
Grant date	Feb 21, 2017

Abstract

A method for processing seismic data includes receiving seismic data and a velocity model (c(x)) for a plurality of locations (x), scaling a dimension of the seismic data according to the velocity model (c(x)) to provide a velocity normalized seismic data, and generating a final image (S(x)) of the subsurface using the velocity normalized seismic data. The velocity normalized seismic data may be a reverse-time migration image (I(x,ξ)) corresponding to the plurality of locations (x) and a plurality of propagation distance offsets (ξ). The method may also include transforming the reverse-time migration image (I(x,ξ)) for the plurality of selected positions (x) to a wavenumber domain to provide velocity normalized wavenumber data (I(k,ψ)) and suppressing data components corresponding to non-physical or undefined reflection angles to provide enhanced wavenumber data (I′(k,ψ)) and using the enhanced wavenumber data (I′(k,ψ)) to generate the final image (S(x)). A corresponding apparatus is also disclosed herein.

First claim

What is claimed is:

1. A method for processing seismic data corresponding to a subsurface, the method comprising:
receiving seismic image data corresponding to a plurality of locations (x);
receiving a velocity model (c(x)) for the plurality of locations (x);
scaling with a computing device a time-shift dimension (τ) associated with the seismic image data according to the velocity model (c(x)) to provide velocity normalized seismic data; and
generating a final image (S(x)) of the subsurface using the velocity normalized seismic data,
wherein the step of scaling converts the time-shift dimension (τ) to a propagation distance offset (ξ) based on the velocity model (c(x)).

2. The method of claim 1 , wherein the seismic image data corresponds to a reverse time migration, beam migration, Kirchhoff migration, one way migration, wave equation, or other migration.

3. The method of claim 1 , wherein the seismic image data is in a time or depth domain.

4. The method of claim 1 , wherein the velocity normalized seismic data is a migration image (I(x,ξ)) obtained by correlation and corresponding to the plurality of locations (x) and a plurality of propagation distance offsets (ξ).

5. The method of claim 4 , wherein the migration image (I(x,ξ)) is generated by post-correlation scaling of the time-shift dimension (τ) according to the velocity model (c(x)).

6. The method of claim 4 , wherein the migration image (I(x,ξ)) is generated by pre-correlation scaling of the time-shift dimension (τ) according to the velocity model (c(x)).

7. The method of claim 4 , wherein generating the final image (S(x)) of the subsurface comprises transforming the migration image (I(x,ξ)) for the plurality of selected positions (x) to a wavenumber domain to provide velocity normalized wavenumber data (I(k,ψ)) and suppressing data components corresponding to non-physical or undefined reflection angles to provide enhanced wavenumber data (I′(k,ψ)).

8. The method of claim 7 , comprising inverse transforming the enhanced wavenumber data (I′(k,ψ)) to provide enhanced offset distance correlation data (I′(x,ξ)).

9. The method of claim 8 , comprising selecting zero offset data for the offset distance correlation data (I′(x,ξ)).

10. The method of claim 7 , comprising inverse transforming the enhanced wavenumber data (I′(k,ψ)) to provide wavenumber offset data (I′(k,ξ)) and inverse phase shifting the wavenumber offset data (I′(k,ξ)) to provide de-propagated wavenumber data (J(k,ξ)).

11. The method of claim 10 , comprising inverse transforming the de-propagated wavenumber data (J(k,ξ)) to provide de-propagated correlation data (J(x,ξ)).

12. The method of claim 11 , comprising stacking a subset of de-propagated correlation data (J(x,ξ)) for the plurality of locations (x).

13. A computing device configured to process seismic data corresponding to a subsurface, the computing device comprising:
an interface configured to receive seismic image data corresponding to a plurality of locations (x);
a processor configured to scale a time-shift dimension (τ) associated with the seismic image data, or an image derived therefrom, according to a velocity model (c(x)) to provide velocity normalized seismic data; and
the processor configured to generate a final image (S(x)) of the subsurface using the velocity normalized seismic data,
wherein the processor scales the time-shift dimension (τ) by converting the time-shift dimension (τ) to a propagation distance offset (ξ) based on the velocity model (c(x)).

14. The apparatus of claim 13 , wherein the seismic image data corresponds to a reverse time migration, beam migration, Kirchhoff migration, one way migration, or wave equation.

15. The apparatus of claim 13 , wherein the seismic image data is in a time or depth domain.

16. The apparatus of claim 13 , wherein the velocity normalized seismic data is a reverse-time migration image (I(x,ξ)) corresponding to the plurality of locations (x) and a plurality of propagation distance offsets (ξ).

17. The apparatus of claim 16 , wherein the processor is configured to generate the final image (S(x)) of the subsurface by transforming the reverse-time migration image (I(x,ξ)) for a plurality of selected positions (x) to a wavenumber domain to provide velocity normalized wavenumber data (I(k,ψ)) and suppressing data components corresponding to non-physical or undefined reflection angles to provide enhanced wavenumber data (I′(k,ψ)).

18. The apparatus of claim 17 , wherein the processor is configured to inverse transform the enhanced wavenumber data (I′(k,ψ)) to provide wavenumber offset data (I′(k,ξ)) and inverse phase shift the wavenumber offset data (I′(k,ξ)) to provide de-propagated wavenumber data (J(k,ξ)).

19. The apparatus of claim 18 , wherein the processor is configured to inverse transform the de-propagated wavenumber data (J(k,ξ)) to provide de-propagated correlation data (J(x,ξ)).

20. A method for processing seismic data corresponding to a subsurface, the method comprising:
correlating a forward propagated source wavefield (D s ) corresponding to a seismic source with a backward propagated receiver wavefield (D r ) corresponding to a seismic receiver as a function of a time shift (τ) to provide time shift correlation data (I(x,τ)) for a selected position (x);
repeating the correlating operation for a plurality of selected positions (x) to provide time shift correlation data (I(x,τ)) for the plurality of selected positions (x);
scaling with a computing system the time shift (τ) of the time shift correlation data (I(x,τ)) according to a velocity mapping function to provide offset distance correlation data (I(x,ξ));
transforming the offset distance correlation data (I(x,ξ)) for the plurality of selected positions (x) to a wavenumber domain to provide velocity normalized wavenumber data (I(k,ψ)) and suppressing data components corresponding to non-physical or undefined reflection angles to provide enhanced wavenumber data (I′(k,ψ)); and
generating a final image (S(x)) of the subsurface using the enhanced wavenumber data (I′(k,ψ)),
wherein the step of scaling converts the time shift (τ) to a propagation distance offset (ξ) according to the velocity mapping function.

Assignees

• Cgg Services Sa
• Cgg Services Sas

Inventors

• Khalil Adel
• Sun James
• Zhang Yu
• Poole Gordon

Classifications

• G01V2210/51

  Migration · CPC title

• G01V2210/679

  Reverse-time modeling or coalescence modelling, i.e. starting from receivers · CPC title

• G01V1/282 Primary

  Application of seismic models, synthetic seismograms · CPC title