SYSTEM AND METHOD FOR SEISMIC DATA INVERSION

• US 20120314538A1
• Filed: 06/08/2011
• Published: 12/13/2012
• Est. Priority Date: 06/08/2011
• Status: Active Grant

• Current Assignee: Chevron USA Incorporated (Chevron Corp.)
• Original Assignee: Chevron USA Incorporated (Chevron Corp.)
• Inventors: Washbourne, John Kenneth, Shah, Nikhil Koolesh, Bube, Kenneth Paul
• Granted Patent: US 9,075,159 B2
• Time in Patent Office: Days
• US Class Current: 367/49
• CPC Class Codes: G01V 1/30 Analysis G01V1/50 takes precedence
• CPC Class Codes: G01V 2210/45 F-x or F-xy domain
• CPC Class Codes: G01V 2210/614 Synthetically generated data

Abstract

A system and computer-implemented method for determining properties of a subsurface region of interest from seismic data is disclosed. In an embodiment, the method includes obtaining seismic data representative of the subsurface region and an initial earth property model for the subsurface region, performing forward modeling using the initial earth property model to create modeled seismic data with similar acquisition specifications as the actual seismic data, transforming the modeled and actual seismic data to a temporal Fourier frequency domain to create frequency domain modeled and actual seismic data wherein the frequency domain modeled and actual seismic data include an amplitude portion and a phase portion, measuring the misfit between the frequency domain modeled seismic data and frequency domain actual seismic data to produce frequency domain residual seismic data, performing phase unwrapping of the phase portion of certain observed frequency components of the frequency domain residual seismic data to create an unwrapped residual phase portion, and inverting the unwrapped residual phase portion to determine desired properties of the subsurface region of interest, wherein the inverting minimizes an objective function defined to measure the misfit. The method may also include phase extrapolation. Additionally, the method may include a second inverting step using the result of the first inverting step as a starting model. The system includes a data source, a user interface, and a processor configured to execute computer modules designed to perform the method.

Claims

1. A computer-implemented method for determining properties of a subsurface region of interest from actual seismic data comprising:

a. obtaining actual seismic data representative of the subsurface region and an initial earth property model for the subsurface region;

b. performing forward modeling, using the initial earth property model, to create modeled seismic data with similar acquisition specifications as the actual seismic data;

c. transforming the modeled and actual seismic data to a temporal Fourier frequency domain to create frequency domain modeled and actual seismic data wherein the frequency domain modeled and actual seismic data includes an amplitude portion and a phase portion;

d. measuring the misfit between the frequency domain modeled seismic data and frequency domain actual seismic data to produce frequency domain residual seismic data;

e. performing phase unwrapping of the phase portion of certain observed frequency components of the frequency domain residual seismic data to create an unwrapped residual phase portion; and

f. inverting the unwrapped residual phase portion to determine desired properties of the subsurface region of interest, wherein the inverting minimizes an objective function defined to measure the misfit, and wherein the performing forward modeling, transforming, measuring, performing phase unwrapping, and inverting steps are performed by a computer processor.

2. The method of claim 1 further comprising extrapolating the unwrapped phase prior to the inverting step.

3. The method of claim 1 further comprising a second inverting step, the second inverting step inverting the amplitude and phase portions of certain observed frequency components of the frequency domain residual seismic data, wherein the second inverting step uses, as a new initial model, the desired properties determined by the inverting the unwrapped residual phase portion.

4. The method of claim 1 wherein the inverting uses a continuation approach to regularize successive iterations.

5. The method of claim 1 wherein the inverting solves a linear inversion problem for each successive linearization.

6. The method of claim 5 wherein the linear inversion problem is solved using a Gauss-Newton method.

7. A system for determining properties of a subsurface region of interest from seismic data comprising:

a. a data source containing computer-readable data including actual seismic data representative of the subsurface region of interest;

b. a processor configured to execute computer-readable code from computer modules, the computer modules comprising;

i. an initial model module configured to receive an initial property model from the data source or generate the initial property model;

ii. a domain transformation module to transform the actual seismic data into a frequency domain to generate frequency domain seismic data;

iii. a data modeling module to generate modeled seismic data from the initial property model or an updated property model;

iv. a residual calculation module to calculate a residual wavefield;

v. a phase preparation module to phase unwrap the frequency domain seismic data;

vi. an inversion module to compute an objective function, a search direction, and update the initial property model or a previously updated property model.

c. an user interface.

8. The system of claim 7, wherein the data source also contains the initial property model.

9. The system of claim 7, wherein the processor repeatedly executes the data modeling module, the residual calculation module, the phase preparation module, and the inversion module to minimize the objective function.

10. The system of claim 7 wherein the computer modules further comprise a phase extrapolation module.

11. The system of claim 7 wherein the inversion module performs a phase-only full waveform inversion until a phase-only objective function has been minimized or a user-specified number of iterations is reached to generate an updated phase-only inverted property model.

12. The system of claim 11 wherein after the updated phase-only inverted property model is generated, the inversion model performs a conventional full waveform inversion using the updated phase-only inverted property model as the initial property model.

13. The system of claim 7 wherein the inversion module uses a continuation approach to regularize successive iterations.

14. The system of claim 7 wherein the inversion module solves a linear inversion problem for each successive linearization.

15. The system of claim 14 wherein the linear inversion problem is solved using a Gauss-Newton method.

16. An article of manufacture comprising a computer readable medium having a computer readable code embodied therein, the computer readable program code adapted to be executed to implement a method for determining properties of a subsurface region of interest from actual seismic data, the method comprising:

a. performing forward modeling, via a computer, using an initial earth property model, to create modeled seismic data with the same or similar acquisition specifications as the actual seismic data;

b. transforming, via a computer, the modeled and the recorded seismic data to a temporal Fourier frequency domain to create frequency domain modeled and frequency domain recorded seismic data wherein the frequency domain modeled seismic data and frequency domain recorded seismic data include an amplitude portion and a phase portion;

c. measuring, via a computer, the misfit between the frequency domain modeled seismic data and frequency domain recorded seismic data to produce frequency domain residual seismic data;

d. performing, via a computer, phase unwrapping of the phase portion of certain observed frequency components of the frequency domain residual seismic data to create an unwrapped residual phase portion;

e. inverting, via a computer, the unwrapped residual phase portion to determine desired properties of the subsurface region of interest, wherein the inverting minimizes an objective function defined to measure the misfit.