Method, device and processing algorithm for multiple and noise elimination from marine seismic data

1. A method for simultaneously denoising and attenuating multiples in seismic data recorded with seismic receivers, the method comprising:

receiving the seismic data, wherein the seismic data includes a pressure P component and a vertical Z component;

separating the seismic data into up-going U wave-fields and down-going D wave-fields;

calculating an up-down deconvolution R based on the up-going U and the down-going D wave-fields;

generating a multiple model M based on the up-down deconvolution R; and

adaptively subtracting the multiple model M from the pressure P component to obtain a corrected pressure P free component from which simultaneously the noise and the multiples were removed,

wherein a multiple is a wave-field that propagates from a source to a receiver and experiences more than one reflection.

2. The method of claim 1 , wherein the pressure P component is measured with a hydrophone and the vertical Z component is measured with a geophone or accelerometer and both components include noise and multiples.

3. The method of claim 1 , wherein the step of calculating the up-down deconvolution R comprises:

deconvolving the up-going U and the down-going D wave-fields such that R=U/D.

4. The method of claim 1 , wherein the step of generating the multiple model M comprises:

defining the multiple model M, in a frequency-wavenumber domain, as a product of the up-down deconvolution R and the pressure P component.

5. The method of claim 4 , further comprising:

redatuming the product of the up-down deconvolution R and the pressure P component with an exponential that includes a wavenumber k z and a vertical distance Δz, where the vertical distance defines a distance between an ocean-bottom and a point where the multiple model M is calculated.

6. The method of claim 1 , further comprising:

applying a debubbling operator to the pressure P component and the vertical Z component prior to separating the seismic data.

7. The method of claim 1 , further comprising:

applying directional designature of the pressure P component before calculating the up-down deconvolution R.

8. The method of claim 7 , wherein the directional designature comprises:

deconvolving the pressure P component with a muted version of the down-going D component.

9. The method of claim 8 , further comprising:

adaptively subtracting the multiple model M from the designatured pressure component to obtain the corrected pressure P free .

10. A computing device for simultaneously denoising and attenuating multiples in seismic data recorded with seismic receivers, the computing device comprising:

an interface configured to receive the seismic data, wherein the seismic data includes a pressure P component and a vertical Z component; and

a processor connected to the interface and configured to, separate the seismic data into up-going U wave-fields and down-going D wave-fields,

calculate an up-down deconvolution R based on the up-going U and the down-going D wave-fields,

generate a multiple model M based on the up-down deconvolution R, and

adaptively subtract the multiple model M from the pressure P component to obtain a corrected pressure P free component from which simultaneously the noise and the multiples were removed,

wherein a multiple is a wave-field that propagates from a source to a receiver and experiences more than one reflection.

11. The computing device of claim 10 , wherein the pressure P component is measured with a hydrophone and the vertical Z component is measured with a geophone or accelerometer and both components include noise and multiples.

12. The computing device of claim 10 , wherein the processor is further configured to:

deconvolve the up-going U and the down-going D wave-fields such that R=U/D.

13. The computing device of claim 10 , wherein the processor is further configured to:

define the multiple model M, in a frequency-wavenumber domain, as a product of the up-down deconvolution R and the pressure P component.

14. The computing device of claim 13 , wherein the processor is further configured to:

redatum the product of the up-down deconvolution R and the pressure P component with an exponential that includes a wavenumber k z and a vertical distance Δz, where the vertical distance defines a distance between an ocean-bottom and a point where the multiple model M is calculated.

15. The computing device of claim 10 , wherein the processor is further configured to:

apply a debubbling operator to the pressure P component and the vertical Z component prior to separating the seismic data.

16. The computing device of claim 10 , wherein the processor is further configured to:

apply directional designature of the pressure P component before calculating the up-down deconvolution R.

17. The computing device of claim 16 , wherein the processor is further configured to:

deconvolve the pressure P component with a muted version of the down-going D component.

18. The computing device of claim 17 , wherein the processor is further configured to:

adaptively subtract the multiple model M from the designatured pressure component to obtain the corrected pressure P free .

19. A computer readable medium configured to permanently store computer instructions such that, when a computer executes the instructions, to simultaneously denoise and attenuate multiples in seismic data recorded with seismic receivers ( 20 ), the instructions comprising:

receiving the seismic data, wherein the seismic data includes a pressure P component and a vertical Z component;

separating the seismic data into up-going U wave-fields and down-going D wave-fields;

calculating an up-down deconvolution R based on the up-going U and the down-going D wave-fields;

generating a multiple model M based on the up-down deconvolution R; and

adaptively subtracting the multiple model M from the pressure P component to obtain a corrected pressure P free component from which simultaneously the noise and the multiples were removed,

wherein a multiple is a wave-field that propagates from a source to a receiver and experiences more than one reflection.

20. The medium of claim 19 , wherein the pressure P component is measured with a hydrophone and the vertical Z component is measured with a geophone and both components include noise and multiples.