Full wavefield inversion with an image-gather-flatness constraint

Patent metadata

Field	Value
Publication number	US-11041971-B2
Application number	US-201916420774-A
Country	US
Kind code	B2
Filing date	May 23, 2019
Priority date	Jul 2, 2018
Publication date	Jun 22, 2021
Grant date	Jun 22, 2021

Abstract

Methods for updating a physical properties model of a subsurface region that combine advantages of FWI and tomography into a joint inversion scheme are provided. One method comprises obtaining measured seismic data; generating simulated seismic data using an initial model; computing at least one of an FWI gradient and a tomography gradient; minimizing a joint objective function E, wherein the objective function E is based on a combination of the FWI gradient and the tomography gradient or preconditioning of the FWI gradient or the tomography gradient; generating a final model based on the minimized joint objective function E; and using the final model to generate a subsurface image. The joint objective function E may be a function of one or both a FWI objective function CFWI and a tomography objective function CTomo. The joint objective function E may be defined as a weighted sum of CFWI and CTomo.

First claim

What is claimed is:

1. A method for updating a physical properties model of a subsurface region in an iterative full wavefield inversion (FWI) of seismic data, said method comprising:
obtaining measured seismic data;
generating, using a computer, simulated seismic data using an initial model;
computing, with the computer, an FWI gradient and a tomography gradient;
minimizing, with the computer, a joint objective function E, wherein the joint objective function E is based on a combination of the MT gradient and the tomography gradient or preconditioning of the FWI gradient or the tomography gradient, wherein the joint objective function E is a function of at least one of a FWI objective function C FWI and a tomography objective function C Tomo and wherein the joint objective function E is defined by:
E =(1−λ) C FWI +λC Tomo

where λ is a weighting hyper-parameter;
generating, with the computer, a final model based on the minimized joint objective function E; and
using the final model to generate an image of the subsurface.

2. The method of claim 1 , wherein minimizing the joint objective function E comprises defining a gradient of the joint objective function F as a function of a gradient of the FWI objective function C FWI , a tomography operator, and a residual depth error.

3. The method of claim 1 , wherein minimizing the joint objective function E comprises defining a gradient of the joint objective function E as a function of the product of a gradient of the FWI objective function C FWI and an approximate inverse of a tomographic Hessian matrix.

4. The method of claim 1 , wherein minimizing the joint objective function E comprises defining a gradient of the joint objective function E as a function of the product of a gradient of the tomography objective function C Tomo and an inverse of a FWI Hessian matrix.

5. The method of claim 4 , wherein the inverse of the FWI Hessian matrix is approximated using point-spread functions.

6. The method of claim 1 , wherein minimizing the joint objective function F comprises:
finding a search direction using the FWI gradient or the tomography gradient;
iteratively using the search direction to perform a line search to find an optimal step size to update a current model; and
iteratively updating the current model based on the optimal step size and the search direction;
wherein in the first iteration the current model is the initial model.

7. The method of claim 1 , wherein computing the tomography gradient comprises measuring residual moveout in image gathers generated from the seismic data and building tomographic equations to obtain the tomography gradient.

8. The method of claim 1 , further comprising using the image of the subsurface to manage hydrocarbons.

9. The method of claim 1 , wherein the final model is a velocity model.

Assignees

• Tang Yaxun
• Lu Rongrong
• He Yaofeng
• Exxonmobil Upstream Res Co

Inventors

• Tang Yaxun
• Lu Rongrong
• He Yaofeng

Classifications

• G01V2210/51

  Migration · CPC title

• G01V1/30 Primary

  Analysis (G01V1/50 takes precedence) · CPC title

• G01V1/282

  Application of seismic models, synthetic seismograms · CPC title

• G01V1/303 Primary

  for determining velocity profiles or travel times · CPC title

• G01V2210/67

  Wave propagation modeling · CPC title