Method of using cascaded sweeps for source coding and harmonic cancellation

• US 6,865,488 B2
• Filed: 01/16/2004
• Issued: 03/08/2005
• Est. Priority Date: 10/17/2000
• Status: Expired due to Fees

• Current Assignee: WesternGeco LLC (Schlumberger NV)
• Original Assignee: WesternGeco LLC (Schlumberger NV)
• Inventors: Sitton, Gary, Nyland, David Lee, Moerig, Rainer, Barr, Frederick J.
• Primary Examiner(s): Barlow, John
• Assistant Examiner(s): LAU, TUNG S
• Application Number:
  US10/758,778
• Publication Number:
  US 20040148104A1
• Time in Patent Office: 417 Days
• Field of Search: 702/14, 702/2, 702/18, 702/17, 367/75, 367/40, 367/48, 367/76, 367/56, 181/111, 736/49
• US Class Current: 702/14
• CPC Class Codes: G01V 1/005 with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
• CPC Class Codes: G01V 1/375 Correlating received seismic signals with the emitted source signal
• CPC Class Codes: G01V 2210/3246 Coherent noise, e.g. spatially coherent or predictable

Abstract

A method of seismic surveying using one or more vibrational seismic energy sources activated by sweep signals. The highest order harmonic that has sufficient strength to cause significant harmonic distortion of a sweep segment is determined. A number of sweep segments in excess of the number of sources is selected. Initial phase angles are selected for each sweep segment of each seismic energy source so that substantially all harmonics up to and including the highest order harmonic are suppressed. Using cascaded sweeps, seismic data are acquired and processed to substantially attenuate harmonics upto the selected order.

Claims

1. A method for recording and processing vibratory source seismic data, the method comprising:

(a) generating a cascaded sweep sequence comprising N sweep segments that are either concatenated or overlapping sequentially, where N is equal to or greater than 2, said N sweep segments being substantially identical, except that the initial phase angles of said N sweep segments are progressively rotated by a constant phase increment of about 2 m π

/N radians where m is an integer and |m|>

1;

(b) using said cascaded sweep sequence to drive at least one vibratory source thereby propagating a seismic wave into the earth at a selected location;

(c) recording a groundforce signal associated with said seismic wave;

(d) recording at least one reflection signal from a location within the earth responsive to said seismic wave; and

(e) using said recorded groundforce signal for processing the signal recorded to produce a processed signal.

2. The method of claim 1 wherein processing the signal recorded comprises cross-correlating the recorded signal with a signal comprising said ground force signal and a first and a last sweep segment of said ground force signal.

3. The method of claim 1, wherein said sweep segments further comprise a sinusoidal wavetrain having a frequency that either increases monotonically with time or decreases monotonically with time.

4. The method of claim 3 wherein said increase or decrease of said frequency is linear with time.

5. The method of claim 1, wherein said sweep segments further comprise a psuedo-random sweep series.

6. The method of claim 1 wherein said at least one vibratory source comprises a plurality of vibratory sources, the method further comprising:

i) selecting a number, n, of said vibratory sources activated by n said cascaded sweep sequence signals;

ii) selecting a number of a highest order harmonic that has sufficient strength to cause significant harmonic distortion of a sweep segment that is part of a sweep sequence;

iii) defining a reference signal associated with each of said sources, each reference signal comprising a number N of sweep segments, N being greater than n, said sweep segments being either concatenated or overlapping; and

iv) selecting initial phase angles for each sweep segment of each reference signal so that substantially all harmonics up to and including said highest order harmonic are suppressed and driving each vibratory source with its sweep sequence and driving each of said sources with its associated reference signal.

7. The method of claim 6 wherein selecting said initial phase angles for generating a sweep sequence associated with each energy source further comprises:

(A) constructing a table with N rows and N columns, where N is the number of sweep segments, each entry e_hj being calculated from the formula e_hj=hj modulo N, where h is the row number and j is the column number;

(B) selecting a number of columns equal to the number of sources such that the number in the first row of each selected column does not appear again in any selected column before the row H=F+1, where F is the highest order harmonic determined to cause significant crossfeed; and

(C) assigning a selected column to each seismic energy source, the initial phase angle for the h-th sweep of that seismic energy source being represented by the number in the h-th row of that column multiplied by 2π

/N.

8. A method for recording and processing vibratory source seismic data, the method comprising:

(a) simultaneously transmitting, with n vibratory sources where n≧

1, associated ground force signals into the earth, each said ground force signal including a fundamental signal and harmonics thereof, said fundamental signal having an initial phase;

(b) recording said n ground force signals;

(c) recording at least one reflection signal from a location within the earth responsive to n ground force signals to give a first recorded signal;

(d) repeating steps (a)-(c) m−

1 times wherein said initial phases are shifted by 2π

/m radians to give a total of m recorded signals and n×

m ground force signals; and

(e) using the n×

m ground force signals for processing the m recorded signals to give processed signals associated with said n vibratory sources.

9. The method of claim 8 wherein processing the m recorded signals further comprises:

i) correlating the first recorded signal with each associated ground force signal to give n intermediate signals;

ii) repeating step (i) m−

1 times to give a total of n×

m intermediate processed signals; and

iii) summing the m intermediate processed signals associated with each vibratory source to give a processed signal associated with each vibratory source.

10. The method of claim 8, wherein said sweep segments further comprise a pseudo-random sweep series.

11. The method of claim 9, wherein said sweep segments further comprise a sinusoidal wave train having a frequency that either increases monotonically with time or decreases monotonically with time.

12. The method of claim 10, wherein said pseudo-random sweep series are phase-rotated by selected increments.

13. The method of claim 11 wherein said increase or decrease of said frequency is linear with time.

14. A method for recording and processing vibratory source seismic data, the method comprising:

(a) selecting a number of seismic energy sources to be used for surveying;

(b) selecting a number of a highest order harmonic that has sufficient strength to cause significant harmonic distortion of a sweep segment that is part of a sweep sequence;

(c) defining a reference signal associated with each of said sources, each reference signal comprising a number N of sweep segments, N being greater than the number of seismic sources, said sweep segments being either concatenated or overlapping;

(d) selecting initial phase angles for each sweep segment of each reference signal so that substantially all harmonics up to and including said highest order harmonic are suppressed and driving each vibratory source with its sweep sequence and driving each of said sources with its associated reference signal, thereby generating seismic waves propagating into the earth;

(e) recording ground force signals associated with each reference signal;

(f) recording at least one reflection signal from a location within the earth responsive to said seismic waves; and

(g) using said ground force signals for processing said at least one recorded reflection signal.

15. The method of claim 14 wherein processing said at least one recorded signal comprises correlating the at least one recorded signal with a plurality of processing signals, each said processing signal comprising said reference signal combined with at least one additional segment of said reference signal.

16. The method of claim 14 wherein selecting initial phase angles further comprises:

i) constructing a table with N rows and N columns, each entry e_hm being calculated from the formula e_hm−

hm modulo (N) where h is the row number and m is the column number;

ii) selecting a number of columns equal to the number of sources such that none of the numbers of each selected column appears again in any other selected column before the row H=F+1; and

iii) assigning a selected column to each seismic energy source, the initial phase angle for the h-th sweep segment of that seismic energy source being represented by the number in the h-th row of that column multiplied by 2π

/N.

17. The method of claim 14, wherein said sweep segments further comprise a sinusoidal wavetrain having a frequency that either increases monotonically with time or decreases monotonically with time.

18. The method of claim 17 wherein said increase or decrease of said frequency is linear with time.

19. The method of claim 14, wherein said sweep segments further comprise a pseudo-random sweep series.

20. The method of claim 19, wherein said pseudo-random sweep series are phase-rotated by selected increments.