Determination of horizontal constraints in subsoil

1. A method of determining a plurality of horizontal constraints in a subsoil, wherein the method comprises:

receiving well data of the subsoil;

dividing the well data in a plurality of contiguous sets of data i respectively associated with a plurality of layers having a constant or variable vertical depth increment, a set of data i being above a set of data i+1;

for each set of data i in the plurality of sets of data, determining, based on the set of data i, at least parameters:

P i , a pressure in subsoil corresponding to the set of data i;

b i , a Biot coefficient corresponding to the set of data i;

ν i , a Poisson's ratio corresponding to the set of data i;

σ ν i , a vertical constraints in subsoil corresponding to the set of data i;

E i , a Young's modulus corresponding to the set of data i;

α i , a rock thermal expansion coefficient corresponding to the set of data i;

T i , a temperature in subsoil corresponding to the set of data i;

for each set of data i in the plurality of sets of data, computing and outputting the horizontal constraints based on an equation:

σ

h

i

=

σ

h

i

-

1

+

b

i

⁡

(

P

i

-

P

i

-

1

)

+

ν

i

1

-

ν

i

⁢

(

σ

v

i

-

σ

v

i

-

1

-

b

i

⁡

(

P

i

-

P

i

-

1

)

)

+

E

i

⁢

α

i

1

-

ν

i

⁢

(

T

i

-

T

i

-

1

)

;

wherein the horizontal constraints are computed according to an ascending order of i assuming that a history of sedimentation and compaction of layer i at present-day at a depth can be represented numerically by a present-day's porosity and density evolution from seabed to the depth, even if layers have different mineral composition compared to that of layer i; and

using the horizontal constraints in estimating a high of a hydrocarbon column, to assess a geomechanical hazards link to reservoir depletion, and for a design of well architecture or well integrity control for petroleum engineering.

2. The method according to claim 1 , wherein determining the parameters comprises:

receiving a porosity value Φ in subsoil corresponding to the set of data i, and a proportion of clay V cl in subsoil corresponding to the set of data i;

computing aspect ratio of clay mineral

d

l

=

λ

⁡

(

Φ

κ

⁢

⁢

V

cl

)

n

 where

d and l are respectively a thickness and a diameter of spheroidal inclusion;

λ, n and κ are three parameters varying respectively in a range of 0.05˜0.06 for λ, 2˜2.5 for n and 0.15 to 0.25 for κ, according to a geological context and applicable for clay content V ct varying in an interval 0.05 to 0.5.

3. 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 a data-processing device.