Processing of seismic data using the S-transform

• US 7,617,053 B2
• Filed: 05/14/2007
• Issued: 11/10/2009
• Est. Priority Date: 05/12/2006
• Status: Active Grant

• Current Assignee: Calgary Scientific Inc.
• Original Assignee: Calgary Scientific Inc.
• Inventors: Pinnegar, C. Robert, McGrory, Robert J., Stewart, Robert R.
• Primary Examiner(s): Dunn, Drew A
• Assistant Examiner(s): VO, HIEN XUAN
• Application Number:
  US11/798,371
• Publication Number:
  US 20070299613A1
• Time in Patent Office: 911 Days
• Field of Search: 702/14, 702/18, 367/7, 367/9, 367/14
• US Class Current: 702/18
• CPC Class Codes: G01V 1/28 Processing seismic data, e.g. for interpretation or for event detection G01V1/48 takes precedence

Abstract

The present invention relates to a method and system for processing seismic data. Seismic data indicative of at least a time-dependent seismic trace are transformed into time-frequency domain using the S-transform. The transformed seismic data are processed for determining at predetermined time instances an instantaneous amplitude and a dominant instantaneous frequency. Based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency seismic data indicative of a likelihood of a presence of hydrocarbons are then determined which are then graphically displayed. The method and system for processing seismic data provides a powerful tool for determining the likelihood of a presence of hydrocarbons based on a double peak hydrocarbon signature in the quotient of the instantaneous amplitude to the dominant instantaneous frequency.

Claims

1. A method for processing seismic data comprising:

a processor in communication with the input port for sensing seismic data;

a) receiving the seismic data, the seismic data being indicative of at least a time-dependent seismic trace;

b) determining second seismic data by transforming the seismic data into time-frequency domain using a time-frequency transform;

c) processing the second seismic data for determining at predetermined time instances an instantaneous amplitude and a dominant instantaneous frequency;

d) determining at the predetermined time instances third seismic data indicative of a likelihood of a presence of hydrocarbons based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency; and

,

e) providing the third seismic data.

2. A method for processing seismic data as defined in claim 1 wherein c) comprises:

determining at each predetermined time instance a maximum amplitude of the modulus of the time-frequency transform and associating the instantaneous amplitude therewith; and

,

determining at each predetermined time instance the dominant instantaneous frequency as the frequency where the maximum amplitude of the modulus of the time-frequency transform occurs.

3. A method for processing seismic data as defined in claim 1 wherein the seismic data are transformed using the Stockwell transform.

4. A method for processing seismic data as defined in claim 1 wherein a) the seismic data are received with the seismic data being indicative of a plurality of time-dependent seismic traces of one of a two-dimensional migrated section and a three-dimensional migrated section.

5. A method for processing seismic data as defined in claim 4 wherein b) comprises transforming each of a plurality of seismic traces into time-frequency domain.

6. A method for processing seismic data as defined in claim 5 comprising processing the second seismic data for averaging a modulus of the time-frequency transform of seismic traces corresponding to at least a predetermined portion of the one of a two-dimensional migrated section and a three-dimensional migrated section.

7. A method for processing seismic data as defined in claim 1 comprising determining plotting data for plotting the third seismic data as a function of space and time.

8. A method for processing seismic data as defined in claim 7 comprising displaying the plotting data in a human comprehensible fashion.

9. A method for processing seismic data as defined in claim 8 wherein the plotting data comprise control data for enabling user interaction.

10. A method for processing seismic data as defined in claim 1 comprising:

processing the third seismic data for detecting first occurrences where the quotient of the instantaneous amplitude to the dominant instantaneous frequency is above a predetermined first threshold;

processing the third seismic data of seismic traces where first occurrences have been detected for detecting second occurrences where the quotient of the instantaneous amplitude to the dominant instantaneous frequency is above a predetermined second threshold;

determining fourth seismic data indicative of seismic traces comprising the first and the second occurrences; and

,

providing the fourth seismic data.

11. A method for processing seismic data as defined in claim 10 comprising:

determining the likelihood of a presence of hydrocarbons in dependence upon the fourth seismic data.

12. A method for processing seismic data as defined in claim 11 comprising:

processing the third and fourth seismic data for determining an approximate depth of a hydrocarbon bearing layer based on time instances of the first occurrences.

13. A method for processing seismic data as defined in claim 12 comprising:

processing the third and fourth seismic data for determining an approximate thickness of a hydrocarbon bearing layer based on time lapses between the time instances of the first occurrences and time instances of the second occurrences.

14. A method for processing seismic data as defined in claim 10 comprising

determining average amplitudes of the modulus of the time-frequency transform amplitudes of seismic traces corresponding to predetermined portions of the one of a two-dimensional migrated section and a three-dimensional migrated section.

15. A method for processing seismic data comprising:

a processor in communication with the input port for sensing seismic data;

receiving the seismic data, the seismic data being indicative of a plurality of time-dependent seismic traces of one of a two-dimensional migrated section and a three-dimensional migrated section;

determining second seismic data by transforming each of the plurality of seismic traces into time-frequency domain using a Stockwell transform;

processing the second seismic data for averaging a modulus of the S-transform of seismic traces corresponding to at least a predetermined portion of the one of a two-dimensional migrated section and a three-dimensional migrated section;

determining at each of a plurality of predetermined time instances a maximum amplitude of the averaged modulus of the Stockwell transform and associating the instantaneous amplitude therewith;

determining at each of the plurality of predetermined time instances the dominant instantaneous frequency as the frequency where the maximum amplitude of the averaged modulus of the Stockwell transform occurs;

determining at each of the plurality of predetermined time instances third seismic data indicative of a likelihood of a presence of hydrocarbons based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency; and

,

providing the third seismic data.

16. A storage medium having stored therein executable commands for execution on a processor, the processor when executing the commands performing:

a) receiving the seismic data, the seismic data being indicative of at least a time-dependent seismic trace;

b) determining second seismic data by transforming the seismic data into time-frequency domain using a time-frequency transform;

c) processing the second seismic data for determining at predetermined time instances an instantaneous amplitude and a dominant instantaneous frequency;

d) determining at the predetermined time instances third seismic data indicative of a likelihood of a presence of hydrocarbons based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency; and

,

e) providing the third seismic data.

17. A storage medium as defined in claim 16 having stored therein executable commands for execution on a processor, the processor when executing the commands performing:

determining at each predetermined time instance a maximum amplitude of the modulus of the time-frequency transform and associating the instantaneous amplitude therewith; and

,

determining at each predetermined time instance the dominant instantaneous frequency as the frequency where the maximum amplitude of the modulus of the time-frequency transform occurs.

18. A storage medium as defined in claim 16 having stored therein executable commands for execution on a processor, the processor when executing the commands performing:

processing the second seismic data for averaging a modulus of the time-frequency transform of seismic traces corresponding to at least a predetermined portion of one of a two-dimensional migrated section and a three-dimensional migrated section.

19. A storage medium as defined in claim 16 having stored therein executable commands for execution on a processor, the processor when executing the commands performing:

processing the third seismic data for detecting first occurrences where the quotient of the instantaneous amplitude to the dominant instantaneous frequency is above a predetermined first threshold;

processing the third seismic data of seismic traces where first occurrences have been detected for detecting second occurrences where the quotient of the instantaneous amplitude to the dominant instantaneous frequency is above a predetermined second threshold;

determining fourth seismic data indicative of seismic traces comprising the first and the second occurrences; and

,

providing the fourth seismic data.

20. A storage medium having stored therein executable commands for execution on a processor, the processor when executing the commands performing:

receiving the seismic data, the seismic data being indicative of a plurality of time-dependent seismic traces of one of a two-dimensional migrated section and a three-dimensional migrated section;

determining second seismic data by transforming each of the plurality of seismic traces into time-frequency domain using a Stockwell transform;

processing the second seismic data for averaging a modulus of the Stockwell transform of seismic traces corresponding to at least a predetermined portion of the one of a two-dimensional migrated section and a three-dimensional migrated section;

determining at each of a plurality of predetermined time instances a maximum amplitude of the averaged modulus of the Stockwell transform and associating the instantaneous amplitude therewith;

determining at each of the plurality of predetermined time instances the dominant instantaneous frequency as the frequency where the maximum amplitude of the averaged modulus of the Stockwell transform occurs;

determining at each of the plurality of predetermined time instances third seismic data indicative of a likelihood of a presence of hydrocarbons based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency; and

,

providing the third seismic data.

21. A system for processing seismic data comprising:

an input port for receiving the seismic data, the seismic data being indicative of at least a time-dependent seismic trace;

a processor in communication with the input port for;

a) determining second seismic data by transforming the seismic data into time-frequency domain using a time-frequency transform;

b) processing the second seismic data for determining at predetermined time instances an instantaneous amplitude and a dominant instantaneous frequency;

c) determining at the predetermined time instances third seismic data indicative of a likelihood of a presence of hydrocarbons based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency; and

,

a storage medium in communication with the processor having stored therein executable commands for execution on the processor, the processor when executing the commands performing at least a portion of a) to c); and

,

an output port in communication with the processor for providing the third seismic data.

22. A system for processing seismic data as defined in claim 21 wherein the processor comprises electronic circuitry designed for performing at least a portion of a) to c).

23. A system for processing seismic data as defined in claim 21 comprising a control port in communication with the processor for receiving control commands for controlling at least a portion of a) to c).

24. A system for processing seismic data as defined in claim 23 comprising a graphical display in communication with the processor for displaying at least one of the seismic data, the second seismic data and the third seismic data in a human comprehensible fashion.

25. A system for processing seismic data as defined in claim 24 wherein the graphical display comprises a graphical user interface.

26. A system for processing seismic data comprising:

an input port for receiving the seismic data, the seismic data being indicative of a plurality of time-dependent seismic traces of one of a two-dimensional migrated section and a three-dimensional migrated section;

a processor in communication with the input port for;

a) determining second seismic data by transforming each of the plurality of seismic traces into time-frequency domain using a stockwell-transform;

b) processing the second seismic data for averaging a modulus of the S-transform of seismic traces corresponding to at least a predetermined portion of the one of a two-dimensional migrated section and a three-dimensional migrated section;

c) determining at each of a plurality of predetermined time instances a maximum amplitude of the averaged modulus of the stockwell-transform and associating the instantaneous amplitude therewith;

d) determining at each of the plurality of predetermined time instances the dominant instantaneous frequency as the frequency where the maximum amplitude of the averaged modulus of the stockwell-transform occurs;

e) determining at each of the plurality of predetermined time instances third seismic data indicative of a likelihood of a presence of hydrocarbons based on a quotient of the instantaneous amplitude to the dominant instantaneous frequency; and

,

a storage medium in communication with the processor having stored therein executable commands for execution on the processor, the processor when executing the commands performing at least a portion of a) to e); and

,

an output port in communication with the processor for providing the third seismic data.