Method for attenuation of multiple reflections in seismic data

• US 7,675,812 B2
• Filed: 06/30/2008
• Issued: 03/09/2010
• Est. Priority Date: 06/30/2008
• Status: Expired due to Fees

• Current Assignee: PGS Geophysical As (PGS ASA)
• Original Assignee: PGS Geophysical As (PGS ASA)
• Inventors: Ferris, Edward James
• Primary Examiner(s): Keith; Jack W.
• Assistant Examiner(s): Saito; Krystine
• Application Number:
  US12/215,862
• Publication Number:
  US 20090323470A1
• Time in Patent Office: 617 Days
• Field of Search: 367 21- 24, 367/38, 367/49
• US Class Current: 367/24
• CPC Class Codes: G01V 1/364 Seismic filtering G01V1/37 takes precedence
• CPC Class Codes: G01V 1/3808 Seismic data acquisition, e.g. survey design
• CPC Class Codes: G01V 2210/21 Frequency-domain filtering, e.g. band pass
• CPC Class Codes: G01V 2210/56 De-ghosting; Reverberation compensation

Abstract

A method for attenuating multiple reflections in seismic signals includes coherency filtering the seismic signals. The seismic signals are low pass filtered and high pass filtered seismic signals at a frequency selected such that substantially only multiple reflections exist in the seismic signals above the frequency. A mask is generated having unity value where the high pass filtered signal amplitude exceeds a selected threshold. The mask has unity zero value for all other filtered signals. The mask value is multiplied by the low pass filtered signals. The multiplied signals are added to the high pass filtered signals and a coherency filter is applied to remove remaining primary energy to generate a model of the multiple reflections. The multiple reflection model is subtracted from the coherency filtered seismic signals.

Claims

1. A method for separating multiple reflections from marine seismic signals, comprising on a computer:

coherency filtering the seismic signals, the seismic signals acquired from a plurality of seismic sensors deployed in a body of water in response to actuation of a seismic energy source in the water;

low pass filtering the seismic signals at a selected frequency;

high pass filtering the seismic signals at the selected frequency, the frequency selected such that substantially only multiple reflections exist in the seismic signals above the frequency

generating a mask having unity value where the high pass filtered signal amplitude exceeds a selected threshold, the mask having zero value for all other high pass filtered signals;

multiplying the mask value by the low pass filtered signals;

adding the multiplied signals to the high pass filtered signals to generate a model of the multiple reflections;

subtracting the multiple reflection model from the coherency filtered seismic signals; and

at least one of storing and displaying a result of the subtracting.

2. The method of claim 1 further comprising sorting the seismic signals into common trace attribute gathers prior to the coherency filtering, low pass filtering and high pass filtering.

3. The method of claim 1 wherein the coherency filtering comprises F-X deconvolution.

4. The method of claim 1 further comprising coherency filtering a result of the subtracting.

5. The method of claim 4 wherein the coherency filtering the subtraction result comprises F-X deconvolution.

6. The method of claim 1 wherein the subtracting comprises adaptive subtracting.

7. The method of claim 1 wherein the frequency is about 40 Hz.

8. The method of claim 1 wherein a combined amplitude spectrum of a roll off portion of the low pass filtering and a roll off portion the high pass filter is substantially unity.

9. The method of claim 1 wherein the low pass filtering and the high pass filtering are substantially zero phase.

10. A computer readable medium comprising a computer program, the program having logic operable to cause a programmable computer to perform acts, comprising:

coherency filtering the seismic signals;

low pass filtering the seismic signals at a selected frequency;

high pass filtering the seismic signals at the selected frequency, the frequency selected such that substantially only multiple reflections exist in the seismic signals above the frequency

generating a mask having unity value where the high pass filtered signal amplitude exceeds a selected threshold, the mask having zero value for all other high pass filtered signals;

multiplying the mask value by the low pass filtered signals;

adding the multiplied signals to the high pass filtered signals to generate a model of the multiple reflections;

subtracting the multiple reflection model from the coherency filtered seismic signals; and

at least one of storing and displaying a result of the subtracting to generate multiple reflection attenuated seismic data.

11. The computer readable medium of claim 10 further comprising sorting the seismic signals into common trace attribute gathers prior to the coherency filtering, low pass filtering and high pass filtering.

12. The computer readable medium of claim 10 wherein the coherency filtering comprises F-X deconvolution.

13. The computer readable medium of claim 10 further comprising coherency filtering a result of the subtracting.

14. The computer readable medium of claim 10 wherein the coherency filtering the subtraction result comprises F-X deconvolution.

15. The computer readable medium of claim 10 wherein the subtracting comprises adaptive subtracting.

16. The computer readable medium of claim 10 wherein the frequency is about 40 Hz.

17. The computer readable medium of claim 10 wherein a combined amplitude spectrum of a roll off portion of the low pass filtering and a roll off portion the high pass filter is substantially unity.

18. The computer readable medium of claim 10 wherein the low pass filtering and the high pass filtering are substantially zero phase.

19. A method for seismic surveying, comprising:

deploying a plurality of seismic sensors at spaced apart locations in a body of water;

at selected times, actuating a seismic energy source in the water;

recording seismic signals detected by the seismic sensors;

coherency filtering the seismic signals;

low pass filtering the seismic signals at a selected frequency;

high pass filtering the seismic signals at the selected frequency, the frequency selected such that substantially only multiple reflections exist in the seismic signals above the frequency

generating a mask having unity value where the high pass filtered signal amplitude exceeds a selected threshold, the mask having zero value for all other high pass filtered signals;

multiplying the mask value by the low pass filtered signals;

adding the multiplied signals to the high pass filtered signals to generate a model of the multiple reflections;

subtracting the multiple reflection model from the coherency filtered seismic signals; and

at least one of storing and displaying a result of the subtracting.

20. The method of claim 19 further comprising sorting the seismic signals into common trace attribute gathers prior to the coherency filtering, low pass filtering and high pass filtering.

21. The method of claim 19 wherein the coherency filtering comprises F-X deconvolution.

22. The method of claim 19 further comprising coherency filtering a result of the subtracting.

23. The method of claim 22 wherein the coherency filtering the subtraction result comprises F-X deconvolution.

24. The method of claim 19 wherein the subtracting comprises adaptive subtracting.

25. The method of claim 19 wherein the frequency is about 40 Hz.

26. The method of claim 19 wherein a combined amplitude spectrum of a roll off portion of the low pass filtering and a roll off portion the high pass filter is substantially unity.

27. The method of claim 19 wherein the low pass filtering and the high pass filtering are substantially zero phase.