Granted Patent B2

US 12,624,615 · App. 17/848,457 · Granted May 12, 2026

Inventor: Didier Yu Ding (Rueil-Malmaison Cedex, FR)

Assignee: IFP ENERGIES NOUVELLES

E21B41/00E21B2200/20

Abstract

Claims (595)

1 . A computer-implemented method of simulating fluid flows in an underground formation comprising a fracture network for exploiting the fluid of the underground formation wherein, from measured properties relative to the formation, a grid representation of the formation is constructed and at least one statistical parameter relative to the fracture network is determined, comprising steps of:

A) constructing from the at least one statistical parameter relative to the fracture network, a porosity model for the underground formation comprising the fracture network, the porosity model comprising a first medium representative of an unfractured matrix of the formation, a second medium representative of fractures of the formation oriented in a first direction, and a third medium representative of fractures of the formation oriented in a second direction, the first and second directions being orthogonal to one another;

B) determining from at least the measurements of properties relative to the formation, from the at least one statistical parameter relative to the fracture network and from the porosity model flow parameters in each cell of the grid representation, the flow parameters comprising:

determining between two neighboring cells in the second direction, a convection transmissibility and a conduction transmissibility of the fluid in the first and second directions in the second medium, the convection and conduction transmissibilities of the fluid being zero in the second direction in the second medium;

determining between two neighboring cells in the first direction, a convection transmissibility and a conduction transmissibility of the fluid in the first and second directions in the third medium, the convection and conduction transmissibilities of the fluid being zero in the first direction in the third medium; and

determining within a cell of the grid representation, mass exchanges by convection, energy exchanges by conduction between the second and third mediums, between the first and second mediums, and between the first medium and the third medium; and

C) simulating from the grid representation and the flow parameters in each cell of the grid, the flows of the fluid in the underground formation comprising the fracture network by use of a flow simulator implementing the porosity model.

2 . A method as claimed in claim 1 , wherein the mass exchanges

F

w

,

i

ff

by convection between the second and third mediums in one of the cells i of the grid are determined with a formula:

F

w

,

i

ff

=

λ

i

f

⁢

T

i

ff

(

P

i

fX

-

P

i

fY

)

wherein

T

i

ff

is a convection transmissibility between the second and third mediums of the cell i,

P

i

fX

and

P

i

fY

correspond to a pressure in the second and third mediums of the cell i respectively, and

λ

i

f

is mobility of the fluid between the second and third mediums of the cell i.

3 . A method as claimed in claim 1 , wherein the energy exchanges

F

H

,

i

ff

by convection between the second and third mediums in one of the cells i of the grid are determined with a formula:

F

H

,

i

ff

=

λ

i

f

⁢

H

i

f

⁢

T

i

ff

(

P

i

fX

-

P

i

fY

)

wherein

H

i

f

is a enthalpy of the fluid between the second and third mediums of the cell i,

T

i

ff

is the convection transmissibility between the second and third mediums of the cell i,

P

i

fX

and

P

i

fY

correspond to the pressure in the second and third mediums of the cell i respectively, and

λ

i

f

is the mobility of the fluid between the second and third mediums of the cell i.

4 . A method as claimed in claim 2 , wherein energy exchanges

F

H

,

i

ff

by convection between the second and third mediums in one of the cells i of the grid are determined with a formula:

F

H

,

i

ff

=

λ

i

f

⁢

H

i

f

⁢

T

i

ff

(

P

i

fX

-

P

i

fY

)

wherein

H

i

f

is a enthalpy of the fluid between the second and third mediums of the cell i,

,

T

i

ff

is the convection transmissibility between the second and third mediums of the cell i,

P

i

fX

and

P

i

fY

correspond to the pressure in the second and third mediums of the cell i respectively, and

λ

i

f

is the mobility of the fluid between the second and third mediums of the cell i.

5 . A method as claimed in claim 2 , wherein the convection transmissibility between the second and third mediums of the cell i is determined with a formula:

T

i

ff

=

α

⁢

∑

j

∈

Ω

i

⁢

(

T

x

,

ij

fX

+

T

y

,

ij

fY

)

where α is a multiplier equal to at least 100, Ω i is all cells next to cell i,

T

x

,

ij

fX

corresponds to the convection transmissibility in the first direction in the second medium for the cell i, and

T

y

,

ij

fY

corresponds to the convection transmissibility in the second direction in the third medium for the cell i.

6 . A method as claimed in claim 3 , wherein the convection transmissibility between the second and third mediums of the cell i is determined with a formula:

T

i

ff

=

α

⁢

∑

j

∈

Ω

i

⁢

(

T

x

,

ij

fX

+

T

y

,

ij

fY

)

where α is a multiplier equal to at least 100, Ω i is all cells next to cell i,

T

x

,

ij

fX

corresponds to the convection transmissibility in the first direction in the second medium for the cell i, and

T

y

,

ij

fY

corresponds to the convection transmissibility in the second direction in the third medium for the cell i.

7 . A method as claimed in claim 1 , wherein energy exchanges

F

D

,

i

ff

by conduction between the second and third mediums in the cell i are determined with a formula:

F

D

,

i

ff

=

G

i

ff

(

Θ

i

f

⁢

X

-

Θ

i

f

⁢

Y

)

where

G

i

ff

is conduction transmissibility between the second and third mediums for the cell i,

Θ

i

fX

and

Θ

i

f

⁢

Y

correspond to the temperature in the second and third mediums respectively.

8 . A method as claimed in claim 2 , wherein energy exchanges

F

D

,

i

ff

by conduction between the second and third mediums in the cell i are determined with a formula:

F

D

,

i

ff

=

G

i

ff

(

Θ

i

f

⁢

X

-

Θ

i

f

⁢

Y

)

where

Θ

j

ff

is conduction transmissibility between the second and third mediums for the cell i,

Θ

i

f

⁢

X

and

Θ

i

f

⁢

Y

correspond to the temperature in the second and third mediums respectively.

9 . A method as claimed in claim 3 , wherein energy exchanges

F

D

,

i

ff

by conduction between the second and third mediums in the cell i are determined with a formula:

F

D

,

i

ff

=

G

i

ff

(

Θ

i

fX

-

Θ

i

fY

)

where

G

i

ff

is conduction transmissibility between the second and third mediums for the cell i,

Θ

i

f

⁢

X

and

Θ

i

f

⁢

Y

respectively correspond to the temperature in the second and third mediums.

10 . A method as claimed in claim 4 , wherein energy exchanges

F

D

,

i

ff

by conduction between the second and third mediums in the cell i are determined with a formula:

F

D

,

i

ff

=

G

i

ff

(

Θ

i

fX

-

Θ

i

fY

)

where

G

i

ff

is conduction transmissibility between the second and third mediums for the cell i,

Θ

i

f

⁢

X

⁢

and

⁢

Θ

i

f

⁢

Y

respectively correspond to the temperature in the second and third mediums.

11 . A method as claimed in claim 5 , wherein energy exchanges

F

D

,

i

ff

by conduction between the second and third mediums in the cell i are determined with a formula:

F

D

,

i

ff

=

G

i

ff

(

Θ

i

fX

-

Θ

i

fY

)

where

G

i

ff

is conduction transmissibility between the second and third mediums for the cell i,

Θ

i

fX

respectively correspond to the temperature in the second and third mediums.

12 . A method as claimed in claim 6 , wherein energy exchanges

F

D

,

i

ff

by conduction between the second and third mediums in the cell i are determined with a formula:

F

D

,

i

ff

=

G

i

ff

(

Θ

i

fX

-

Θ

i

fY

)

where

G

i

ff

is conduction transmissibility between the second and third mediums for the cell i,

Θ

i

fX

and

Θ

i

fY

respectively correspond to the temperature in the second and third mediums.

13 . A method as claimed in claim 12 , wherein conduction transmissibility between the second and third mediums for the cell i is determined with a formula:

G

i

ff

=

Λ

i

ff

⁢

C

i

ff

where

Λ

i

ff

is an arithmetic mean of an effective thermal conductivity of the second and third mediums,

C

i

ff

is a geometric coefficient depending on the dimensions of the cell i, the dimensions of one of the matrix blocks into which the first medium is broken, and the opening of the fractures of the second and third mediums.

14 . A method as claimed in claim 8 , wherein the conduction transmissibility between the second and third mediums for the cell i is determined with a formula:

G

i

ff

=

Λ

i

ff

⁢

C

i

ff

where

Λ

i

ff

is an arithmetic mean of an effective thermal conductivity of the second and third mediums,

C

i

ff

is a geometric coefficient depending on the dimensions of the cell i, the dimensions of one of the matrix blocks into which the first medium is broken, and the opening of the fractures of the second and third mediums.

15 . A method as claimed in claim 9 , wherein conduction transmissibility between the second and third mediums for the cell i is determined with a formula:

G

i

ff

=

Λ

i

ff

⁢

C

i

ff

where

Λ

i

ff

is an arithmetic mean of an effective thermal conductivity of the second and third mediums,

C

i

ff

is a geometric coefficient depending on the dimensions of the cell i, the dimensions of one of the matrix blocks into which the first medium is broken, and the opening of the fractures of the second and third mediums.

16 . A method as claimed in claim 10 , wherein conduction transmissibility between the second and third mediums for the cell i is determined with a formula:

G

i

ff

=

Λ

i

ff

⁢

C

i

ff

where

Λ

i

ff

is an arithmetic mean of an effective thermal conductivity of the second and third mediums

C

i

ff

is a geometric coefficient depending on the dimensions of the cell i, the dimensions of one of the matrix blocks into which the first medium is broken, and the opening of the fractures of the second and third mediums.

17 . A method as claimed in claim 11 , wherein the conduction transmissibility between the second and third mediums for the cell i is determined with a formula:

G

i

ff

=

Λ

i

ff

⁢

C

i

ff

where

Λ

i

ff

is an arithmetic mean of an effective thermal conductivity of the second and third mediums,

C

i

ff

is a geometric coefficient depending on the dimensions of the cell i, the dimensions of one of the matrix blocks into which the first medium is broken, and the opening of the fractures of the second and third mediums.

18 . A method of exploiting a fluid of an underground formation comprising a fracture network, wherein the method as claimed in claim 1 is performed and, from at least simulation of the flows in the underground formation, an exploitation scheme comprising at least one site for at least one of an injection well and at least one production well is determined for the fluid, and the fluid of the underground formation is exploited at least by drilling the wells at the site and by providing wells with exploitation infrastructures.

19 . A non-transitory computer-readable medium storing a computer program product which, when executed on a computer, implements the method as claimed in claim 1 .

Assignments (1)

Priority Claims (1)

Continuity (1)

References Cited (13)