Petroleum exploration and prediction apparatus and method

• US 6,952,649 B2
• Filed: 10/01/2003
• Issued: 10/04/2005
• Est. Priority Date: 10/04/2002
• Status: Expired due to Fees

• Current Assignee: Daniel R. Cook, John D. Walther
• Original Assignee: Daniel R. Cook, John D. Walther
• Inventors: Cook, Daniel R., Walther, John D.
• Primary Examiner(s): MCELHENY JR, DONALD E
• Application Number:
  US10/676,945
• Publication Number:
  US 20040068375A1
• Time in Patent Office: 734 Days
• Field of Search: 702 10- 14, 702/16, 702/17, 702/18, 367/73
• US Class Current: 702/14
• CPC Class Codes: A61B 5/245 specially adapted for magnetoencephalographic [MEG] signals
• CPC Class Codes: A61B 5/316 Modalities, i.e. specific diagnostic methods
• CPC Class Codes: A61B 5/369 Electroencephalography [EEG] devices for psychotechnics A61B5/16
• CPC Class Codes: G01V 1/28 Processing seismic data, e.g. for interpretation or for event detection G01V1/48 takes precedence
• CPC Class Codes: G01V 1/282 Application of seismic models, synthetic seismograms
• CPC Class Codes: G01V 1/30 Analysis G01V1/50 takes precedence
• CPC Class Codes: G01V 2210/51 Migration
• CPC Class Codes: G01V 2210/512 Pre-stack
• CPC Class Codes: G06F 18/40 Software arrangements specially adapted for pattern recognition, e.g. user interfaces or toolboxes therefor
• CPC Class Codes: G06F 2218/12 Classification; Matching

Abstract

A method for predicting the state of a geological formation. The method may include generating a separation key effective to extract a first feature from a signal or signals corresponding to a first state and a second feature, distinct from the first feature, from a signal or signals corresponding to a second state. The separation key may list at least one feature operator and a weighting table. The at least one feature operator may expand a test signal collected from a geological formation of unknown state in at least one of frequency space and time space to generate a plurality of feature segments. A weighting table may weight the plurality of feature segments. The weighted plurality of feature segments may be superimposed to form a third feature. The geological formation may be classified as having one of the first state and second state based on the correspondence of the third feature to one of the first feature and second feature.

Claims

1. A method comprising:

providing a first signal segment corresponding to a first geological formation having a first characteristic;

providing a second signal segment corresponding to a second geological formation having a second characteristic, distinct from the first characteristic;

expanding the first and second signal segments in at least one of frequency space and time space by applying at least one feature operator to the first and second signal segments to generate a plurality of first feature segments corresponding to the first signal segment and a plurality of second feature segments corresponding to the second signal segment;

weighting the pluralities of first and second feature segments by applying a weight table to the pluralities of first and second feature segments to generate a weighted plurality of first feature segments and a weighted plurality of second feature segments, respectively; and

collapsing the weighted plurality of first feature segments to form a first feature comprising a non-random first pattern and the weighted plurality of second feature segments to form a second feature having a second pattern distinct from the non-random first pattern.

2. The method of claim 1, further comprising generating a separation key listing the at least one feature operator used to expand the first and second signal segments and the weighting table applied to the pluralities of first and second feature segments.

3. The method of claim 2, further comprising providing a third signal segment corresponding to a third geological formation.

4. The method of claim 3, further comprising expanding, weighting, and collapsing the third signal segment in accordance with the at least one feature operator and weighting table listed in the separation key to generate a third feature.

5. The method of claim 4, further comprising classifying the third geological formation as having one of the first and second characteristics based on the correspondence of the third feature to one of the first feature and the second feature.

6. The method of claim 5, wherein the first characteristic is selected from the group consisting of porosity above a first threshold value, density above a second threshold value, pressure above a third threshold value, shale in a concentration above a fourth threshold value, hydrocarbon production above a fifth threshold value, temperature above a sixth threshold value, and sand in a concentration above a seventh threshold value.

7. The method of claim 6, wherein the second characteristic is selected from the group consisting of porosity below the first threshold value, density below the second threshold value, pressure below the third threshold value, shale in a concentration below the fourth threshold value, hydrocarbon production below the fifth threshold value, temperature below the sixth threshold value, and sand in a concentration below the seventh threshold value.

8. The method of claim 1, wherein the first signal segment is a portion of a seismic trace.

9. The method of claim 8, wherein a dominant waveform exists in the first signal segment and is an apparently substantially random waveform.

10. The method of claim 9, wherein the dominant waveform of the first signal segment effectively obscures other waveforms.

11. The method of claim 10, wherein weighting comprises applying weights selected based on their effectiveness to neutralize the random contribution of the dominant waveform.

12. The method of claim 1, wherein expanding comprises dividing the first and second signal segments into selected frequency bands.

13. The method of claim 12, wherein a first feature operator of the at least one feature operator provides selected central frequencies and corresponding frequency widths to define the boundaries of the selected frequency bands.

14. The method of claim 12, wherein expanding comprises dividing the first and second signal segments into selected time segments and wherein a second feature operator of the at least one feature operator provides selected central times and corresponding times widths to define the boundaries of the selected time segments.

15. The method of claim 14, wherein collapsing comprises superimposing the weighted plurality of first feature segments to form a first feature comprising a non-random first pattern and the weighted plurality of second feature segments to form a second feature having a second pattern distinct from the non-random first pattern.

16. A method comprising:

providing a first signal corresponding to a first geological formation tapped by a hydrocarbon well producing hydrocarbons above a threshold rate;

providing a second signal corresponding to a second geological formation tapped by a hydrocarbon well producing hydrocarbons below the threshold rate;

expanding the first and second signal segments in at least one of frequency space and time space by applying at least one feature operator to the first and second signal segments to generate a plurality of first feature segments corresponding to the first signal segment and a plurality of second feature segments corresponding to the second signal segment;

weighting the pluralities of first and second feature segments by applying a weight table to the pluralities of first and second feature segments to generate a weighted plurality of first feature segments and a weighted plurality of second feature segments, respectively;

superimposing the weighted plurality of first feature segments to form a first feature comprising a non-random first pattern and the weighted plurality of second feature segments to form a second feature having a second pattern distinct from the non-random first pattern; and

generating a separation key listing the at least one feature operator used to expand the first and second signal segments and the weighting table applied to the pluralities of first and second feature segments.

17. The method of claim 16, wherein the first signal is a pre-stack, migrated seismic trace.

18. The method of claim 17, wherein the second seismic trace is a pre-stack, migrated seismic trace.

19. The method of claim 16, wherein the first seismic trace is a post-stack, migrated seismic trace.

20. The method of claim 19, wherein the second seismic trace is a post-stack, migrated seismic trace.

21. The method of claim 16, wherein the hydrocarbon well is an oil well.

22. The method of claim 16, wherein the hydrocarbon well is a gas well.

23. The method of claim 16, wherein expanding comprises dividing the first and second signal segments into selected frequency bands and wherein a first feature operator of the at least one feature operator provides selected central frequencies and corresponding frequency widths to define the boundaries of the selected frequency bands.

24. The method of claim 23, wherein the first feature operator imposes a Gaussian weighting at centered at each of the selected central frequencies to define the selected frequencies bands.

25. The method of claim 23, wherein expanding comprises dividing the first and second signal segments into selected time segments and wherein a second feature operator of the at least one feature operator provides selected central times and corresponding times widths to define the boundaries of the selected time segments.

26. The method of claim 25, wherein the second feature operator imposes a Gaussian weighting centered at each of the selected central times to define the selected time segments.

27. The method of claim 23, wherein collapsing comprises superimposing the weighted plurality of first feature segments to form a first feature comprising a non-random first pattern and the weighted plurality of second feature segments to form a second feature having a second pattern distinct from the non-random first pattern.

28. The method of claim 27, further comprising providing a third signal corresponding to a third geological formation containing a prospective hydrocarbon well location.

29. The method of claim 28, further comprising processing third signal in accordance with the at least one feature operator and weighting table listed in the separation key to produce a third feature.

30. The method of claim 29, further comprising classifying the prospective hydrocarbon well location as one of hydrocarbon well producing above the threshold rate and a hydrocarbon well producing below the threshold rate based on the correspondence of the third feature to one of the first feature and the second feature.

31. A method for predicting the state of a geological formation, the method comprising:

providing a separation key effective to extract a first feature from a signal corresponding to a first state and a second feature, distinct from the first feature, from a signal corresponding to a second state, the separation key listing at least one feature operator and a weighting table;

providing a first test signal collected from the geological formation;

applying the at least one feature operator to expand the first test signal in at least one of frequency space and time space to generate a plurality of feature segments;

generating a weighted plurality of feature segments by applying the weighting table to the plurality of feature segments;

collapsing the weighted plurality of feature segments to generate a third feature; and

classifying the geological formation as having one of the first state and second state based on the correspondence of the third feature to one of the first feature and second feature.

32. The method of claim 31, further comprising providing a second test signal collected from the geological formation at a time after the first test signal was collected.

33. The method of claim 32, further comprising processing the second test signal by applying the at least one feature operator, applying the weight table, and collapsing to generate a forth feature.

34. The method of claim 33, further comprising quantifying changes in the geological formation between the time the first test signal was collected and the time the second test signal was collected by comparing the fourth feature to the first feature.

35. The method of claim 31, further comprising selecting a color spectrum comprising at least two colors and having a first extreme and a second extreme, opposite the first extreme, and assigning the color of the first extreme to the first feature and the color of the second extreme to the second feature.

36. The method of claim 35, further comprising assigning the third feature a color selected from the color spectrum based on the relative correspondence of the third feature to the first and second features.

37. The method of claim 36, further comprising providing a plurality of test signals, each collected from a different location on the geological formation.

38. The method of claim 37, further comprising processing the plurality of test signals by applying the at least one feature operator, applying the weight table, and collapsing to generate a plurality of features.

39. The method of claim 38, further comprising assigning each of the plurality of features a color selected from the color spectrum based on the relative correspondence thereof to the first and second features.

40. The method of claim 39, further comprising generating a two-dimensional image of the geological formation by coloring portions of the image in accordance with the color assigned to the feature of the plurality of features corresponding to the portions.

41. A method for predicting the presence of extractable hydrocarbons in a geological formation, the method comprising:

providing a separation key effective to extract a first feature from signal collected from a first geological formation tapped by a hydrocarbon well producing above a threshold rate and a second feature from signal collected from a second geological formation tapped by a hydrocarbon well producing below the threshold rate, the separation key listing at least one feature operator and a weighting table;

providing a test signal collected from a third geological formation;

applying the at least one feature operator to expand the test signal in at least one of frequency space and time space to generate a plurality of feature segments;

generating a weighted plurality of feature segments by applying the weighting table to the plurality of feature segments;

superimposing the weighted plurality of feature segments to generate a third feature; and

classifying the third geological formation as one of producing hydrocarbons above the threshold rate and producing hydrocarbons below the threshold rate based on the correspondence of the third feature to one of the first feature and second feature.