Method for separating seismic signals from two or more distinct sources

• US 6,882,938 B2
• Filed: 07/30/2003
• Issued: 04/19/2005
• Est. Priority Date: 07/30/2003
• Status: Active Grant

• Current Assignee: PGS Americas, Inc. (PGS ASA)
• Original Assignee: PGS Americas, Inc. (PGS ASA)
• Inventors: Brittan, John, Martinez, Ruben D., Vaage, Svein Torleif
• Primary Examiner(s): Barlow, John
• Assistant Examiner(s): Walling, Meagan S
• Application Number:
  US10/630,385
• Publication Number:
  US 20050027454A1
• Time in Patent Office: 629 Days
• Field of Search: 702/14, 702/15, 702/16, 702/17, 702/18, 367/38, 367/43
• US Class Current: 702/17
• CPC Class Codes: G01V 1/364 Seismic filtering G01V1/37 takes precedence
• CPC Class Codes: G01V 2210/20 Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out

Abstract

A method is disclosed for separating energy resulting from actuating at least two different seismic energy sources from seismic signals. The sources are actuated to provide a variable time delay between successive actuations of a first one and a second one of the sources. The method includes sorting the seismic signals such that events therein resulting from actuations of the first source are substantially coherent in all spatial directions, coherency filtering the first source coherency sorted signals, sorting the seismic signals such that events therein resulting from actuations of the second source are substantially coherent in all spatial directions, and coherency filtering the second source coherency sorted signals.

Claims

1. A method for separating energy resulting from actuating at least two different seismic energy sources from seismic signals, the method comprising:

sorting the seismic signals such that events therein resulting from actuations of the first source are substantially coherent in all spatial directions, the seismic signals resulting from a first one and a second one of the at least two sources actuated such that there is a variable time delay between successive actuations thereof;

coherency filtering the first source coherency sorted signals;

sorting the seismic signals such that events therein resulting from actuations of the second source are substantially coherent in all spatial directions; and

coherency filtering the second source coherency sorted signals.

2. The method of claim 1 further comprising:

subtracting the coherency filtered first source coherency sorted signals and the coherency filtered second source coherency sorted signals from the seismic signals to generate residual seismic signals;

sorting the residual seismic signals such that events therein resulting from actuations of the first source are substantially coherent in all spatial directions;

coherency filtering the first source coherency sorted residual signals;

sorting the residual seismic signals such that events therein resulting from actuations of the second source are substantially coherent in all spatial directions; and

coherency filtering the second source coherency sorted residual signals.

3. The method of claim 2 further comprising:

determining whether a value of a parameter in the coherency filtered, first source and second source sorted residual signals exceeds a selected threshold, the parameter related to an amount of energy in a seismic signal;

subtracting the coherency filtered, first source sorted residual signals and the coherency filtered second source sorted residual signals from the residual signals; and

repeating the sorting and coherency filtering until the value drops below the selected threshold.

4. The method of claim 3 wherein the parameter comprises peak trace amplitude.

5. The method of claim 3 wherein the parameter comprises trace energy.

6. The method of claim 3 wherein the parameter comprises average trace amplitude.

7. The method of claim 2 further comprising adding the coherency filtered, first source sorted residual signals to the coherency filtered, first source sorted seismic signals.

8. The method of claim 2 further comprising adding the coherency filtered, second source sorted residual signals to the coherency filtered, second source sorted seismic signals.

9. The method of claim 1 wherein the sorting the seismic signals comprises generating a common channel plane gather with respect to the first source.

10. The method of claim 1 wherein the sorting the seismic signals comprises generating a common channel plane gather with respect to the second source.

11. The method of claim 1 wherein the coherency filtering comprises slant stack processing.

12. The method of claim 11 wherein the slant stack filtering comprises transforming traces into the tau-p domain, excluding portions of the transformed traces corresponding to energy other than a coherency reference, and inverse transforming portions of the traces having the excluded energy into the time-space domain.

13. The method of claim 12 wherein the transforming into the tau-p domain comprises performing a Radon transform.

14. The method of claim 1 wherein the sorting comprises sorting seismic traces into the common channel/common shot domain.

15. The method of claim 1 wherein the sorting comprises sorting seismic traces into the common depth point/common offset domain.

16. A method for seismic surveying, comprising:

towing a first seismic energy source and at least one seismic sensor system;

towing a second seismic energy source at a selected distance from the first seismic energy source;

actuating the first seismic energy source and the second seismic energy source in a plurality of firing sequences, each of the firing sequences including firing of the first source and the second source and recording signals generated by the at least one seismic sensor system, a time interval between firing the first source and the second source varied between successive ones of the firing sequences;

sorting the seismic signals such that events therein resulting from actuations of the first source are substantially coherent in all spatial directions;

coherency filtering the first source coherency sorted signals;

sorting the seismic signals such that events therein resulting from actuations of the second source are substantially coherent in all spatial directions; and

coherency filtering the second source coherency sorted signals.

17. The method of claim 16 further comprising:

subtracting the coherency filtered first source coherency sorted signals and the coherency filtered second source coherency sorted signals from the seismic signals to generate residual seismic signals;

sorting the residual seismic signals such that events therein resulting from actuations of the first source are substantially coherent in all spatial directions;

coherency filtering the first source coherency sorted residual signals;

sorting the residual seismic signals such that events therein resulting from actuations of the second source are substantially coherent in all spatial directions; and

coherency filtering the second source coherency sorted residual signals.

18. The method of claim 16 further comprising:

determining whether a value of a parameter in the coherency filtered, first source and second source sorted residual signals exceeds a selected threshold, the parameter related to an amount of energy in a seismic signal;

subtracting the coherency filtered, first source sorted residual signals and the coherency filtered second source sorted residual signals from the residual signals; and

repeating the sorting and coherency filtering until the value drops below the selected threshold.

19. The method of claim 18 wherein the parameter comprises peak trace amplitude.

20. The method of claim 18 wherein the parameter comprises trace energy.

21. The method of claim 18 wherein the parameter comprises average trace amplitude.

22. The method of claim 17 further comprising adding the coherency filtered, first source sorted residual signals to the coherency filtered, first source sorted seismic signals.

23. The method of claim 17 further comprising adding the coherency filtered, second source sorted residual signals to the coherency filtered, second source sorted seismic signals.

24. The method of claim 16 wherein the sorting the seismic signals comprises generating a common channel plane gather with respect to the first source.

25. The method of claim 16 wherein the sorting the seismic signals comprises generating a common channel plane gather with respect to the second source.

26. The method of claim 16 wherein the coherency filtering comprises slant stack processing.

27. The method of claim 26 wherein the slant stack processing comprises transforming traces into the tau-p domain, excluding portions of the transformed traces corresponding to energy other than a coherency reference, and inverse transforming portions of the traces having the excluded energy into the time-space domain.

28. The method of claim 27 wherein the transforming into the tau-p domain comprises performing a Radon transform.

29. The method of claim 16, wherein the time interval is varied systematically.

30. The method of claim 16 wherein the time interval is varied quasi-randomly.

31. The method of claim 16, wherein the time interval varied is randomly.

32. The method of claim 16, wherein the time interval is varied in steps of about 100 milliseconds.

33. The method as defined in claim 16, wherein the time interval is at least as long as a wavelet time of the first source.

34. The method of claim 16 wherein the sorting comprises sorting seismic traces into the common channel/common shot domain.

35. The method of claim 16 wherein the sorting comprises sorting seismic traces into the common depth point/common offset domain.

36. A program stored in a computer readable medium, the program including logic operable to cause a programmable computer to perform steps comprising:

sorting seismic signals resulting from actuating at least two different seismic energy sources, the sources actuated to provide a variable time delay between successive actuations of a first one and a second one of the sources, the sorting performed such that events therein resulting from actuations of the first source are substantially coherent in all spatial directions;

coherency filtering the first source coherency sorted signals;

sorting the seismic signals such that events therein resulting from actuations of the second source are substantially coherent in all spatial directions; and

coherency filtering the second source coherency sorted signals.

37. The program of claim 36 further comprising logic operable to cause the computer to perform:

subtracting the coherency filtered first source coherency sorted signals and the coherency filtered second source coherency sorted signals from the seismic signals to generate residual seismic signals;

sorting the residual seismic signals such that events therein resulting from actuations of the first source are substantially coherent in all spatial directions;

coherency filtering the first source coherency sorted residual signals;

sorting the residual seismic signals such that events therein resulting from actuations of the second source are substantially coherent in all spatial directions; and

coherency filtering the second source coherency sorted residual signals.

38. The program of claim 37 further comprising logic operable to cause the computer to perform:

determining whether a value of a parameter in the coherency filtered, first source and second source sorted residual signals exceeds a selected threshold, the parameter related to an amount of energy in a seismic signal;

subtracting the coherency filtered, first source sorted residual signals and the coherency filtered second source sorted residual signals from the residual signals; and

repeating the sorting and coherency filtering until the value drops below the selected threshold.

39. The program of claim 36 wherein the parameter comprises peak trace amplitude.

40. The program of claim 36 wherein the parameter comprises trace energy.

41. The program of claim 36 wherein the parameter comprises average trace amplitude.

42. The program of claim 37 further comprising logic operable to cause the computer to perform adding the coherency filtered, first source sorted residual signals to the coherency filtered, first source sorted seismic signals.

43. The program of claim 37 further comprising logic operable to cause the computer to perform adding the coherency filtered, second source sorted residual signals to the coherency filtered, second source sorted seismic signals.

44. The program of claim 36 wherein the sorting the seismic signals comprises generating a common channel plane gather with respect to the first source.

45. The program of claim 36 wherein the sorting the seismic signals comprises generating a common channel plane gather with respect to the second source.

46. The program of claim 36 wherein the coherency filtering comprises slant stack processing.

47. The program of claim 46 wherein the slant stack filtering comprises transforming traces into the tau-p domain, excluding portions of the transformed traces corresponding to energy other than a coherency reference, and inverse transforming portions of the traces having the excluded energy into the time-space domain.

48. The program of claim 47 wherein the transforming into the tau-p domain comprises performing a Radon transform.

49. The program of claim 36 wherein the sorting comprises sorting seismic traces into the common channel/common shot domain.

50. The program of claim 36 wherein the sorting comprises sorting seismic traces into the common depth point/common offset domain.