Method of simulating the sedimentary deposition in a basin respecting the thicknesses of the sedimentary sequences

• US 7,337,069 B2
• Filed: 06/24/2005
• Issued: 02/26/2008
• Est. Priority Date: 06/30/2004
• Status: Active Grant

• Current Assignee: Institut Francais Du Petrole
• Original Assignee: Institut Francais Du Petrole
• Inventors: Cacas, Marie-Christine, Dobranszky, Gabriela, Masson, Roland
• Primary Examiner(s): Barlow; John
• Assistant Examiner(s): Le; Toan M.
• Application Number:
  US11/165,264
• Publication Number:
  US 20060015260A1
• Time in Patent Office: 977 Days
• Field of Search: 702/14, 702/2, 702/5, 702 11- 13, 702/17, 703/6, 703/10, 367/73
• US Class Current: 702/14
• CPC Class Codes: G01V 1/282 Application of seismic models, synthetic seismograms
• CPC Class Codes: G01V 11/00 Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
• CPC Class Codes: G01V 2210/661 Model from sedimentation process modeling, e.g. from first principles

Abstract

The method of the invention allows, by means of an iterative inversion algorithm, prediction of the spatial distribution of the lithologic composition of sediments deposited in a sedimentary basin during a geologic time interval, and the temporal evolution of the depositional profile throughout filling of the basin, while respecting exactly the thicknesses of the sedimentary sequences measured otherwise which is useful for locating hydrocarbon reservoirs. Input data of thickness maps of the sedimentary layer, location and composition of the sediment supply at the boundaries of the sedimentary basin, and physical parameters characterizing transport of the sediments during the period of consideration are applied to an iterative inversion loop initialized by the accommodation provided by a stationary multilithologic diffusive model which works as a fixed-point algorithm correcting the accommodations by means of a preconditioning of the residue on the sequence thicknesses obtained by an instationary model.

Claims

1. A method for forming a model of a spatial distribution of a lithologic composition of sediments deposited in a sedimentary basin having sedimentary sequences, and a temporal evolution of bathymetry, comprising:

estimating, from interpretations of measurements resulting from seismic surveys and well data, input data comprising;

maps giving current thicknesses of the sedimentary sequences;

a map of initial topography of a basement of the sedimentary basin;

a distribution of a sedimentary supply in the basin, comprising at least descriptions of location in space, lithologic composition, and evolution of supply volume in time during filling of the basin; and

diffusion coefficient values characterizing transport of the sediments to the basin; and

establishing the three-dimensional map by applying the input data to a first iterative inversion procedure, comprising;

calculating a thickness of each of the sequences by solving instationary diffusive transport equations;

updating an accommodation at an end of each of the sequences by applying a second iterative procedure, comprising calculating a residue between current thicknesses of the sedimentary sequences and calculated thicknesses, and correcting the accommodation by preconditioning the residue until a sufficiently low residue is obtained on the thicknesses; and

locating hydrocarbon reservoirs from the three-dimensional map.

2. A method as claimed in claim 1, wherein the accommodation at the end of a given sequence is initialized by the accommodation at the end of the sequence provided by a stationary diffusive model.

3. A method as claimed in claim 1, wherein preconditioning is obtained by application to the residue, on the thicknesses, of a tangent stationary diffusive model.

4. A method as claimed in claim 1, wherein an initial accommodation of a sequence is determined by a final accommodation calculated for the sequence located immediately below a sequence which is considered.

5. A method as claimed in claim 1, wherein mean values in time of the sedimentary supply for a sequence are considered in stationary and tangent stationary models.

6. A method as claimed in claim 1 wherein, when close to convergence, updating the accommodation at the end of each of the sequences is obtained by calculation of a Newton step of an inverse problem.

7. A method as claimed in claim 2, wherein preconditioning is obtained by application to the residue, on the thicknesses, of a tangent stationary diffusive model.

8. A method as claimed in claim 7, wherein an initial accommodation of a sequence is determined by a final accommodation calculated for the sequence located immediately below a sequence which is considered.

9. A method as claimed in claim 8, wherein mean values in time of the sedimentary supply for a sequence are considered in stationary and tangent stationary models.

10. A method as claimed in claim 7, wherein mean values in time of the sedimentary supply for a sequence are considered in stationary and tangent stationary models.

11. A method as claimed in claim 2, wherein an initial accommodation of a sequence is determined by a final accommodation calculated for the sequence located immediately below a sequence which is considered.

12. A method as claimed in claim 11, wherein mean values in time of the sedimentary supply for a sequence are considered in stationary and tangent stationary models.

13. A method as claimed in claim 2, wherein mean values in time of the sedimentary supply for a sequence are considered in the stationary and tangent stationary models.

14. A method as claimed in claim 2 wherein, when close to convergence, updating the accommodation at the end of each of the sequences is obtained by calculation of Newton step of an inverse problem.

15. A method as claimed in claim 3, wherein an initial accommodation of a sequence is determined by a final accommodation calculated for the sequence located immediately below a sequence which is considered.

16. A method as claimed in claim 15, wherein mean values in time of the sedimentary supply for a sequence are considered in stationary and tangent stationary models.

17. A method as claimed in claim 3, wherein mean values in time of the sedimentary supply for a sequence are considered in stationary and tangent stationary models.

18. A method as claimed in claim 3 wherein, when close to convergence, updating the accommodation at the end of each of the sequences is obtained by calculation of Newton step of an inverse problem.

19. A method as claimed in claim 4, wherein mean values in time of the sedimentary supply for a sequence are considered in stationary and tangent stationary models.

20. A method as claimed in claim 6, wherein a Newton step of the inverse problem is calculated by means of a generalized minimal residual algorithm preconditioned by a tangent stationary model.