SPE 79698 Coupled Analysis of Geomechanics and Fluid Flow in Reservoir Simulation
SPE 79698 Coupled Analysis of Geomechanics and Fluid Flow in Reservoir Simulation
SPE 79698 Coupled Analysis of Geomechanics and Fluid Flow in Reservoir Simulation
Abstract
Generally, in classical reservoir studies, the geomechanical
behavior of the porous media is taken into account by the rock
compressibility. Inside the reservoir simulator, the rock
compressibility is assumed to be constant or to vary with the
pressure of the oil phase. It induces some changes in the
porosity field.
During the depletion phase or the cold water injection of HPHT reservoirs, the stress state in and around a reservoir can
change dramatically. This process might result in rock
movements such as compaction, induced fracturing,
enhancement of natural fractures and/or fault activation, which
continuously modify the reservoir properties such as the
permeabilities and the fault transmissibilities.
Modifications of such parameters strongly affects the flow
pattern in the reservoir and ultimately the recovery factor.
To capture the link between flow and in situ stresses, it
becomes essential to conduct coupled reservoirgeomechanical simulations.
This paper compares the use of 5 types of approach for the
reservoir simulations:
A classical approach with rock compressibility using only
a reservoir simulator,
A loose coupled approach between a reservoir simulator
(finite volumes) and a geomechanical simulator (finite
elements). At given user-defined steps, the hydrocarbons
pressures calculated by the reservoir simulator are
transmitted to the geomechanical tool which computes
the actual stresses and feeds back the modifications of the
petrophysical properties (porosities and permeabilities) to
the reservoir simulator.
A one way coupling: this approach is a simplification of
the loose coupled approach, the modifications are not fed
back to the reservoir simulator.
SPE 79698
ij = Dijkl kl pl
kl bPpij 3K d Tij
m
1
pl
= Pp + b kk pl
kk + kk flT
0fl
(1)
(2)
(3)
(4)
+ P2 2PPc = 0
Pc = Pco * exp pl
kk
(5)
(6)
Compaction drive
In standard reservoir simulation, the equations of flow are
independent of stress and the contribution of rock
compressibility to production is expressed as a function of
pore pressure20,21. The assumptions underlying this separation
between stress and flow are twofold:
The incremental effective stress ratio, defined as the ratio
of horizontal to vertical effective stress increments, is
uniform in the reservoir.
One of the stresses (the overburden stress normally)
remains constant or changes uniformly everywhere in
the reservoir.
(7)
h = (K 1) bPp + 3K d T
(8)
The multiplying factor (K 1) b is generally referred to
as the reservoir stress path. Eqs. 7 and 8 are sufficient to
calculate the volumetric strain as a function of changes of
pressure and temperature. For an elastic rock for instance, the
volumetric strain is given by :
bPp
elkk = (1 + 2K )
+ T
3
K
d
(9)
1 Vp 1 b
b2
(1 + 2K )
=
+
Vp Pp K s
3K d
(10)
1 Vp b
= (1 + 2K )
Vp T
(11)
and
Cel
PT =
1
Vp
Vp
= Fc * C pp
lab
Pp field
{ }
(12)
SPE 79698
X 2
V P (P ) = V P Pref .1 + X +
where
( )
X = C. P Pref
( )
SPE 79698
Layer 1 & 2
Layer 3
Porosity
40%
20%
Layer 4 & 5
38%
Kh
100mD
100mD
100mD
Kv
25mD
0,25mD
25mD
SPE 79698
SPE 79698
CPU
1
C=25 C=350
E-7
E-7
1,13
1,16
1,39
5
No
WW
With
WW
235
200
SPE 79698
Young's Poisson's
Rock layers modulus
ratio
(GPa)
Overburden 1
Overburden 2
C
Rese
B
rvoir
A
Sides
Pentland
Trias
Basement 3
12.5
10
12.5
12.5
12.5
12.5
20
20
20
0.19
0.12
0.19
0.19
0.19
0.19
0.19
0.19
0.19
24
20
24
24
24
24
24
24
24
0.2
0.1
0.15
0.2
0.1
0.2
0.1
0.1
0.1
Faults F3 F2 F5
Fault F8
Closing/opening
displacement
Opening (in
reservoir) and
closing elsewhere
Closing
Opening
Shear
displacement
direction
and directions
direction
SPE 79698
SPE 79698
.....Cam-Clay
hardening parameter
....Tensor
of total strain
Plastic multiplier
..Biot Modulus
Poissons ratio
Stress tensor
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
10
SPE 79698
1,05
-1
C = 35,0E-06 PSI
1,00
porosity multiplier
-1
C = 2,5E-06 PSI
0,95
0,90
1000
2000
3000
4000
5000
6000
7000
pressure P
(PSI)
Porosity multiplier
1,00
0,98
0,96
0,94
20%
0,92
38%
0,90
40%
0,88
0
1000
2000
3000
4000
5000
6000
7000
porosity multiplier
0,95
0,90
0,85
20 %
38 %
40 %
0,80
0,75
0
0,1
0,2
0,3
0,4
0,5
0,6
water sauration
0,7
0,8
0,9
SPE 79698
11
Seafloor
Layer
Porosity
Youngs modulus of rock
Overburden
bar
bar/m
Reservoir
Side
burden
Underburden
dimensionless
1&2
20%
38%
40%
133842
42754
37662
0,148
0,115
0,111
0,33
0,24
0,23
dimensionless
dimensionless
0,75
dimensionless
0,75
degrees
bar
Cohesion
Friction angle
Radius of ellipse Pc
Saturated Radius of Ellipse
factor
Critical saturation
Hardening parameter H
H =
90
357
96
83
degrees
69
38
35
bar
43
10
4&5
dimensionless
bar
1
1085
dimensionless
125
98
0,75
dimensionless
0,40
dimensionless
1 + e0
= Compressibility coefficient
= Swelling coefficient
e0 = void ratio
15,6
20,2
20,8
2845 m
0.204bar/m
580 bar
24131 bar
0.43
93333 ft
0.9 ft/psi
8400 psi
350000 psi
12
SPE 79698
Billions
Billions
10,1
Constant compresibility (C=25E-7)
Constant compresibility (C=350E-7)
Pseudo coupling (no water weakening)
Pseudo coupling (with water weakening)
Iterative coupling (no water weakening)
Iterative coupling (with water weakening)
FGPT (MSCF)
9,9
FRPV (RB)
9,7
2,5
1,5
9,5
0,5
9,3
0
0
9,1
10
15
20
25
30
years
35
8,9
0
10
15
20
25
30 years 35
1,01
7500
FPR (PSI)
6500
0,99
0,98
Pseudo coupling 1
Iterative coupling (1,1)
0,97
0,96
5500
0,95
4500
0,94
3500
4000
4500
5000
5500
6000
6500
7000
7500
3500
0
10
15
20
25
30
years
35
1,01
800
700
600
FOPT (STB)
Millions
0,99
0,98
Pseudo coupling 1
Iterative coupling (1,1)
0,97
0,96
500
400
0,95
300
0,94
3500
4000
4500
5000
5500
6000
6500
7000
7500
200
100
0
0
10
15
20
25
30
years
35
13
1,01
700
Pseudo coupling (with water weakening)
600
500
0,99
FOPT (STB)
1,00
Millions
SPE 79698
0,98
400
300
0,97
200
0,96
100
0,95
0
0
10
15
20
25
30
0,94
0,1
0,2
0,3
0,4
0,5
(Water Saturation)current - (Water Saturation)initial
0,6
0,7
FGPT (MSCF)
7500
Pseudo coupling (with water weakening)
7000
FPR (PSI)
6500
Billions
2,5
1,5
1
Pseudo coupling (with water weakening)
0,5
6000
5500
0
5000
4500
10
15
20
25
4000
3500
Billions
10
15
20
25
30
years
35
9,9
Pseudo coupling (with water weakening)
Iterative coupling (with water weakening)
9,8
9,7
FRPV (RB)
years
9,6
9,5
9,4
9,3
0
10
15
20
25
30 years 35
30
years
35
35
14
SPE 79698
Fig. 19 -- Interface normal (a) and tangential (b) stresses along the
faults.
Fig. 16 -- Principal stresses directions in reservoir zone.
(Pa)
8.E+07
6.E+07
4.E+07
2.E+07
p'
(Pa)
0.E+00
-2.5E+08 -2.0E+08 -1.5E+08 -1.0E+08 -5.0E+07 0.0E+00
SPE 79698
15
1,1
1,1
constant compressibility
pseudo coupling 2 with condensate banking
pseudo coupling 1
one way coupling
Multiplier
0,9
0,8
Porosity
0,7
K/Ki
0,9
0,6
0
200
400
600
800
Pressure
1000
1200
1400
0,8
0
bars
10
15
20
25
30
35
years
1,1
1
constant compressibility
0,95
0,9
Multiplier
0,8
Pseudo coupling 1
0,7
0,9
Porosity
0,6
K/Ki
0,5
0,85
0,4
10
15
20
25
30
35
years
0,3
0
500
1000
1500
Pressure - bars
500
250000
constant compressibility
Pseudo coupling 2 with condensate banking
Pseudo coupling 1
450
200000
FPR (bars)
400
350
150000
constant compressibility
Pseudo coupling 2 with condensate banking
100000
Pseudo coupling 1
One way coupling
50000
300
0
10
15
20
25
30
35
years
0
0
10
15
20
25
30 years 35
16
SPE 79698
1200
1000
800
600
400
200
0
0
0,5
1
1,5
2
2,5
Max. vertical displacement (m)