Efficient solutions of inverse problems

Patent metadata

Field	Value
Publication number	US-11143788-B2
Application number	US-201515572444-A
Country	US
Kind code	B2
Filing date	May 8, 2015
Priority date	May 8, 2015
Publication date	Oct 12, 2021
Grant date	Oct 12, 2021

Abstract

A method of solving a geophysical inverse problem for estimating a physical parameter, the method comprising providing a model vector representing the physical parameter, transforming the model vector by a first operator to provide a first transformed model vector, solve the inverse problem for the first transformed model vector to provide a first solution, transforming the model vector by a second operator to produce a second transformed model vector, solve the inverse problem for the second transformed model vector to provide a second solution, calculating a weighted sum of the first solution and the second solution.

First claim

The invention claimed is:

1. A computer-implemented method of estimating a physical parameter of a subsurface of the Earth by solving a geophysical inverse problem using a set of measured geophysical data and simulations of the set of measured geophysical data, wherein the simulations are dependent on a model vector representing the physical parameter, the method comprising:
providing a model vector representing the physical parameter,
transforming the model vector by a first operator to provide a first transformed model vector, wherein transforming the model vector by the first operator comprises reducing a number of free parameters in the inverse problem according to a first compression factor,
solving the inverse problem for the first transformed model vector to provide a first solution,
transforming the same model vector by a second operator to produce a second transformed model vector, wherein the step of transforming the model vector by the second operator comprises reducing a number of free parameters according to a second compression factor,
solving the inverse problem for the second transformed model vector to provide a second solution,
calculating a weighted sum of the first solution and the second solution,
providing an estimation of the physical parameter based on the calculated weighted sum of the first solution and the second solution, wherein the weighted sum of the first solution and the second solution minimizes a distance between the set of measured geophysical data and the simulation of the set of measured geophysical data,
determining a presence of hydrocarbons or water using the estimation of the physical parameter, and
exploring hydrocarbons or water in the subsurface of the Earth based on the determined presence of hydrocarbons or water.

2. The computer-implemented method of claim 1 , wherein the physical parameter is a conductivity of the subsurface of the Earth.

3. The computer-implemented method according to claim 1 , wherein transforming the model vector by a first operator comprises applying a first Fourier transform to the model vector and wherein transforming the model vector by a second operator comprises applying a second Fourier transform to the model vector.

4. The computer-implemented method according to claim 3 , wherein transforming the model vector by the first operator or the second operator further comprises applying a filter.

5. The computer-implemented method according to claim 1 , wherein reducing the number of free parameters in the inverse problem comprises interpolating the model vector defined on a fine grid to provide the first transformed model vector defined on a first coarse grid, wherein the resolution of the first coarse grid is lower than the resolution of the fine grid, and interpolating the model vector defined on a fine grid to provide the second transformed model vector defined on a second coarse grid, wherein the resolution of the second coarse grid is lower than the resolution of the fine grid.

6. The computer-implemented method according to claim 5 , wherein the coarse grid is uniform or non-uniform.

7. The computer-implemented method according to claim 5 , wherein the coarse grid depends on prior knowledge of the physical structure of a formation.

8. The computer-implemented method of claim 5 , wherein the first coarse grid is shifted with respect to the second coarse grid.

9. The computer-implemented method according to claim 1 , wherein the set of measured geophysical data are one of: a set of seismic, electric, magnetic or gravimetric measurements.

10. The computer-implemented method according to claim 1 , wherein the model vector comprises a model update in an iterative inversion method.

11. The computer-implemented method of claim 10 , wherein the minimal distance is estimated with a line search.

12. The computer-implemented method of claim 1 , wherein in the inversion problem is a linear or a non-linear inversion problem.

13. A computer system arranged to carry out the method of claim 1 .

14. Computer software installed on a non-transitory computer readable medium of a computer system and arranged to cause the computer system to carry out the method of claim 1 .

15. The computer-implemented method of claim 1 , further comprising the step of constructing a 3D model of the subsurface of the Earth using the estimation of the physical parameter.

16. A computer-implemented method, comprising:
estimating a physical parameter of a subsurface of the Earth by solving a geophysical inverse problem using a set of measured geophysical data and simulations of the set of measured geophysical data, wherein the simulations are dependent on a model vector representing the physical parameter, comprising:
providing a model vector representing a physical parameter of a subsurface of the Earth;
transforming the model vector by a first operator to provide a first transformed model vector by reducing a number of free parameters in an inverse problem according to a first compression factor;
solving the inverse problem for the first transformed model vector to provide a first solution;
transforming the same model vector by a second operator to produce a second transformed model vector by reducing the number of free parameters according to a second compression factor;
solving the inverse problem for the second transformed model vector to provide a second solution,
calculating a weighted sum of the first solution and the second solution, and
providing an estimation of the physical parameter based on the calculated weighted sum of the first solution and the second solution, wherein the weighted sum of the first solution and the second solution minimizes a distance between the set of measured geophysical data and the simulation of the set of measured geophysical data;

determining a presence of hydrocarbons or water using the estimation of the physical parameter; and
exploring hydrocarbons or water in the subsurface of the Earth based on the determined presence of hydrocarbons or water,
wherein the set of measured geophysical data are one of: a set of seismic, electric, magnetic or gravimetric measurements.

17. The computer-implemented method of claim 16 , further comprising the step of constructing a 3D model of the subsurface of the Earth using the estimation of the physical parameter.

18. A computer-implemented method, comprising:
acquiring a set of geophysical data from a subsurface of the Earth;
estimating a physical parameter of the subsurface of the Earth by solving a geophysical inverse problem using the set of measured geophysical data and simulations of the set of measured geophysical data, wherein the simulations are dependent on a model vector representing the physical parameter, comprising:
providing a model vector representing a physical parameter of a subsurface of the Earth;
transforming the model vector by a first operator to provide a first transformed model vector by reducing a number of free parameters in an inverse problem according to a first compression factor;
solving the inverse problem for the first transformed model vector to provide a first solution;
transforming the same model vector by a second operator to produce a second transformed model vector by reducing the number of free parameters according to a second compression factor;
solving the inverse problem for the second transformed model vector to provide a second solution,
calculating a weighted sum of the first solution and the second solution, and
providing an estimation of the physical parameter based on t

Assignees

• Statoil Petroleum As

Inventors

• Causse Emmanuel
• Nguyen Anh Kiet
• Nordskag Janniche Iren
• Maaø Frank Agnar

Classifications

• G06F17/16

  Matrix or vector computation {, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization (matrix transposition G06F7/78)} · CPC title

• G01V2003/086

  Processing · CPC title

• G06F2111/10

  Numerical modelling · CPC title

• G01V3/165

  operating with magnetic or electric fields produced or modified by the object or by the detecting device (with electromagnetic waves G01V3/17) · CPC title

• G01V3/083

  Controlled source electromagnetic [CSEM] surveying · CPC title