Method for estimating elastic parameters through inversion of 4D seismic measures

Patent metadata

Field	Value
Publication number	US-9690001-B2
Application number	US-201113811378-A
Country	US
Kind code	B2
Filing date	Jul 19, 2011
Priority date	Jul 21, 2010
Publication date	Jun 27, 2017
Grant date	Jun 27, 2017

Abstract

A 4D seismic technique, where a base seismic trace is measured at a first time in a region of the subsoil, and then a monitor seismic trace corresponding to the base seismic trace is measured at a second time. To interpret the 4D measurements, assumptions are made about the variation of elastic parameters in permeable layers at predefined positions in one direction between the first and the second time. Elastic parameters include the density (ρ) and the speed of propagation of the pressure waves (V P ) in the permeable layers. Numerical evaluation is performed for a capability of each assumption about the variation of elastic parameters to give an account of a change between the measured base seismic trace and the measured monitor seismic trace, and the variation of the elastic parameters is estimated in accordance with an assumption of optimum capability.

First claim

The invention claimed is:

1. A method for estimating elastic parameters of a region of the subsoil, comprising:
measuring at least one base seismic trace in a first phase;
measuring, in a second phase, at least one monitor seismic trace corresponding to the base seismic trace;
making assumptions of variation of elastic parameters in permeable layers of predefined positions along a direction between the first and the second phases, the elastic parameters including a density (φ and a speed of propagation of pressure waves (V P ) in each of said permeable layers;
numerically evaluating a capacity of each assumption of variation of elastic parameters to account for an evolution between the measured base seismic trace and the measured monitor seismic trace; and
estimating the variation of the elastic parameters in accordance with an assumption having an optimum evaluated capacity.

2. The method as claimed in claim 1 , wherein the base and monitor seismic traces are measured by sending seismic waves at normal incidence toward layers following one another along said direction and by collecting the seismic waves reflected by interfaces between said layers.

3. The method as claimed in claim 1 , wherein the base and monitor seismic traces are measured by sending seismic waves at non-normal incidence toward layers following one another along said direction and by collecting the seismic waves reflected by the interfaces between said layers, and wherein the elastic parameters also include a speed of propagation of shear waves (V S ) in the permeable layers along said direction.

4. The method as claimed in claim 1 , wherein the elastic parameters further include the position, along said direction, of at least one interface delimiting one of said permeable layers.

5. The method as claimed in claim 1 , wherein the variations of elastic parameters are taken into account in permeable layers along a well drilled in the subsoil.

6. The method as claimed in claim 5 , wherein the positions of the permeable layers along the well are defined on the basis of measurements performed in the well.

7. The method as claimed in claim 5 , wherein the positions of the permeable layers along the well are defined on the basis of positions of perforations made in a casing of the well.

8. The method as claimed in claim 5 , wherein the numerical evaluation of the capacity of an assumption of variation of elastic parameters to account for the evolution between the measured base seismic trace and the measured monitor seismic trace comprises:
computing a simulated base seismic trace on the basis of a wavelet representative of an incident seismic signal and of values of the elastic parameters measured in the well in the first phase;
obtaining values of the elastic parameters in the well for the second phase on the basis of the values measured in the well in the first phase and of said assumption of variation;
computing a simulated monitor seismic trace on the basis of said wavelet and of the values of the elastic parameters in the well obtained for the second phase; and
comparing a difference between the measured monitor seismic trace and the measured base seismic trace with a difference between the simulated monitor seismic trace and the simulated base seismic trace.

9. The method as claimed in claim 5 , wherein the numerical evaluation of the capacity of an assumption of variation of elastic parameters to account for the evolution between the measured base seismic trace and the measured monitor seismic trace comprises:
obtaining values of the elastic parameters in the well for the second phase on the basis of values of the elastic parameters measured in the well in the first phase and of said assumption of variation;
estimating a perturbation in amplitude of the measured base seismic trace resulting from a change of the elastic parameters in the well from the values measured in the first phase to the values obtained for the second phase;
computing a seismic pseudo-trace by combining one of the measured seismic traces with the estimated amplitude perturbation; and
comparing, on a same timescale, the other one of the measured seismic trace with the computed seismic pseudo-trace.

10. The method as claimed in claim 9 , wherein the estimation of the amplitude perturbation comprises:
computing a simulated base seismic trace on the basis of a wavelet representative of an incident seismic signal and of values of the elastic parameters measured in the well in the first phase;
computing a simulated monitor seismic trace on the basis of said wavelet and of values of the elastic parameters in the well obtained for the second phase, the simulated base and monitor seismic traces being computed with a same depth-time conversion law; and
subtracting the simulated base seismic trace from the simulated monitor seismic trace to obtain the estimated amplitude perturbation.

11. The method as claimed in claim 9 , 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 the timescale of the measured monitor seismic trace to take into account a variation of speed of propagation of the seismic waves in accordance with the assumption of variation of elastic parameters; and
subtracting the estimated amplitude perturbation from the modified measured monitor seismic trace.

12. The method as claimed in claim 1 , wherein the numerical evaluation of the capacity of an assumption of variation of the elastic parameters to account for the evolution between the measured base seismic trace and the measured monitor seismic trace comprises:
estimating a perturbation in amplitude of the measured base seismic trace resulting from the variation of the elastic parameters;
computing a seismic pseudo-trace by combining one of the measured seismic traces with the estimated amplitude perturbation; and
comparing, on a same timescale, the other one of the measured seismic traces with the computed seismic pseudo-trace.

13. The method as claimed in claim 12 , wherein the amplitude perturbation 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.

14. The method as claimed in claim 12 , 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 the timescale of the measured monitor seismic trace to take into account a variation of speed of propagation of the seismic waves in accordance with the assumption of variation of elastic parameters; and
subtracting the estimated amplitude perturbation from the modified measured monitor seismic trace.

15. The method as claimed in claim 1 , wherein the positions of the permeable layers along said direction are defined on the basis of a reservoir grid.

Assignees

• Thore Pierre
• Hubans Christian
• Total Sa

Inventors

• Thore Pierre
• Hubans Christian

Classifications

• G01V1/308 Primary

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

• G01V2210/612

  Previously recorded data, e.g. time-lapse or 4D · CPC title

• G01V2210/6242

  Elastic parameters, e.g. Young, Lamé or Poisson · CPC title

• G01V1/306 Primary

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

• G01V1/30 Primary

  Analysis (G01V1/50 takes precedence) · CPC title