Method for determining a lithology map

1. A method for determining a lithology map, wherein the method comprises:

receiving a geophysical attribute image, the image comprising a plurality of points, each point in the plurality of points being associated with n geophysical attributes, each attribute of the n geophysical attributes having a respective value, n being an integer greater than 0, wherein for each point in the plurality of points, the geophysical attributes having an order index;

receiving first information data representing a plurality of zones in the image, each point of the image being contained in a zone;

for each integer j in [1; n]:

for at least one current point (x) of the image, computing a uniform value (F i (Y(x))) associated with the current point (x) based on:

a j th geophysical attribute value (Y(x)), wherein the j th geophysical attribute value comprises at least one value in a group comprising a seismic value, an acoustic impedance value and a pseudo V-Clay value, and

distributions values (f i (y)) of the j th geophysical attribute value associated with points of the zone containing the current point;

receiving second information data representing geological probability values for a plurality of facies, said geological probability values being associated with the points of the image;

for at least one current point (x) of the image where a uniform value is computed, determining a determined facies in the plurality of facies based on:

n computed uniform values computed for each geophysical attribute and associated with the current point, and

the geological probability values associated with said current point;

providing a final facies repartition image based on said determined facies;

estimating an oil/gas reserve in a subsoil from said final facies repartition image; and wherein, the zones represent at least one zone in a group comprising an environment of deposition zone (ECM), an architectural elements zone (AE), and a large geological body.

2. The method according to claim 1 , wherein the computation of uniform values is performed for all points of the images.

3. The method according to claim 1 , wherein the distribution value of a value of the m th geophysical attribute in a zone, m being an integer in [1; n], is computed based on:

a number of points in said zone associated with said value of a m th geophysical attribute; and

a number of points in said zone.

4. The method according to claim 1 , wherein the distribution value of a first value of a p th geophysical attribute in a zone, p being an integer in [1; n], is computed based on:

a number of points in said zone associated with a value of the p th geophysical attribute, a distance between the value of the p th geophysical attribute and the first value being lower than a predetermined threshold value; and

a number of points in said zone.

5. The method according to claim 1 , wherein computing a uniform value associated with a point comprises:

determining a first sum, the first sum being a finite sum of distribution values of geophysical attribute values equal or below the geophysical attribute value associated with said point in a zone containing said point; and

determining a second sum, the second sum being a finite sum of distribution values of a geophysical attribute values in a zone containing said point.

6. The method according to claim 5 , wherein the uniform value is function of the division of the first sum by the second sum.

7. The method according to claim 1 , wherein the method further comprises:

receiving well data representing geological information for a plurality of facies associated with well points of the image;

for each facies, determining modified geological probability values associated with points of the image based on the well data for said facies and the received second information data representing geological probability values for said facies;

and wherein, for at least one point of the image, determining a facies is based on the computed uniform value for each j in [1; n], and the modified geological probabilities associated with said point.

8. The method according to claim 7 , wherein determining modified geological probability values for a facies comprises performing an exact interpolation process based on the well data for said facies and the received second information data representing geological probability values for said facies.

9. The method according to claim 7 , wherein determining a facies uses a truncation diagram.

10. The method according to claim 1 , wherein n is equal to 2.

11. 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 the steps of claim 1 when the computer program is run by the data-processing device.

12. A device for determining a lithology map, wherein the device comprises:

an interface for receiving a geophysical attribute image, the image comprising a plurality of points, each point in the plurality of points being associated with n geophysical attributes, each attribute in the n geophysical attributes having a respective value, n being an integer greater than 0, wherein for each point in the plurality of points, the geophysical attributes are ordered;

an interface for receiving first information data representing a plurality of zones in the image, each point of the image being contained in a zone;

for each j in [1; n]:

for at least one current point (x) of the image, a circuit for computing a uniform value (F i (Y(x))) associated with the current point (x) based on:

a j th geophysical attribute value (Y(x)), and wherein the j th geophysical attribute value comprises at least one value in a group comprising a seismic value, an acoustic impedance value and a pseudo V-Clay value, and

distributions values (f i (y)) of the j th geophysical attribute value associated with points of the zone containing the current point;

an interface for receiving second information data representing geological probability values for a plurality of facies, said geological probability values being associated with the points of the image;

for at least one current point (x) of the image where a uniform value is computed, a circuit for determining a determined facies in the plurality of facies based on:

n computed uniform values computed for each geophysical attribute and associated with the current point, and

the geological probability values associated with said current point;

an interface to provide a final facies repartition image based on said determined facies

a circuit to estimate oil/gas reserves in a subsoil from said final facies repartition image,

and wherein, the zones represents at least one zone in a group comprising an environment of deposition zone (EOD), an architectural elements zone (AE) and large geological body.