Method for estimating elastic parameters by inverting 4D seismic measurements

Patent metadata

Field	Value
Publication number	US-9377549-B2
Application number	US-201113825287-A
Country	US
Kind code	B2
Filing date	Sep 16, 2011
Priority date	Sep 20, 2010
Publication date	Jun 28, 2016
Grant date	Jun 28, 2016

Abstract

According to the invention, elastic parameters, including density, pressure wave propagation speed, and/or shear in pervious layers located in a dense underground area along an array of horizontal positions, are estimated by inverting 4D seismic data. Firstly, an estimate of variations in the elastic parameters in one or more starting positions of the array, which can be located on bored wells within the area in question, is obtained. Then, a propagation algorithm is used in order to gradually carry out the 4D data inversion on the basis of the starting positions. The inversion takes into account the previously estimated parameter variations. A spatial variation in the depth and/or the thickness of the pervious lavers in question can also be taken into account. Propagation is based on positions that are consecutively selected as providing optimal values for a cost function assessed in order to invert the 4D data.

First claim

The invention claimed is:

1. A method for estimating elastic parameters of a region of the subsoil meshed according to an array of horizontal positions, the method comprising:
(a) obtaining seismic traces associated with said horizontal positions;
(b) obtaining an estimation of variations of elastic parameters in layers of the subsoil in at least one starting position of the array;
(c) determining a boundary consisting of the positions of the array which are adjacent to a starting position;
(d) obtaining, at each position of the boundary, estimations of variations of the elastic parameters in the layers by minimizing a cost function derived from a seismic wave propagation model accounting for the seismic traces associated with said position of the boundary, the minimization taking into account the estimations obtained at a starting position adjacent to said position of the boundary, and marking each position of the boundary as having been processed;
(e) selecting at least one propagation position in the boundary, where the minimized cost function has, within the boundary, a minimal value, then removing the selected propagation position from the boundary;
(f) adding to the boundary each position of the network not marked as processed and adjacent to a propagation position selected in the preceding step (e);
(g) obtaining, at each position added to the boundary in the preceding step (f), estimations of variations of the elastic parameters in the layers by minimizing the cost function derived from the seismic wave propagation model accounting for the seismic traces associated with said position added to the boundary, the minimization taking into account the estimations obtained at the selected propagation position, and marking each position added to the boundary in the preceding step (f) as having been processed;
(h) repeating from step (e).

2. The method as claimed in claim 1 , wherein the seismic traces obtained in step (a) include seismic traces measured by sending seismic waves at normal incidence toward layers of the subsoil and by collecting the seismic waves reflected by interfaces between said layers.

3. The method as claimed in claim 2 , wherein the elastic parameters include a density and a speed of propagation of pressure waves in at least some of said layers.

4. The method as claimed in claim 1 , wherein the seismic traces obtained in step (a) include seismic traces measured by sending seismic waves at non-normal incidence toward layers of the subsoil and by collecting the seismic waves reflected by the interfaces between said layers.

5. The method as claimed in claim 4 , wherein the elastic parameters include a speed of propagation of shear waves in at least some of said layers along said direction.

6. The method as claimed in claim 1 , wherein at least one starting position of the array is given by the position of a well drilled in the subsoil.

7. The method as claimed in claim 1 , wherein each of steps (d) and (g) is followed by a step of removing from the boundary any position where the minimized cost function exceeds a limit value.

8. The method as claimed in claim 1 , wherein the repeating step (h) is performed if the boundary contains at least one position adjacent to a position of the array not marked as processed.

9. The method as claimed in claim 1 , wherein step (a) of obtaining seismic traces associated with the horizontal positions comprises:
measuring, in a first phase, base seismic traces associated with said horizontal positions; and
measuring, in a second phase, monitor seismic traces associated with said horizontal positions and corresponding to the base seismic traces,
wherein the estimations of variations of elastic parameters obtained in steps (b), (d) and (g) comprise estimations of variations of the elastic parameters between the first and the second phases in layers of the subsoil,
and wherein said propagation model accounts for an evolution between the base and monitor seismic traces associated with a position of the boundary.

10. The method as claimed in claim 9 , wherein the elastic parameters having variations estimated in steps (d) and (g) are elastic parameters in permeable layers of the subsoil.

11. The method as claimed in claim 10 , wherein steps (d) and (g) further comprise obtaining estimations of variations of depths of said permeable layers.

12. The method as claimed in claim 11 , wherein the permeable layers are located on the basis of a reservoir grid, and wherein the geometry of the reservoir grid is adjusted according to depths of layers estimated in step (d) or (g).

13. The method as claimed in claim 10 , wherein steps (d) and (g) further comprise obtaining estimations of variations of thicknesses of said permeable layers.

14. The method as claimed in claim 13 , wherein the permeable layers are located on the basis of a reservoir grid, and wherein the geometry of the reservoir grid is adjusted according to thicknesses of layers estimated in step (d) or (g).

15. The method as claimed in claim 9 , wherein a number of assumptions as to the variation of the elastic parameters in permeable layers of the subsoil are defined at a position of the boundary based on the variations estimated at the starting or propagation position adjacent to said position of the boundary, and wherein the minimization of the cost function in step (d) or (g) of estimation, at said position of the boundary, of the variations of the elastic parameters in the permeable layers comprises:
for each of said assumptions:
estimating a perturbation in amplitude of the measured base seismic trace associated with said position of the boundary following the variation of the elastic parameters conforming to said assumption;
computing a seismic pseudo-trace by combining one of the measured seismic traces associated with the position of the boundary with the estimated amplitude perturbation; and
comparing, on a same timescale, the other measured seismic trace associated with the position of the boundary with the computed seismic pseudo-trace in order to evaluate the cost function;
estimating the variation of the elastic parameters in said permeable layers in accordance with an assumption for which the evaluated cost function is minimal.

16. The method as claimed in claim 15 , wherein the perturbation in amplitude is estimated as a function of impedance variations in said permeable layers, deduced from the assumption of variation of the elastic parameters, and from a wavelet representative of an incident seismic signal.

17. The method as claimed in claim 15 , wherein the measured seismic trace combined with the estimated amplitude perturbation for the computation of the seismic pseudo-trace is the measured monitor seismic trace, the combination comprising:
modifying a timescale of the measured monitor seismic trace to take into account the variation in speed of propagation of the seismic waves according to the assumption of variation of elastic parameters; and
subtracting the estimated amplitude perturbation from the modified measured monitor seismic trace.

18. The method as claimed in claim 15 , wherein the assumptions defined at a position of the boundary based on the variations estimated at the starting or propagation position adjacent to said position of the boundary further comprise assumptions as to the variation of depths of said permeable layers, and wherein the minimization of the cost function in step (d) or (g) further comprises, for each assumption, estimating a variation of the depths of said permeable layers in accordance

Assignees

• Thore Pierre
• Hubans Christian
• Total Sa

Inventors

• Thore Pierre
• Hubans Christian

Classifications

• G01V1/306 Primary

  for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles · CPC title

• G01V1/308

  Time lapse or 4D effects, e.g. production related effects to the formation (fluid flow per se E21B47/00) · CPC title

• G01V1/28

  Processing seismic data, e.g. for interpretation or for event detection (G01V1/48 takes precedence) · CPC title