01-Principe Des Essais de Puits
01-Principe Des Essais de Puits
01-Principe Des Essais de Puits
Le réservoir est un objet dans lequel les variables de différents types (porosité et
perméabilité de la roche, viscosité de fluide…etc.) interagissent et produisent des signaux
observables (débit, pression). Comme, le réservoir est matériellement inaccessible, il est
nécessaire de procéder à son identification.
e
Dans ce chapitre, le problème que nous étudierons consiste à déterminer le modèle
représentatif du réservoir à partir de la mesure de la pression et du débit en fonction du
temps sur une courte durée de temps.
Durant des décennies, les modèles de réservoir considérés par l’interprétation des
essais de puits étaient limités à de géométries simples se comportant comme un réservoir
homogène. Dans les années 1970, les chercheurs ont tenté d’élargir les conditions
d’interprétation en prenant en compte par exemple la capacité d’expansion ou encore en
traitant les réservoirs de type particuliers (fissurés, fracturés…etc.).
Dans la littérature, nous pouvons classer ces méthodes en deux grandes familles :
1. Méthode conventionnelles
Elles ont été mises au point à partir des années 30. Elles étaient les seules
disponibles jusqu’aux années 70. Elles consistent à repérer sur l’évolution de pression les
différentes périodes d’écoulement caractéristiques qui se succèdent.
Ces méthodes sont apparues dans les années 70 mais n’ont pas été diffusées et ni
pris toutes leurs extension que dans les années 80.
Les méthodes des courbes types ont été fortement améliorées par l’utilisation
simultanée de la dérivée de la pression à partir de l’année 83. Elles ont en commun
d’interpréter d’un seul coup la globalité de l’évolution de pression enregistrée au
cours d’un essai de puits. Cette propriété permet à l’interprétateur de déterminer la
succession des écoulements visibles dans l’essai. Il peut ainsi porter un diagnostic sur
son puits et son réservoir.
À partir de 1990, une nouvelle méthode moderne a été mise au point par
Economides et al. qui consiste à unifier une méthodologie d’interprétation des essais de
puits. Elle repose principalement sur le calcul de la réponse impulsionnelle du réservoir
et de la régression non linéaire entre la pression enregistrée au fond du puits et la
pression calculée à partir d’un modèle du puits implanté sur un micro-ordinateur. Cette
méthode a permis de remédier à une bonne partie des limitations des courbes types.
Depuis que les méthodes modernes ont été introduites dans les essais de puits, le
déroulement d’une interprétation a été sensiblement modifié et amélioré.
Le problème d’identification du réservoir peut être subdivisé en quatre problèmes
élémentaires :
Régression
Identification du
Préparation
Déconvolution nonlinéaire
modèle initial
des données
1. seules les informations présentes sur les points de capture sont enregistrées ;
2. tout le reste est perdu.
q(t ) Réservoi p (t )
r
sortie p (t ) ( pi p(t ) ) : perte de charge dans le réservoir) sont reliés par un produit de
t
dpu ( )
p(t ) q(t ) d .................................... Error! No text of specified style in document.-1
0
d
p (t ) : La chute de pression.
ti
d
p(ti ) q(ti ) p( )d .................................... Error! No text of specified style in document.-2
0
d
i 1
d
p(ti ) q(ti )
t j 1
pu ( )d .................................... Error! No text of specified style in
j 0
tj d
document.-3
i 1
d
p(ti ) q(ti t j )
t j 1
pu ( )d .................................... Error! No text of specified style in
j 0
tj d
document.-4
j 1
d
Comme d p ( )d p (t
j
u u j 1 ) pu (t j ) alors
.................................... Error! No text of specified style in
i 1
p(ti ) q(ti t j ) pu (t j 1 ) pu (t j )
j 0
document.-5
p(t0 ) 0
p(t2 ) q(t2 t0 )pu (t1 ) pu (t0 ) q(t 2 t1 )pu (t 2 ) pu (t1 )
p(t3 ) q(t3 t0 )pu (t1 ) pu (t0 ) q(t3 t1 )pu (t2 ) pu (t1 ) q(t3 t2 )pu (t3 ) pu (t2 )
p(ti ) q(ti t0 )pu (t1 ) pu (t0 ) q(ti t1 )pu (t 2 ) pu (t1 ) ..........
q(ti t j ) pu (t j 1 ) pu (t j ) ........
q(ti ti 1 )pu (ti ) pu (ti 1 )
p(t n1 ) q(t n1 t0 )pu (t n1 ) pu (t0 ) q(t n1 t1 )pu (t 2 ) pu (t1 ) ..........
q(t n1 t j )pu (t n1 ) pu (t n ) ........
q(t n t n1 )pu (t n ) pu (t n1 )
Avec :
p(t1 )
p(t 2 )
.
.
p .................................... Error! No text of specified style in document.-7
p (t i )
.
.
p(t )
n 1
Et
pu (t1 )
pu (t 2 )
.
.
pu .................................... Error! No text of specified style in document.-8
pu (ti )
.
.
p (t )
u n1
Avec : q(ti t j ) 0 si ti t j 0
Figure Error! No text of specified style in document.-2 : calcul de la réponse impulsionnelle à partir d’un test
multi-rate
Alors le calcul de la réponse impulsionnelle revient à un calcul inverse.
pu Q 1p
I.1.3.1 Diagnostic
Le diagnostic sert à déterminer la succession des écoulements visibles au cours de
l’essai. Le repérage des ces écoulements détermine la configuration du réservoir puits qui
sera ensuite réutilisée dans l’interprétation. Le diagnostic est réalisé à l’aide de la dérivée
de la pression.
dp
p '
d ln t
La majorité des écoulements visibles au cours d’un essai de puits se traduit par une
évolution de la pression soit linéairement en fonction du logarithme du temps, soit
linéaire d’une puissance du temps, voir la figure suivante
I.1.3.2 Interprétation
L’interprétation vise à quantifier les paramètres de la configuration réservoir/puits.
Elle est réalisée avec les courbes types, la dérivée de la pression et les méthodes
conventionnelles.
Considérons les m couples d’observation (t1, p1), (t2, p2), (t3, p3),…(tm, pm), on
cherche les paramètres (k, w, λ, …etc.) de sorte que la somme des carrés des déviations
i m 2
S (k , w, , pi ...) pmésuré pcalculé (ti k , w, , pi ,...) ................... Error! No text of specified style in
i 1
document.-10
soit minimale.
Le minimum de la somme des carrés des résidus S est atteint lorsque le gradient
s’annule (condition nécessaire).
i m
pcalculé (ti k , w, , pi ,...)
p i 1
mésuré (ti ) pcalculé (ti k , w, , pi ,...)*
k
0
i m
pcalculé (ti k , w, , pi ,...)
p i 1
mésuré (ti ) pcalculé (ti k , w, , pi ,...)*
w
0
i m
pcalculé (ti k , w, , pi ,...)
p i 1
mésuré (ti ) pcalculé (ti k , w, , pi ,...)*
0
n
i m
pcalculé (ti k , w, , pi ,...)
pi 1
mésuré (ti ) pcalculé (ti k , w, , pi ,...) *
pi
0
équations
Écrivons ce système d’équation non linéaire sous forme abrégée. L’ensemble des
arguments k, w, λ,… etc. peut être considéré comme un vecteur de dimension n.
k
w
.................................... Error! No text of specified style in document.-11
.
.
.
k , w, ,... du modèle et que la pression pcalculé (ti ) soit continûment dérivable dans
un certain domaine convexe qui contient et k , le développement de Taylor autour de
k donne:
j 1 j
specified style in document.-15
pmésuré (t1 ) pcalculé (t1 ) j 1 j
pmésuré (t 2 ) pcalculé (t 2 ) j n p
calculé t 3
k
pmésuré (t3 ) pcalculé (t3 ) k
k
pmésuré (t3 ) pcalculé (t3 ) j 1 j
F ( )
. .
. .
pmésuré (t m ) pcalculé (t m ) j n p
calculé t m
k
pmésuré (t m ) pcalculé (t m )
k
k
j
j 1
.................................... Error! No text of specified style in document.-16
j n pcalculé t1 k
k
j 1 j
j n pcalculé t 2
k
pmésuré (t1 ) pcalculé (t1 )
k
k
j 1 j
pmésuré (t 2 ) pcalculé (t 2 ) j n
k
pcalculé t3
k
k
pmésuré (t3 ) pcalculé (t3 k ) j
j 1
. .
. .
pmésuré (t m ) pcalculé (t m k ) j n pcalculé t m k
k
j
j 1
………………………….....................................
Error! No text of specified style in document.-17
F ( ) F ( k ) J ( ) * k .................................... Error! No text of specified style in document.-18
Les formules (II-13) et (II-17) donnent :
J ( )T * F (m) J ( )T * F ( k ) J ( ) * k 0 .................................... Error! No text of
specified style in document.-19
J ( )T * F ( k ) J ( )T J ( ) * k 0 .................................... Error! No text of specified style in
document.-20
Par conséquent,
document.-21
(Méthode de Newton).
Algorithme de newton :
4. résoudre J ( k )T * J ( k ) * J ( k )T F ( k )
5. calculer k 1 , k 1 k
k 1 k
6. si :
k 1
convergence atteinte
écrire la solution k 1
arrêt
7. si le nombre maximal d’itération N est atteint :
convergence non atteinte en N itérations
arrêt
8. retour à l’étape 4.