Dynamic resolution full-waveform inversion method and system

Patent metadata

Field	Value
Publication number	US-12560732-B2
Application number	US-202318302119-A
Country	US
Kind code	B2
Filing date	Apr 18, 2023
Priority date	Dec 22, 2022
Publication date	Feb 24, 2026
Grant date	Feb 24, 2026

Abstract

A method for raw seismic data inversion includes receiving raw seismic data acquired over an underground formation, receiving an initial velocity model of the underground formation, performing a dynamic resolution full waveform inversion, DR-FWI, on the raw seismic data so that a tomography component of a velocity gradient G is compensated to generate a compensated tomographic component while a diving wave component and a migration component are preserved, outputting an updated velocity model based on an illumination compensated velocity gradient G′, which is calculated based on the compensated tomographic component), and locating natural resources in the underground formation using the updated velocity model updated by the DR-FWI.

First claim

What is claimed is:

1 . A method for raw seismic data inversion, the method comprising:
receiving raw seismic data acquired over an underground formation; receiving an initial velocity model of the underground formation; performing a dynamic resolution full waveform inversion, DR-FWI, on the raw seismic data so that a tomography component of a velocity gradient G is compensated to generate a compensated tomographic component while a diving wave component and a migration component are preserved; outputting an updated velocity model based on an illumination compensated velocity gradient G′, which is calculated based on the compensated tomographic component; and locating natural resources in the underground formation using the updated velocity model updated by the DR-FWI.

2 . The method of claim 1 , wherein the step of performing comprises:
calculating the velocity gradient G from full traces of the seismic data; and separating the velocity gradient G into (1) a diving wave component and (2) the migration component and the tomography component by applying a diving wave mask.

3 . The method of claim 2 , further comprising:
separating the tomography component from the migration component; and compensating only the tomography component to obtain the compensated tomographic component.

4 . The method of claim 3 , wherein the step of compensating comprises:
multiplying the tomography component with a varying scalar to obtain the compensated tomographic component.

5 . The method of claim 4 , wherein the varying scalar is equal to a square root of a ratio between (1) an energy of a forward propagated wavefield, and (2) an energy of a scattered wavefield.

6 . The method of claim 4 , wherein the varying scalar is equal to an inverse of a reflectivity calculated from a current model at each iteration.

7 . The method of claim 4 , wherein the varying scalar is equal to a root mean square of a ratio of the migration component and the tomography component.

8 . The method of claim 4 , further comprising:
combining the compensated tomographic component with the unchanged migration component and the diving component to obtain the illumination compensated velocity gradient G′.

9 . The method of claim 8 , further comprising:
estimating a step length or Hessian matrix for the DR-FWI based on the illumination compensated velocity gradient G′; determining whether a cost function reached a predefined threshold; and returning to the step of performing if the cost function did not reach the predefined threshold.

10 . A computing device for raw seismic data inversion, the device comprising:
an interface for receiving raw seismic data acquired over an underground formation, and for receiving an initial velocity model of the underground formation; and a processor connected to the interface and configured to: perform a dynamic resolution full waveform inversion, DR-FWI, on the raw seismic data so that a tomography component of a velocity gradient G is compensated to generate a compensated tomographic component while a diving wave component and a migration component are preserved; output an updated velocity model based on an illumination compensated velocity gradient G′, which is calculated based on the compensated tomographic component; and locate natural resources in the underground formation using the updated velocity model updated by the DR-FWI.

11 . The device of claim 10 , wherein the processor is further configured to:
calculate the velocity gradient G from full traces of the seismic data; and separate the velocity gradient G into (1) a diving wave component and (2) the migration component and the tomography component by applying a diving wave mask.

12 . The device of claim 11 , wherein the processor is further configured to:
separate the tomography component from the migration component; and compensate only the tomography component to obtain the compensated tomographic component.

13 . The device of claim 12 , wherein the processor is further configured to:
multiply the tomography component with a varying scalar to obtain the compensated tomographic component.

14 . The device of claim 13 , wherein the varying scalar is equal to a square root of a ratio between (1) an energy forward propagated wavefield, and (2) an energy of a scattered wavefield.

15 . The device of claim 13 , wherein the varying scalar is equal to an inverse of a reflectivity calculated from a current model at each iteration.

16 . The device of claim 13 , wherein the varying scalar is equal to a root mean square of a ratio of the migration component and the tomography component.

17 . The device of claim 13 , wherein the processor is further configured to:
combine the compensated tomographic component with the unchanged migration component and the diving component to obtain the illumination compensated velocity gradient G′.

18 . The device of claim 17 , the processor is further configured to:
estimate a step length or Hessian matrix for the DR-FWI based on the illumination compensated velocity gradient G′; determine whether a cost function reached a predefined threshold; and return to the step of performing if the cost function did not reach the predefined threshold.

19 . A non-transitory computer readable recording medium storing executable codes that when executed by a computer make the computer perform a method for raw seismic data inversion, the method comprising:
receiving raw seismic data acquired over an underground formation; receiving an initial velocity model of the underground formation; performing a dynamic resolution full waveform inversion, DR-FWI, on the raw seismic data so that a tomography component of a velocity gradient G is compensated to generate a compensated tomographic component while a diving wave component and a migration component are preserved; outputting an updated velocity model based on an illumination compensated velocity gradient G′, which is calculated based on the compensated tomographic component; and locating natural resources in the underground formation using the updated velocity model updated by the DR-FWI.

20 . The medium of claim 19 , wherein the step of performing comprises:
calculating the velocity gradient G from full traces of the seismic data; separating the velocity gradient G into (1) a diving wave component and (2) the migration component and the tomography component by applying a diving wave mask, separating the tomography component from the migration component; and compensating only the tomography component to obtain the compensated tomographic component.

Assignees

• Cgg Services Sas

Inventors

• Wang Min
• Xie Yi

Classifications

• G01V2210/675

  Wave equation; Green's functions · CPC title

• G01V2210/6222

  Velocity; travel time · CPC title

• G01V2210/512

  Pre-stack · CPC title

• G01V1/303 Primary

  for determining velocity profiles or travel times · CPC title

• G01V2210/51

  Migration · CPC title