METHOD FOR SEPARATING INDEPENDENT SIMULTANEOUS SOURCES

• US 20100039894A1
• Filed: 08/17/2009
• Published: 02/18/2010
• Est. Priority Date: 08/15/2008
• Status: Active Grant

• Current Assignee: BP Corporation North America Incorporated (BP PLC)
• Original Assignee: BP Corporation North America Incorporated (BP PLC)
• Inventors: Abma, Raymond L.
• Granted Patent: US 8,295,124 B2
• Time in Patent Office: Days
• US Class Current: 367/52
• 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/362 Effecting static or dynamic corrections; Stacking

Abstract

This is a method of separating simultaneous sources that uses an inversion-type approach. Each source will preferably activated at a random time with respect to the others. These random delays tend to make the interference between sources incoherent while the reflections create coherent events within a series of shots. The shot separation is performed via a numerical inversion process that utilizes the sweeps for each shot, the start times of each shot, and the coherence of reflection events between nearby shots. This method will allow seismic surveys to be acquired faster and cheaper.

Claims

1. A method of seismic exploration above a region of the subsurface containing structural or stratigraphic features conducive to the presence, migration, or accumulation of hydrocarbons, wherein is provided a blended source seismic survey containing at least two interfering seismic source excitations therein, comprising the steps of:

(a) within a computer, accessing at least a portion of said blended seismic source survey;

(b) selecting one of said at least two interfering seismic source excitations;

(c) separating said selected source excitation from said at least two interfering seismic source excitations using sparse inversion, thereby producing a single source gather; and

(d) using said single source gather to explore for hydrocarbons within said region of the subsurface.

2. The method according to claim 1, wherein step (c) comprises the steps of:

(c1) selecting a threshold value,

(c2) accessing a plurality of seismic traces that contain reflection signals from said selected source excitation, and

(c3) using at least said selected threshold value to apply a coherency constraint to said accessed plurality of seismic traces, thereby producing a single source gather.

3. The method according to claim 2, wherein step (c3) comprises the steps of:

(1) determining an activation time for said selected seismic source excitation,

(2) using at least said activation time to shift said accessed plurality of seismic traces to a zero time,

(3) Fourier transforming said time shifted accessed plurality of seismic traces,

(4) thresholding said Fourier transformed time shifted plurality of seismic traces,

(5) inverse Fourier transforming said thresholded Fourier transformed time shifted plurality of seismic traces, thereby producing a single source gather.

4. The method of claim 1, wherein step (d) comprises the steps of:

(d1) performing at least steps (b) and (c) a plurality of times for a plurality of different source excitations, thereby producing a plurality of single source gathers, and,

(d2) using said plurality of single source gathers to explore for hydrocarbons within said region of the subsurface.

5. The method of claim 1, wherein step (c) comprises the steps of:

(c1) separating said selected source excitation from said at least two interfering seismic source excitations by solving the equation;

d=Γ

S m,

for the matrix m thereby producing a single source gather, where

m is said single source gather,

d is a matrix representation of the recorded seismic data,

S is a matrix or operator that describes the similarity between nearby shots, and

Γ

is a matrix that defines the blending or mixing of the individual sources.

6. The method of claim 1, wherein step (c) comprises the steps of:

(c1) separating said selected source excitation from said at least two interfering seismic source excitations by solving the equation;

W d=WΓ

S m,

for the matrix m thereby producing a single source gather, where

m is said single source gather,

d is a matrix representation of the recorded seismic data,

S is a matrix or operator that describes the similarity between nearby shots,

Γ

is a matrix that defines the blending or mixing of the individual sources, and,

W is a known matrix of weights.

7. The method of claim 1, wherein each of said at least two interfering seismic source activations are separated in time by a random period of time.

8. A method of seismic exploration above a region of the subsurface containing structural or stratigraphic features conducive to the presence, migration, or accumulation of hydrocarbons, wherein is provided a blended source seismic survey containing at least two interfering seismic source excitations from at least one seismic source, comprising the steps of:

(a) within a computer, accessing at least a portion of said blended seismic source survey;

(b) selecting one of said at least one seismic sources;

(c) identifying at least one seismic source excitation from said selected seismic source;

(d) separating at least one identified seismic source excitation from said at least two interfering seismic source excitations using sparse inversion, thereby producing at least one single source gather; and

(e) using said at least one single source gather to explore for hydrocarbons within said region of the subsurface.

9. The method according to claim 8, wherein step (d) comprises the steps of:

(d1) selecting a threshold value,

(d2) accessing a plurality of seismic traces that contain reflection signals from at least one of said at least one identified seismic source excitations, and

(d3) using at least said selected threshold value to apply a coherency constraint to said accessed plurality of seismic traces, thereby producing at least one single source gather.

10. The method according to claim 9, wherein step (d3) comprises the steps of:

(1) determining an activation time for each of said at least one identified seismic source excitations,

(2) accessing a plurality of seismic traces that contain reflection signals from said at least one identified seismic source excitations,

(3) using at least said activation time to shift said accessed plurality of seismic traces that contain reflection signals from said selected source excitation to a zero time,

(4) Fourier transforming said time shifted accessed plurality of seismic traces that contain reflection signals from said selected source excitation,

(5) thresholding said transformed time shifted plurality of seismic traces that contain reflections signals from said selected source excitation,

(6) inverse Fourier transforming said thresholded transformed time shifted plurality of seismic traces that contain reflection signals from said selected source excitation, thereby producing a single source gather.

11. The method of claim 8, wherein step (d) comprises the steps of:

(d1) separating at least one of said identified seismic source excitations from said at least two interfering seismic source excitations by solving the equation;

d=Γ

S m,

for the matrix in thereby producing at least one single source gather, where

m is said at least one single source gather,

d is a matrix representation of the recorded seismic data,

S is a matrix or operator that describes the similarity between nearby shots, and

Γ

is a matrix that defines the blending or mixing of the individual sources.

12. The method of claim 8, wherein step (d) comprises the steps of:

(d1) separating at least one of said identified seismic source excitations from said at least two interfering seismic source excitations by solving the equation;

W d=WΓ

S m,

for the matrix m thereby producing a single source gather, where

m is said single source gather,

d is a matrix representation of the recorded seismic data,

S is a matrix or operator that describes the similarity between nearby shots,

Γ

is a matrix that defines the blending or mixing of the individual sources, and,

W is a known matrix of weights.

13. The method of claim 8, wherein each of said at least two interfering seismic source activations are separated in time by a random time period.