Joint inversion with unknown lithology

Patent metadata

Field	Value
Publication number	US-9453929-B2
Application number	US-201214111519-A
Country	US
Kind code	B2
Filing date	Mar 9, 2012
Priority date	Jun 2, 2011
Publication date	Sep 27, 2016
Grant date	Sep 27, 2016

Abstract

Method for joint inversion of geophysical data to obtain 3-D models of geological parameters for subsurface regions of unknown lithology. Two or more data sets of independent geophysical data types are obtained, e.g. seismic and electromagnetic. Then they are jointly inverted, using structural coupling, to infer geophysical parameter volumes, e.g. acoustic velocity and resistivity. Regions of common lithology are next identified based on similar combinations of geophysical parameters. Then a joint inversion of the multiple data types is performed in which rock physics relations vary spatially in accordance with the now-known lithology, and 3-D models of geological properties such as shale content and fracture density are inferred. The computational grid for the last inversion may be defined by the lithology regions, resulting in average geological properties over such regions, which may then be perturbed to determine uncertainty in lithologic boundaries.

First claim

The invention claimed is:

1. A computer-implemented method for joint inversion of two or more sets of geophysical data of different types, measured in surveys of a subsurface region, to obtain a model of at least one geological property for the subsurface region, said method comprising:
using a computer to jointly invert the sets of geophysical data, using structural coupling between the different data types, to obtain models of geophysical properties corresponding to the sets of geophysical data, said structural coupling being determined from the geophysical data or from a priori knowledge of the subsurface region;
partitioning the subsurface region into sub-regions based on similar combinations of geophysical parameters in the geophysical property models that correlate to particular lithologies, thereby defining lithology sub-regions;
determining mathematical rock physics relationships appropriate for each lithology, said rock physics relationships relating geological properties to geophysical properties; and
using a computer to jointly invert the sets of geophysical data, using the rock physics relationships according to lithology sub-region, to obtain a model of one or more of the geological properties for the subsurface region,
wherein at least one of the sets of geophysical data are filtered using a low-pass filter before the jointly inverting with structural coupling, which inverting uses a computational grid with coarseness of scale consistent with frequency content passed by the low-pass filter.

2. The method of claim 1 , wherein the two or more different types of geophysical data are two or more of a group consisting of: seismic reflection data, seismic refraction data, controlled source electromagnetic data, magnetotelluric data, gravity data, and magnetic data.

3. The method of claim 1 , wherein assigning a lithology to each sub-region is implemented using at least one of a group consisting of: a database correlating lithology to geophysical properties; a neural network; statistical cluster analysis; and pattern recognition.

4. The method of claim 1 , wherein an initial geological properties model for the joint inversion using rock physics relationships is developed using the models of geophysical properties obtained from the joint inversion with structural coupling, substituted into the rock physics relationships.

5. The method of claim 1 , wherein a geological properties model to constrain inversion for the joint inversion using rock physics relationships is developed using the models of geophysical properties obtained from the joint inversion with structural coupling, substituted into the rock physics relationships.

6. The method of claim 5 , further comprising using the models of geophysical properties obtained from the joint inversion with structural coupling as additional constraint on the joint inversion using rock physics relationships.

7. The method of claim 1 , wherein in the joint inverting using rock physics relationships, each lithology sub-region is treated as a single cell in a computational grid, and the geological properties obtained are average properties over each lithology sub-region.

8. The method of claim 7 , further comprising generating a plurality of perturbations of the average properties and then for each perturbation performing a joint inversion of the sets of geophysical data using geometrical coupling between the different data types and also using said average properties, said geometrical coupling comprising shape and size of the lithology sub-regions, then determining uncertainty in boundaries of the lithology sub-regions.

9. A method for producing hydrocarbons from a subsurface region, comprising:
obtaining two or more data sets of different geophysical data types produced by surveys of the subsurface region;
jointly inverting the two or more data sets to obtain a model of at least one geological property for the subsurface region, using a method of claim 1 ;
using the model in determining where to drill a well into the subsurface region; and
drilling the well and producing hydrocarbons from it.

10. A computer-implemented method for joint inversion of two or more sets of geophysical data of different types, measured in surveys of a subsurface region, to obtain a model of at least one geological property for the subsurface region, said method comprising:
using a computer to jointly invert the sets of geophysical data, using structural coupling between the different data types, to obtain models of geophysical properties corresponding to the sets of geophysical data, said structural coupling being determined from the geophysical data or from a priori knowledge of the subsurface region;
partitioning the subsurface region into sub-regions based on similar combinations of geophysical parameters in the geophysical property models that correlate to particular lithologies, thereby defining lithology sub-regions;
determining mathematical rock physics relationships appropriate for each lithology, said rock physics relationships relating geological properties to geophysical properties; and
using a computer to jointly invert the sets of geophysical data, using the rock physics relationships according to lithology sub-region, to obtain a model of one or more of the geological properties for the subsurface region, wherein the using the computer to jointly invert includes
simulating synthetic data for each data type using an initial model of a corresponding geophysical property, then determining a misfit between synthetic data and corresponding measured data, then using the misfit to adjust the initial model, then iterating beginning with the simulating synthetic data until the misfit is less than a predetermined tolerance or other stopping condition is reached,
wherein the using the misfit to adjust the initial model comprises selecting a mathematical objective function that includes a data misfit term for each data type multiplied by a weight for that data type and a term expressing the structural coupling.

11. The method of claim 10 , wherein the term expressing the structural coupling comprises a cross-gradient constraint.

12. The method of claim 10 , wherein the jointly inverting the sets of geophysical data using structural coupling further comprises in at least one iteration changing at least one of a group consisting of: type of structural coupling; frequency content of one or more data types; offsets for data types using multiple offsets, where “offset” is source-receiver spacing; time window used for time-domain data; spacing of discrete cells in a grid on which a model is expressed for purposes of inversion; data types used; and weights assigned to each data type.

13. The method of claim 10 , wherein the jointly inverting the sets of geophysical data using the rock physics relationships further comprises in at least one iteration changing at least one of a group consisting of: frequency content of one or more data types; offsets for data types using multiple offsets, where “offset” is source-receiver spacing; time window used for time-domain data; spacing of discrete cells in a grid on which a model is expressed for purposes of inversion; data types used; and weights assigned to each data type.

14. A computer-implemented method for joint inversion of two or more sets of geophysical data of different types, measured in surveys of a subsurface region, to obtain a model of at least one geological property for the subsurface region, said method comprising:
using a computer to jointly invert the sets of geophysical data, using structural coupling between the different data types, to obtain models of

Assignees

• Schmedes Jan
• Dicaprio Christopher J
• Jing Charlie
• Leahy Garrett M
• Mullur Anoop A
• Saltzer Rebecca L
• Exxonmobil Upstream Res Co

Inventors

• Schmedes Jan
• Dicaprio Christopher J
• Jing Charlie
• Leahy Garrett M
• Mullur Anoop A
• Saltzer Rebecca L

Classifications

• G01V11/00 Primary

  Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00 · CPC title

• G01V1/282 Primary

  Application of seismic models, synthetic seismograms · CPC title

• G01V3/38

  Processing data, e.g. for analysis, for interpretation, for correction · CPC title