Elastic impedance estimation for inversion of far offset seismic sections

• US 6,058,073 A
• Filed: 03/30/1999
• Issued: 05/02/2000
• Est. Priority Date: 03/30/1999
• Status: Expired due to Term

• Current Assignee: Atlantic Richfield Company Incorporated (BP PLC)
• Original Assignee: Atlantic Richfield Company Incorporated (BP PLC)
• Inventors: VerWest, Bruce J.
• Primary Examiner(s): Oda, Christine K.
• Assistant Examiner(s): JOLLY, ANTHONY S
• Application Number:
  US09/281,378
• Time in Patent Office: 399 Days
• Field of Search: 367/31, 367/83, 367/38, 367/47, 367/68
• US Class Current: 367/31
• CPC Class Codes: G01V 1/28 Processing seismic data, e.g. for interpretation or for event detection G01V1/48 takes precedence
• CPC Class Codes: G01V 2210/66 Subsurface modeling

Abstract

A computer-implemented method, and system implementing the method, of performing impedance inversion of seismic survey data are disclosed. Sonic and density well logs are used to generate an elastic impedance model, at selected angles of incidence or offset range, using an expression for elastic impedance that depends upon the measured or estimated compressional velocity, shear velocity, and density of the subsurface layers in the survey region. The elastic impedance expression also includes a reference density value, and is dependent upon the ray parameter (or angle of incidence). This elastic impedance model is used, for example by way of a pseudo-density log, in estimating the input wavelet, and in carrying out an impedance inversion of the seismic survey over the entire survey region, for at least two angles of incidence. Low-frequency components of the elastic impedance, at the selected angles of incidence, are merged with the corresponding impedance inversion sections, to produce a merged impedance model at the selected angles of incidence. Rock properties, such as Poisson'"'"'s ratio, are then retrieved from a comparison, such as a ratio, of the impedance traces at the varying angles of incidence for common midpoints.

Claims

1. A method of operating a computer to perform an inversion of seismic survey traces associated with a plurality of midpoints in a survey region, and obtained over a plurality of offset distances, comprising the steps of:

retrieving well log data including compressional velocity, shear velocity, and density data, over depth for at least one well location in the survey region;

selecting a first angle of incidence;

generating an elastic impedance model using the compressional velocity, shear velocity, and density log data and dependent upon the first angle of incidence;

retrieving, from memory, seismic survey data for a plurality of midpoints in the survey region having offset distances at or near the first angle of incidence;

stacking the retrieved seismic survey data by midpoint;

performing an inversion of the stacked retrieved seismic survey data according to the generated elastic impedance model, to generate a first impedance section for the survey region associated with the first angle of incidence;

combining the first impedance section with a low-frequency impedance model associated with the first angle of incidence, to generate a first merged impedance model for the survey region associated with the first angle of incidence;

selecting a second angle of incidence;

repeating the generating, retrieving, stacking, performing, and combining steps to generate a second merged impedance model for the survey region associated with the second angle of incidence; and

producing a survey of rock properties in the survey region using a comparison of the first and second merged impedance models at corresponding depth points.

2. The method of claim 1, further comprising:

deriving a reference density for the survey region from the retrieved density data;

wherein the steps of generating an elastic impedance model also uses the reference density.

3. The method of claim 2, wherein the steps of generating an elastic impedance model comprise evaluating the expression:

##EQU27## for each of a plurality of layers in the survey region, where α

is the compressional velocity for the layer, where β

is the shear velocity for the layer, where ρ

is the density for the layer, where ρ

₀ is the reference density, and where p is the ray parameter corresponding to the selected angle of incidence.

4. The method of claim 1, further comprising:

before the step of performing an inversion, retrieving seismic survey data for at least one midpoint respectively corresponding to the at least one well location in the survey region; and

tying the well log data to the retrieved seismic survey data at the corresponding at least one midpoint and at least one well location.

5. The method of claim 1, further comprising:

before the step of performing an inversion, retrieving seismic survey data for at least one midpoint respectively corresponding to the at least one well location in the survey region; and

estimating an input wavelet by applying the elastic impedance model to the retrieved seismic survey data for the at least one midpoint.

6. The method of claim 5, wherein the step of performing an inversion comprises:

applying the estimated input wavelet to the stacked retrieved seismic data for each of the plurality of midpoints, to recover a reflectivity series for each of the plurality of midpoints; and

determining an elastic impedance value for a plurality of depths for each of the plurality of midpoints, from the reflectivity series.

7. The method of claim 1, wherein the step of producing a survey of rock properties in the survey region comprises:

retrieving the first and second merged impedance models for a first one of the plurality of midpoints; and

for each of a plurality of depths associated with the first midpoint, determining the value of a selected rock property from a ratio of the impedance values of the first and second merged impedance models at the depth.

8. The method of claim 7, further comprising:

deriving a reference density for the survey region from the retrieved density data;

wherein the steps of generating an elastic impedance model comprise evaluating the expression;

##EQU28## for each of a plurality of layers in the survey region, where α

is the compressional velocity for the layer, where β

is the shear velocity for the layer, where ρ

is the density for the layer, where ρ

₀ is the reference density, and where p is the ray parameter corresponding to the selected angle of incidence;

and wherein the step of determining the value of a selected rock property comprises determining, from the ratio of the impedance values of the first and second merged impedance models, a ratio of shear velocity to compressional velocity for each of the plurality of depths.

9. The method of claim 8, wherein the step of determining the value of a selected rock property further comprises determining Poisson'"'"'s ratio from the ratio of shear velocity to compressional velocity, for each of the plurality of depths.

10. A computer-readable memory having a storage medium configured so that, when read and used by a computer, the computer is directed to perform an inversion of seismic survey traces associated with a plurality of midpoints in a survey region, and obtained over a plurality of offset distances, said computer directed by a plurality of operations comprising:

retrieving, from a memory device, well log data including compressional velocity, shear velocity, and density data, over depth for at least one well location in the survey region;

selecting a first angle of incidence;

generating an elastic impedance model using the compressional velocity, shear velocity, and density log data and dependent upon the first angle of incidence;

retrieving, from a memory device, seismic survey data for a plurality of midpoints in the survey region having offset distances at or near the first angle of incidence;

stacking the retrieved seismic survey data by midpoint;

performing an inversion of the stacked retrieved seismic survey data according to the generated elastic impedance model, to generate a first impedance section for the survey region associated with the first angle of incidence;

combining the first impedance section with a low-frequency impedance model associated with the first angle of incidence, to generate a first merged impedance model for the survey region associated with the first angle of incidence;

selecting a second angle of incidence;

repeating the generating, retrieving, stacking, performing, and combining steps to generate a second merged impedance model for the survey region associated with the second angle of incidence; and

producing and storing, in a memory device, a survey of rock properties in the survey region using a comparison of the first and second merged impedance models at corresponding depth points.

11. The computer-readable memory of claim 10, wherein the plurality of operations further comprises:

deriving a reference density for the survey region from the retrieved density data;

and wherein the operations of generating an elastic impedance model comprise evaluating the expression;

##EQU29## for each of a plurality of layers in the survey region, where α

is the compressional velocity for the layer, where β

is the shear velocity for the layer, where ρ

is the density for the layer, where ρ

₀ is the reference density, and where p is the ray parameter corresponding to the selected angle of incidence.

12. The computer-readable memory of claim 10, wherein the plurality of operations further comprises:

before the step of performing an inversion, retrieving seismic survey data for at least one midpoint respectively corresponding to the at least one well location in the survey region; and

estimating an input wavelet by applying the elastic impedance model to the retrieved seismic survey data for the at least one midpoint;

and wherein the operation of performing an inversion comprises;

applying the estimated input wavelet to the stacked retrieved seismic data for each of the plurality of midpoints, to recover a reflectivity series for each of the plurality of midpoints; and

determining an elastic impedance value for a plurality of depths for each of the plurality of midpoints, from the reflectivity series.

13. The computer-readable memory of claim 10, wherein the plurality of operations further comprises:

deriving a reference density for the survey region from the retrieved density data;

wherein the operations of generating an elastic impedance model comprise evaluating the expression;

##EQU30## for each of a plurality of layers in the survey region, where α

is the compressional velocity for the layer, where β

is the shear velocity for the layer, where ρ

is the density for the layer, where ρ

₀ is the reference density, and where p is the ray parameter corresponding to the selected angle of incidence;

and wherein the operation of producing a survey of rock properties in the survey region comprises;

retrieving the first and second merged impedance models for a first one of the plurality of midpoints; and

determining, from a ratio of the impedance values of the first and second merged impedance models at each of a plurality of depths associated with the first midpoint, a ratio of shear velocity to compressional velocity.

14. The computer-readable memory of claim 13, wherein the operation of producing a survey of rock properties in the survey region further comprises determining Poisson'"'"'s ratio from the ratio of shear velocity to compressional velocity, for each of the plurality of depths.

15. A digital computing system for performing an inversion of seismic survey traces associated with a plurality of midpoints in a survey region, and obtained over a plurality of offset distances, comprising:

a memory resource, for storing data corresponding to well log data including compressional velocity, shear velocity, and density data, over depth for at least one well location in the survey region, data corresponding to seismic survey data for a plurality of midpoints in the survey region, the seismic survey data corresponding to varying offset distances, and horizon data corresponding to the location of boundaries between adjacent formations in the survey region;

a computer system output device; and

a system computer, coupled to the memory resource and to the output device, for performing the operations of;

retrieving well log data from the memory resource for at least one well location in the survey region;

selecting a first angle of incidence;

generating an elastic impedance model using the compressional velocity, shear velocity, and density log data and dependent upon the first angle of incidence;

retrieving, from the memory resource, seismic survey data for a plurality of midpoints in the survey region having offset distances at or near the first angle of incidence;

stacking the retrieved seismic survey data by midpoint;

performing an inversion of the stacked retrieved seismic survey data according to the generated elastic impedance model, to generate a first impedance section for the survey region associated with the first angle of incidence;

combining the first impedance section with a low-frequency impedance model associated with the first angle of incidence, to generate a first merged impedance model for the survey region associated with the first angle of incidence;

selecting a second angle of incidence;

repeating the generating, retrieving, stacking, performing, and combining steps to generate a second merged impedance model for the survey region associated with the second angle of incidence; and

producing and forwarding, to the output device, a survey of rock properties in the survey region using a comparison of the first and second merged impedance models at corresponding depth points.

16. The digital computing system of claim 15, wherein the system computer is further programmed to perform the operations of:

deriving a reference density for the survey region from the retrieved density data;

and wherein the operations of generating an elastic impedance model comprise evaluating the expression;

##EQU31## for each of a plurality of layers in the survey region, where α

is the compressional velocity for the layer, where β

is the shear velocity for the layer, where ρ

is the density for the layer, where ρ

₀ is the reference density, and where p is the ray parameter corresponding to the selected angle of incidence.

17. The digital computing system of claim 15, wherein the system computer is further programmed to perform the operations of:

before the step of performing an inversion, retrieving seismic survey data for at least one midpoint respectively corresponding to the at least one well location in the survey region; and

estimating an input wavelet by applying the elastic impedance model to the retrieved seismic survey data for the at least one midpoint;

and wherein the operation of performing an inversion comprises;

applying the estimated input wavelet to the stacked retrieved seismic data for each of the plurality of midpoints, to recover a reflectivity series for each of the plurality of midpoints; and

determining an elastic impedance value for a plurality of depths for each of the plurality of midpoints, from the reflectivity series.

18. The digital computing system of claim 15, wherein the system computer is further programmed to perform the operations of:

deriving a reference density for the survey region from the retrieved density data;

and wherein the operations of generating an elastic impedance model comprise evaluating the expression;

##EQU32## for each of a plurality of layers in the survey region, where α

is the compressional velocity for the layer, where β

is the shear velocity for the layer, where ρ

is the density for the layer, where ρ

₀ is the reference density, and where p is the ray parameter corresponding to the selected angle of incidence;

and wherein the operation of producing a survey of rock properties in the survey region comprises;

retrieving the first and second merged impedance models for a first one of the plurality of midpoints; and

determining, from a ratio of the impedance values of the first and second merged impedance models at each of a plurality of depths associated with the first midpoint, a ratio of shear velocity to compressional velocity.

19. The digital computing system of claim 18, wherein the operation of producing a survey of rock properties in the survey region further comprises determining Poisson'"'"'s ratio from the ratio of shear velocity to compressional velocity, for each of the plurality of depths.