Granted Patent B2

US 12,360,268 · App. 17/395,100 · Granted Jul 15, 2025

Inventors: Weiguang He (Beijing, CN); Yubing Li (Beijing, CN); Yue Ma (Beijing, CN); Yi Luo (Dhahran, SA)

Assignee: SAUDI ARABIAN OIL COMPANY

G01V1/282G01V1/345G01V2210/624G01V2210/66

Abstract

Claims (68)

1. A method, comprising:

obtaining, by a seismic acquisition system, an observed seismic dataset for a subterranean region of interest;

receiving, by a seismic processing system, the observed seismic dataset from the seismic acquisition system;

receiving, by the seismic processing system, a seismic velocity model for the subterranean region of interest;

using the seismic processing system:

generating a simulated seismic dataset based, at least in part, on the seismic velocity model and a geometry of the observed seismic dataset,

determining a transformed observed seismic dataset, wherein determining the transformed observed seismic dataset comprises applying a nonlinear amplitude transform to the observed seismic dataset, wherein the nonlinear amplitude transform is selected from the group consisting of a hyperbolic tangent, an arc tangent, and an error function,

determining a transformed simulated seismic dataset, wherein determining the transformed simulated seismic dataset comprises applying the nonlinear amplitude transform to the simulated seismic dataset,

forming an objective function based, at least in part, on the transformed observed seismic dataset and on the transformed simulated seismic dataset,

determining an updated seismic velocity model based, at least in part, upon finding an extremum of the objective function,

forming a seismic image of the subterranean region of interest based, at least in part, on the updated seismic velocity model, and

identifying a portion of the subterranean region of interest with a likelihood of containing oil and gas based, at least in part, on the seismic image of the subterranean region of interest;

determining, using a computer processor, a well path through the portion of the subterranean region of interest based, at least in part, on the identified portion of the subterranean region of interest; and

drilling, using a drilling system, a well guided by the well path.

2. The method of claim 1 , wherein the nonlinear amplitude transform magnifies a first value of a seismic amplitude more than a second value of the seismic amplitude, wherein the first value is smaller than the second value.

3. The method of claim 1 , wherein the seismic velocity model comprises an elastic seismic velocity model.

4. The method of claim 1 , wherein the seismic velocity model comprises an anisotropic seismic velocity model.

5. The method of claim 1 , wherein determining the updated seismic velocity model comprises:

determining an adjoint source for each of a plurality of seismic receiver locations;

simulating, for each of the plurality of seismic receiver locations, propagation of a seismic wave generated by the adjoint source backward in time;

simulating, for each of a plurality of seismic source locations, propagation of a seismic wave forward in time;

determining a seismic velocity gradient using an imaging condition based, at least in part, on the simulated propagation of the seismic waves propagated forward in time and the simulated propagation of the seismic waves propagated backward in time;

determining a seismic velocity increment model by scaling the seismic velocity gradient based, at least in part, on the extremum of the objective function; and

updating the seismic velocity model based, at least in part, on the seismic velocity increment model.

6. A non-transitory computer readable medium storing instructions executable by a seismic processing system, wherein the instructions, when executed, cause the seismic processing system to perform steps comprising:

receiving, from a seismic acquisition system, an observed seismic dataset for a subterranean region of interest;

receiving a seismic velocity model for the subterranean region of interest;

generating a simulated seismic dataset based, at least in part, on the seismic velocity model and a geometry of the observed seismic dataset;

determining a transformed observed seismic dataset, wherein determining the transformed observed seismic dataset comprises applying a nonlinear amplitude transform to the observed seismic dataset, wherein the nonlinear amplitude transform is selected from the group consisting of a hyperbolic tangent, an arc tangent, and an error function;

determining a transformed simulated seismic dataset, wherein determining the transformed simulated seismic dataset comprises applying the nonlinear amplitude transform to the simulated seismic dataset;

forming an objective function based, at least in part, on the transformed observed seismic dataset and on the transformed simulated seismic dataset;

determining an updated seismic velocity model based, at least in part, upon finding an extremum of the objective function;

forming a seismic image of the subterranean region of interest based, at least in part, on the updated seismic velocity model;

identifying a portion of the subterranean region of interest with a likelihood of containing oil and gas based, at least in part, on the seismic image; and

determining a well path through the portion of the subterranean region of interest based, at least in part, on the identified portion of the subterranean region of interest.

7. The non-transitory computer readable medium of claim 6 , wherein the nonlinear amplitude transform magnifies a first value of a seismic amplitude more than a second value of the seismic amplitude, wherein the first value is smaller than the second value.

8. The non-transitory computer readable medium of claim 6 , wherein the seismic velocity model comprises an elastic seismic velocity model.

9. The non-transitory computer readable medium of claim 6 , wherein the seismic velocity model comprises an anisotropic seismic velocity model.

10. The non-transitory computer readable medium of claim 6 , wherein determining the updated seismic velocity model comprises:

determining an adjoint source for each of a plurality of seismic receiver locations based, at least in part, on the observed seismic dataset and the simulated seismic dataset;

simulating, for each of a plurality of seismic source locations and the plurality of seismic receiver locations, propagation of a seismic wave generated by the adjoint source backward in time;

simulating, for each of the plurality of seismic source locations, propagation of a seismic wave forward in time;

determining a seismic velocity gradient using an imaging condition based, at least in part, on the simulated propagation of the seismic waves propagated forward in time and the simulated propagation of the seismic waves propagated backward in time;

determining a seismic velocity increment model by scaling the seismic velocity gradient based, at least in part, on the extremum of the objective function; and

updating the seismic velocity model based, at least in part, on the seismic velocity increment model.

11. A system, comprising:

a seismic acquisition system, configured to acquire an observed seismic dataset for a subterranean region of interest; and

a seismic processing system configured to:

receive the observed seismic dataset from the seismic acquisition system,

receive a seismic velocity model for the subterranean region of interest,

generate a simulated seismic dataset based, at least in part, on the seismic velocity model and a geometry of the observed seismic dataset,

determine a transformed observed seismic dataset, wherein determining the transformed observed seismic dataset comprises applying a nonlinear amplitude transform to the observed seismic dataset, wherein the nonlinear amplitude transform is selected from the group consisting of a hyperbolic tangent, an arc tangent, and an error function,

determine a transformed simulated seismic dataset, wherein determining the transformed simulated seismic dataset comprises applying the nonlinear amplitude transform to the simulated seismic dataset,

form an objective function based, at least in part, on the transformed observed seismic dataset and on the transformed simulated seismic dataset,

determine an updated seismic velocity model based, at least in part, upon finding an extremum of the objective function,

form a seismic image of the subterranean region of interest based, at least in part, on the updated seismic velocity model, and

identify a portion of the subterranean region of interest with a likelihood of containing oil and gas based, at least in part, on the seismic image; and

a drilling system configured to drill a wellbore to penetrate the portion.

12. The system of claim 11 , wherein the nonlinear amplitude transform magnifies a first value of a seismic amplitude more than a second value of the seismic amplitude, wherein the first value is smaller than the second value.

13. The system of claim 11 , wherein the seismic velocity model comprises an elastic seismic velocity model.

14. The system of claim 11 , wherein the seismic velocity model comprises an anisotropic seismic velocity model.

15. The system of claim 11 , wherein determining the updated seismic velocity model comprises:

determining an adjoint source for each of a plurality of seismic receiver locations based, at least in part, on the observed seismic dataset and the simulated seismic dataset;

simulating, for each of a plurality of seismic source locations and the plurality of seismic receiver locations, propagation of a seismic wave generated by the adjoint source backward in time;

simulating, for each of the plurality of seismic source locations, propagation of a seismic wave forward in time;

determining a seismic velocity gradient using an imaging condition based, at least in part, on the simulated propagation of the seismic waves propagated forward in time and the simulated propagation of the seismic waves propagated backward in time;

determining a seismic velocity increment model by scaling the seismic velocity gradient based, at least in part, on the extremum of the objective function; and

updating the seismic velocity model based, at least in part, on the seismic velocity increment model.

Assignments (3)

Continuity (1)

References Cited (27)