Method for enhancing the resolution of a seismic image

Patent metadata

Field	Value
Publication number	US-2016124102-A1
Application number	US-201414896618-A
Country	US
Kind code	A1
Filing date	Jun 6, 2014
Priority date	Jun 7, 2013
Publication date	May 5, 2016
Grant date	—

Abstract

A method for enhancing a physical parameter map in a zone of a seismic image. The dip of points of the image is obtained. For one of these points, called second point, a correction factor of a physical parameter is obtained with a residual move-out algorithm from a common image gather. A first point is selected on a line substantially perpendicular to the dip at the second point. The selection, involves at least one parameter among whether the difference between the dip at the second point and the dip at the first point is below a first preset value; and the spacing between the first and the second point is below a second preset value. An inversion algorithm gives a corrected interval value of the physical parameter to update the physical parameter map.

First claim

1 . A method for enhancing a physical parameter map in a zone of a seismic image, said method comprising:
/a/ obtaining the zone of the seismic image, said zone comprising a set of points; /b/ obtaining a dip at each point in the set of points; /c/ selecting at least a first point in the set of points, and a second point in the set of points, said selection being based on at least one parameter among:
a difference between a dip at said second point and a dip at the first point is below a first preset value; and
a spacing between the first point and the second point is below a second preset value;
/d/ obtaining a common image gather associated with the second point in said set of points; /e/ applying a residual move-out analysis algorithm to the common image gather to obtain a correction factor of a physical parameter at said second point; /f/ applying an inversion algorithm to points between the first point and the second point in the set of points and along a line joining said first point and second point, to obtain an interval value of said physical parameter to update the physical parameter.

2 . The method of claim 1 wherein the first preset value is based on a scalar product between a unitary vector perpendicular to a dip of the first point and a unitary vector perpendicular to a dip of the second point.

3 . The method of claim 1 , wherein the second preset value is a value above which a ray normal to the first point misses the second point by more than 5% of the distance separating the first point from the second point.

4 . The method of claim 1 s, wherein the selection of a second point in step /f/ is further based on a parameter expressed as:
for a distance between the first point and points of said set of points smaller than a third preset value, a dip gradient is below a fourth preset value.

5 . The method of claim 1 , wherein a common image gather associated with the first point in said set of points is obtained and a correction factor of the physical parameter at said second point is obtained by applying a residual move-out analysis algorithm to the common image gather associated with the first point, step /f/ using a correction factor of a physical parameter at said first point and the correction factor of the physical parameter at said second point to obtain an interval value of said physical parameter to update the physical parameter map.

6 . The method of claim 1 , wherein each point of said set of points is associated to x, y, z coordinates and to a physical parameter v(x,y,z) at said point of coordinates x, y, z, the method further comprising introducing a weighting function W defining a weighting factor w(x,y,z) for each point of the set of points, a physical parameter profile Vp(x,y,z) being obtained at each point of coordinates x, y, z, wherein:
Vp(x,y,z)=V(x,y,z)*W(x,y,z)/ΣW(x,y,z)

where ΣW(x,y,z) is a sum of all weighting factors of the set of points.

7 . The method of claim 1 , wherein an enhanced migration velocity analysis comprising a predetermined number of iterations is implemented, said enhanced migration velocity analysis comprising:
obtaining multi-offset seismic data; and for each current iteration of said predetermined number of iteration; /a′/transforming multi-offset seismic data into a depth-migrated seismic image using a physical parameter map; implementing steps /a/ to /f/.

8 . The method of claim 7 , wherein a of the first point as a function of offset is determined prior to a current iteration, and the enhanced migration velocity analysis is continued until said depth as a function of offset differs from a depth at zero-offset by less than a fifth preset value.

9 . The method of claim 1 , wherein a first physical parameter map is estimated based on a migration velocity analysis applied to multi-offset seismic data.

10 . The method of claim 1 , wherein a seismic ray path is a zero-offset seismic ray path.

11 . The method of claim 1 , wherein the inversion algorithm is based on a Dix inversion algorithm.

12 . The method of claim 1 , wherein the physical parameter is a seismic wave velocity.

13 . The method of claim 1 , wherein the physical parameter is an elasticity coefficient of a medium through which seismic waves propagate.

14 . A non-transitory computer readable storage medium, having stored thereon a computer program comprising program instructions, the computer program being loadable into a data-processing unit and adapted to cause the data-processing unit to carry out a method for enhancing a physical parameter map in a zone of a seismic image, said method comprising:
/a/ obtaining the zone of the seismic image, said zone comprising a set of points; /b/ obtaining a dip at each point in the set of points: /c/ selecting at least a first point in the set of points, and a second point in the set of points, said selection being based on at least one parameter among:
a difference between a dip at said second point and a dip at the first point is below a first preset value; and
a spacing between the first point and the second point is below a second preset value;
/d/ obtaining a common gather associated with the second point in said set of points; /e/ applying a residual move-out analysis algorithm to the common image gather to obtain a correction factor of a physical at said secondpoint; /f/ applying an inversion algorithm to points between the first point and the second point in the set of points and along a line joining said first point and second point to obtain an interval value of said physical parameter to update the physical parameter map, when the computer program is run by the data-processing device.

Assignees

• Total Sa

Inventors

• Hardy Pierre

Classifications

• G01V2210/522

  Dip move-out [DMO] · CPC title

• G01V1/325

  Transforming one representation into another · CPC title

• G01V1/362 Primary

  Effecting static or dynamic corrections; Stacking · CPC title

• G01V2210/52

  Move-out correction · CPC title

• G01V2210/48

  Other transforms · CPC title