Uncertainty estimation for large-scale nonlinear inverse problems using geometric sampling and covariance-free model compression

Patent metadata

Field	Value
Publication number	US-9619590-B2
Application number	US-201113634522-A
Country	US
Kind code	B2
Filing date	Mar 14, 2011
Priority date	Mar 19, 2010
Publication date	Apr 11, 2017
Grant date	Apr 11, 2017

Abstract

A method for uncertainty estimation for nonlinear inverse problems includes obtaining an inverse model of spatial distribution of a physical property of subsurface formations. A set of possible models of spatial distribution is obtained based on the measurements. A set of model parameters is obtained. The number of model parameters is reduced by covariance free compression transform. Upper and lower limits of a value of the physical property are mapped to orthogonal space. A model polytope including a geometric region of feasible models is defined. At least one of random and geometric sampling of the model polytope is performed in a reduced-dimensional space to generate an equi-feasible ensemble of models. The reduced-dimensional space includes an approximated hypercube. Probable model samples are evaluated based on data misfits from among an equi-feasible model ensemble determined by forward numerical simulation. Final uncertainties are determined from the equivalent model ensemble and the final uncertainties are displayed in at least one map.

First claim

What is claimed is:

1. A method comprising:
obtaining an inverse model of a spatial distribution of a physical property of subsurface formations made using time-varying measurements acquired by sensors deployed proximate the subsurface formations, wherein the spatial distribution of the physical property varies in at least in at least two dimensions;
obtaining a set of possible models of the spatial distribution based on the time-varying measurements;
obtaining a set of model parameters;
reducing a number of model parameters in the set by a covariance-free compression transform;
mapping upper and lower limits of a physical bound value to an orthogonal space defining a model polytope including a geometric region of feasible models;
performing at least one of random and geometric sampling of the model polytope in a reduced-dimensional space to generate an equi-feasible ensemble of models, the reduced-dimensional space consisting of a hypercube approximation;
evaluating probable model samples based on misfits from among the equi-feasible ensemble of models, the misfits being determined by forward simulation for time-varying measurements and at least in two-dimensions, and data misfit rejection;
computing final uncertainties based on the equi-feasible model ensemble; and
displaying the final uncertainties in at least one map.

2. The method of claim 1 wherein the at least one map comprises a probability cut-off map.

3. The method of claim 1 wherein the at least one map comprises a parameter value cut-off map.

4. The method of claim 1 wherein the time-varying measurements comprise controlled source electromagnetic field measurements.

5. The method of claim 1 wherein the covariance free compression transform comprises singular value decomposition.

6. The method of claim 1 wherein the covariance free compression transform comprises a discrete cosine transform.

Assignees

• Tompkins Michael J
• Fernandez-Martinez Juan Luis
• Schlumberger Technology Corp

Inventors

• Tompkins Michael J
• Fernandez-Martinez Juan Luis

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

• G06F17/11

  for solving equations {, e.g. nonlinear equations, general mathematical optimization problems (optimization specially adapted for a specific administrative, business or logistic context G06Q10/04)} · CPC title

• G06F30/20

  Design optimisation, verification or simulation (optimisation, verification or simulation of circuit designs G06F30/30) · CPC title

• G01V20/00 Primary

  Geomodelling in general · CPC title

• G06F17/5009 Primary

  Physics · mapped topic