METHOD FOR GENERATION OF IMAGES RELATED TO A SUBSURFACE REGION OF INTEREST

• US 20100118651A1
• Filed: 09/24/2009
• Published: 05/13/2010
• Est. Priority Date: 11/10/2008
• Status: Abandoned Application

• Current Assignee: Chevron USA Incorporated (Chevron Corp.)
• Original Assignee: Chevron USA Incorporated (Chevron Corp.)
• Inventors: Ergas, Raymond, Liu, Wei, Nemeth, Tamas, Pell, Oliver
• Application Number:
  US12/566,252
• Publication Number:
  US 20100118651A1
• Time in Patent Office: Days
• US Class Current: 367/50
• CPC Class Codes: G01V 1/282 Application of seismic models, synthetic seismograms
• CPC Class Codes: G01V 1/34 Displaying seismic recordings or visualisation of seismic data or attributes
• CPC Class Codes: G01V 2210/614 Synthetically generated data
• CPC Class Codes: G01V 2210/679 Reverse-time modeling or coalescence modelling, i.e. starting from receivers

Abstract

A method and system for generating images of a subsurface region of interest. In general, one embodiment of the present invention includes establishing boundary conditions utilizing seismic data and initial conditions which include excitation from source locations in an earth model. Source wavefields are then propagated forward through the earth model to a maximum time, and saved at a plurality of checkpoints sparsely in time and also corresponding boundary values of the source wavefields at each time step are saved. Source wavefields are also propagated backward through the earth model from the maximum time utilizing the plurality of checkpoints when available and the saved boundary values at each time step. Receiver wavefields are propagated backward concurrently through the earth model from the maximum time. Imaging conditions are applied at selected time steps to both the backward propagated source wavefields and receiver wavefields and those wavefields are utilized to generate images related to the subsurface region.

Claims

1. A computer-implemented method of generating images related to a subsurface region;

the method comprising;

obtaining seismic data and an earth model related to the subsurface region, wherein both the seismic data and the earth model are stored on electronic media;

utilizing at least one processor, configured to communicate with the electronic media and arranged to execute machine executable instructions stored in a processor accessible memory for performing steps comprising;

establishing boundary conditions utilizing the seismic data and initial conditions which include excitation from source locations in the earth model;

propagating forward source wavefields through the earth model to a maximum time;

saving source wavefields at a plurality of checkpoints sparsely in time and saving corresponding boundary values of the source wavefields at each time step;

propagating backward the source wavefields through the earth model from the maximum time, utilizing the plurality of checkpoints when available and the saved boundary values at each time step, and concurrently propagating backward receiver wavefields or seismic data through the earth model from the maximum time; and

applying imaging conditions at selected time steps to both the backward propagated source wavefields and receiver wavefields from seismic data, wherein the backward propagated source wavefields and receiver wavefields are utilized to generate images related to the subsurface region.

2. The method of claim 1 wherein the computer-implemented method is utilized for reverse time migration, waveform inversion, or an application requiring reverse order access of data components in order to operate with other data components.

3. The method of claim 1 wherein the earth model is extended with a slow-velocity cortex or a boundary condition of properties to accommodate source wavefields that propagated outside of the original boundary.

4. The method of claim 1 wherein the wavefields are propagated utilizing a numerical solver.

5. The method of claim 4 wherein the numerical solver includes comprises at least one reverse time migration, Gaussian beam migration, Kirchhoff migration or waveform inversion.

6. The method of claim 4 wherein the numerical solver includes a wave-equation based migration.

7. The method of claim 1 wherein the wavefield propagation is performed in time, frequency or wavelet domains.

8. A system configured to generate images related to a subsurface region of interest, the system comprising:

at least one data storage device having seismic data and an earth model related to the subsurface region of interest;

at least one processor, configured and arranged to execute machine executable instructions stored in a processor accessible memory for performing a method comprising;

establishing boundary conditions utilizing the seismic data and initial conditions which include excitation from source locations in the earth model;

propagating forward source wavefields through the earth model to a maximum time;

utilizing the data storage device to save source wavefields at a plurality of checkpoints sparsely in time and saving corresponding boundary values of the source wavefields at each time step;

propagating backward the source wavefields through the earth model from the maximum time utilizing the plurality of checkpoints when available and the saved boundary values at each time step, and concurrently propagating backward receiver wavefields or seismic data through the earth model from the maximum time; and

applying imaging conditions at selected time steps to both the backward propagated source wavefields and receiver wavefields, wherein the backward propagated source wavefields and receiver wavefields are utilized to generate images related to the subsurface region on a display device.

9. The system of claim 8 wherein the method is utilized for reverse time migration, waveform inversion, or an application requiring reverse order access of data components in order to operate with other data components.

10. The system of claim 8 wherein the earth model is extended with a slow-velocity cortex or a boundary condition of properties to accommodate source wavefields that propagated outside of the original boundary.

11. The system of claim 8 wherein the wavefields are propagated utilizing a numerical solver.

12. The system of claim 11 wherein the numerical solver comprises at least one of reverse time migration, Gaussian beam migration, Kirchhoff migration or waveform inversion.

13. The system of claim 11 wherein the numerical solver includes a wave-equation based migration.

14. The system of claim 8 wherein the wavefield propagation is performed in time, frequency or wavelet domains.