Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data

• US 7,373,251 B2
• Filed: 12/22/2004
• Issued: 05/13/2008
• Est. Priority Date: 12/22/2004
• Status: Expired due to Fees

• Current Assignee: Cgg Services Sa (CGG)
• Original Assignee: Compagnie Generale De Geophysique, Marathon Oil Company (Marathon Oil Corporation)
• Inventors: Crozat, Thierry, Lafet, Yves, Deschizeaux, Bernard, Allo, Fabien, Molle, Amelie Rodrigue, Lanfranchi, Pierre, Bornard, Raphael, Hamman, Jeffry G., Caldwell, Donald H., Coleou, Thierry
• Primary Examiner(s): Barlow; John E.
• Assistant Examiner(s): Le; Toan M.
• Application Number:
  US11/022,597
• Publication Number:
  US 20060136162A1
• Time in Patent Office: 1,238 Days
• Field of Search: 702/14, 702/11, 702/12, 702/13, 367/18, 367/73, 703/9, 703/10, 175/50
• US Class Current: 702/14
• CPC Class Codes: G01V 1/306 for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles

Abstract

A method is provided for predicting a value of a designated rock or fluid property in a subterranean geologic volume. A first predicted value of the designated rock or fluid property is also assigned to a volume of a multi-dimensional, multi-scale model. A first predicted value of seismic response for the model volume is calculated from a response model using the first predicted value of the designated rock or fluid property, wherein the response model is responsive to changes in predicted values of the designated rock or fluid property. A synthetic trace is generated and iteratively compared to the corresponding trace obtained from one or more sets of actual seismic data to determine a difference while consistency is maintained between the types, scales and dimensions of values and data. The first predicted value of the designated rock or fluid property is adjusted in response to the difference to create a second predicted value of the designated rock or fluid property, wherein the second predicted value generally reduces the difference.

Claims

1. A method for determining a value of a designated rock or fluid property in a subterranean geologic volume, the method comprising:

assigning a first predicted value of a rock or fluid property to a cell within a multi-dimensional layered model that is correlated to a subterranean geologic volume;

calculating first predicted values of seismic response for said cell from a response model using said first predicted value of said rock or fluid property and based upon compressional velocity, shear velocity and bulk density, wherein said response model is responsive to changes in predicted values of said rock or fluid property;

comparing synthetic traces that are generated from said first predicted values of seismic response to the corresponding traces obtained from sets of actual seismic data corresponding to the subterranean geologic volume to ascertain a difference between said synthetic traces and said corresponding traces, said step of comparing involving maintaining consistency between the types, scales and dimensions of said first predicted value of a designated rock or fluid property and said sets of actual seismic data;

adjusting said first predicted value of said designated rock or fluid property in response to said difference to create a second predicted value of said designated rock or fluid property, wherein said second predicted value generally reduces said difference; and

storing said second predicted value for later use.

2. The method of claim 1 wherein said multi-dimensional layered model has X, Y and multiple Z axes to store and manipulate data obtained from different domains, the process further comprising:

assigning first data on one of the Z axis based in a given domain to a volume within the multi-dimensional layered model;

assigning second data on another of the Z axis based in one or more other domains to the volume within the multi-dimensional layered model;

changing a variable in one of the first data or the second data; and

updating the other of the first data or the second data in response to a rule of change so as to provide consistency between the first data and the second data.

3. A method for determining an error-minimized value of a designated rock or fluid property at a location in a subterranean geologic volume which is characterized by a multi-dimensional layered model volume having a plurality of model subvolumes, said method comprising:

assigning a first predicted value of a designated rock or fluid property to each of said plurality of model subvolumes;

calculating first predicted values of a seismic response for the interface between at least two groups of subvolumes based upon compressional velocity, shear velocity and bulk density;

comparing synthetic traces that are generated from said first predicted values of seismic response to the corresponding traces obtained from a plurality of response volumes of actual seismic data corresponding to the subterranean geologic volume to determine a first difference between said synthetic traces and said corresponding traces, said step of comparing involves maintaining consistency between the types, scales and dimensions of said first predicted value of a designated rock or fluid property and said plurality of response volumes of actual seismic data;

adjusting said first predicted value of said designated rock or fluid property in response to said difference to create a second predicted value of said designated rock or fluid property, wherein said second predicted value generally reduces said first difference; and

storing said second predicted value for later use.

4. The method of claim 3 further comprising calculating a second predicted values of seismic response for said interface between said at least two groups using said second predicted value of said designated rock or fluid property.

5. The method of claim 4 further comprising comparing synthetic traces that are generated from said second predicted values of seismic response to the corresponding traces obtained from said plurality or response volumes of actual seismic data corresponding to the subterranean geologic volume to determine a second difference, wherein said second difference is generally less than said first difference.

6. The method of claim 3 further comprising iteratively repeating said steps of adjusting said predicted value of said designated rock or fluid property, calculating said predicted values of said seismic response using said adjusted predicted value of said designated rock or fluid property, and comparing said synthetic traces that are generated from said calculated predicted values of seismic response to the corresponding traces obtained from said plurality of response volumes of actual seismic data corresponding the subterranean geologic volume to determine said difference between said synthetic traces and said corresponding traces.

7. The method of claim 6 wherein said iterative repetition is terminated when said difference does not exceed a predetermined tolerance or other suitable criteria.

8. The method of claim 7 wherein said predicted value of said designated rock or fluid property upon termination of said iterative repetition is an error minimized value of said designated rock or fluid property.

9. The method of claim 3 wherein said response model has a plurality of rock and fluid properties as model parameters including said designated rock or fluid property.

10. The method of claim 9 wherein said plurality of rock and fluid properties include elastic rock and fluid moduli.

11. The method of claim 10 further comprising assigning estimated values of said plurality of rock and fluid properties to said model subvolume in addition to said first predicted value of said designated rock or fluid property.

12. The method of claim 3 further comprising providing consistency of the said model subvolume rock and fluid properties with constraints and/or relationships at a location in said geologic volume.

13. The method of claim 3 wherein said designated rock or fluid property is selected from the group consisting of mineral or the matrix modulus, water modulus, oil modulus, gas modulus, frame modulus, frame shear modulus, mineral density, water density, oil density, gas density, porosity, water saturation, oil saturation, and gas saturation.

14. The method of claim 3 wherein said first predicted values of seismic response for said cell are simultaneously calculated from said response model.

15. The method of claim 14 wherein said synthetic traces that are generated from said first predicted values of seismic response are simultaneously compared to the corresponding traces obtained from sets of actual seismic data corresponding to the subterranean geologic volume.

16. The method of claim 3 wherein said plurality of response volumes correspond to different values of angle of incidence.

17. The method of claim 3 wherein said plurality of response volumes correspond to different time of acquisition.

18. The method of claim 3 wherein said plurality of response volumes correspond to different shooting direction.

19. The method of claim 3 wherein said plurality of response volumes correspond to different mode of propagation.

20. A method of characterizing a subterranean geologic volume, the method comprising:

assigning a first predicted value of a rock or fluid property to a cell within a multi-dimensional layered model that is correlated to a subterranean geologic volume;

simultaneously calculating first predicted values of seismic response for multiple angles of incidence for said cell from a response model using said first predicted value of said rock or fluid property and based upon compressional velocity, shear velocity and bulk density;

simultaneously comparing synthetic traces that are generated from said first predicted values of seismic response for multiple angles of incidence to corresponding traces obtained from corresponding multiple angles of incidence of actual seismic data for the subterranean geologic volume to ascertain a difference between said synthetic traces and said corresponding traces;

adjusting said first predicted value of said designated rock or fluid property in response to said difference to create a second predicted value of said designated rock or fluid property, wherein said second predicted value generally reduces said difference; and

storing said second predicted value for later use.

21. The method of claim 20 wherein said step of simultaneously comparing involves maintaining consistency between the types, scales and dimensions of said first predicted values of a designated rock or fluid property and said corresponding multiple angles of incidence of actual seismic data.