02-Performance de Réservoir
02-Performance de Réservoir
02-Performance de Réservoir
La relation entre le débit et la chute de pression dans le milieu poreux peut être très complexe et
dépend de plusieurs paramètres, tels que les propriétés pétrophysiques de la roche et les propriétés
des fluides, régime d’écoulement, la saturation des roches en fluide, la compressibilité des fluides, la
formation endommagée ou stimulée….
Plusieurs méthodes ont été élaborées pour construire la courbe IPR des puits à huile et à gaz, parmi
ces méthodes on a :
Qliq
IP=
p r −p wf ........................ 1-1
Exp :
Un puit produit 282 STB /day avec 1765 psi de pression de fond en écoulement. La pression du
réservoir est de2085 psi , quel est l’IP de ce puits.
En intégrant cette équation le long du ray de drainage du puits L’index de productivité peut être
exprimé par la formule suivante :
PR
0 , 00708⋅h k0
J= ⋅∫ dP
( )
e r P wf µ0⋅B 0
( P R−P wf )⋅ln 0 , 472⋅r
w ........................ 1-2
k0
f ( p )=
Si la fonction µ0⋅B 0 est indépendante de la pression où qu’elle peut être évaluée à la
0 , 00708⋅k 0⋅h
J=
( )
r
µ0⋅B 0⋅ln 0 , 472⋅ e
rW
........................ 1-3
Si toutes les variables dans les équations d’inflow pourraient être calculées, les équations résultantes
de l’intégration de la loi de Darcy, pourraient être employées pour construire les IPR.
Malheureusement, les informations sont rarement existantes d’une manière suffisante pour
appliquer cette équation. Donc des méthodes empiriques doivent être employées pour prévoir le
débit d’inflow d’un puits. La plupart de ces méthodes exigent au moins un test de stabilisation sur un
puits, et certains exigent plusieurs essais dans lesquels Pwf et q0 étaient mesurés.
Une méthode généralisée pour estimer la productivité d’un réservoir a été proposée par J.V. Vogel. Il
a écrit l’équation suivante :
( ) ( )
2
q0 PWf P
=1− 0,2⋅ − 0,8⋅ Wf
q0 ( max ) PR PR
q
Où max est le débit maximum avec une pression nulle au fond en écoulement.
En utilisant la courbe d’IPR de Vogel, Il est possible de calculer le débit maximum que le réservoir
peut produire. Il correspond à une pression de fond d’écoulement nulle (Pwf = 0). A ce point, le
‘’draw-down’’ sur la formation est maximum et est égal à la pression statique du réservoir Pr. Sur la
courbe, c’est le point (x=1, y=0). Cette valeur est théorique et ne peut jamais être atteinte.
Exp
1. Réservoir saturé :(
P R ≤Pb )
:
( ) ( )
2
q0 Pwf P wf
= 1−0 , 2⋅ −0 , 8⋅
q0 ( max) PR PR
[ ( ) ( )]
2
J⋅Pb P wf Pwf
q 0 =J⋅( P R − Pb ) + ⋅ 1−0 , 2⋅ −0 ,8⋅
1,8 PR PR P ≤P b
pour wf
q0
J=
P R−Pwf P ≥P b
pour wf
1.3 MODIFICATION DE STANDING
La méthode présentée par Vogel pour générer les IPR ne considère pas le changement de la
perméabilité absolue du réservoir. Standing propose une procédure pour modifier la méthode de
Vogel pour qu’elle tienne compte de l’endommagement ou la stimulation du réservoir.
q
drawdown ideal P R − P ' wf J' J
FE= = = =
drawdown actuel P R − Pwf q J'
J ........................ 1-4
( ) ( )
2
q0 P ' wf P ' wf
=1−0 , 2⋅ −0 , 8⋅
q FE=1 PR PR
0 ( max )
P , P' wf et FE peut être obtenue par la résolution de l’équation (1-4) pour P ' wf :
Une relation reliant wf
( )
P ' wf P
=1−FE+ FE⋅ wf
PR PR
( ) ( )
2
q0 P wf P
⇒ =1 ,8⋅( FE )⋅ 1− −0 , 8⋅( FE)2⋅ 1− wf
FE=1
q 0 (max ) PR PR
(
Pwf ≥P R⋅¿ 1−
1
FE)¿
FE≤1
Ou : cette contrainte est toujours satisfaite si
Exp :
Pour tenir compte des cas ou FE >1 Harrison a introduit la relation suivante :
( )
P wf
1, 792⋅
q0 PR
=1 , 2−0 , 2⋅e
q0FE=1
( max )
Dans le cas d’un réservoir sous saturé L'équation de Vogel peut être modifiée pour (FE#1) on obtient:
[ ( ) ( )]
2
J⋅Pb P wf Pwf
q 0 =J⋅( P R −Pb ) + ⋅ 1 , 8⋅ 1− −0 ,8⋅( FE)⋅ 1−
1,8 Pb Pb
( ) ( )
2
qo p p
=1−0 . 52 wf −0 . 48 wf
q o,max pr pr
( ) ( )
2
qw p p
=1−0 . 52 wf −0. 48 wf
q w,max pr pr
n
q 0 =C⋅( P2R−P2wf )
2
P2 R − P Wf
= A⋅q0 + B⋅q20 cas gaz
N horz =
√ k moy
k horz Coefficient de perméabilité horizontal
N vert=
√ k moy
k vert Coefficient de perméabilité vertical
2. Composant horizontal de drainage - plan horizontal
[ ] √ ( )
2
L eq
1+ 1−
2⋅r eq Leq
Dhorz =ln ×
Leq heq
2⋅r eq
D vert =ln
[ heq
2⋅π⋅r w ]
3. Terme de stockage de réservoir
W s= Dhorz + D vert + S
4. Générer les IPR en utilisant l’équation
( ( ) ( ))
2 n
qo p wf pwf
= 1−α −(1−α )
q o,max pr pr
( )( )
q0 PWf P Wf
= 1− ⋅ 1+ 0 , 8⋅
q0 ( max ) PR PR
( )
q o ,max p wf
J= 1+0 . 8
pr pr
¿ 1. 8 q o ,max
J = lim J =
p wf → pr pr
¿
J pr
q o ,max=
D’où 1.8
¿ 0 , 00708⋅k⋅h
J= ⋅f ( P R )
( )
re
ln 0 , 472⋅
rw
En utilisant la définition de on peut prédire le future index de
productivité du réservoir en utilisant l’index de productivité actuel
¿
JF f ( P RF )
¿ =
JP f ( PRP )
Avec
k ro
f ( p r )=
µ o Bo
D’où
q 0 ( max ) F =q 0 (max ) P⋅
[ P RF⋅f ( P RF )
PRP⋅f ( P RP ) ]
1.8.2 Méthode de Fetkovich
Cette méthode consiste à ajuster le coefficient d’écoulement C dans l’équation de Fetkovitch. Il a
f ( pr ) pr , et C peut être déterminé à partir de :
supposé que est une fonction linéaire de
( )
pr , f
C f =C p
pr , p
D’où
( )
Pr ,f n
q 0 , f =C p⋅ ⋅( p2r , f −P2wf )
Pr , p
2
p − p wf = Aqo + Bq o
Huile : r
2 2 2
p − pwf = Aq g + Bq g
Gaz : r
Où
1. calcul de
A R et BR
A R=
141 , 2⋅μo⋅Bo
K oR⋅h [(
⋅ ln
0 , 472⋅r e
rw )
+S d
]
2 , 3⋅10−14⋅β R⋅Bo⋅ρ o
BR =
h2⋅r w
2 , 33⋅1010
βR =
K 1R,2
( ) ( )
KR rd
Sd = −1 ⋅ln
Kd rw
Avec:
r w : Rayon du puits.
A B
2. Calcul des coefficients p et p
141 , 2⋅μo⋅Bo
A p= ⋅( S p + S dp )
K oR⋅h
2 ,3⋅10−14⋅β dp⋅B2o⋅ρo
B p=
r p⋅L2p⋅N 2
)[ ( ( )) ]
0,5
S p=
( h
hp
−1 ⋅ ln
h KR
⋅
rw K v
−2
( )( ) ( )
h KR KR r dp
Sdp = ⋅ − ⋅ln
L p⋅N K dp K d rp
2 ,33⋅1010
β dp=
K 1,
dp
2
Avec :
h : La hauteur totale de la formation
L P : La longueur de perforation.
r p : Rayon de perforations.
1 , 47⋅107
β G=
K 0G, 55
Avec :
K G : Perméabilité du gravier.