System and method for 3D frequency domain waveform inversion based on 3D time-domain forward modeling

• US 7,725,266 B2
• Filed: 05/31/2007
• Issued: 05/25/2010
• Est. Priority Date: 05/31/2006
• Status: Active Grant

• Current Assignee: BP Corporation North America Incorporated (BP PLC)
• Original Assignee: BP Corporation North America Incorporated (BP PLC)
• Inventors: Sirgue, Laurent, Etgen, John T., Albertin, Uwe, Brandsberg-Dahl, Sverre
• Primary Examiner(s): Cosimano; Edward R
• Application Number:
  US11/756,384
• Publication Number:
  US 20070282535A1
• Time in Patent Office: 1,090 Days
• Field of Search: 731/520.1, 731/520.2, 735/70, 735/84, 735/87, 735/94, 166/244.1, 166/250.01, 166/254.1, 175/40, 175/50, 181/101, 181/108, 324/323, 367/14, 367/15, 367/25, 367/36, 367/37, 367/38, 367/73, 702/1, 702/2, 702/6, 702/11, 702/14, 702/16, 702/17, 702/18, 702/127, 702/187, 702/189, 702/190, 702/191, 703/6, 703/9, 703/10
• US Class Current: 702/17
• 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 2210/67 Wave propagation modeling
• CPC Class Codes: G01V 3/083 Controlled source electromagnetic [CSEM] surveying

Abstract

According to a first preferred aspect of the instant invention, there is provided an efficient method of computing a 3D frequency domain waveform inversion based on 3D time domain modeling. In the preferred arrangement, 3D frequency domain wavefields are computed using 3D time-domain modeling and a discrete Fourier transformation that is preferably computed “on the fly” instead of solving the large systems of linear equations that have traditionally been required by direct frequency domain modeling. The instant invention makes use of the theory of gradient-based waveform inversion that estimates model parameters (for example velocities) by matching modeled data to field data sets. Preferably the modeled data are calculated using a forward modeling algorithm.

Claims

1. A method of geophysical exploration within a predetermined volume of the earth containing subsurface structural and stratigraphic features conducive to the generation, migration, accumulation, or presence of hydrocarbons, comprising the steps of:

a. selecting a 3D seismic data set, wherein at least a portion of said 3D seismic data set images said at least a portion of said predetermined volume of the earth, and wherein said 3D seismic data set comprises a plurality of seismic traces;

b. computing a partial Fourier transform of at least a portion of said seismic traces in said 3D seismic data set, thereby producing a transformed 3D seismic data set;

c. choosing an initial subsurface model;

d. using at least said initial subsurface model to calculate a 3D frequency domain forward propagated wave field by time domain forward modeling;

e. calculating a residual volume between said 3D frequency domain wave field and said transformed 3D seismic data set;

f. calculating a 3D frequency domain back propagated wave field by 3D time domain forward modeling using at least said back propagated wave field;

g. using at least said 3D frequency domain back propagated wave field and said 3D frequency domain forward propagated wave field to update said initial subsurface model, thereby producing an updated subsurface model; and

,

h. using said at least a portion of said updated subsurface model to explore for said subsurface structural and stratigraphic features conducive to the generation, migration, accumulation or presence of hydrocarbons within said predetermined volume of the earth.

2. A method of geophysical exploration according to claim 1, wherein step (h) comprises the steps of:

(h1) performing at least steps (d) through (g) with said initial subsurface model replaced by said updated subsurface model, thereby producing a further updated subsurface model, and,

(h2) using said at least a portion of said further updated subsurface model to explore for said subsurface structural and stratigraphic features conducive to the generation, migration, accumulation or presence of hydrocarbons within said predetermined volume of the earth.

3. A method of geophysical exploration according to claim 1, wherein step (g) comprises the steps of:

(g1) using at least said 3D frequency domain back propagated wave field and said 3D frequency domain forward propagated wave field to calculate a gradient volume,

(g2) calculating a step length, and,

(g3) using at least said gradient volume and said step length to update said initial subsurface model, thereby producing an updated subsurface model.

4. A method of geophysical exploration according to claim 1, wherein step (h) comprises the steps of

(h1) storing at least a portion of said updated subsurface model on a computer readable medium;

and,

(h2) using said stored at least a portion of said updated subsurface model to explore for said subsurface structural and stratigraphic features conducive to the generation, migration, accumulation or presence of hydrocarbons within said predetermined volume of the earth.

5. A method according to claim 4, wherein said computer readable medium is selected from the group consisting of computer RAM, computer ROM, a PROM chip, flash RAM, a ROM card, a RAM card, a floppy disk, a magnetic disk, a magnetic tape, a magneto-optical disk, an optical disk, a CD-ROM disk, or a DVD disk.

6. A method according to claim 1, wherein step (h) comprises the steps of:

(h1) displaying at least a portion of said updated subsurface model on a generally flat medium; and

,

(h2) using at least said displayed portion of said updated subsurface model to explore for said subsurface structural and stratigraphic features conducive to the generation, migration, accumulation or presence of hydrocarbons.

7. A method according to claim 6, wherein said generally flat medium is selected from a group consisting of a computer monitor and a piece of paper.

8. A method of geophysical exploration for natural resources within a predetermined volume of the earth, comprising the steps of:

a. selecting a 3D seismic data set, wherein at least a portion of said 3D seismic data set images said at least a portion of said predetermined volume of the earth, and wherein said 3D seismic data set comprises a plurality of seismic traces;

b. computing a partial Fourier transform of at least a portion of said seismic traces in said 3D seismic data set, thereby producing a transformed 3D seismic data set;

c. choosing a current subsurface model;

d. using at least said current subsurface model to calculate a 3D frequency domain forward propagated wave field by time domain forward modeling;

e. calculating a residual volume between said 3D frequency domain wave field and said transformed 3D seismic data set;

f. calculating a 3D frequency domain back propagated wave field by 3D time domain forward modeling using at least said back propagated wave field;

g. using at least said 3D frequency domain back propagated wave field and said 3D frequency domain forward propagated wave field to update said current subsurface model, thereby producing an updated subsurface model;

h. choosing said current subsurface model to be said updated subsurface model;

i. performing steps (d) through (g), thereby producing a further updated subsurface model; and

,

j. using said at least a portion of said further updated subsurface model to explore for said natural resources within said predetermined volume of the earth.

9. A method of geophysical exploration according to claim 8, wherein step (b) comprises the steps of:

(b1) computing a full Fourier transform of at least a portion of said seismic traces in said 3D seismic data set, thereby producing a transformed 3D seismic data set.

10. A method of geophysical exploration according to claim 8, wherein step (g) comprises the steps of:

(g1) using at least said 3D frequency domain back propagated wave field and said 3D frequency domain forward propagated wave field to calculate a gradient volume,

(g2) calculating a step length, and,

(g3) using at least said gradient volume and said step length to update said current subsurface model, thereby producing an updated subsurface model.

11. A method of geophysical exploration according to claim 8, wherein step (h) comprises the steps of

(h1) storing at least a portion of said updated subsurface model on a computer readable medium;

and,

(h2) using said stored at least a portion of said updated subsurface model to explore for said subsurface structural and stratigraphic features conducive to the generation, migration, accumulation or presence of hydrocarbons within said predetermined volume of the earth.

12. A method according to claim 11, wherein said computer readable medium is selected from the group consisting of computer RAM, computer ROM, a PROM chip, flash RAM, a ROM card, a RAM card, a floppy disk, a magnetic disk, a magnetic tape, a magneto-optical disk, an optical disk, a CD-ROM disk, or a DVD disk.

13. A method of geophysical exploration for natural resources within a predetermined volume of the earth, comprising the steps of:

a. selecting a 3D seismic data set, wherein at least a portion of said 3D seismic data set images said at least a portion of said predetermined volume of the earth, and wherein said 3D seismic data set comprises a plurality of seismic traces;

b. selecting at least one frequency;

c. computing a partial Fourier transform of at least a portion of said seismic traces in said 3D seismic data set with respect to said at least one frequency, thereby producing a transformed 3D seismic data set;

d. choosing a current subsurface model;

e. using at least said current subsurface model to calculate a 3D frequency domain forward propagated wave field by time domain forward modeling;

f. calculating a residual volume between said 3D frequency domain wave field and said transformed 3D seismic data set;

g. calculating a 3D frequency domain back propagated wave field by 3D time domain forward modeling using at least said back propagated wave field;

h. using at least said 3D frequency domain back propagated wave field and said 3D frequency domain forward propagated wave field to update said current subsurface model, thereby producing an updated subsurface model;

i. performing at least steps (d) through (g) using said updated subsurface model in place of said current subsurface model, thereby producing a further updated subsurface model; and

j. using said at least a portion of said further updated subsurface model to explore for said natural resources within said predetermined volume of the earth.

14. A method of geophysical exploration according to claim 13, wherein step (h) comprises the steps of:

(h1) using at least said 3D frequency domain back propagated wave field and said 3D frequency domain forward propagated wave field to calculate a gradient volume,

(h2) calculating a step length, and,

(h3) using at least said gradient volume and said step length to update said current subsurface model, thereby producing an updated subsurface model.

15. A method of geophysical exploration according to claim 13, wherein step (j) comprises the steps of

(j1) storing at least a portion of said further updated subsurface model on a computer readable medium;

and,

(j2) using said stored at least a portion of said further updated subsurface model to explore for said natural resources within said predetermined volume of the earth.

16. A method according to claim 15, wherein said computer readable medium is selected from the group consisting of computer RAM, computer ROM, a PROM chip, flash RAM, a ROM card, a RAM card, a floppy disk, a magnetic disk, a magnetic tape, a magneto-optical disk, an optical disk, a CD-ROM disk, or a DVD disk.