Abonnez-vous à DeepL Pro pour traduire des fichiers plus volumineux.

Visitez www.DeepL.com/pro pour en savoir plus.

CHAPITRE 1

QUELQUES CONCEPTS DE BASE DE L'INGENIERIE DES RESERVOIRS 1.1

INTRODUCTION

En illustrant les principales fonctions d'un ingénieur de réservoir, à savoir l'estimation des
hydrocarbures en place, le calcul d'un facteur de récupération et la fixation d'une échelle de
temps pour la récupération, ce chapitre introduit de nombreux concepts fondamentaux de
l'ingénierie des réservoirs.

La description du calcul du pétrole en place se concentre principalement sur la détermination

des régimes de pression des fluides et sur le problème de la localisation des contacts entre
les fluides dans le réservoir. La récupération primaire est décrite en termes généraux en
considérant l'importance des compressibilités isothermes des fluides du réservoir, tandis que
la détermination du facteur de récupération et la fixation d'une échelle de temps sont
illustrées par la description de l'ingénierie des réservoirs de gaz volumétriques. Le chapitre
se termine par une brève description quantitative du comportement des phases des
systèmes d'hydrocarbures à composants multiples.

1.2 CALCUL DES VOLUMES D'HYDROCARBURES

Considérons un réservoir initialement rempli d'huile liquide. Le volume d'huile dans le

réservoir (huile en place) est de

OIP = Vφ(1- Swc )(vol.rés.) (1.1)

où V = le volume apparent net de la roche

réservoir φ laporosité, ou fraction de volume de la roche qui est poreuse etSwc

= la saturation en eau connate ou irréductible, exprimée en

fraction du volume des pores.

Le produit Vφ est appelé volume de pore (PV) et représente le volume total du réservoir qui
peut être occupé par des fluides. De même, le produit Vφ (1-Swc) est appelé volume de pores
d'hydrocarbures (HCPV) et représente le volume total du réservoir qui peut être rempli
d'hydrocarbures, qu'il s'agisse de pétrole, de gaz ou des deux.

L'existence de la saturation de l'eau connate, qui est normalement de 10-25% (PV), est un
exemple de phénomène naturel fondamental pour l'écoulement des fluides dans les milieux
poreux. En effet, lorsqu'un fluide en déplace un autre dans un milieu poreux, la saturation du
fluide déplacé ne peut jamais être réduite à zéro. Cela s'applique à condition que les fluides
ne soient pas miscibles (ne se mélangent pas), ce qui implique qu'il existe une tension
superficielle finie à l'interface entre eux.

Ainsi, le pétrole, qui est généré dans une roche mère profonde, en migrant dans un piège
rempli d'eau, déplace une partie de l'eau, mais pas toute, ce qui entraîne la présence d'une
saturation en eau contiguë. Comme l'eau est immobile, sa seule influence dans le réservoir
est la saturation en eau.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 2
RÉSERVOIRS
Les calculs d'ingénierie visent à réduire le volume du réservoir qui peut être occupé par les
hydrocarbures.

Le volume de pétrole calculé à l'aide de l'équ. (1.1) est exprimé en volume de réservoir.
Étant donné que tous les pétroles, aux pressions et températures élevées qui prévalent dans
les réservoirs, contiennent différentes quantités de gaz dissous par unité de volume, il est
plus judicieux d'exprimer les volumes de pétrole aux conditions du réservoir de stockage
(surface), au moment où le pétrole et le gaz se seront séparés. Ainsi, le volume de pétrole
initialement en place dans le réservoir est de

STOIIP = n = vφ(1- Swc )/Boi (1.2)

où Boi est le facteur de volume de formation du pétrole, dans les conditions initiales, et a pour
unité le volume du réservoir/le volume du stock, généralement le baril du réservoir/le baril du
stock (rb/stb). Ainsi, un volume de Boi rb de pétrole produira un stb de pétrole à la surface,
ainsi que le volume de gaz qui était à l'origine dissous dans le pétrole dans le réservoir. La
détermination du facteur de volume de formation du pétrole et son application générale dans
l'ingénierie des réservoirs seront décrites en détail au chapitre 2.

Dans l'équ. (1.2), les paramètres φ et Swc sont normalement déterminés par l'analyse
pétrophysique et la manière dont ils sont évalués ne sera pas décrite dans ce texte1 . Le
volume net en vrac, V, est obtenu à partir de l'analyse géologique et de l'analyse de la
pression des fluides.

Le géologue fournit des cartes en courbes de niveau du sommet et de la base du réservoir,

comme le montre la figure 1.1. Ces cartes comportent des courbes de niveau tous les 15 mètres
d'altitude environ

Fig. 1.1 (a) Carte en courbes de niveau de la partie supérieure du réservoir, et (b) coupe
transversale du réservoir, le long de la ligne X-Y

et le problème consiste à déterminer le niveau auquel se situe le contact entre l'huile et l'eau
(OWC). La mesure du volume de la roche encaissée du réservoir au-dessus de ce niveau
donnera alors le volume net en vrac V. Dans la situation décrite à la figure 1.1 (b), il ne serait
pas possible de déterminer ce contact en inspectant les diagraphies effectuées dans le puits,
puisque seule la zone pétrolière a été pénétrée. Une telle technique pourrait toutefois être
appliquée si l'OWC se trouvait un peu plus haut dans le réservoir.

La manière dont le contact entre l'huile et l'eau, ou les contacts entre les fluides en général,
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 3
peut être localisé nécessite une connaissance des régimes de pression des fluides dans le
RÉSERVOIRS
réservoir, qui sera décrite dans la section suivante.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 4
RÉSERVOIRS
1.3 RÉGIMES DE PRESSION DES FLUIDES

La pression totale à n'importe quelle profondeur, résultant du poids combiné de la roche de

formation et des fluides, qu'il s'agisse d'eau, de pétrole ou de gaz, est connue sous le nom de
pression de recouvrement. Dans la majorité des bassins sédimentaires, la pression de
recouvrement augmente linéairement avec la profondeur et présente généralement un
gradient de pression de 1 psi/pied (fig. 1.2).

14, 7Pression (psia)

Profond
eur
PC GP
(pieds)

pression des morts-

surpression terrains
(OP )

sous-pression
hydrostatique
normal
pression

Fig. 1.2 Régimes de pression des morts-terrains et de pression hydrostatique (FP =

pression des fluides ; GP = pression des grains)

À une profondeur donnée, la pression des morts-terrains peut être assimilée à la somme de la
pression des fluides (FP) et de la pression des grains ou de la matrice (GP) agissant entre les
différentes particules de roche, à savoir

OP = FP + GP (1.3)

et, en particulier, puisque la pression des terrains de couverture reste constante à toute
profondeur donnée, alors

d FP() = -d GP( )
(1.4)

En d'autres termes, une réduction de la pression du fluide entraînera une augmentation

correspondante de la pression du grain, et vice versa.

Les régimes de pression des fluides dans les colonnes d'hydrocarbures sont dictés par la
pression de l'eau qui prévaut à proximité du réservoir. Dans un cas parfaitement normal, la
pression de l'eau à n'importe quelle profondeur peut être calculée comme suit
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 5
dDdp RÉSERVOIRS
pw eau × D + 14,7 (psia) (1.5)
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 6
RÉSERVOIRS
où dp/dD, le gradient de pression de l'eau, dépend de la composition chimique (salinité) et,
pour l'eau pure, a une valeur de 0,4335 psi/ft.

L'ajout de la pression de surface d'une atmosphère (14,7 psia) permet d'exprimer la pression
en unités absolues plutôt qu'en unités manométriques (psig), qui sont mesurées par rapport à
la pression atmosphérique. Dans de nombreux cas, l'ingénierie des réservoirs se préoccupe
principalement des différences de pression, qui sont les mêmes que l'on utilise des pressions
absolues ou des pressions manométriques, et sont désignées simplement par psi.

L'équation (1.5) suppose qu'il y a continuité de la pression de l'eau jusqu'à la surface et que la
salinité ne varie pas avec la profondeur. La première hypothèse est valable dans la majorité
des cas, même si les sables aquifères sont généralement entrecoupés de schistes
imperméables, puisque toute rupture de la continuité surfacique de ces joints apparents
entraînera l'établissement d'une continuité de la pression hydrostatique jusqu'à la surface.
Cette dernière hypothèse est toutefois assez naïve, car la salinité peut varier
considérablement avec la profondeur. Néanmoins, pour le moment, nous supposerons un
gradient de pression hydrostatique constant, à des fins d'illustration. Comme nous le verrons
plus loin, ce qui importe vraiment à l'ingénieur, c'est la définition du régime de pression
hydrostatique à proximité des sables contenant des hydrocarbures.

Contrairement à cette situation normale, on rencontre des pressions hydrostatiques

anormales qui peuvent être définies par l'équation suivante

pw = dDdp eau × D +14,7 + C(psia) (1.6)

où C est une constante qui est positive si l'eau est en surpression et négative si elle est en
dépression.

Pour que l'eau contenue dans un sable soit soumise à une pression anormale, le sable doit
être efficacement isolé des strates environnantes, de sorte que la continuité de la pression
hydrostatique jusqu'à la surface ne puisse pas être établie. Bradley2 a dressé une liste de
diverses conditions susceptibles de provoquer des pressions de fluide anormales dans les
sables aquifères fermés :

– changement de température ; une augmentation de température d'un degré

Fahrenheit peut provoquer une augmentation de pression de 125 psi dans un système d'eau
douce scellé.

– des changements géologiques tels que le soulèvement du réservoir ou son

équivalent, l'érosion de surface, qui entraînent tous deux une pression de l'eau dans le
sable du réservoir trop élevée pour sa profondeur d'enfouissement ; l'effet inverse se produit
dans un réservoir enfoncé, dans lequel une pression de fluide anormalement basse peut se
produire.

– l'osmose entre des eaux de salinité différente, le schiste de scellement jouant le rôle
de membrane semi-perméable dans cet échange ionique ; si l'eau à l'intérieur du scellement
est plus salée que l'eau environnante, l'osmose provoquera une pression anormalement
élevée et vice-versa.

Certaines de ces causes de pression anormale sont interactives. Par exemple, si un bloc de
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 7
réservoir est soulevé, la surpression
RÉSERVOIRS qui en résulte est partiellement atténuée par une baisse
de la température du réservoir.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 8
RÉSERVOIRS
Le manuel de géologie de Chapman3 fournit une description complète de la mécanique de la
surpression. Les ingénieurs de réservoir ont toutefois tendance à être plus pragmatiques que
les géologues en ce qui concerne les pressions anormales, les principales questions étant les
suivantes : les sables contenant de l'eau sont-ils soumis à une pression anormale et, dans
l'affirmative, quel est l'effet de cette pression sur l'étendue des accumulations d'hydrocarbures
?

Jusqu'à présent, seules les pressions hydrostatiques ont été prises en compte. Les régimes
de pression des hydrocarbures sont différents en ce sens que les densités du pétrole et du
gaz sont inférieures à celle de l'eau et que, par conséquent, les gradients de pression sont
plus faibles, les chiffres typiques étant les suivants

dp
= 0,45 psi/pied
dDwater

dp
= 0,35 psi/pied
dDoil

dp
= 0,08 psi/pied
dDgas

Ainsi, pour le réservoir contenant à la fois du pétrole et un gascon libre, illustré à la figure 1.3,
l'utilisation des gradients ci-dessus donnerait la distribution de pression illustrée sur le côté
gauche du diagramme.

Au contact entre le pétrole et l'eau, à 5500 ft, la pression dans le pétrole et l'eau doit être égale,
sinon il n'y aurait pas d'interface statique. La pression dans l'eau peut être déterminée à l'aide de
l'équation (1.5), arrondie au psi le plus proche, comme suit (1.5), arrondie au psi le plus proche,
comme suit

pw = 0,45 D + 15 (psia) (1.7)

Pression (psia)
2250 2375 2500 Exploration
Bien
5000 2265 2369

Profo GAZ
ndeur
(pieds GOC : p o = p w = 2385 GOC 5200'
) (EN
ANGLAIS)
RÉSULTATS DES
5250
PÉ TESTS
TR à 5250 ft
OL po = 2402 psia
E OWC 5500'
dp
5500 OWC : p dD = 0.35psi/ft
o= pw = 2490
EAU
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 9
RÉSERVOIRS

Fig. 1.3Régimes de pression dans le pétrole et le gaz pour une accumulation typique
d'hydrocarbures

qui suppose un régime de pression hydrostatique normal. Par conséquent, au contact entre
l'huile et l'eau
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 10
RÉSERVOIRS
po = pw = .45 × 5500 + 15 = 2490 (psia)

L'équation linéaire de la pression d'huile, au-dessus du contact entre l'huile et l'eau,

est alors po = 0,35D + constante

et puisque po = 2490 psia à D = 5500 ft, la constante peut être évaluée pour donner l'équation
suivante

po = 0,35D + 565 (psia) (1.8)

Au contact gaz-pétrole à 5200 ft, la pression dans les deux fluides doit être égale et peut être
calculée, en utilisant l'équation (1.8), comme étant de 2385 psia. (1.8), à 2385 psia.
L'équation de la ligne de pression du gaz peut alors être déterminée comme suit

po = 0,08D +1969 (psia) (1.9)

Enfin, en utilisant cette dernière équation, la pression du gaz au sommet de la structure, à

5000 pieds, peut être calculée comme étant de 2369 psia. Les lignes de pression dans la
colonne d'hydrocarbures sont tracées dans le diagramme de profondeur de pression, fig.
1.3, qui montre qu'au sommet de la structure, la pression du gaz dépasse la pression
hydrostatique normale de 104 psi. Ainsi, dans un puits foré à travers un schiste d'étanchéité
sur la crête de la structure, il y aura une forte augmentation de la pression de 2265 psi à
2369 psia lors de la première pénétration dans le réservoir à 5000 pieds. L'ampleur de la
discontinuité de pression lors du forage dans un réservoir d'hydrocarbures dépend de la
distance verticale entre le point de pénétration du puits et le contact entre l'eau et les
hydrocarbures et, pour une valeur donnée de cette distance, sera beaucoup plus grande si le
réservoir ne contient que du gaz.

Lors du forage d'un puits d'exploration et de la découverte d'un nouveau réservoir, l'un des
principaux objectifs est de déterminer la position des contacts fluides qui, comme décrit dans
la section précédente, facilitera le calcul du pétrole en place.

Considérons le puits d'exploration illustré à la figure 1.3, qui pénètre dans le réservoir près du
sommet de la colonne de pétrole. Le contact entre le gaz et le pétrole dans le réservoir sera
clairement "vu", à 5200 pieds, sur les diagraphies effectuées dans le puits. En revanche, le
contact entre le pétrole et l'eau ne sera pas visible puisqu'il se situe à environ 225 pieds en
dessous du point où le puits pénètre dans la base du réservoir. La position du contact ne peut
être déduite que du résultat d'un essai de puits, tel qu'une tige de forage4 ou un essai de
formation par câble5,6 , au cours duquel la pression et la température sont mesurées et un
échantillon de pétrole récupéré. L'analyse de l'échantillon permet de calculer la densité du
pétrole dans les conditions du réservoir et donc le gradient de pression du pétrole (voir
exercice 1, chapitre 2). Ensemble, la mesure de la pression et le gradient de pression
suffisent à définir la ligne droite qui est la relation pression-profondeur dans la colonne de
pétrole. Si un tel test était effectué à une profondeur de 5250 pieds, dans le puits de la fig.
1.3, l a pression mesurée serait de 2402 psia et le gradient de pression calculé de 0,35
psi/pi, ce qui est suffisant pour définir la ligne de pression du pétrole comme suit

po = 0,35D + 565 (psia) (1.8)

et l'extrapolation de cette ligne pour atteindre la ligne de pression hydrostatique normale

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 11
situera le contact entre l'huile et l'eau à 5500 ft.
RÉSERVOIRS

Ce type d'analyse repose essentiellement sur la connaissance du régime de pression

hydrostatique. Si, par exemple, l'eau est surpressée de seulement 20 psi, le contact entre l'huile
et l'eau sera à
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 12
RÉSERVOIRS
Ce fait peut être vérifié par une inspection visuelle de la figure 1.3 ou en exprimant l'équation de
la ligne d'eau en surpression, équ. (1.6) comme pw = 0.45D + 35 (psia)

et en résolvant simultanément l'équation (1.8) pour la condition que pw = po au contact

huile-eau. (1.8) pour la condition que pw = po au contact entre le pétrole et l'eau. La différence
de 200 pieds dans la position du contact peut faire une énorme différence dans le calcul du
pétrole en place, en particulier si l'étendue du réservoir est importante.

C'est pour cette raison que les ingénieurs de réservoir sont prêts à consacrer beaucoup de
temps (et donc d'argent) à la définition du régime de pression hydrostatique dans un nouveau
champ. Un moyen simple d'y parvenir consiste à effectuer une série de tests de formation par
câble5,6 dans le puits d'exploration, généralement après la diagraphie et avant la pose du
tubage, au cours desquels les pressions sont délibérément mesurées dans les sables
contenant de l'eau, à la fois au-dessus et en dessous du ou des réservoirs d'hydrocarbures.
La série de mesures de pression à différentes profondeurs permet de déterminer la ligne de
pression hydrostatique, équ. (1.6), d'être définie avec précision à proximité de l'accumulation
d'hydrocarbures, que le régime de pression soit normal ou anormal.

Ces tests sont répétés dans les premiers puits forés dans un nouveau champ ou une
nouvelle zone jusqu'à ce que les ingénieurs soient convaincus de l'uniformité de la pression
hydrostatique. Le non-respect de cette procédure peut entraîner une erreur importante dans
l'estimation des hydrocarbures en place, ce qui peut conduire à la formulation de plans de
développement de gisements terriblement imprécis.
Pression (psia)
2250 2375 2500 Exploration
bien
5000

GDT GAZ
Profo 5150 '
ndeur DPGWC
(pieds 5281 '
) COLONNE DE
PÉTROLE
Fig. 1.3
RÉSULTATS DES
5500 TESTS
à 5100 ft
pg = 2377 psia
DPOWC
5640 ' dpg
dD = .08 psi / ft

Fig. 1.4 Illustration de l'incertitude dans l'estimation de l'étendue possible d'une colonne de
pétrole, résultant de tests de puits dans la calotte gazeuse

La figure 1.4 illustre un autre type d'incertitude associé à la détermination des contacts fluides
à partir des mesures de pression. Le réservoir est le même que celui de la figure 1.3, mais
dans ce cas, le puits d'exploration n'a pénétré que dans le gascon. Un essai de puits est
effectué à une profondeur de 5100 pieds, ce qui permet de déterminer que la pression du gaz
est de 2377 psia et, d'après l'analyse d'un échantillon prélevé (voir exercice 1.1), que le
gradient de gaz dans le réservoir est de
0,08 psi/ft. À partir de ces données, l'équation de la conduite de pression de gaz peut être
définie comme suit
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 13
RÉSERVOIRS

po = 0,08D +1969 (psia) (1.9)

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 14
RÉSERVOIRS
N'ayant pas vu de pétrole dans le puits, l'ingénieur peut penser qu'il a pénétré dans un
réservoir de gaz uniquement, et extrapoler l'équ. (1.9) pour atteindre la ligne de pression
hydrostatique normale pw = 0,45D +15 (psia) (1.7)

à une profondeur de 5281 pieds, à laquelle pw = pg. Ce niveau est indiqué sur la figure 1.4
comme étant le contact gaz-eau le plus profond possible (DPGWC), en supposant qu'il n'y
ait pas de pétrole.

Par ailleurs, étant donné que le point le plus profond où du gaz a été observé dans le puits
est de 5150 pieds (GDT - gas down to), il n'y a aucune raison physique pour qu'une colonne
de pétrole ne s'étende pas immédiatement en dessous de ce point. La pression du pétrole au
sommet d'une telle colonne serait égale à la pression du gaz, qui peut être calculée à l'aide
de l'équ. (1.9) à 2381 psia. Par conséquent, l'équation de la conduite de pression d'huile, en
supposant le gradient d'huile utilisé précédemment de 0,35 psi/pied, serait po =
0,35D + 579 (psia)
et la résoudre simultanément avec l'équ. (1.7), pour la condition que po = pw, donne le contact
pétrole-eau à une profondeur de 5640 pieds. Cette profondeur est indiquée sur la figure 1.4
comme le contact pétrole-eau le plus profond possible (DPOWC) et correspond à la colonne
de pétrole maximale possible. Par conséquent, malgré le fait que le puits ait été
soigneusement testé, il subsiste un degré élevé d'incertitude quant à l'étendue d'une
éventuelle colonne de pétrole. Elle pourrait en effet être nulle (DPGWC - 5281 ft) ou, dans le
cas le plus défavorable, être de l'ordre d'un million d'euros.
Dans le cas le plus optimiste, elle pourrait s'étendre sur 490 pieds (DPOWC - 5640 pieds)
ou, au contraire, prendre n'importe quelle valeur entre ces deux limites. La figure 1.4 montre
également la colonne de pétrole actuelle de la figure 1.3.

C'est pourquoi la question se pose toujours, lorsqu'on pénètre dans un réservoir ne

contenant que du gaz, de savoir s'il existe une colonne de pétrole significative, ou un bord de
pétrole, en aval-pendage, qui pourrait être exploité. Le seul moyen sûr de le savoir est de
forer un autre puits plus en aval de la structure ou, si cela est mécaniquement possible, de
reboucher et de dévier le trou d'origine. Lors de la planification du forage d'un puits
d'exploration, il n'est donc pas toujours opportun d'orienter le puits vers le point le plus élevé
de la structure. Cela tend à maximiser les chances de trouver des hydrocarbures, mais va à
l'encontre de l'un des principaux objectifs du forage de puits d'exploration, qui est d'obtenir
autant d'informations que possible sur les réservoirs et leur contenu.

Après avoir déterminé les contacts fluides dans le réservoir, à l'aide des méthodes décrites
dans cette section, l'ingénieur est alors en mesure de calculer le volume net en vrac V
nécessaire pour calculer les hydrocarbures en place. Dans la figure 1.1 (a), par exemple,
cela peut être fait en planifiant les contours au-dessus de l'OWC7,8 .

Enfin, en ce qui concerne l'application de l'équ. (1.2), le chiffre correct du STOIIP ne sera
obtenu que si tous les paramètres de l'équation sont réellement représentatifs de leurs
valeurs moyennes dans l'ensemble du réservoir. Comme il est impossible d'obtenir de tels
chiffres, il est plus courant de représenter chaque paramètre de l'équation STOIIP par une
distribution de probabilité plutôt que par une valeur déterminée. Par exemple, il peut y avoir
plusieurs interprétations géologiques différentes de la structure donnant une dispersion des
valeurs du volume net en vrac V, qui pourrait être exprimée comme une distribution de
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 15
probabilité de la valeur deRÉSERVOIRS
ce paramètre.

L'équation STOIIP est ensuite évaluée à l'aide d'une procédure de calcul statistique,
proportionnelle à la qualité des données d'entrée, et les résultats sont exprimés en termes de
distribution de probabilité de la STOIIP. L'avantage de cette méthode est que, si une valeur
moyenne de la STOIIP peut être extraite de la distribution finale, les résultats peuvent
également être formulés en termes d'incertitude liée à ce chiffre, exprimée, par exemple, sous
la forme d'une norme.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 16
RÉSERVOIRS
écart par rapport à la moyenne9,10 . Si l'incertitude est très grande, il peut être nécessaire de
forer un ou plusieurs puits supplémentaires pour réduire la fourchette avant de procéder à
l'exploitation du gisement.

1.4 RÉCUPÉRATION DE PÉTROLE : FACTEUR DE RÉCUPÉRATION

L'équation (1.2), pour le STOIIP, peut être convertie en une équation pour le calcul de la
récupération ultime de pétrole simplement en la multipliant par le facteur de récupération
(FR), qui est un nombre entre zéro et l'unité représentant la fraction de pétrole récupérable,
comme suit

Récupération ultime (RU) = (V (1 Sφ -wc )/B ) RFoi × (1.10)

S'il est facile de dire "simplement en multipliant par le facteur de récupération", il est beaucoup
moins facile de déterminer quel devrait être le facteur de récupération pour un réservoir donné
et, en fait, c'est la détermination de ce chiffre qui est la tâche la plus importante de l'ingénieur
de réservoir.

Tout d'abord, on peut clairement distinguer deux types de facteurs de récupération. Le

premier est régi par les circonstances économiques actuelles et, de plus en plus, par des
considérations environnementales et écologiques, tandis que le second peut être classé
comme un facteur de récupération purement technique dépendant de la physique du système
réservoir-fluide. Malheureusement, le premier, bien qu'il soit peut-être le plus intéressant,
n'est pas le sujet de ce livre.

Les deux principales catégories de récupération d'hydrocarbures sont appelées récupération

primaire et récupération supplémentaire. La récupération primaire est le volume
d'hydrocarbures qui peut être produit grâce à l'utilisation de l'énergie naturelle disponible dans
le réservoir et l'aquifère adjacent. En revanche, la récupération supplémentaire est le pétrole
obtenu en ajoutant de l'énergie au système réservoir-fluide. Le type le plus courant de
récupération supplémentaire est l'injection d'eau, qui consiste à injecter de l'eau dans le
réservoir et à déplacer le pétrole vers les puits de production, augmentant ainsi l'énergie
naturelle du système. La mécanique de la récupération supplémentaire sera décrite plus loin,
au chapitre 4, sec. 9 et au chapitre 10 ; pour l'instant, seule la récupération primaire sera prise
en compte.

Toute la mécanique de la récupération primaire repose sur l'expansion des fluides dans le
réservoir et peut être mieux comprise en considérant la définition de la compressibilité
isotherme.

1 ∂V
c=- (1.11)
V ∂pT

La compressibilité isotherme est couramment appliquée dans la majorité des calculs

d'ingénierie des réservoirs, car elle est considérée comme une approximation raisonnable
selon laquelle, lorsque des fluides sont produits et qu'ils retirent de la chaleur du réservoir par
convection, la roche couverture et la roche de base, qui sont supposées agir comme des
sources de chaleur d'étendue infinie, remplacent immédiatement cette chaleur par
conduction, de sorte que la température du réservoir reste constante. Par conséquent, la
compressibilité, lorsqu'elle est mentionnée dans ce texte, doit toujours être interprétée comme
la compressibilité isotherme.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 17
RÉSERVOIRS
Le signe négatif est nécessaire dans l'équ. (1.11) parce que la compressibilité est définie
comme un nombre positif, alors que le différentiel, ∂V/∂p, est négatif, puisque les fluides se
dilatent lorsque leur pression de confinement diminue. Lorsque l'on utilise la définition de la
compressibilité de manière isolée, pour décrire l'épuisement du réservoir, il est plus illustratif
de l'exprimer sous la forme suivante

dV = cV p∆ (1.12)
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 18
RÉSERVOIRS
où dV est une expansion et ∆p une chute de pression, toutes deux positives. Il s'agit de
l'équation de base qui sous-tend toutes les formes de mécanismes de récupération primaire.
Dans le réservoir, si ∆p est considéré comme la chute de pression depuis la pression initiale
jusqu'à une pression inférieure,

pi - p, alors dV sera l'expansion correspondante du fluide, qui se manifeste par une production.

L'art de concevoir un facteur de récupération primaire élevé, en utilisant l'énergie naturelle du

réservoir, consiste à s'assurer que le dV, c'est-à-dire la production, est le fluide le plus
précieux du point de vue commercial dans le réservoir, à savoir le pétrole. La manière dont
cela peut être réalisé est illustrée schématiquement dans la fig.
1.5.
dV tot = production de
pétrole= dV o+ dV dV
+w g

aquifèr gascap
hu
e Vw dV W ile dVg Vg

Vo

Fig. 1.5 Récupération primaire de pétrole résultant de l'expansion du pétrole, de l'eau et du

gaz

Le diagramme illustre le fait assez évident que pour produire un réservoir de pétrole, il faut
forer des puits dans la zone pétrolière. Si le réservoir est en contact avec un gascon et un
aquifère, la production de pétrole due à une chute de pression uniforme, ∆p, dans
l'ensemble du système, aura des composantes dues à l'expansion séparée du gaz pétrolier
et de l'eau, d'où dVTOT = Production de pétrole = dVo + dVw + dVg

dans laquelle le bilan est exprimé en volumes de fluide aux conditions du réservoir. En

appliquant l'équ. (1.12), on peut l'exprimer comme suit : dVTOT = co Vo ∆p + cw Vw ∆p + cg Vg

∆p

En considérant les figures suivantes comme typiques des compressibilités des trois
composants à une pression de 2000 psia :

co = 15 × 10 /psi−6

cw = 3 × 10 /psi−6

1
≈ ; voir sec.1.5
cg = 500 × 10 /psi−6 p

il est évident que la contribution au dVTOT fournie par l'expansion du pétrole et de l'eau ne sera
significative que si Vo et Vw, les volumes initiaux de pétrole et d'eau, sont tous deux
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 19
importants. En revanche, en
RÉSERVOIRS
raison de sa très grande compressibilité, même un volume
relativement faible de gaz gascon contribuera de manière significative à la production de
pétrole.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 20
RÉSERVOIRS
Par conséquent, s'il est évident qu'il ne faut pas produire un aquifère, mais plutôt laisser l'eau
se dilater et déplacer le pétrole, il en va de même pour le gaz contenu dans le gascon, qui,
bien qu'ayant une valeur commerciale, est souvent conservé dans le réservoir et autorisé à
jouer son rôle très important en contribuant à la récupération primaire grâce à son expansion.
Les mécanismes de la récupération primaire du pétrole seront examinés plus en détail au
chapitre 3.

1.5 INGÉNIERIE DES RÉSERVOIRS DE GAZ VOLUMÉTRIQUES

L'ingénierie des réservoirs de gaz volumétriques est introduite à ce stade précoce du livre en
raison de la relative simplicité du sujet. Elle sera donc utilisée pour illustrer la manière dont un
facteur de récupération peut être déterminé et une échelle de temps attachée à la
récupération.

Cette simplicité s'explique par le fait que le gaz est l'une des rares substances dont l'état,
défini par la pression, le volume et la température (PVT), peut être décrit par une relation
simple impliquant les trois paramètres. Une autre substance de ce type est la vapeur saturée,
mais pour l'huile contenant du gaz dissous, par exemple, une telle relation n'existe pas et,
comme le montre le chapitre 2, les paramètres PVT doivent être dérivés de manière
empirique afin de définir l'état du mélange.

L'équation d'état d'un gaz idéal, c'est-à-dire un gaz dans lequel les attractions
intermoléculaires et le volume occupé par les molécules sont tous deux négligeables, est la
suivante

pV = nRT (1.13)

dans laquelle, pour les unités conventionnelles utilisées dans l'industrie

p = pression (psia) ; V = volume (pi3)

T = température absolue - degrés Rankine (°R=460+°F)

n= lenombre de moles lb, où une mole lb est le poids moléculaire du gaz exprimé
en livres.
et= laconstante universelle des gaz qui, pour les unités ci-dessus, a la
valeur de 10,732 psia.cu.ft/lb. mole.°R.

Cette équation résulte des efforts combinés de Boyle, Charles, Avogadro et Gay Lussac, et
n'est applicable qu'à des pressions proches de la pression atmosphérique, pour laquelle elle a
été dérivée expérimentalement, et à laquelle les gaz se comportent comme des gaz idéaux.

De nombreuses tentatives ont été faites dans le passé pour expliquer les écarts d'un gaz réel
par rapport à l'équation d'état d'un gaz idéal, dans des conditions extrêmes. L'une des plus
célèbres est l'équation de van der Waals qui, pour une mole Ib.d'un gaz, peut être exprimée
comme suit

Va
(p + 2 ) (V -b) = RT (1.14)

En utilisant cette équation, on affirme que la pression p, mesurée à la paroi d'un récipient
contenant un gaz réel, est inférieure à ce qu'elle serait si le gaz était idéal. En effet, la quantité
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 21
de mouvement d'une molécule de gaz sur le point de heurter la paroi est réduite par les
RÉSERVOIRS
attractions intermoléculaires ; la pression, qui est proportionnelle au taux de variation de la
quantité de mouvement, est donc réduite. Pour corriger ce phénomène, le terme a/V2 doit être
ajouté à la pression observée, où a est a
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 22
RÉSERVOIRS
constante en fonction de la nature du gaz. De même, le volume V, mesuré en supposant que
les molécules occupent un espace négligeable, doit être réduit pour un gaz réel par le facteur
b qui, lui aussi, dépend de la nature du gaz.

Le principal inconvénient de l'utilisation de l'équ. (1.14) pour décrire le comportement

des gaz réels rencontrés dans les réservoirs est que la pression maximale pour laquelle
l'équation est applicable est encore bien inférieure à la gamme normale des pressions
des réservoirs.

Des équations d'état plus récentes et plus performantes ont été dérivées, par exemple les
équations de Beattie-Bridgeman et de Benedict-Webb-Rubin (qui ont été résumées de
manière pratique au chapitre 3 de la référence 18), mais l'équation la plus couramment
utilisée dans la pratique par l'industrie est la suivante

pV = ZnRT (1.15)

où les unités sont les mêmes que pour l'équation (1.13) et où Z, sans dimension, est
appelé facteur Z. (1.13) et Z, qui est sans dimension, est appelé le facteur Z. En
exprimant l'équation sous la forme

P V =nRT Z

le facteur Z peut être interprété comme un terme par lequel la pression doit être corrigée pour
tenir compte de l'écart par rapport à l'équation des gaz idéaux.

Le facteur Z est fonction à la fois de la pression et de la température absolue mais, pour les
besoins de l'ingénierie des réservoirs, l'intérêt principal réside dans la détermination de Z, en
tant que fonction de la pression, à température constante du réservoir. La relation Z(p)
obtenue est alors appropriée pour la description de l'épuisement isotherme des réservoirs.
Trois façons de déterminer cette relation sont décrites ci-dessous. a) Détermination
expérimentale

Une quantité de n moles de gaz est chargée dans un récipient cylindrique dont le volume peut
être modifié par le mouvement d'un piston. Le récipient est maintenu à la température du
réservoir, T, pendant toute la durée de l'expérience. Si Vo est le volume de gaz à la pression
atmosphérique, l'application de la loi des gaz réels, équ. (1.15),

14,7 Vo = nRT

puisque Z=1 à la pression atmosphérique. À toute pression supérieure p, pour

laquelle le volume correspondant du gaz est V, alors pV = ZnRT

et en divisant ces équations, on obtient

Z= pV
14.7Vo

En faisant varier p et en mesurant V, la fonction isotherme Z(p) peut être facilement obtenue.
Il s'agit de la méthode la plus satisfaisante pour déterminer la fonction, mais dans la majorité
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 23
des cas, le temps et les dépenses nécessaires ne sont pas justifiés, car il existe des
RÉSERVOIRS
méthodes fiables de calcul direct, décrites ci-dessous.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 24
RÉSERVOIRS
b) La corrélation du facteur Z de Standing et Katz

Cette corrélation nécessite une connaissance de la composition du gaz ou, au moins, de la

gravité du gaz. Les hydrocarbures naturels sont principalement composés de membres de la
série des paraffines (CnH2n+2) auxquels s'ajoutent des impuretés non hydrocarbures telles que
le dioxyde de carbone, l'azote et le sulfure d'hydrogène. Le gaz naturel diffère du pétrole en
ce sens qu'il se compose principalement des éléments les plus légers de la série des
paraffines, le méthane et l'éthane, qui représentent généralement plus de 90 % du volume. La
composition typique d'un gaz est indiquée dans le tableau 1.1.

Pour utiliser la corrélation de Standing-Katz11 , il est d'abord nécessaire, à partir de la

connaissance de la composition du gaz, de déterminer la pression et la température pseudo-
critiques du mélange comme suit

ppc = n pi ci (1.16)
i

et

Tpc = nTi ci (1.17)

i

où la somme porte sur tous les composants présents dans le gaz. Les paramètres pci et Tci
sont la pression et la température critiques du composant ith , énumérés dans le tableau 1.1,
tandis que les ni sont les fractions de volume ou, pour un gaz, les fractions molaires de
chaque composant (loi d'Avogrado). L'étape suivante consiste à calculer ce que l'on appelle
la pression et la température pseudo-réduites

ppr = pppc (1.18)

et

Tpr = TTpc (1.19)

où p et T sont la pression et la température auxquelles il est nécessaire de déterminer Z.

Dans la majorité des problèmes d'ingénierie des réservoirs, qui sont isothermes, Tpr est
constante et ppr variable.

Avec ces deux paramètres, le tableau de corrélation de Standing-Katz, fig. 1.6, qui consiste
en un ensemble d'isothermes donnant Z en fonction de la pression pseudo-réduite, peut être
utilisé pour déterminer le facteur Z.

Par exemple, pour la composition du gaz indiquée dans le tableau 1.1, à une pression de 2000

psia et à une température de 180° F, le lecteur peut vérifier que ppc = 663,3 et Tpc = 374,1.

donner
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 25
RÉSERVOIRS
32000
ppr = 663, = 3,02 et Tpr = 374,1640= 1,71

à partir de laquelle, en utilisant la figure 1.6, le facteur Z peut être obtenu comme 0,865.

Composant Poids Pression Constante

moléculaire critique s Temp. Composition
(oR) typique (volume ou
668 mole)
fraction, ni)
CH4
Le méthane 16.04 343 .8470

C2 H6
Éthane 30.07 708 550 .0586

C3 H8
Propane 44.10 616 666 .0220

i-C4 H10 Isobutane 58.12 529 735 .0035

n-C4 H10 Butane normal 58.12 551 765 .0058

i-C5 H12 Isopentane 72.15 490 829 .0027

n-C5 H12 Pentane normal 72.15 489 845 .0025

n-C6 H14 Hexane normal 86.18 437 913 .0028

n-C7 H14 Heptane normal 100.20 397 972 .0028

n-C8 H18
Octane normal 114.23 361 1024 .0015

n-C9 H20 Nonane normal 128.26 332 1070 .0018

n-C10 H22 Décane normal 142.29 304 1112 .0015

CO2
Dioxyde de carbone 44.01 1071 548 .0130

H2S Sulfure d'hydrogène 34.08 1306 672 .0000

N2
Azote 28.01 493 227 .0345

TABLEAU 1.1 Constantes physiques des constituants communs des hydrocarbures

gazeux12 , et composition typique d'un gaz

Par convention, la composition des gaz naturels est indiquée en termes de composants
individuels jusqu'à l'hexane, l'heptane et les composants plus lourds étant regroupés sous le
nom de C7+ (heptanes-plus). Lors de l'analyse en laboratoire, le poids moléculaire et la
densité de ce groupe sont mesurés, ce qui permet de déterminer la pression et la
température pseudo-critiques du C7+ à partir de corrélations standard8,13 . Cela facilite à son
tour le calcul du facteur Z à l'aide de la méthode décrite ci-dessus.

Pour le calcul du facteur Z, on a supposé que les composants autres que les hydrocarbures,
à savoir le dioxyde de carbone, le sulfure d'hydrogène et l'azote, pouvaient être inclus dans
les sommations (équations (1.16) et (1.17)) afin d'obtenir la pression et la température
pseudo critiques. (1.16) et (1.17), pour obtenir la pression et la température pseudo-critiques.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 26
RÉSERVOIRS
Cette approche n'est valable que si les fractions volumiques des composants autres que les
hydrocarbures sont faibles, par exemple inférieures à 5 % en volume. Pour les quantités plus
importantes, des corrections aux procédures de calcul ci-dessus figurent dans le manuel
d'Amyx, Bass et Whiting8 . En revanche, si les fractions volumiques des composants autres
que les hydrocarbures sont très importantes (la teneur en dioxyde de carbone du champ de
Kapuni, en Nouvelle-Zélande, par exemple, est de 45 % en volume), il est préférable de
déterminer le facteur Z de manière expérimentale, comme décrit au point a) ci-dessus.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 27
RÉSERVOIRS

Fig. 1.6 Le tableau de corrélation du facteur Z de Standing et Katz11 (reproduit avec l'aimable
autorisation de la SPE de l'AIME)

Si la composition du gaz n'est pas disponible, la corrélation de Standing-Katz peut encore être
utilisée à condition que la gravité du gaz, basée sur l'échelle air = 1, à la pression
atmosphérique et à 60°F, soit connue (voir sec. 1.6). Dans ce cas, la fig. 1.7 est utilisée pour
obtenir la pression et la température pseudo-critiques ; ensuite, les équations (1.18) et (1.18)
sont utilisées. (1.18) et (1.19) peuvent être appliqués pour calculer les pseudo paramètres
réduits nécessaires pour obtenir le facteur Z de la fig. 1.6.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 28
RÉSERVOIRS

Fig. 1.7 Propriétés pseudo-critiques de divers gaz naturels et fluides de puits de condensats19

c) Calcul direct des facteurs Z

La corrélation Standing-Katz est très fiable et a été utilisée avec confiance par l'industrie au
cours des trente-cinq dernières années. Cependant, avec l'avènement des ordinateurs, il est
apparu nécessaire de trouver une technique pratique pour calculer les facteurs Z, afin de les
utiliser dans les réservoirs de gaz.
plutôt que d'introduire l'ensemble du tableau de corrélation à partir duquel les facteurs Z
pourraient être extraits par consultation d'une table. Takacs14 a comparé huit méthodes
différentes de calcul des facteurs Z qui ont été développées au fil des ans. Elles se répartissent
en deux catégories principales : celles qui tentent d'ajuster analytiquement les courbes des
isothermes de Standing-Katz et celles qui calculent les facteurs Z à l'aide d'une équation d'état.
Parmi ces dernières, la méthode de Hall-Yarborough15 mérite d'être mentionnée car elle est à la
fois extrêmement précise et très simple à programmer, même pour les petites calculatrices de
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 29
bureau, puisqu'elle ne nécessite que cinq registres de stockage.
RÉSERVOIRS
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 30
RÉSERVOIRS
Les équations de Hall-Yarborough, élaborées à partir de l'équation d'état de Starling-Carnahan,
sont les suivantes

Z = 0,06125ppr t e -1,2(1- t)2 (1.20)

y

où ppr = la pression pseudo-réduite t = la réciproque de la température

pseudo-réduite (Tpc/T)

et y = la densité "réduite" qui peut être obtenue comme solution de

l'équation.

-0,06125 ppr t e-1,2(1 t)- 2 + y + y(12 -+y)y33 - y4 -(14,76t - 9,76t2 + 4,58t )3 y2

(2,18 + 2,82t)
+ (90,7t - 242,2t2 + 42,4t )3 y =0 (1.21)

Cette équation non linéaire peut être facilement résolue pour y en utilisant la technique itérative
simple de Newton Raphson. Les étapes de l'application de cette technique sont les suivantes :

1) faire une première estimation de yk , où k est un compteur d'itérations (qui dans

ce cas est l'unité, par exemple y1 = 0,001)

2) substituer cette valeur dans l'équ. (1.21) ; à moins que la valeur correcte de y
n'ait été initialement sélectionnée, l'équ. (1.21) aura une petite valeur non nulle
Fk .

3) en utilisant l'expansion de la série de Taylor du premier ordre, une

meilleure estimation de y peut être déterminée comme suit

yk 1+ = yk -F kdFk (1.22)
dy

où l'expression générale de dF/dk peut être obtenue comme la dérivée de l'équ.

(1.21), c'est-à-dire

dFdy = 1+ 4y +(14y-2y)-44y3 + y4- (29.52t -19.52t2 + 9.16t )3 y

(1,18 + 2,82t)
+ (2,18 + 2,82t) (90,7t - 242,2t2 + 42,4t )3 y (1.23)

4) itérer, en utilisant les équations (1.21) et (1.22), jusqu'à obtenir une

convergence satisfaisante (F (1.21) et (1.22), jusqu'à l'obtention d'une
convergence satisfaisante (Fk ≈ 0).

5) la substitution de la valeur correcte de y dans l'équ. (1.20) donne le facteur Z.

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 31
(N.B. il semble y avoir uneRÉSERVOIRS
erreur typographique dans l'article original de Hall-Yarborough15 ,
dans la mesure où les équations présentées pour F (équ. 8) et dF/dy (équ. 11), contiennent
1-y3 et 1-y4 dans l'équation de F (équ. 12), et 1-y dans l'équation de F (équ. 13).
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 32
RÉSERVOIRS
dénominateurs du deuxième et du premier terme, respectivement, au lieu de (1-y)3 et (1-y)4
comme dans les équations (1.21) et (1.23) de ce texte. (1.21) et (1.23) de ce texte).

Takacs14 a déterminé que la différence moyenne entre le tableau de corrélation de Standing-

Katz et la méthode analytique de Hall-Yarborough est de - 0,158 % et la différence absolue
moyenne de 0,518 %. La figure 1.8 montre une relation isotherme entre le facteur Z et la
pression, obtenue à l'aide de la méthode Hall-Yarborough, pour un gaz d'une densité de 0,85
et à une température de réservoir de 200°F. Le tracé coïncide, à une épaisseur de crayon
près, avec la relation similaire suivante
obtenu par l'application de la méthode décrite au point b) ci-dessus.

Le graphique montre qu'il y a un écart significatif par rapport au comportement du gaz idéal qui
est particulièrement visible dans la plage de pression intermédiaire à environ 2500 psia. À
cette pression, l'utilisation de l'équation du gaz idéal (1.13) entraînerait une erreur de près de
25 % dans les volumes de gaz calculés.

1.6 APPLICATION DE L'ÉQUATION D'ÉTAT DES GAZ RÉELS

La détermination du facteur Z en fonction de la pression et de la température facilite

l'utilisation de l'équation simple suivante

pV = ZnRT (1.15)

pour définir complètement l'état d'un gaz réel. Cette équation est une relation PVT et dans
l'ingénierie des réservoirs, en général, la principale utilisation de ces fonctions est de relier
les volumes de surface aux volumes de réservoir d'hydrocarbures. Pour un gaz réel, en
particulier, cette relation est exprimée par le facteur d'expansion du gaz E, où

Z - Facto r

0.95

0.85

0.75
0 1000 2000 3000 4000 5000
Pression ( psia )
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 33
RÉSERVOIRS
Fig. 1.8 Facteur Z isotherme en fonction de la pression (gravité du gaz = 0,85 ; température = 200°
F)

E = Vsc = volume de nmoles de gaz dans des conditions normales

Vvolume de n moles de gaz aux conditions du réservoir

et en appliquant l'équ. (1.15) dans les conditions normales et dans les conditions du réservoir,
cela donne

E = Vsc = p Tsc Zsc (1.24)

Vpsc T Z

Pour les unités de champ définies en relation avec l'équ. (1.13), et pour les conditions standard
de psc =
14,7 psia, Tsc = (460+60) = 520°R et Zsc = 1, l'équ. (1.24) peut être réduite à

p
E = 35,37 (vol/vol) (1.25)
ZT
A une pression de 2000 psia et à une température de réservoir de 180°F, le gaz dont la
composition est détaillée dans le tableau 1.1 a un facteur Z de 0,865, comme déjà déterminé
à l'article 1.5(b). 1.5(b). Par conséquent, le facteur d'expansion du gaz correspondant est

E= = 127,8 (vol/vol)

En particulier, le gaz initialement en place (GIIP) dans un réservoir peut être calculé à l'aide
d'une équation similaire à l'équ. (1.2) pour le pétrole, à savoir

G = V (1 Sφ - wc )Ei ( 1.26) où Ei est évalué à la pression initiale.

Les autres paramètres importants qui peuvent être exprimés de manière pratique à l'aide de
l'équation d'état sont la densité réelle du gaz, la gravité et la compressibilité isotherme.

Puisque la masse de n moles de gaz est nM, où M est le poids moléculaire, la densité est la
suivante

ρ= nM = nM = Mp (1.27)
V ZnRT/p ZRT

La comparaison de la densité d'un gaz, quelles que soient la pression et la température, à la

densité de l'air dans les mêmes conditions donne

ρgas = (M/Z)gas
ρair (M/Z)air

et, en particulier, aux conditions normales

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 34
RÉSERVOIRS
M 97M
ρρgasair = γg = Mgasair = 28.
(1.28)

où γg est la gravité du gaz par rapport à l'air dans des conditions normales et est
conventionnellement exprimée, par exemple, comme γg = 0,8 (air = 1).

Par conséquent, si la gravité du gaz est connue, M peut être calculé à l'aide de l'équ. (1.28) et
substituée à l'équ. (1.27) pour obtenir la densité à n'importe quelle pression et température. Par
ailleurs, si la composition du gaz est connue, M peut être calculée comme suit

M = nMi i (1.29)
i

et de nouveau substituée dans l'équ. (1.27). Les masses moléculaires des différents
composants du gaz, Mi, sont indiquées dans le tableau 1.1. Il est également utile de se
rappeler la densité de l'air dans des conditions normales (quel que soit l'ensemble d'unités
utilisé par le lecteur). Pour les unités indiquées, ce chiffre est

(ρair) sc = 0,0763 ib/cu.ft, ce qui permet d'évaluer la densité du gaz dans des

conditions normales comme suit

ρsc = 0,0763 g(γ lbs/cu.ft) (1.30)

La dernière application de l'équation d'état consiste à dériver une expression pour la

compressibilité isotherme d'un gaz réel. En résolvant l'équ. (1.15) pour V donne

V = ZnRT
p

dont la dérivée par rapport à la pression est

∂V = - ZnRT 1 - 1 ∂Z
∂p p p Z ∂p

et en substituant ces deux expressions dans la définition de la compressibilité

isotherme, équ. (1.11), on obtient

cg = - V1 ∂∂Vp = p1 -Z1∂∂Zp (1.31)

Dans la figure 1.9, un graphique de la compressibilité du gaz définie par l'équation (1.31) est
comparé à l'expression approximative. (1.31) est comparé à l'expression approximative.

p1
cg = (1.32)
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 35
RÉSERVOIRS
pour le gaz de gravité 0,85 dont le facteur Z isotherme est représenté à la fig. 1.8 à 200°F.
Comme on peut le voir, l'approximation, équ. (1.32), est valable dans la plage de pression
intermédiaire entre 2000-2750 psia où ∂Z/∂p est petit, mais elle est moins acceptable à
des pressions très élevées ou très basses.

EXERCICE 1.1 GRADIENT DE PRESSION DE GAZ DANS LE RÉSERVOIR

1) Calculer la densité du gaz, dans des conditions normales, dont la

composition est indiquée dans le tableau 1.1.
2) Quel est le gradient de pression du gaz dans le réservoir à 2000 psia et 180°
F (Z = 0,865).

EXERCICE 1.1 SOLUTION

1) Le poids moléculaire du gaz peut être calculé comme suit

M = nMi i = 19,91
i

et donc, en utilisant l'équ. (1.28), la gravité est

COMPRESSIBILITÉ
DES GAZ

Fig. 1.9 Compressibilité isotherme d'un gaz en fonction de la pression (gravité du gaz = 0,85 ;
température

= 200° F) γg = MMair = 28.9719.91 = 0.687 (air = 1)

La densité dans les conditions normales peut être évaluée à l'aide de l'équation
(1.27) comme suit (1.27) comme suit

Mp
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 36
RÉSERVOIRS
sc = Z RTscscsc= 1 10.732×19.91 14.7× × 520 = 0.0524 Ib/cu.ft

ou, alternativement, en utilisant l'équ. (1.30) comme

ρ
sc= 0,0763γg = 0,0524 lb/cu.ft

2) La densité du gaz dans le réservoir peut être calculée directement à l'aide de

l'équ. (1.27), ou en considérant la conservation de la masse d'une quantité donnée de
gaz comme suit

(ρV)sc = (ρV)res ou ρres = ρscE

qui, en utilisant l'équ. (1.25), peut être évaluée à 2000 psia et 180o F comme ρres =

35.37ZTρscp = 35.370.865××0.0524(180 +×460)2000 = 6.696 Ib/cu.ft

Pour convertir ce nombre en un gradient de pression en psi/pied, il faut procéder à la

manipulation suivante.

dDdp = 6,696 ftIb2 × ft1 × inchft2 2 = 6,696144 = 0,0465 psi/ft.

gaz

1.7 BILAN DES MATIÈRES GAZEUSES : FACTEUR DE RÉCUPÉRATION

L'équation du bilan matière, pour tout système d'hydrocarbures, est simplement un bilan de
volume qui assimile la production totale à la différence entre le volume initial d'hydrocarbures
dans l e réservoir et le volume actuel. Dans l'ingénierie des réservoirs de gaz, l'équation est
très simple et nous l'examinerons maintenant pour les cas distincts où il n'y a pas d'afflux
d'eau dans le réservoir et pour ceux où l'afflux est important.

a) Réservoirs à épuisement volumétrique

Le terme d'épuisement volumétrique, ou simplement d'épuisement, appliqué à la

performance d'un réservoir signifie qu'à mesure que la pression diminue, en raison de la
production, il y a une quantité insignifiante d'afflux d'eau dans le réservoir à partir de
l'aquifère adjacent. Cela implique que l'aquifère doit être de petite taille (voir section 1.4).
Ainsi, le volume du réservoir occupé par les hydrocarbures (HCPV) ne diminuera pas au
cours de l'épuisement. L'équation (1.26) permet d'obtenir une expression du volume des
pores d'hydrocarbures sous la forme suivante (1.26) comme suit

HCPV = Vφ (1-Swc) = G/Ei

où G est le gaz initial en place exprimé aux conditions normales. Le bilan matière, également
exprimé aux conditions normales, pour un volume de production Gp donné, et une baisse
conséquente de la pression moyenne du réservoir ∆p = pi-p est alors,

Production =GIIP - Gaz non produit

 (sc) (sc) - (sc)
(1.33)
-
Gp =G (HCPV)E
QUELQUES BASESC CON CEPTS DANS L'INGENIERIE DES RESERVOIRS 24
GE
Gp =G
Ei

qui peut être exprimée

comme suit
(1.34)
Gp =1- E
G Ei

ou, en utilisant l'équ. (1.25),

comme
p = pi 1 - Gp
(1.35)
Z Zi G

Le rapport Gp/G est la récupération fractionnelle du gaz à n'importe quel stade de l'épuisement
et, si le facteur d'expansion du gaz E, dans l'équ. (1.34), est évalué à la pression d'abandon
proposée, la valeur correspondante de Gp/G est le facteur de récupération du gaz.

Avant de décrire comment l'équation du bilan matière est utilisée en pratique, il convient de
reconsidérer plus en détail le bilan exprimé par l'équation (1.33). (1.33). L'équation suppose
implicitement que l'afflux d'eau étant négligeable, le volume des pores d'hydrocarbures reste
constant pendant la déplétion. Cette hypothèse néglige toutefois deux phénomènes
physiques liés à la baisse de pression. Premièrement, l'eau connate dans le réservoir se
dilate et deuxièmement, à mesure que la pression du gaz (fluide) diminue, la pression du
grain augmente conformément à l'équation (1.4). (1.4).

En conséquence, les particules de roche se tasseront plus étroitement et le volume des

pores sera réduit. Ces deux effets peuvent être combinés pour donner la variation totale du
volume des pores d'hydrocarbures comme suit

d(HCPV) = - dVw +dVf (1.36)

où Vw et Vf représentent respectivement le volume initial de l'eau connate et le volume des

pores (PV). Le signe négatif est nécessaire puisqu'une expansion de l'eau connate entraîne
une réduction du HCPV. Ces changements de volume peuvent être exprimés, en utilisant
l'équ. (1.11), en termes de compressibilités de l'eau et des pores, cette dernière étant définie
comme suit

cf =- V1f ∂(∂GPVf ) (1.4)

où GP est la pression du grain qui est liée à la pression du fluide par d(FP) = -

d(GP)

par conséquent

cf =- V1f∂(∂FPVf ) = V1 f∂∂Vpf (1.37)

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 25
RÉSERVOIRS
où p est la pression du fluide. L'équation (1.36) peut maintenant être exprimée comme

d(HCPV) = cw Vw dp + cf Vf dp

ou comme une réduction du volume des pores d'hydrocarbures en tant que

d(HCPV) = - (cw Vw + cf Vf ) ∆p (1.38)

où ∆p = pi - p, la chute de pression du fluide (gaz). Enfin, en exprimant les volumes d'eau

interstitielle et d'eau connexe comme suit

Vf = PV = (1HCPV- Swc ) = E(1 Si -G wc ) et

Vw = PV × Swc = E(1 Si GS- wcwc )

la réduction du volume des pores d'hydrocarbures, équ. (1.38), peut être incluse dans l'équ.

(1.33), ce qui donne Gp = 1- 1 - (c Sw wc + c ) pf ∆ E (1.39) G 1- Swc Ei

comme bilan matière modifié. En insérant les valeurs typiques de cw = 3 × 10−6 /psi, cf = 10
× 10−6 /psi et Swc = 0,2 dans cette équation, et en considérant une chute de pression
importante de ∆p = 1000 psi, le terme entre parenthèses devient

1- ×10-6 ×103 = 1 0,013-

En d'autres termes, l'inclusion du terme représentant la réduction du volume des pores

d'hydrocarbures, due à l'expansion de l'eau et à la réduction du volume des pores, ne modifie
le bilan matière que de 1,3 % et est donc souvent négligée. La raison de cette omission est
que les compressibilités de l'eau et des pores sont généralement, mais pas toujours,
insignifiantes par rapport à la compressibilité du gaz, cette dernière étant définie dans la
section 1.6 comme étant approximativement la réciproque de la compressibilité de l'eau et de
celle des pores. 1.6 comme étant approximativement l'inverse de la pression. Comme décrit
au chapitre 3, sec. 8, la compressibilité interstitielle peut parfois être très importante dans les
réservoirs non consolidés peu profonds et les valeurs
des valeurs supérieures à 100 × 10−6 / psi ont été mesurées, par exemple, dans les champs
de la côte Bolivar au Venezuela. Dans de tels réservoirs, il serait inadmissible d'omettre la
compressibilité des pores dans le bilan matière du gaz. Dans un réservoir qui ne contient que
du pétrole liquide, sans gaz libre, la prise en compte des effets de l'eau connate et de la
compressibilité des pores doit être incluse dans le bilan matière, car ces compressibilités ont
le même ordre de grandeur que le pétrole liquide lui-même (voir chapitre 3, section 5).

Dans la majorité des cas, le bilan matière d'un réservoir de gaz de type déplétion peut être
décrit de manière adéquate à l'aide de l'équation (1.35). (1.35). Cette équation indique qu'il
existe une relation linéaire entre p/Z et la récupération fractionnelle Gp/G, ou la production
cumulée Gp, comme le montrent les figures 1.10(a) et (b), respectivement. Ces diagrammes
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 26
illustrent l'une des techniques de base de l'ingénierie des réservoirs, qui consiste à essayer
RÉSERVOIRS
de réduire toute équation, aussi complexe soit-elle, à l'équation d'une ligne droite.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 27
RÉSERVOIRS
extrapolées, alors que les fonctions non linéaires, en général, ne peuvent pas l'être. Ainsi, un
tracé de p en fonction de Gp/G ou de Gp aurait moins d'utilité que les représentations illustrées
à la figure 1.10, puisque les deux seraient non linéaires.

pi (a) pi (b)
Zi Zi

p p
Z Z

p
Z ab

G p/ G RF (RF) comp Gp G

Fig. 1.10 Représentations graphiquesdu bilan matière pour un réservoir de gaz à

déplétion volumétrique ; équ. (1.35)

La figure 1.10(a) montre comment le facteur de récupération (FR) peut être déterminé en
entrant l'ordonnée à la valeur de (p/Z)ab correspondant à la pression d'abandon. Cette
pression est dictée en grande partie par la nature du contrat gazier, qui spécifie généralement
que le gaz doit être vendu à un certain taux constant et à une pression de surface constante,
cette dernière étant la pression au point de livraison, le gazoduc. Une fois que la pression
dans le réservoir est tombée à un niveau inférieur à la somme des pertes de charge
nécessaires pour transporter le gaz du réservoir au gazoduc, le taux de production du plateau
ne peut plus être maintenu. Ces pertes de charge comprennent le rabattement de pression
dans chaque puits, qui est la différence entre les pressions moyennes d'écoulement du
réservoir et du fond du trou, provoquant l'écoulement du gaz dans le puits de forage ; la perte
de charge nécessaire à l'écoulement vertical vers la surface, et la perte de charge dans le
traitement et le transport du gaz jusqu'au point de livraison. Par conséquent, les réservoirs de
gaz sont souvent abandonnés à des pressions assez élevées. Il est toutefois possible
d'augmenter la récupération en produisant le gaz à des pressions beaucoup plus basses et
en le comprimant à la surface pour obtenir la récupération (RF)comp, comme le montre la figure
1.10(a). Dans ce cas, le coût d'investissement des compresseurs et leurs coûts d'exploitation
doivent être compensés par l'augmentation de la récupération du gaz.

La figure 1.10(b) illustre également les techniques importantes de l'ingénierie des réservoirs,
à savoir la "correspondance de l'historique" et la "prédiction". Les points encerclés dans le
diagramme, reliés par la ligne continue, représentent l'historique observé du réservoir. En
d'autres termes, pour les valeurs enregistrées de la production cumulative de gaz, les
pressions ont été mesurées dans les puits de production et une pression moyenne du
réservoir a été déterminée, comme décrit en détail dans les chapitres 7 et 8.

Comme les valeurs tracées de p/Z en fonction de Gp forment une ligne droite, l'ingénieur peut
être enclin à penser que le réservoir est un type d'épuisement et procéder à l'extrapolation
de la tendance linéaire pour prédire les performances futures. Dans ce cas, la prédiction
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 28
RÉSERVOIRS
porterait sur la manière dont la pression diminue en fonction de la production et, si le taux de
marché est fixe, en fonction du temps. En particulier, l'extrapolation à l'abscisse donnerait la
valeur du GIIP qui peut être vérifiée par rapport à l'estimation volumétrique obtenue comme
décrit dans les sections 1.2 et 1.3. 1.2 et 1.3. Cette technique de
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 29
RÉSERVOIRS
Faire correspondre l'historique de la pression de production observée en construisant un
modèle mathématique approprié, bien qu'en l'occurrence très simple, équ. (1.35), et en
utilisant le modèle pour prédire les performances futures est un élément fondamental du sujet
de l'ingénierie des réservoirs. b)
Réser
voirs d'entraînement de l'eau

Si la réduction de la pression dans le réservoir entraîne une expansion de l'eau de l'aquifère

adjacent et, par conséquent, un afflux dans le réservoir, l'équation du bilan matière doit être la
suivante
est alors modifiée comme suit

Production =GIIP - Gaz non produit

(sc) (sc) (sc) (1.40)

Gp =G - EGi - We E

où, dans ce cas, le volume des pores d'hydrocarbures à la pression inférieure est
réduit de la quantité We, qui est la quantité cumulative d'afflux d'eau résultant de la
chute de pression. L'équation suppose qu'il n'y a pas de différence entre les volumes
d'eau en surface et dans le réservoir et néglige à nouveau les effets de l'expansion de
l'eau connexe et de la réduction du volume des pores.

Si une partie de l'afflux d'eau a été produite, on peut en tenir compte en soustrayant
ce volume, Wp, de l'afflux, We, dans la partie droite de l'équation. Avec une légère
manipulation algébrique, l'équ. (1.40) peut être exprimée comme suit

p = pi 1- G p 1- W Ee i (1.41) Z Zi G G

où We Ei /G représente la fraction du volume initial des pores d'hydrocarbures inondée

par l'eau et est, par conséquent, toujours inférieure à l'unité. Si l'on compare ce
résultat au bilan matière de l'épuisement, équ. (1.35), on constate que l'effet de l'afflux
d'eau est de maintenir la pression du réservoir à un niveau plus élevé pour une
production cumulée de gaz donnée. En outre, l'équ. (1.41) est non linéaire,
contrairement à l'équ. (1.35), ce qui complique à la fois la comparaison des données
historiques et les prévisions. Des tracés typiques de cette équation, pour différentes
forces aquifères, sont présentés à la fig. 1.11.

Au cours de la phase de comparaison des historiques, une partie distincte du modèle

mathématique doit être conçue pour calculer l'afflux d'eau cumulé correspondant à une
chute de pression totale donnée dans le réservoir ; cette partie de la comparaison des
historiques est décrite comme "l'ajustement de l'aquifère". Pour un aquifère dont les
dimensions sont du même ordre de grandeur que le réservoir lui-même, le modèle
simple suivant peut être utilisé
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 30
RÉSERVOIRS

ligne d'épuisement
p
Z
C

p B
Z ab

Gp
G

Fig. 1.11 Représentation graphiquede l'équation du bilan de matière pour un

réservoir de gaz d'entraînement de l'eau, pour différentes forces aquifères ; équ. (1.41)

Nous = cW p∆

où c = compressibilité totale de l'aquifère (cw + cf)

W = le volume total d'eau, qui dépend principalement de la géométrie de l'aquifère

et∆p = la chute de pression à la limite initiale entre le réservoir et l'aquifère.

Ce modèle suppose que, l'aquifère étant relativement petit, une chute de pression
dans le réservoir est instantanément transmise dans l'ensemble du système réservoir-
aquifère. Dans ce cas, le bilan matière serait tel que représenté par le graphe A de la
figure 1.11, qui n'est pas significativement différent de la ligne d'épuisement.

Pour obtenir la réponse à la pression indiquée par les lignes B et C, le volume de

l'aquifère doit être considérablement plus important que celui du réservoir et il est
alors inadmissible de supposer une transmission instantanée de la pression dans
l'ensemble du système. Il y aura alors un décalage entre une perturbation de la
pression dans le réservoir et la réponse complète de l'aquifère. La construction d'un
modèle d'aquifère, incluant cette dépendance temporelle, est assez complexe et le
sujet sera reporté au chapitre 9, où l'utilisation d'un tel modèle pour la correspondance
historique et la prédiction sera décrite en détail.

L'un des aspects malheureux du retard dans la réponse de l'aquifère est que,
initialement, tous les tracés de bilan matière de la fig. 1.11 semblent linéaires et, si
l'historique de la production et de la pression est insuffisant pour montrer l'écart par
rapport à la linéarité, on peut être tenté d'extrapoler les premières tendances, en
supposant un réservoir de type épuisement, ce qui aboutirait à la détermination d'une
valeur trop importante du GIIP. Dans ce cas, une grande différence entre cette valeur
et l'estimation volumétrique du GIIP peut être diagnostique pour décider de l'existence
ou non d'un aquifère. Il s'ensuit également que la tentative de construire un modèle
mathématique pour décrire la performance du réservoir sur la base de données
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 31
RÉSERVOIRS
historiques insuffisantes peut produire des résultats erronés lorsqu'elle est utilisée
pour prédire la performance future du réservoir.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 32
RÉSERVOIRS

Si l'on dispose d'un historique de la pression de production, il est possible d'estimer le

GIIP, dans un réservoir à eau, en utilisant la méthode suivante, décrite par Bruns et
al.16 .
Le bilan matière de l'épuisement, équ. (1.34), est d'abord résolu pour déterminer le
gaz apparent en place comme suit

Gp (1.42)
Ga = 1 - E/Ei

S'il existe une conduite d'eau active, la valeur de Ga calculée à l'aide de cette
équation, pour des valeurs connues de E et de Gp, ne sera pas unique. Les valeurs
calculées successives de Ga augmenteront à mesure que l'écart de p/Z au-dessus de
la ligne d'équilibre des matériaux d'épuisement s'accroît, en raison du maintien de la
pression assuré par l'aquifère. La valeur correcte du gaz en place peut toutefois être
obtenue à partir de l'équ. (1.40) comme suit

E
G = Gp - We (1.43)
1-E/Ei

où We est l'afflux d'eau cumulé calculé à l'aide d'un modèle mathématique de l'aquifère,
au moment où E et Gp ont été mesurés.

Nous - trop petits

Ga Nous - modèle correct

Nous - trop grands

Nous E/ (1 -E / E i)

Fig. 1.12 du GIIP dans un réservoir de gaz d'entraînement de l'eau Leslignes

courbes en pointillésrésultent du choix d'un modèle d'aquifère incorrect, dépendant du
temps (voir chapitre 9).

En soustrayant l'équ. (1.43) de l'équ. (1.42), on obtient

Ga = G + 1 -W EeE/Ei (1.44)

Si les valeurs calculées de Ga, équ. (1.42), sont tracées en fonction de WeE/(1-E/Ei), le
résultat devrait être une ligne droite, à condition que le bon modèle d'aquifère ait été
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 33
RÉSERVOIRS

sélectionné,
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 34
RÉSERVOIRS

comme le montre la figure 1.12, et la valeur correcte de G peut être déterminée par
extrapolation linéaire à l'ordonnée. La sélection du bon modèle d'aquifère (ajustement
de l'aquifère) est un processus d'essais et d'erreurs qui se poursuit jusqu'à l'obtention
d'une ligne droite.

Une autre caractéristique intéressante de la figure 1.11 est que la récupération

maximale possible de gaz, indiquée par les points encerclés, dépend du degré de
maintien de la pression offert par l'aquifère, étant plus faible pour les aquifères les plus
réactifs. La raison en a déjà été mentionnée au sec. 1.2 ; dans le déplacement non
miscible d'un fluide par un autre, tout le fluide déplacé ne peut pas être retiré de
chaque espace poreux.
Ainsi, lorsque l'eau progresse dans le réservoir, une saturation résiduelle en gaz est
piégée derrière le front. Cette saturation en gaz, Sgr, est assez élevée, de l'ordre de 30
à 50 % du volume des pores7,17 , et est largement indépendante de la pression à
laquelle le gaz est piégé. Dans ces conditions, l'application de l'équation d'état, équ.
(1.15), au gaz piégé par pied cube de volume poreux derrière le front d'inondation,
donne

pZ Sgr = nRT

et, puisque Sgr est indépendant de la pression, pour une déplétion

isotherme n
ce qui indique qu'une plus grande quantité de gaz est piégée à haute pression qu'à basse
pression.

La récupération ultime du gaz dépend à la fois de la nature de l'aquifère et de la

pression d'abandon. Pour la valeur de (p/Z)ab indiquée à la figure 1.11, l'aquifère
donnant la réponse de pression correspondant à la ligne B est le plus favorable. Si le
choix de la pression d'abandon est sous le contrôle de l'ingénieur, le choix de
l'aquifère ne l'est malheureusement pas. Il est donc extrêmement important de
mesurer avec précision les pressions et la production de gaz pour permettre la
construction d'un modèle aquifère fiable qui, à son tour, peut être utilisé pour les
prévisions de performance.

L'un des aspects les plus aventureux de l'ingénierie des réservoirs de gaz est que les
contrats de vente de gaz, spécifiant le taux du marché et la pression du gazoduc, sont
généralement conclus entre l'exploitant et l'acheteur très tôt dans la vie du gisement,
lorsque la quantité de données historiques est minimale. L'exploitant est alors
contraint de prendre des décisions importantes sur la durée pendant laquelle il sera
en mesure de répondre à la demande du marché, sur la base de données plutôt
maigres. Des études de sensibilité sont généralement menées à ce stade, en utilisant
les équations simples de bilan matière présentées dans ce chapitre et en faisant
varier les principaux paramètres, à savoir

– le GIIP

– le modèle d'aquifère, basé sur les configurations géométriques possibles de

l'aquifère
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 35
– RÉSERVOIRS
la pression d'abandon ; l'application ou non d'une compression

en surface - l e nombre de puits de production et leur conception

mécanique.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 36
RÉSERVOIRS

(ce dernier point n'a pas été abordé jusqu'à présent, car il nécessite le développement
d'équations d'entrée de puits, qui seront décrites au chapitre 8). Les résultats d'une
telle étude peuvent donner une première orientation sur la meilleure façon de
développer la réserve de gaz.

EXERCICE 1.2 BILAN MATIÈRE DU GAZ

Les données suivantes sont disponibles pour un réservoir de gaz récemment découvert :
GWC = 9700 ft
Profondeur du = 9537 ft
centroïde
Volume net en vrac = 1,776 × 1010 cu.ft
(V)
φ= 0.19
Swc = 0 20 γg = 0,85

Bien qu'un échantillon de gaz ait été prélevé au cours d'un bref essai de
production, la pression du réservoir n'a pas été enregistrée en raison d'une
défaillance de l'outil. On sait cependant que le régime de pression de l'eau dans
la localité est pw = 0,441D + 31 psia.

et que le gradient de température est de 1,258°F/100 ft, avec une température de

surface ambiante de 80° F.

1) Calculer le volume du GIIP.

2) Il est prévu de conclure un contrat de vente de gaz dans lequel les points
suivants ont été stipulés par l'acheteur.

a) au cours des deux premières années, une augmentation du taux de

production à partir de zéro
-100 MMscf/d (millions) doivent être atteints pendant le développement du champ
b) le taux de plateau doit être maintenu pendant 15 ans à une pression de
livraison au point de vente qui correspond à une pression de réservoir
minimale de 1200 psia. Cette dernière exigence peut-elle être satisfaite ?
(Supposons que l'aquifère soit petit, de sorte que l'équation du bilan
matière de l'épuisement puisse être utilisée).

3) Lorsque les besoins du marché ne pourront plus être satisfaits, le taux de

production diminuera de manière exponentielle de 20 % par an jusqu'à ce qu'il
soit réduit à 20 millions de pieds cubes par jour (ce gaz sera utilisé comme
combustible dans les opérations de la société ou comprimé pour répondre à une
partie des besoins actuels du marché). (Ce gaz sera soit utilisé comme
combustible dans les opérations de la société, soit comprimé pour répondre à
une partie des besoins actuels du marché).

Quel sera le facteur de récupération total pour le réservoir et quelle sera la durée de
vie du projet ?
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 37
RÉSERVOIRS
EXERCICE 1.2 SOLUTION

a) Pour déterminer le GIIP, il faut d'abord calculer la pression initiale du gaz à

la profondeur du centroïde du réservoir. C'est-à-dire la profondeur à laquelle il y a

Comme il y a autant de gaz en haut qu'en bas, les pressions à utiliser dans l'équation du
bilan de matière seront toujours évaluées à cette profondeur.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 38
RÉSERVOIRS

Pour ce faire, la pression de l'eau au niveau du contact gaz-eau doit d'abord être calculée
comme suit : pw =

0,441 × 9700 + 31 = 4309 psia = pgGWC

et la température comme

T = (1,258× ) + 80 + 460 = 662oR

pour ce gaz de gravité 0,85, le tracé du facteur Z isotherme à 200°F (660°R), fig. 1.8,
peut être utilisé pour déterminer le facteur Z au GWC, avec une erreur négligeable.
Ainsi

ZGWC = 0,888

et EGWC = 35,37 p = 35,37 × 4309 = 259,3 ZT 0,888 × 662

Le gradient de pression dans le gaz, au niveau du GWC, peut maintenant être calculé,
comme décrit dans l'exercice 1.1, comme suit

dp = ρscE = 0,0763× 0,85× 259,3 = 0,117 psi/pied dD 144

144

La pression du gaz au centroïde est donc

p = pgGWC - dDdp GWC × ∆D (1.45)

p = 4309 - 0,117 × (9700 - 9537) = 4290 psia

et la température absolue au centroïde est

T = (1,258× ) + 80 + 460 = 660oR

On pourrait améliorer cette estimation en réévaluant le gradient de gaz au centroïde,

pour p = 4290 psia et T = 660°R, et en faisant la moyenne de cette valeur avec la
valeur originale au GWC pour obtenir un gradient de gaz plus fiable à utiliser dans
l'équation (1.45). (1.45). Toutefois, les gradients de gaz sont généralement si faibles
que cette correction est rarement nécessaire. Le lecteur peut vérifier que, dans le cas
présent, la correction ne modifierait la pression centroïde que de moins de la moitié
d'un psi.

Pour une pression centroïde et une température de 4290 psia et 660°R, le GIIP peut être
estimé comme suit

GIIP = G = Vφ (1-Swc) Ei (1.26)

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 39
RÉSERVOIRS

= 1,776x10 0,19 0.835. 374290 =

699.70 10× 9 scf

2) Le calendrier de production global peut être divisé en trois parties, les

périodes de montée en puissance, de plateau de production et de déclin, comme le
montre la figure 1.13.
Gp3
Taux G G
Q p1 p2
(MMscf /d )
Qo = 100 (MMscf/d)

50 (MMscf /d )
20 (MMscf /d)

t1 t2 t3

Fig. 1.13Calendrier de développement d'un champ gazier (exercice 1.2)

Il faut d'abord déterminer GP2 , c'est-à-dire la production cumulée lorsque la pression du

réservoir est tombée à 1200 psia et que le taux de plateau ne peut plus être maintenu.
Lorsque p = 1 200 psia, Z = 0,832 (fig. 1.8) et en utilisant le bilan matière de
l'épuisement, équ. (1.35),

p 1200

Gp2 = G(1- ZpZii 1200) = 699,70×10 (19 =

0.8870.8324290 )

GP2 = 491,04 10×9 scf

Étant donné que la production cumulée pendant la période

d'accumulation de deux ans est GP1 = Qavg × 2 × 365 = 50 × 106 × 2 × 365

= 36,5 × 109 scf de gaz.

au taux plateau de 100 millions de pieds cubes par jour.

Gp2 - Gp1 = (491,04 - 36,50) × 109 = 454,54 × 109 scf

et le temps pendant lequel ce taux peut être maintenu est t2 =

p2Q-oGp1 = 100 10454.54× 6××103659 = 12.45 ans

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 40
RÉSERVOIRS

Par conséquent, la durée pendant laquelle le taux de plateau peut être maintenu sera
inférieure à l'exigence d'environ 2,5 ans.

3) Pendant la période de décroissance exponentielle, le taux à tout moment

après le début de la décroissance peut être calculé comme suit

Q = Qo e -bt

où Qo est le taux à t = 0, c'est-à-dire 100 MMscf/d, et b est le facteur de décroissance

exponentielle de
0,2 p.a. Par conséquent, le temps nécessaire pour que le taux tombe à 20 MMscf/d sera
de

t = 1ln Qo = 1 ln100 = 8,05 ans b Q 0.2

20

Si gp est la production cumulative de gaz à l'instant t, mesurée à partir du début de la

baisse, alors
t t

-bt
gp = Qdt = Q eo dt
o o

c'est-à-dire

gp = Qbo (1- e-bt)

et lorsque t = 8,05 ans.

gp(8.05) = (1- e-0.2 × 8.05) =146.02 10×9 scf

Par conséquent, la récupération totale cumulée au moment de l'abandon sera de

Gp3 = Gp2 + gp(8.05) = (491.04 +146.02) 10× 9 = 637.06 10× 9 scf

et le facteur de
récupération G
RF = Gp3= 637.06 10699.70 10××99 = 0.91 ou 91% GIIP

qui sera récupéré après une période totale de t1 + t2 +

 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 41
t3 = RÉSERVOIRS
2 + 12,45 + 8,05 = 22,5 ans.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 42
RÉSERVOIRS

Cet exercice simple couvre le spectre de l'activité d'ingénierie des réservoirs, à savoir
l'estimation des hydrocarbures en place, le calcul d'un facteur de récupération et
l'établissement d'une échelle de temps pour la récupération. Cette dernière est
imposée par le taux de marché global requis pour le champ, c'est-à-dire

temps = production cumulée

taux
d'occup
ation

Plus loin dans le livre, aux chapitres 4, 6 et 8, la méthode de calcul des taux de puits
individuels est décrite, ce qui signifie que l'échelle de temps peut être fixée par le type
d'expression le plus habituel.

∆m p( )nD = 3,161 × 10− 12

βγΤ g

Q2 = FQ2 µwh rp
w2
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 43
RÉSERVOIRS

1.8 COMPORTEMENT DE LA PHASE HYDROCARBURE

Ce sujet a été largement traité dans des ouvrages spécialisés8,13,18 et est décrit ici de
manière quelque peu superficielle, simplement pour fournir une compréhension qualitative
de la différence entre les divers systèmes d'hydrocarbures tels qu'ils existent dans le
réservoir.

Considérons tout d'abord l'expérience simple dans laquelle un cylindre contenant l'un des
membres les plus légers de la série des paraffines, le C2 H6-éthane, est soumis à une pression
continuellement croissante à température constante. À une certaine pression unique (la
pression de vapeur) au cours de cette expérience, le C2 H6, qui était totalement en phase
gazeuse à basse pression, se condense en liquide. Si cette expérience était répétée à une
série de températures différentes, le diagramme de phase résultant, qui est la relation
pression-température, pourrait être dessiné comme indiqué à la figure 1.14(a).
CP 50 % - C 2 H 6
100 % - C 2 H 6 100 % - C 7 H 16
LIQUID 50 % - C 7 H 16
E
P CP P P LIQUID
LIQUIDE E
+
LIQUID CP GAZ
GAZ
E GAZ GAZ

T T T
(a) (b) (c)

Fig. 1.14 Diagrammes de phasepour (a) l'éthane pur ; (b) l'heptane pur et (c) pour
un mélange 50-50 des deux

La ligne définissant les pressions auxquelles se produit le passage du gaz au liquide, à

différentes températures, est connue sous le nom de ligne de pression de vapeur. Elle se
termine au point critique (PC) où il n'est plus possible de distinguer si le fluide est liquide ou
gazeux, les propriétés intensives des deux phases étant identiques. Au-dessus de la ligne
de pression de vapeur, le fluide est entièrement liquide, tandis qu'en dessous, il est à l'état
gazeux.

Si l'expérience ci-dessus était répétée pour un membre plus lourd de la série des paraffines,
par exemple C7 H16 - heptane, les résultats seraient ceux de la figure 1.14(b). Une différence
évidente entre
(a) et (b) est qu'à des températures et des pressions plus basses, la tendance à
l'augmentation de la température et de la pression est plus grande.
hydrocarbure plus lourd, C7 H16, à l'état liquide.

Pour un système à deux composants, le diagramme de phase d'un mélange composé à 50

% de C2 H6 et à 50 % de C7 H16 est illustré à la figure 1.14 (c). Dans ce cas, s'il existe des
régions où le mélange de fluides est soit entièrement gazeux, soit entièrement liquide, il
existe également une région clairement définie dans laquelle les états gazeux et liquides
peuvent coexister ; c'est ce que l'on appelle la région biphasée. La forme de l'enveloppe
définissant la région à deux phases dépend de la composition du mélange, étant plus
inclinée verticalement si le C2H6 est le composant prédominant et plus horizontale s'il s'agit
du C7 H16.

Les hydrocarbures naturels sont plus complexes que le système illustré à la figure 1.14, car
ils contiennent un grand nombre de membres de la série des paraffines et généralement
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 44
RÉSERVOIRS

quelques impuretés non hydrocarbures. Néanmoins, un diagramme de phase peut être

défini de la même manière pour les mélanges complexes et un tel diagramme pour un gaz
naturel typique est illustré à la figure 1.15(a).
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 45
RÉSERVOIRS

Les lignes définissant la région à deux phases sont décrites comme la ligne du point de
bulle, qui sépare le liquide de la région à deux phases, et la ligne du point de rosée, qui
sépare le gaz de la région à deux phases. En d'autres termes, lorsque l'on franchit la ligne
du point de bulle de
CP C A
LIQUID LIQUID
E D E
A

CP
P P
CT B

90%
E 70%
90%
70 % x 50 %
50 %
30% 30%
10% GAZ 10 % GAZ
B
T T
(a) (b)

Fig. 1.15 Diagrammes de phase schématiqueshydrocarbures à plusieurs

composants ; (a) pour un gaz naturel ; (b) pour un pétrole

Si l'on traverse la ligne du point de rosée à partir du gaz, les premières bulles de gaz
apparaissent, tandis que si l'on traverse la ligne du point de rosée à partir du gaz, les
premières gouttes de liquide (rosée) apparaissent. Les lignes à l'intérieur de la région
biphasée représentent des saturations de liquide constantes.

Pour un champ de gaz, comme décrit dans les sec. 1.5 - 1.8, la température du réservoir
doit être telle qu'elle dépasse ce que l'on appelle le cricondentherme (CT), qui est la
température maximale à laquelle les deux phases peuvent coexister pour le mélange
d'hydrocarbures en question. Si la pression et la température initiales du réservoir sont telles
qu'elles coïncident avec le point A de la figure 1.15(a), alors pour une déplétion isotherme du
réservoir, ce qui est généralement supposé, la pression diminuera de A vers le point B et la
ligne du point de rosée ne sera jamais franchie. Cela signifie qu'il n'y a que du gaz sec dans
le réservoir, quelle que soit la pression. Cependant, lors de la production du gaz à la surface,
la pression et la température diminuent et l'état final se situe en un point X de l'enveloppe à
deux phases, la position du point dépendant des conditions de séparation de la surface.

Les équations de bilan matière présentées dans ce chapitre, équs. (1.35) et (1.41),
supposent qu'un volume de gaz dans le réservoir est produit sous forme de gaz à la
surface. Si, en raison de la séparation en surface, de petites quantités d'hydrocarbures
liquides sont produites, le volume cumulé de liquide doit être converti en un volume de gaz
équivalent et ajouté à la production cumulée de gaz pour obtenir la valeur correcte de Gp à
utiliser dans l'équation du bilan matière.

Ainsi, si n livres-moles de liquide de poids moléculaire M ont été produites, la masse totale
de liquide est nM = γoρw × (volume de liquide).

où γo est la gravité du pétrole (eau = 1) et ρw est la densité de l'eau (62,43 Ib/cu.ft). Comme
les volumes d'hydrocarbures liquides sont généralement mesurés en barils de stockage
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 46
RÉSERVOIRS

(1 bbl = 5,615 cu.ft), alors le nombre de livres-moles d'hydrocarbures liquides produites

dans Np stb est de
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 47
RÉSERVOIRS

n = 350,5 γoNp
M

En exprimant ce nombre de moles d'hydrocarbures en volume de gaz équivalent dans les

conditions normales, on obtient

Vsc = nRTpscsc= 350,5γoMNp × 10,73214,7× 520

ouVsc = 1,33×105 γoMNp

La correction apportée par l'ajout du volume de gaz équivalent à la production cumulative de

gaz est généralement assez faible, de l'ordre d'un pour cent ou moins, et est parfois
négligée.

Si la pression et la température initiales du réservoir sont telles que le gaz se trouve au

point C (fig. 1.15(a)), le liquide commencera à se condenser dans le réservoir au cours de
la déplétion isotherme lorsque la pression sera tombée en dessous du point de rosée au
point D.

La saturation maximale du liquide déposé dans le réservoir, lorsque la pression se situe

entre les points D et E dans la région à deux phases, est généralement assez faible et
souvent inférieure à la saturation critique qui doit être dépassée pour que le liquide devienne
mobile. Ce phénomène est analogue aux saturations résiduelles, évoquées précédemment,
pour lesquelles l'écoulement cesse. Par conséquent, les hydrocarbures liquides déposés
dans le réservoir, que l'on appelle condensats liquides rétrogrades, ne sont pas récupérés
et, étant donné que les composants les plus lourds ont tendance à se condenser en premier,
cela représente une perte de la partie la plus précieuse du mélange d'hydrocarbures. On
peut imaginer que la poursuite de l'abaissement de la pression en dessous du point de
rosée à E entraînerait une nouvelle vaporisation du condensat liquide. Cela ne se produit
pas, cependant, parce qu'une fois que la pression tombe en dessous du point D, le poids
moléculaire global des hydrocarbures restant dans le réservoir augmente, étant donné
qu'une partie des paraffines les plus lourdes est laissée dans le réservoir sous forme de
condensat rétrograde. Par conséquent, l'enveloppe de phase composite des fluides du
réservoir a tendance à se déplacer vers le bas et vers la droite, ce qui empêche la
revaporisation.

Il est parfois économiquement viable de produire un champ de condensats de gaz par le

processus de recyclage du gaz sec. C'est-à-dire, à partir du début de la production au point
C, fig. 1.15(a), séparer le condensat liquide du gaz sec à la surface et réinjecter ce dernier
dans le réservoir de manière à ce que le gaz sec déplace le gaz humide vers les puits de
production. Étant donné que seule une quantité relativement faible de fluide est extraite du
réservoir au cours de ce processus, la chute de pression est faible et, pour un projet réussi,
l'objectif devrait être de fonctionner à une pression supérieure au point de rosée jusqu'à ce
qu'une percée de gaz sec se produise dans les puits de production. Ensuite, l'injection est
interrompue et le gaz sec restant est produit.

Les équations du bilan matière du gaz sec peuvent également être appliquées aux
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 48
RÉSERVOIRS
réservoirs de condensats de gaz si le facteur Z monophasé est remplacé par le facteur Z
biphasé. Ce facteur doit être déterminé expérimentalement en laboratoire en effectuant une
expérience d'épuisement à volume constant.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 49
RÉSERVOIRS

Un volume de gaz, G scf, est chargé dans une cellule PVT à une pression initiale pi, qui est
supérieure au point de rosée, et à la température du réservoir. La pression est réduite par
étapes au fur et à mesure que le gaz est retiré de la cellule, et mesurée comme Gp′ scf, sans
modifier ni le volume ni la température de la cellule. Cela simule la production du réservoir
dans des conditions d'épuisement volumétrique et, par conséquent, l'application de
l'équation du bilan matière de type épuisement, (1.35), et la résolution explicite de Z
donnent

Z2-phase = p p G'p (1.46)

i 1-
Zi G

Tant que la pression n'a pas chuté au point de rosée, le facteur Z mesuré dans cette
expérience est identique au facteur Z obtenu en utilisant la technique décrite à la section
1.5 (a). 1.5(a). En revanche, en dessous du point de rosée, les deux techniques donnent
des résultats différents.

Cette dernière expérience, qui permet de déterminer le facteur Z monophasé, suppose

implicitement qu'un volume de fluides du réservoir, inférieur à la pression du point de rosée,
est produit dans sa totalité jusqu'à la surface. Dans l'expérience d'épuisement à volume
constant, cependant, il est tenu compte du fait qu'une partie du fluide reste dans le réservoir
sous forme de condensat liquide, ce volume étant également enregistré en fonction de la
pression au cours de l'expérience. Par conséquent, si un échantillon de condensat gazeux
est analysé à l'aide des deux techniques expérimentales, le facteur Z biphasé déterminé
lors de l'épuisement à volume constant sera inférieur au facteur Z monophasé. Cela
s'explique par le fait que le condensat liquide rétrograde n'est pas inclus dans le facteur Z
cumulatif.
la production de gaz Gp′ dans l'équ. (1.46), qui est donc inférieure à ce qu'elle serait en
supposant que tous les fluides sont produits jusqu'à la surface, comme dans l'expérience à
phase unique.

La figure 1.15(b) montre un diagramme de phase typique pour le pétrole. Comme nous
l'avons déjà noté, étant donné que le pétrole contient une plus grande proportion des
éléments les plus lourds de la série des paraffines, l'enveloppe à deux phases est plus
inclinée horizontalement que dans le cas du gaz.

Si la température et la pression initiales sont telles que l'huile du réservoir se trouve au point
A du diagramme, il n'y aura qu'une seule phase dans le réservoir, à savoir l'huile liquide
contenant du gaz dissous. Une réduction isotherme de la pression amènera finalement
l'huile au point de bulle B. Ensuite, une nouvelle réduction de la pression produira un
système à deux phases dans le réservoir : l'huile liquide, contenant une quantité de gaz
dissous proportionnelle à la pression, et un volume de gaz libéré. Malheureusement, lorsque
l'huile liquide et le gaz sont soumis à la même différence de pression dans le réservoir, le
gaz, plus mobile, se déplace à une vitesse beaucoup plus grande que l'huile. Cela entraîne
un certain chaos dans le réservoir et complique considérablement la description de
l'écoulement des fluides.
 QUELQUES CONCEPTS DE BASE DE L'INGÉNIERIE DES 50
RÉSERVOIRS

De ce point de vue, il est préférable de produire le réservoir à une pression proche (ou
supérieure) du point de bulle, ce qui simplifie grandement la description mathématique. De
plus, comme nous le verrons au chapitre 3, un tel mode de fonctionnement permet
d'améliorer la récupération du pétrole. La manière dont la pression du réservoir peut être
maintenue au niveau ou au-dessus du point de bulle est traditionnellement l'injection d'eau,
dont une description détaillée est présentée au chapitre 10.