Amplitude contrast seismic attribute

Patent metadata

Field	Value
Publication number	US-9341729-B2
Application number	US-201213438649-A
Country	US
Kind code	B2
Filing date	Apr 3, 2012
Priority date	Apr 6, 2011
Publication date	May 17, 2016
Grant date	May 17, 2016

Abstract

A method for visualizing seismic data of a subterranean formation, including obtaining an estimated dip field of the subterranean formation, wherein the estimated dip field represents a measure of deviation of a stratigraphic layer from flat, extracting a matrix data item surrounding a voxel of the seismic data, wherein the matrix data item is extracted from the seismic data based on a value of the estimated dip field surrounding the voxel, generating modified seismic data by at least applying a matrix operator to the seismic data, wherein the matrix operator calculates a partial derivative of the seismic data using the matrix data item, and displaying the modified seismic data.

First claim

What is claimed is:

1. A method for visualizing of a subterranean formation, comprising:
obtaining an estimated dip field of the subterranean formation, wherein the estimated dip field represents a measure of deviation of a stratigraphic layer from flat;
selecting a voxel from a seismic volume;
extracting, by a computer processor and from the seismic volume, a matrix data item, wherein the matrix data item comprises seismic data surrounding the voxel in the seismic volume, wherein the matrix data item is extracted from the seismic volume based on a value of the estimated dip field surrounding the voxel;
generating, by the computer processor and using a matrix operator, a plurality of amplitude gradients from the seismic data in the matrix data item;
normalizing, using the plurality of amplitude gradients, an amplitude value for the voxel to produce a normalized amplitude value for the voxel; and
generating, using the normalized amplitude value, an image representing a portion of the subterranean formation.

2. The method of claim 1 ,
wherein the matrix operator comprises a three dimensional (3D) operator,
wherein the matrix data item comprises a 3 by 3 by 3 matrix represented by Data(x, y, z) for x, y, z in {−1, 0, 1},
wherein the voxel comprises the seismic data, represented by Input(i, j, k), at position (i, j, k) of a voxel space,
wherein the matrix data item is extracted based on Data(x, y, z)=Input(i+x, j+y, k+z+x*dip IL (i, j, k)+y*dip XL (i, j, k)),
wherein dip IL (i, j, k) represents an inline dip of the estimated dip field at the position (i, j, k) of the voxel space, and
wherein dip XL , (i, j, k) represents a crossline dip of the estimated dip field at the position (i, j, k) of the voxel space.

3. The method of claim 1 , wherein the matrix operator comprises a Sobel operator.

4. The method of claim 1 , wherein normalizing the amplitude value further comprises:
adjusting, based on pre-determined weighting factors, contributions of first, second, and third partial derivatives to the normalized amplitude value along a perpendicular direction of the stratigraphic layer.

5. The method of claim 1 ,
wherein normalized amplitude value is represented by
S =√{square root over ( G X 2 +G Y 2 +W Z 2 G Z 2 )},

wherein G x , G y , and G z are proportional to a first, second, and third partial derivatives, respectively, and
wherein W z is a pre-determined fraction.

6. The method of claim 5 ,
wherein W z is in the range of approximately [0, 0.4].

7. A system for visualizing of a subterranean formation, comprising:
a seismic data processing module executing on a computer processor and configured to:
obtain an estimated dip field of the subterranean formation, wherein the estimated dip field represents a measure of deviation of a stratigraphic layer from flat;
select a voxel from a seismic volume;
extract, from the seismic volume, a matrix data item, wherein the matrix data item comprises seismic data surrounding the voxel in the seismic volume, wherein the matrix data item is extracted from the seismic volume based on a value of the estimated dip field surrounding the voxel; and
generate, using a matrix operator, a plurality of amplitude gradients from the seismic data in the matrix data item;
normalize, using the plurality of amplitude gradients, an amplitude value for the voxel to produce a normalized amplitude value for the voxel; and
generate, using the normalized amplitude value, an image representing a portion of the subterranean formation.

8. The system of claim 7 ,
wherein the matrix operator comprises a three dimensional (3D) operator,
wherein the matrix data item comprises a 3 by 3 by 3 matrix represented by Data(x, y, z) for x, y, z in {−1, 0, 1},
wherein the voxel comprises the seismic data, represented by Input(i, j, k), at position (i, j, k) of a voxel space,
wherein the matrix data item is extracted based on Data(x, y, z)=Input(i+x, j+y, k+z+x*dip IL (i, j, k)+y*dip XL (i, j, k)),
wherein dip IL (i, j, k) represents an inline dip of the estimated dip field at the position (i, j, k) of the voxel space, and
wherein dip XL (i, j, k) represents a crossline dip of the estimated dip field at the position (i, j, k) of the voxel space.

9. The system of claim 7 , wherein the matrix operator comprises a Sobel operator.

10. The system of claim 7 , wherein normalizing the amplitude value further comprises:
adjusting, based on pre-determined weighting factors, contributions of first, second, and third partial derivatives to the normalized amplitude value along a perpendicular direction of the stratigraphic layer,
wherein the matrix operator comprises first, second, and third matrix operators calculate the first, second, and third partial derivatives, respectively, along three orthogonal directions.

11. The system of claim 7 ,
wherein the normalized amplitude value is represented by
S =√{square root over ( G X 2 +G Y 2 +W Z 2 G Z 2 )},

wherein G x , G y , and G z are proportional to the first, second, and third partial derivatives, respectively, and
wherein W z is a pre-determined fraction.

12. The system of claim 11 ,
wherein W z is in the range of approximately [0, 0.4].

13. A non-transitory computer readable medium storing instructions for visualizing a subterranean formation, the instructions when executed causing a processor to:
obtain an estimated dip field of the subterranean formation, wherein the estimated dip field represents a measure of deviation of a stratigraphic layer from flat;
select a voxel from a seismic volume;
extract, from the seismic volume, a matrix data item, wherein the matrix data item comprises seismic data surrounding the voxel in the seismic volume, wherein the matrix data item is extracted from the seismic volume based on a value of the estimated dip field surrounding the voxel;
generate, using a matrix operator, a plurality of amplitude gradients from the seismic data in the matrix data item;
normalize, using the plurality of amplitude gradients, an amplitude value for the voxel to produce a normalized amplitude value for the voxel; and
generate, using the normalized amplitude value, an image representing a portion of the subterranean formation.

14. The non-transitory computer readable medium of claim 13 ,
wherein the matrix operator comprises a three dimensional (3D) operator,
wherein the matrix data item comprises a 3 by 3 by 3 matrix represented by Data(x, y, z) for x, y, z in {−1, 0, 1},
wherein the voxel comprises the seismic data, represented by Input(i, j, k), at position (i, j, k) of a voxel space,
wherein the matrix data item is extracted based on Data(x, y, z)=Input(i+x, j+y, k+z+x*dip IL (i, j, k)+y*dip XL (i, j, k)),
wherein dip IL (i, j, k) represents an inline dip of the estimated dip field at the position (i, j, k) of the voxel space, and
wherein dip XL (i, j, k) represents a crossline dip of the estimated dip field at the position (i, j, k) of the voxel space.

15. The non-transitory computer readable medium of claim 13 , wherein the matrix operator comprises a Sobel operator.

16. The non-transitory computer readable medium of claim 13 , wherein normalizing the amplitude value further comprises:
adjusting, based on pre-determined weighting factors, contributions of first, second, and third partial derivatives to the normalized amplitude value along a perpendicular direction of the stratigraphic layer,
wherein the matrix operator comprises first, second, and third matrix operators calculate the first, s

Assignees

• Aqrawi Ahmed Adnan
• Bø Trond Hellem
• Schlumberger Technology Corp

Inventors

• Aqrawi Ahmed Adnan
• Bø Trond Hellem

Classifications

• G01V1/307

  for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity · CPC title

• G01V1/325 Primary

  Transforming one representation into another · CPC title

• G01V2210/74

  Visualisation of seismic data · CPC title