Sujets

POROSITÉ

5.1 Introduction

Le pétrole brut et le gaz naturel s'accumulent dans les pores infinitésimaux des roches
et non dans des lagoons dans le sous-sol comme on a tendance à penser. La porosité avec
la perméabilité est la propriété physique la plus importante d'un réservoir
hydrocarbure, c'est pourquoi il est nécessaire de décrire correctement cette propriété. La
la porosité du réservoir est l'un des facteurs qui déterminent la stratégie à suivre
pour la gestion du réservoir. De plus, la porosité est essentielle pour créer des modèles
géologiques du réservoir qui caractérisent de manière adéquate le réservoir et de
de cette manière, on estime le pourcentage de pétrole brut ou de gaz sur place qui se
recuperera. Les estimations conventionnelles des réserves sont sensibles au type de
porosité que présente le réservoir.

5.2 Définition de la porosité

La porosité est une propriété de la roche qui décrit les espaces qui se
ils se trouvent dans la matrice de la roche donc la porosité est une indication de
la capacité de stocker des fluides d'une roche. Elle peut également être définie comme la
raisin du volume des pores et du volume total de la roche. Cela est représenté en
l'équation suivante :

Volume des pores VPoro VTotal− VMatrice

φ Total = = =
Volume total VTotal VTotal

Volume des pores VTotal− VMatrice

φ Total % = • 100 = • 100
Volume total VTotal

La valeur de la porosité est généralement multipliée par 100 pour la représenter

comme pourcentage et est représentée par le symbole grec Phi. La porosité varie de 0
un 1% ou si elle est exprimée en pourcentage, elle varie de 0 à 100. La porosité des
Les strates sédimentaires sont fonction de différents facteurs, certains de ces facteurs
fils

/ ÊÀi ÛiÊÌ ÃÊ Ì Vi]ÊÛ Ã Ì\Ê

ÜÜÜ° Vi °V ÉÕ V ° Ì
 • Raison d'enterrement

• Profondeur d'ensevelissement

Problème 5-1. Déterminer la porosité d'un échantillon de noyau qui a les

dimensions suivantes : hauteur 4 cm, diamètre 2,3 cm et poids 0,3024 N. L'échantillon
le noyau est une arénite qui a une densité de 2,65
g/cm3.

Solution

Calculer le volume total de l'échantillon de noyau.

π
VTotal=VCylindre = • D 2• h
4

π
VTotal =VCylindre = • ( 2,3 cm ) 2
• 4 cm = 16,62 cm3
4

Calculer le volume de la matrice ou de la partie solide du noyau.

Poids
masse g Poids
ρ Arenisca = = =
Volume VArenisca g•VArenisca

Déterminer le volume de la matrice (volume de grès)

m 100 cm 1000 g
0,3024 kg • •
Poids s 2
1m 1 kg
VArenisca = = = 11,63 cm3
g• ρ Arenisca
981 cm/s • 2, 65 g/cm
2 3

VTotal− VMatrice 16,62 cm3− 11,63 cm3

φ Total = = = 0,30 = 30,0 %
VTotal 16,62 cm3

5.3 Classification de la porosité

La porosité peut être classée selon différents paramètres, comme par exemple :
morphologie, connectivité et par leur origine. Ci-après, les
classifications mentionnées.
5.3.1 Classification par sa morphologie

La porosité par sa morphologie peut être divisée en trois types qui sont :

1Porosité caternaria– la porosité caternante représente les pores qui

sont connectés à d'autres pores par plus d'une cavité. Les hydrocarbures
que peuvent contenir ce type de pores peuvent être produits à l'étape
primaire par la pression naturelle du réservoir, à l'étape secondaire et
La période ternaire peut également produire des hydrocarbures.

[Link]é de cul-de-sac – cette porosité représente les pores qui sont

connectés à d'autres pores par une seule cavité. Les hydrocarbures
accumulés dans ce type de pores ne peuvent généralement pas être déplacés et
produits avec l'injection de gaz ou d'eau mais peut produire
hydrocarbures par l'expansion du gaz naturel qui se produit en raison de la
réduction de la pression dans le réservoir.

[Link]é fermée– la porosité fermée représente les pores qui ne

ils ont des cavités qui les relient à d'autres pores de la formation, c'est
comparable à la porosité ineffective que l'on a dans la classification par
connectivité.

La porosité caternarie et de Cul-de-sac additionnées forment la porosité effective que

se présente dans le classement par connectivité.

5.3.2 Classification par sa connectivité

La connectivité fait référence aux pores qui sont interconnectés entre eux
par conduits. On peut définir trois types de porosité selon
connectivité des pores dans les roches.

[Link]é effective– la porosité effective se réfère à tous les pores de

la formation qui se trouvent interconnectées par des cavités. Ce
le type de porosité est la principale source qui fournira la production
de pétrole brut ou de gaz naturel aux différentes étapes de production du
réservoir.

VPoro Efectivo
φEfficace =
VTotal

[Link]é inefficace– la porosité inefficace consiste en des pores de la

formation qui sont isolés. Ce type de porosité ne fournira pas
une source de production de pétrole ou de gaz naturel à différentes étapes de
production à moins qu'un travail de fracture ou un autre type ne soit effectué
 stimulation dans le ou les puits. Ces travaux de stimulation permettent
créer des canaux qui relient un pourcentage des pores isolés avec les
poros qui composent la porosité effective et de cette manière on peut
produire un certain pourcentage du pétrole ou du gaz naturel qui se trouve
cumulé dans la porosité ineffective.

VPoro Inefectivo
φInefficace =
VTotal

3. Porosité totale– la porosité totale est une indication de tous les pores
existants dans la formation, cela inclut les pores connectés et non
connectés. En d'autres termes, c'est la somme de la porosité effective et
ineffectif

φTotal= φEfficace+ φInefficace

Figure [Link]ésentation de la porosité effective et ineffective dans un milieu poreux.

Problème 5-2. Déterminer le pourcentage de la porosité inefficace d'un échantillon

de noyau qui est composée principalement de grès dont la densité est de 2,66
g/cm3Les dimensions de l'échantillon sont les suivantes :

4,2 cm
Radio : 1,25 cm.

L'échantillon est d'abord pesé et enregistre un poids de 0,5 kg m/s2. Puis, se

désagrège l'échantillon et le pèse à nouveau en enregistrant une réduction de son poids de 8
pour cent. Supposons que l'échantillon ait été nettoyé et séché avant de procéder aux
mesures.

Solution

π
VTotal =VCylindre = • (2• 1,25 cm)2• 4,2 cm = 20,62 cm3
4

Poids
pâte g Poids Poids
ρ Areniscas = = = ⇒ VArenisca =
Volume VArenisca g• VArenisca g• ρ Areniscas

m 100 cm 1000 g
0,5 kg • •
s 2
1m 1 kg
VArenisca, 1 = = 19,16 cm3
2
981 cm/s • 2,66 g/cm 3

m 100 cm 1000 g
[ 0,5− (0,5• 0,08) kg] • •
s 2
1m 1 kg
VArenisca, 2 = = 17,63 cm3
981 cm/s • 2, 66 g/cm
2 3

VTotal− VMatrice 20,62− 19,16

φ Efficace, 1 (%) = • 100 = • 100 = 7,1 %
VTotal 20,62

VTotal − VMatrice 20,62− 17,63

φ Total, 2 (%) = • 100 = • 100 = 14,5 %
VTotal 20,62

φ Total= φ Efficace+ φ Inefficace ⇒ φ Inefficace= φ Total− φ Efficace

φ Inefficace= 14,5 %− 7,1 % = 7,4 %

5.3.3 Classification par son origine

La porosité d'une roche peut également être classée par son origine ou par les
mécanismes qui donnent origine à la porosité de la roche. Selon cette classification la
La porosité se classe en :
 Primaire

2. Secondaire

Porosité primaire

La porosité primaire est la porosité qui se crée dans la roche lorsqu'elle est
formée osseuse lorsque les sédiments se déposent. Cette porosité peut
sous-classifier dans :

• Inter-particules – ce type de porosité se perd rapidement dans

Arenites carbonatées durant les processus de compactage et
Cimentation, cependant cela est retenu dans des grès silicoclastiques.

• Intra particules– ce type de porosité est formé par les pores qui se
ils se trouvent entre les grains des roches carbonatées.

Porosité Secondaire

La porosité secondaire est la porosité qui se forme par le biais de processus

postérieurs à la création de la roche.

5.4 Mesure de la porosité

L'équation de base pour calculer la porosité est assez simple, mais dans le
l'industrie pétrolière cette équation est principalement utilisée dans les calculs de
laboratoire. L'estimation de la porosité par d'autres méthodes est plus complexe que
l'équation de base de la porosité.

VTotal− VMatrice
φ=
VTotal

La porosité peut être mesurée par trois méthodes :

1. Registres de puits– l'un des principaux usages des profils est le

estimation de la porosité, en général des cartes du sous-sol sont réalisées
que montrent les changements de porosité. Selon la lithologie de
Dans les réservoirs, différents types de profils sont utilisés pour estimer la porosité.
par exemple : profil sonore, profil de neutron et/ou profil de densité.

2.Études sismiquesLes études sismiques fournissent un bon

quantité d'informations pour définir la structure des formations
géologiques du sous-sol, mais il est aussi utile pour caractériser les
 propriétés pétrophysiques du réservoir ; parmi ces propriétés pétrophysiques
il y a la porosité, la lithologie, etc.

[Link] directe des échantillons de noyau- cette méthode généralement

enveloppe l'injection d'un gaz ou d'un liquide dans l'échantillon pour déterminer le
volume de pore et ensuite déterminer indépendamment le volume
total de l'échantillon.

5.5 Porosité et enregistrements de puits

L'évaluation de la formation à l'aide de capteurs mesure les propriétés des roches,

ainsi que des propriétés des fluides contenus dans la formation, par
capteurs qui mesurent différents paramètres de la formation, comme par exemple :
propriétés électriques, acoustiques, radioactives, électromagnétiques, etc. Ceci
l'information aide les analystes (géologues, géophysiciens et ingénieurs pétroliers) dans
le calcul du volume d'hydrocarbures que peut contenir la formation. Les
les géologues utilisent également ces informations pour corréler des formations d'un
puits à un autre et de cette manière, il est possible de créer des cartes du réservoir. Aide aux
géophysiciens à calibrer les informations sismiques obtenues auparavant (convertir la
dimension de temps en profondeur). Aide les ingénieurs à recueillir
informations pour créer des modèles 3D du réservoir et réaliser des travaux de simulation.

La porosité d'une formation peut être déterminée par la mesure d'un seul
un enregistrement ou une combinaison de ceux-ci. La combinaison de différents enregistrements aide
à corriger les effets de la variation de la lithologie qui se présente normalement dans
toutes les formations. Les enregistrements utilisés pour estimer la porosité des
les formations sont :

• Enregistrement sonore

• Enregistrement de la densité

• Enregistrement de neutron

• Enregistrement d'IRM

• Registres de résistivité

Dans cette liste d'enregistrements, les 4 premiers sont connus sous le nom de profils de porosité.
et des cinq nommés, aucun ne mesure directement la porosité de la formation.
Pour obtenir la porosité, d'abord d'autres paramètres sont mesurés.
formation qui permet d'estimer la porosité.
5.5.1 Enregistrement sonore

En 1946, le premier profil sonique a été produit, au départ cette technologie a permis
la meilleure détermination de la profondeur et la meilleure localisation des forages
dans les formations productrices. Les profils soniques travaillent généralement
émettant des signaux de 20 kHz à 30 kHz, mais il est possible de trouver des équipements qui
travaillent à des niveaux supérieurs. Les outils mesurent l'intervalle de temps (Δt) que
requiert une onde de compression sonore pour parcourir un pied de distance et se
mesuré en microsecondes par pied (μseg/ft). En général, l'intervalle de temps est
augmenté par la présence d'hydrocarbures dans la formation.

Si un enregistrement sonique est utilisé pour estimer la porosité de la formation, les

les équations utilisées dans ce processus sont les suivantes :

Δt = φ • Δt f+ (1− φ ) Δt ma
Journal
[ Wyllie ]

5 Δt − Δt ma
φ= •
Journal
[ RHG ]
8 Δt Journal

Où :

Δt intervalle de temps mesuré

Journal
φ porosité
Δt fintervalle de temps de voyage du fluide Δt maintervalle de temps de voyage de la matrice

Pour l'estimation de la porosité avec cette équation, on utilise différents

valeurs pour l'intervalle de temps de voyage de la matrice et le temps de voyage du fluide
selon le type de formation et le fluide que contient la formation qui se
il est en train d'étudier. Ci-dessous se trouve un tableau avec les valeurs de l'intervalle
de temps de voyage pour différentes lithologies et fluides.

Tableau 5-1. Intervalles de temps pour différentes matrices et fluides.16

Intervalle de temps
Litologie/Fluide
(Wyllie) (RHG)
Arenisca 55,5 - 51,0 56,0
Calcaire 47,6 49,0
Dolomie 43,5 44,0
Anhidrite 50,0 ---
Sal 66,7 ---
Fer (tuyau de revêtement) 57,0 ---
Lodo à base d'eau douce 189,0 ---
Lodo à base d'eau salée 185,0 ---

16
Schlumberger, 1972.
Les intervalles de temps présentés dans le tableau ont des unités de micro-
secondes par pied et appartiennent à deux équations présentées par différents auteurs,
ces auteurs sont : M. R. J. Wyllie et Raymer-Hunt-Gardner (RGH). Les valeurs de
Les intervalles de temps pour la matrice et le fluide sont obtenus à partir de tableaux et la valeur
pour l'intervalle de temps de l'enregistrement, cela prend des lectures de l'outil. La
outil de capteurs mesure l'intervalle de temps en microsecondes par pied
(µsec/ft) y généralement esta designado comoDTdans les registres.

Problème 5-3. Déterminer la porosité de la formation (calcaire) enregistrée par le

enregistrement sonore qui a eu les lectures suivantes du puits. Dans un intervalle de 4 pieds
se registro unΔtfmoyenne de 61,1μseg/ft. Le puits a été initialement foré avec
lodo à base d'eau douce, après 400 mètres de profondeur, a été remplacé par
un lodo à base d'eau salée jusqu'à atteindre la formation productrice.

Solution

Δt = φ • Δt f+ (1− φ ) Δt ma
Journal
⇒ Δt = φ Δt f+ Δt ma− φ Δt ma
Journal

Δt − Δt ma
⇒
Journal
Δt Journal
= φ (Δt f− Δt ma)+ Δt ma φ=
Δt f− Δt ma

Δt − Δt ma 61,1 µseg/pi− 47,6 µseg/ft

=0,0983 = 9,83 %
Journal
φWyllie = =
Δt f− Δt ma 185 µseg/pieds− 47,6 µseg/ft

5 Δt − Δt ma 5 61,1 µseg/ft− 47,6 µseg/f

=
Journal
φRHG = • • = 0,1381 = 13,81 %
8 Δt Journal 8 61,1 µseg/ft

5.5.2 Enregistrement de densité

Le profil de densité (density Log en anglais) a été introduit en 1962 et appartient au

type de profils de radiation. Dans ce type de profilage, la formation est bombardée
avec des rayons gamma et un lecteur mesure les réflexions des rayons gamma. Le nombre
de rayons gamma qui sont réfléchis et qui retournent à l'équipement de mesure par rapport aux
Des rayons qui ne retournent pas sont liés à la densité de la formation. Vous pouvez mesurer.
la densité totale (partie solide de la roche et fluides) et la densité de la matrice (partie
solide de la roche).

Pour estimer la porosité de ce type d'enregistrements, on utilise l'équation de

balance de masa que se montre ci-dessous :

ρ b= φ • ρ + (1−
f
φ) ρ ma
Où :

ρ bdensité totale φ porosité

ρ fdensité du fluide dans les pores ρ madensité de la matrice

Les valeurs de la densité du fluide dans les pores et de la densité de la matrice se

toman de tableaux (voir tableau 5-2) et la valeur de la densité totale est obtenue des
lectures que prend le dispositif de mesure. La densité totale dans les enregistrements est
représentée par RHOB et ses unités sont généralement présentées en grammes par
centimètre cube (g/cm)3).

Tableau [Link]és des matrices et des fluides en grammes par centimètre cube.17
Lithologie/Fluide ρmaoρf
Arenisca 2,644
Calcaire 2 710
Dolomie 2,877
Anhidrite 2 960
Salut 2 040
Eau douce 1 000
Eau salée 1 150

En général, l'introduction de valeurs erronées de la densité de la matrice présente

un plus grand effet sur l'estimation de la porosité que l'introduction de valeurs
érronés sur la densité du fluide de pore.

Problème 5-4. Déterminer la porosité du registre de densité qui a les

caractéristiques suivantes : RHOB = 2,73 ; formation de grès ; pression de
formation 2112 psia et boue de forage à base d'eau salée.

Solution

ρ b= φ • ρ + (1−
f
φ) ρ ma
⇒ ρ b= φ • ρ + ρf − φ maρ ma

ρ b− ρ
ρ b= φ (ρ −f ρ ) + ρma ma
⇒ φ= ma

ρ f− ρ ma

ρ b− ρ 2,73− 2 644
φ= ma
= = − 0,0523
ρ f− ρ ma
1,0− 2, 644

La porosité aboutit à un nombre négatif lorsque les paramètres utilisés pour

les calculs sont incorrects. Par exemple : lorsque la densité de la matrice est

17
Asquith, G. et Krygowski, D.:Analyse de base des enregistrements de forage, Méthodes AAPG dans la série d'exploration, n° 16, 2, L'Association américaine des géologues pétroliers, Tulsa, Oklahoma, p.
40, 2004.
prise comme calcaire mais c'est de l'anhydrite, la porosité sera négative. C'est
Il est important de noter qu'en obtenant cette porosité, ce n'est pas une indication de mal.
fonctionnement de l'équipement, même avec la valeur incorrecte de la porosité, on obtient
bonne information qui aide à identifier la présence d'anhydrite ou d'autres
minéraux lourds. Un effet similaire peut se produire à cause de la boue de forage.
Lorsque de la barytine est ajoutée pour augmenter le gradient de pression (exemple : 14
lb/gal) de lodo de perforation, peut entraîner des valeurs négatives de la porosité
dérivée à partir de l'enregistrement de densité.18

5.5.3 Enregistrement du neutron

Le profil de neutron (neutron Log en anglais) a ses débuts dans les années 60.
appartient au type de profils de radiation. Le profil des neutrons mesure également les
rayons gamma qui reviennent de la formation. Cela peut mesurer la concentration de
hydrogène dans la formation, avec cela on peut déterminer la porosité et le
contenu de la formation. La technologie des profils de neutrons a évolué
améliorant et maintenant on a le profil de neutron de pulsation et le profil de neutron
spectroscopique.

Dans les premiers enregistrements des neutrons, les échelles étaient exprimées en comptages, mais dans les
de nouveaux équipements les enregistrements mesurent en unités de porosité apparente par rapport
d'un certain type de minéralogie. Le minéral qui est généralement sélectionné est le
calcite, dans ce cas les valeurs de porosité qui sont calculées avec l'enregistrement de
les neutrons seront des valeurs vraies dans les zones de calcaires. En cas de sélection de la
calcite comme minéral de référence et des mesures sont effectuées dans des zones où il n'y a pas
s'ils ont des calcaires, une correction pour le minéral de la zone doit être effectuée ou combinée
l'enregistrement avec un enregistrement de densité pour estimer la véritable porosité. Sans
embargo, on peut également sélectionner le quartz comme minéral pour les formations
de grès. La désignation de l'enregistrement de neutron est PHIN ou NPHI et les valeurs
de porosidad se peuvent être lus directement à partir du registre.

Combinaison d'enregistrements (enregistrement de neutron et enregistrement de densité)

La combinaison des enregistrements de densité et de neutron est utilisée pour estimer un

porosité qui est pratiquement libre des effets de la lithologie de la formation
étudiée. Chaque enregistrement mesure une porosité apparente qui est vraie uniquement
dans le cas où la lithologie de la zone étudiée concorde avec la lithologie
utilisée par l'ingénieur pour mettre l'échelle de l'enregistrement. La combinaison de ces
deux types d'enregistrements permettent d'annuler les effets de la lithologie de la formation et
estimer la vraie porosité en calculant une moyenne des deux lectures
dans les deux enregistrements ou en appliquant l'équation suivante :

18
Asquith, G. et Krygowski, D. :Analyse de base des relevés de forage, Méthodes de l'AAPG en exploration, Nº 16, 2, L'Association Américaine des Géologues Pétroliers, Tulsa, Oklahoma, pag.
40, 2004.
 φn2+ φ 2d
φ=
2

Où :

φnporosité des neutrons φdporosité de densité

Problème 5-5. Déterminer la porosité par combinaison d'enregistrements de neutrons et

de densité. La formation a été forée avec un boue de forage à base d'eau
salada jusqu'à trouver la formation productrice à une profondeur de 3000 pieds. Se
ils ont réalisé plusieurs évaluations, du registre de neutron on a obtenu une porosité de
19,4 pour cent et l'enregistrement de la densité a donné 2,31 g/cm3pour la densité totale. La
la formation est en roche calcaire.

Solution

ρ b− ρ 2,31− 2, 71
φ= ma
= =0,2564
ρ f− ρ ma
1,15− 2, 71

φn2+ φ 2 d 0,1942+ 0,25642

φ= = = 0,2274
2 2

5.5.4 Registre d'IRM

Le profil de résonance magnétique nucléaire (résonance magnétique nucléaire ou RMN en

Ingles) détecte les fluides dans le voisinage du puits mais la partie solide de la
la formation est pratiquement invisible. Ce type de profil travaille en envoyant des pulsations
que polarisent les hydrogènes contenus dans l'eau, le pétrole et le gaz naturel de
poros. À mesure que les hydrogènes polarisés se réorganisent à leur état naturel
induisent des signaux qui sont mesurés par l'outil de RMN, l'amplitude de ceux-ci
les signaux sont directement liés à la porosité de la formation

En plus de la détermination de la porosité, la RMN permet également de déterminer

permeabilité et saturation initiale de l'eau. En général, la perméabilité
déterminée par des mesures de RMN est plus précise que celles obtenues de
correlations entre porosité et perméabilité. Il est important de se rappeler que ce type
de profil n'utilise pas la radiation pour réaliser ses mesures (l'équipe de RMN
ne contient aucune source de rayonnement).19

19
Buller D.:Avancées dans la technologie de logging par câbleConseil de transfert de technologie pétrolière (PTTC), septembre, 2001.
5.5.5 Enregistrements de résistivité

Les outils pour mesurer la résistivité des formations se divisent en deux

types, these are:
Laterologue
2. Induction
Latérologue

Ce type de profils nécessite un circuit électrique complet pour fonctionner.

le courant électrique passe d'un électrode d'origine à travers la formation et
retourne vers l'équipe de profilage qui reçoit le courant électrique par un
électrode de surface ou un électrode de retour. La résistivité est fonction de
décrément de tension entre l'électrode d'origine et l'électrode de retour. Les
Les laterologues ont besoin d'établir un contact électrique avec la formation à travers une boue.
conductif ou par contact physique direct. Ces profils sont capables d'évaluer
formaciones hautement résistives et sont bons pour évaluer les strates d'épaisseur
maigre.20

Induction

Les équipes d'induction n'ont pas besoin de contacter la formation parce que
Ces équipements transmettent des ondes électromagnétiques qui induisent des courants électriques.
dans la formation. Les courants électriques induits sont fonction de la résistivité,
plus la conductivité de la formation est élevée, plus les
courants induits dans la formation. Les courants induits sont mesurés par
alambres récepteurs et transformés en lectures de résistivité. Ces profils
ils travaillent mieux dans des formations hautement conductrices et peuvent fonctionner avec des boues non
conducteurs.

Première équation d'Archie

Archie a publié en 1942 un document proposant deux équations qui décrivent

le comportement de la résistivité dans les roches réservoir. Une de ces équations
relie la porosité de la formation à la résistivité du fluide qui se
trouve dans les pores. Cette équation est connue sous le nom de première équation de
Archie et sous forme d'équation se représente de la manière suivante :

Ro a
F= =
Rw φm

Bonner S., Fredette M., Lovell J., Montaron B., Rosthal R., Tabanou J., Wu P., Clark B., Mills R. et Williams R.:Résistivité lors du forage – Images de la chaîneRevue des champs pétroliers
20

pag. 6, 1996.
Où :

F: facteur de formation

Rorésistivité attendue de la zone, si elle était saturée

complètement avec de l'eau de formation

Rwrésistivité de l'eau

a : facteur de tortuosité

m : exponente de cimentation

φ porosité

Les valeurs du facteur de cimentation (m) augmentent à mesure que la

les arénites sont mieux consolidées, cependant, la tortuosité dans le système de
les pores affectent également les mesures de résistivité. Pour cela, d'autres auteurs
ils ont ajouté un facteur de tortuosité (a).

En dégageant la porosité de la première équation d'Archie, on obtient :

⎡ R w• un ⎤ m
φ= ⎢ R ⎥
⎣ o ⎦

Ceci et la deuxième équation d'Archie sont appliquées de manière courante dans l'analyse
de registres de puits lorsque les formations ne contiennent pas d'argiles. Dans le cas de
Les formations argileuses utilisent des équations qui sont des adaptations de la
équations d'Archie qui tiennent compte de l'effet des argiles sur les
mesures.

Les valeurs des variables sont obtenues à partir des sources suivantes : la valeur de Rwes
la valeur de la résistivité de l'eau à la température de la formation, les constantes
myunsont déterminées à partir d'études pétrophysiques et de la valeur de Roest pris
des mesures de résistivité des enregistrements.

Les valeurs demgénéralement, cela dépend du degré de consolidation de la roche

mais c'est aussi une question d'expérience. La valeur du facteur de tortuosité
généralement c'est 1, cependant, cela peut varier en fonction du personnel qui réalise
les études. Voici une échelle de valeurs pourmlesquels sont
une bonne approximation.
 Tableau [Link] de l'exposant de cimentation pour différents degrés de consolidation.
Niveau de consolidation Exposant de cimentation
Non consolidées 1,3
Très légèrement cimentées 1,4 - 1,5
Légèrement cimentées 1,6 - 1,7
Modérément cimentées 1,8 - 1,9
Altement cimentées 2,0 - 2,2

5.6 Porosité et études sismiques

Les informations obtenues des études sismiques étaient initialement utilisées uniquement
pour définir les structures géologiques du réservoir. Cependant, ceci
les informations sismiques peuvent également être utilisées pour caractériser des propriétés
petrophysiques du réservoir ; une de ces propriétés est la porosité. Le
La caractérisation du réservoir nécessite la construction d'un modèle spatial
que représente le réservoir. Cela se fait en rassemblant des informations de différents
études réalisées précédemment et en les intégrant. Le modèle obtenu de la
la caractérisation du réservoir est utilisée dans les simulateurs d'écoulement ; les simulateurs
de flux aident à prédire le comportement du réservoir. Le plus grand défi dans la
La caractérisation d'un réservoir est l'intégration de tous les types d'informations.
pour obtenir un modèle qui représente de manière adéquate le réservoir.

Les études sismiques 3D fournissent des informations sur une vaste zone de
réservoir, en revanche, les données provenant d'échantillons de noyau et
Les registres de puits fournissent des informations sur des zones spécifiques et isolées de
réservoir. En raison de cette caractéristique des études sismiques 3D, cela a été
développé des méthodes pour utiliser les données sismiques dans la description des
propriétés du réservoir dans des endroits éloignés des puits.

La grande variété d'origines et d'échelles de l'information utilisée dans la

La caractérisation du réservoir doit être prise en considération. En plus de la
différentes échelles est la quantité d'informations dont on dispose ; ces
les considérations présentent des problèmes complexes lors de la phase d'intégration
les informations. Pour résoudre ce type de problèmes, on utilise des outils
computational tools non conventionnelles, parmi ce type d'outils se trouvent les
systèmes intelligentsqui imitent le mécanisme de fonctionnement de l'esprit
humaine.21

Les études sismiques 3D sont un outil qui fournit une représentation

des éléments stratigraphiques et structuraux du réservoir, aidant à localiser
la position de structures compliquées dans le sous-sol. Ce type d'études
permet également d'extraire des informations supplémentaires sur le réservoir. Tout le
Les informations obtenues grâce aux études 3D sont manipulées de manière numérique, cela
permet d'obtenir de nouvelles informations par des moyens mathématiques à partir des

21
Nikravesh. M. et Aminzadeh, F.:Tendances de caractérisation des réservoirs intelligents : passé, présent et futur, Journal of Petroleum Science and Engineering, Vol. 31, p. 67-79, 2001.
signaux sismiques. Les analyses pour obtenir des informations supplémentaires peuvent être par
régression et/ou géostatistique des données sismiques. La modélisation inverse
(régression) des propriétés du réservoir à partir des données sismiques est
connue sous le nom deinversion sismique(seismic inversion in English). These analyses
ils peuvent établir empiriquement des relations entre les caractéristiques des données sismiques (amplitude,
force de réflexion, force de fréquence et phase instantanée) avec les propriétés
physiques du réservoir (porosité, perméabilité, etc.).22

Pour déterminer la porosité d'une formation à partir de données sismiques, on

nécessite une relation entre les propriétés élastiques de la roche et la porosité.
Les relations entre la porosité et la vitesse peuvent devenir compliquées parce que
la rigidité des roches dépend non seulement de la porosité et de la minéralogie, mais aussi
dépend de la micro-structure ou en d'autres termes l'agencement des
composants solides à l'échelle des pores.

5.6.1 Acquisition et traitement

Les études sismiques sont la méthode la plus utilisée lors de l'exploration de

hydrocarbures car ils montrent une vue structurale de la géologie du sous-sol. Le
la théorie de base de la méthode sismique est basée sur des mouvements de signaux qui se
se déplacent à travers les différents strates du sous-sol. La méthode de réflexion est
le méthode la plus utilisée dans les études sismiques. Après la collecte des données
on poursuit avec l'étape de traitement de l'information, cela pour améliorer les
signaux, minimiser l'interférence (bruits) et améliorer la résolution. Enfin, les
Les images traitées sont compilées pour produire le résultat final qui est un
section sismique.

Tant l'acquisition que le traitement des données sismiques ont un impact

dans l'image sismique finale. L'interprétation des images sismiques doit donner
début de la conception de l'étude sismique. Cela est dû au fait que peu importe l'efficacité
que soit la méthode d'interprétation, si la qualité des données acquises est mauvaise,
L'image sismique ne représentera pas adéquatement les caractéristiques du réservoir.
Lors de la conception de l'étude sismique, il faut tenir compte de :

• Malla marine contre malla terrestre

• Caractéristiques du bruit associé aux sédiments de surface

• Facteurs environnementaux tels que la végétation

• Vie animale

• Constructions rurales et urbaines

22
Technologie sismique 3D : un aperçu, Conseil de Transfert de Technologie Pétrolière (PTTC), 1999.
L'acquisition de données sismiques est généralement déterminée par les
les exigences de l'entreprise mais restreintes par les coûts qu'elle est prête à
couvrir la compagnie.

Le traitement de l'information comporte plusieurs étapes. Tout d'abord, on effectue le

désalojo normal(normal move out o NMO en anglais). Cela consiste à corriger les
effets de séparation entre le récepteur et la source qui affectent les temps de
il est arrivé des réflexions. Deuxièmement, on travaille sur leévacuation de pente(tremper
déplacer DMO en anglais). Cette étape du traitement permet de corriger le
différence dans les temps d'arrivée ou temps de voyage des ondes réfléchies,
en plus de prendre en compte les réflexions inclinées.

La méthode utilisée pour améliorer la résolution est connue sous le nom dedéconvolution
(deconvolution en anglais). Une bonne quantité d'algorithmes est disponible pour
réaliser cette tâche sur les données sismiques. Un autre traitement courant des données
sismiques est lamigrationsismique, la migration sismique compense la distorsion
introduite en raison de la propagation des ondes et de la géométrie des points de
acquisition de données. Le traitement de migration est généralement effectué
après l'empilement (tri) des données sismiques, cependant,
récemment, des méthodes de traitement de la migration ont été développées
applicables avant l'empilement des données. Ce dernier est généralement appliqué à
structures géologiques complexes.23

Trazado de profil sismique vertical (profil sismique vertical ou VSP en anglais)

Le VSP diffère des études sismiques conventionnelles par le positionnement de

les récepteurs des signaux. Dans les études de VSP, les récepteurs sont positionnés
dans un puits au lieu de la surface terrestre comme cela se fait dans les études sismiques
conventionnels. Pendant les études sismiques VSP, les strates du sous-sol
ils agissent comme un filtre, c'est pourquoi, en positionnant les récepteurs à une certaine profondeur, on
réduire la distance que les signaux doivent parcourir depuis leur source jusqu'aux
récepteurs. Cela se traduit par des données avec des fréquences plus élevées (meilleure résolution).24

5.7 Mesure directe de la porosité

Les deux méthodes précédentes, les relevés de puits et les études sismiques, déterminent la
porosité de la formation indirectement. La porosité des formations est une
propriété statique, donc elle peut être mesurée en l'absence de flux. Cependant,
Pour déterminer la porosité effective, il faut avoir un flux pour déterminer ce que
pourcentage des pores sont interconnectés.

Les mesures directes de la porosité sont effectuées dans des laboratoires. Pour les
Les mesures utilisent des échantillons extraits de la formation productrice, ces échantillons

23
Imagerie sismique des jeux structuraux, stratigraphiques et diagenétiques, Conseil de transfert de technologie pétrolière (PTTC), 2003.
24
Gadallah, M.R.:Sismologie des réservoirs : Géophysique en langage non technique, PennWell Books, Tulsa, Oklahoma, 1994.
sont connues commemoyennes de noyauLes études sont généralement réalisées en
petits échantillons cylindriques extraits de l'échantillon de noyau, cependant, si
dispose des équipements et pour certaines études, il est possible de réaliser les mesures dans
l'échantillon de noyau. Lorsqu'on utilise des échantillons cylindriques de l'échantillon de
noyau, des mesures doivent être effectuées sur plusieurs échantillons pour obtenir des résultats
statistiquement représentatifs.

Pour mesurer la porosité, on mesure 2 des 3 paramètres suivants :

• Volume total de l'échantillon (VT)

• Volume de la matrice ou volume de grains (VM)

• Volume de pore (VP)

Volume des pores VP VT− VM

φ Total = = =
Volume total VT VT

5.7.1 Échantillon de noyau

Vers 1920, les premiers échantillons de carottes de puits ont été extraits aux États-Unis.
(Californie, Texas et Colorado). Les outils utilisés coupaient et extrayaient
échantillons de noyau du fond du puits et puis étaient analysés en laboratoire
pour obtenir des propriétés des formations qui étaient forées. Les échantillons de
Le noyau est des sections de roche qui sont extraites du puits arrêtant les travaux de
perforation ou au cours des travaux de forage.

Les tests realizados en le laboratoire sur les échantillons de noyau peuvent être
vous contrôlez avec une grande précision, cependant, le processus de collecte des
échantillons de noyau et la préparation de ceux-ci pour les tests de laboratoire peuvent
affecter considérablement les mesures effectuées au laboratoire. D'autre part
D'autre part, la collecte des échantillons de noyau peut s'avérer coûteuse à
grands champs d'hydrocarbures. La collecte d'échantillons de noyau est
inévitable peu importe si le réservoir est petit ou grand, mais le bon
La représentation des propriétés est complexe en raison des propriétés de
Les réservoirs changent en fonction des coordonnées spatiales.

En raison des coûts engendrés par l'extraction et l'analyse des échantillons de carottage
ces choses ne se réalisent pas dans l'ensemble du puits. Il est donc nécessaire de déterminer les
intervalles où se fera l'extraction de noyaux. Par l'utilisation de
cartes géologiques et coupes transversales de la formation sont déterminées par les
profondeurs d'où l'on souhaite extraire des échantillons de carottes pour une utilisation ultérieure
analyse en laboratoire. L'analyse des échantillons de noyau peut être réalisée dans
le noyau complet ou en petits échantillons cylindriques ou petits noyaux. La
l'utilisation de petits échantillons se fait généralement lorsque la roche est
homogène, sinon il est préférable de faire les mesures sur l'échantillon de
noyau.

5.7.2 Préparation du noyau

Pour obtenir les propriétés en utilisant les échantillons de cœur, on coupe

plusieurs petits échantillons cylindriques de l'échantillon de noyau. Les échantillons
Les cylindres ont environ les dimensions suivantes :

• Diamètre : 1,5 pouces

• Longueur : 2,5 pouces

La mesure des propriétés est réalisée sur tous les échantillons cylindriques, mais
Auparavant, les échantillons sont préparés. Après avoir sélectionné et
Pour couper les échantillons cylindriques, il faut les nettoyer complètement. Il est nécessaire
nettoyer les échantillons car ils contiennent du fluide de forage, des résidus de
pétrole brut, eau, sels, etc. Mais en même temps, il faut faire attention à ne pas
changer les caractéristiques de la roche pendant le processus de nettoyage et de séchage. Pour
Le nettoyage des échantillons utilise différents fluides, parmi lesquels figurent :

• Toluène (C7H8)

• Toluène saturé avec CO2

• Hexano (C6H14)

• Méthanol (CH4O)

• Trichloroéthylène (C2HCl3)

La sélection du fluide nettoyant dépend des caractéristiques de la roche et ce

doit être sélectionné de manière à ne pas altérer les caractéristiques de la roche. Les
les échantillons sont immergés dans le fluide pour enlever les contaminants, ensuite, on
procède à sécher les échantillons cylindriques. Après cette procédure les
Les échantillons sont prêts à obtenir les propriétés de la roche. Les carottes doivent être
exposés à des températures d'environ 100 ºC pour enlever l'humidité et
solvant se trouvant dans les pores.

5.8 Détermination du volume total du noyau

Le volume total du noyau est composé du volume occupé par les grains de
roche également connue sous le nom de matrice ou partie solide et le volume qu'ils occupent
pores et fluides.
 VTotal=VMatrice+VPoro

Pour la détermination du volume total du noyau, on utilise généralement deux

méthodes, ces méthodes sont :

1. Déplacement volumétrique

2. Mesure des dimensions du noyau

5.8.1 Déplacement volumétrique

Le déplacement volumétrique se fait en utilisant un liquide comme du mercure.

paraffine ou tout fluide qui sature ou déplace l'échantillon de noyau. Il y en a deux
méthodes pour ce type de mesures, ces méthodes sont :

• Méthode par gravité

• Méthode volumétrique

Méthode par gravité (méthode d'Archimède)

Cette méthode permet de déterminer le volume total de l'échantillon par le biais de la

mesure du poids de l'échantillon dans différentes conditions. Le volume total de la
L'échantillon par cette méthode est déterminé par l'équation suivante :

WSam. − W Somme.
VTotal =
gjeρ Fluide

Pour mesurer la porosité par cette méthode, il faut suivre les étapes suivantes :
Après avoir obtenu les échantillons, les nettoyer et les sécher, on pèse l'échantillon sec et
nettoyer (pores non saturés). Ensuite, on procède à les plonger dans le fluide
sélectionné pour l'essai. En plongeant les échantillons dans le fluide, ils se saturent.
poros de l'échantillon ; une fois que les échantillons sont saturés, ils sont pesés de manière normale (au
aire libre). Ensuite, l'échantillon est immergé dans un récipient contenant le même fluide et se
le pèse dans le récipient. Une fois ces 3 poids obtenus pour les différents
On utilise l'équation suivante pour déterminer la porosité de chaque échantillon.
échantillon.

WSam. - WSeco m Sam.-mSeco

φTotal = ó =
φTotal
WSam. - WSomme. m Sam.-mSomme.

Lors de la mesure des poids, des erreurs s'introduisent. Par exemple : la plus grande
la source d'erreurs se présente lors de la mesure du poids saturé des échantillons.
Après avoir saturé les échantillons avec le fluide, ils doivent être pesés, l'échantillon peut
avoir un fluide adhéré à sa surface externe, ce qui introduit un certain degré d'erreur
dans la mesure du poids et donc de la porosité. S'il se sèche trop
des échantillons, il est possible d'extraire le liquide des pores, donc l'échantillon ne serait pas
complètement saturée et un erreur serait introduite dans la mesure de la porosité.

Figure 5-2. Étapes pour la mesure de la porosité par la méthode d'Archimède.

Dédaction de l'équation pour calculer la porosité

Pour la déduction, se référer à la figure 5-2; de (1) nous avons :

WMatrice= W Seco

De (2) nous avons :

WFluide
ρ Fluide = ; VFluide=VPoro ; WMatrice= W Seco
giV Fluide

WSam. =WMatrice+WFluide ⇒ WFluide= W −[Link]

WFluide WSam. − WMatrice WSam. − W Seco

VPoro = = =
giρ Fluide
giρ Fluide
giρ Fluide

De (3) nous avons :

WSomme.= WSeco − F Flottaison ; FFlottabilité =WFluide

 WFluide= W Seco
− W Somme.

WFluide WFluide WSeco − W Somme.

ρ Fluide = ; VFluide =VMatrice ⇒ VMatrice = =
gjeVFluide gjeρ Fluide gjeρ Fluide

Alors, le volume total est :

WSam. − W Seco WSeco − W Somme. WSam. − W Somme.

VTotal =VPoro+VMatrice = + =
giρ Fluide
giρ Fluide
giρ Fluide

WSam. − W Seco
VTotal− VMatrice VPoro gjeρ
φTotal = = = Fluide

VTotal VTotal WSam. − W Somme.

gjeρ Fluide

WSam. − W Seco
φTotal =
WSam. − W Total.

En termes de masse

gjem Sam.− gjem Seco m Sam.− m Seco

φTotal = =
gjem Sam.− gjem Somme. m Sam.− m Somme.

Méthode volumétrique

Pour cette méthode, on utilise un picnomètre à mercure, le picnomètre utilise le

déplacement de mercure pour déterminer le volume total de l'échantillon
en se basant sur le fait qu'à pression atmosphérique, le mercure n'entrera pas dans
pores avec des diamètres inférieurs à 15 microns (micron = 1 millionième de mètre).
Cela est dû au fait que le mercure est un fluide pour la plupart des roches, non.
humectant. La tension interfaciale entre le mercure et la matrice de l'échantillon ne
permettent au mercure d'entrer dans les pores ou les fissures jusqu'à ce qu'un certain
une pression externe soit exercée sur le mercure pour entrer dans l'échantillon.

Le picnomètre est un conteneur d'un volume connu; bien que le picnomètre

est utilisé pour déterminer la densité, ce que le picnomètre mesure est le volume.
Tout d'abord, l'équipement est calibré à l'aide du piston ; une pression est exercée avec le piston sur
le mercure (sélectionnant une force de compression de sorte qu'elle puisse être répétée
à chaque test) pour établir une ligne de base. Deuxièmement, l'échantillon est introduit dans
l'équipe avec le mercure le processus de compactage est répété avec le même
force de compression utilisée lors de l'étape de calibration. La différence entre les
Les positions (A) et (B) sont utilisées pour calculer le volume total de l'échantillon
en utilisant l'équation suivante (voir figure 5-3) :

VTotal= π jer 2jeh

Figure 5-3. Détermination du volume d'un échantillon par déplacement de mercure.

Déduction de l'équation pour calculer le volume total (pycnomètre)

Le volume de mercure durant le test de calibration pour établir la ligne

base es

VMercure - 1 = π jer 2jeh 0

Le volume de mercure plus le volume total de l'échantillon est

VMercure + montre = π jer 2jeh +[ 0h ]

Le volume de l'échantillon est obtenu en soustrayant le volume de mercure

( [ 0h
VTotal de l'échantillon = π jer 2jeh + ]) − (π jer jeh )
2
0

En simplifiant l'équation

(
VTotal de l'échantillon = π jer jeh+
2

0
)π jer (jeh
2
) − (π jer jeh )
2
0

VTotal de l'échantillon = π jer 2jeh

Problème 5-6. Déterminer la porosité d'un échantillon de carotte qui a été
immergée dans l'eau salée. L'échantillon a été pesé sec et a enregistré un poids de 21,7
grammes, puis, il s'est saturé d'eau salée et pesait 25,6 grammes. Enfin, le
L'échantillon saturé a été immergé dans de l'eau salée et pesé, enregistrant 13,5 grammes.

Solution

m Sam.− m Seco 25,6− 21,7

=
φTotal = =0,3223 = 32,23 %
m Sam.− m Somme. 25,6− 13,5

5.8.2 Mesure des dimensions du noyau

Une autre méthode pour déterminer le volume total du noyau est par la mesure
de ses dimensions. Généralement, le noyau a une forme cylindrique donc
il suffit de mesurer le diamètre et la hauteur du noyau pour déterminer son
volume total, en utilisant l'équation du volume d'un cylindre.

π
VTotal= π • r 2• h ó VTotal = • d 2• h
4

Où :

r : radio de la muestra

h : hauteur de l'échantillon

d : diamètre de l'échantillon

La détermination du volume total des échantillons par cette méthode a une

exactitude limitée. Cela est dû au fait que les échantillons peuvent contenir des irrégularités
à sa surface, ce qui se produit généralement dans la pratique. Ces irrégularités ne
sont prises en compte dans les calculs du volume total des échantillons, par donc
tant de degré d'erreur est introduit dans les calculs.

5.9 Détermination du volume de la matrice

Le volume de la matrice se réfère à la partie solide de l'échantillon de noyau. Cette

Une partie de l'échantillon est composée de différents matériaux selon le type de
roche. Les types de roche les plus importants dans l'industrie pétrolière sont :

• Arenisca

• Calcaire
 • Dolomite ou dolomie (dolomite en anglais)

Le volume de la matrice d'échantillons de noyau peut être déterminé de trois

formes, celles-ci sont :

1. En utilisant la densité de l'échantillon

2. Expansión de gaz (Lois de Boyle)

3. Méthode de déplacement

5.9.1 En utilisant la densité de la roche

La densité de la roche qui compose le noyau peut être utilisée pour la

détermination du volume de la matrice. La densité peut être estimée à partir de
analyse de la lithologie de l'échantillon. Par exemple, la valeur approximative de la
la densité d'un grès est de 2,644 g/cm3. En utilisant une balance, on peut
déterminer la masse de l'échantillon, avec la masse et la densité du noyau nous pouvons
calculer le volume de la matrice en utilisant l'équation suivante :

Masse du noyau m Noyau

ρ Matrice = =
Volume de la matrice VMatrice

m Noyau
VMatrice =
ρ Matrice

L'exactitude des mesures du volume de la matrice par cette méthode dépend

de la valeur de la densité de l'échantillon. La densité de l'échantillon peut être mesurée
ou peut être assumée. Si la densité de l'échantillon est connue, le degré de
la précision de la méthode est élevée. En revanche, si la densité de l'échantillon est
assumée, la précision de la méthode est limitée. Cela est dû au fait que les échantillons de
Les roches sont constituées de différents composants qui à leur tour ont différents
matrices avec des densités hétérogènes, par conséquent, la densité doit être mesurée.

Tableau [Link]é des roches communes dans l'industrie pétrolière.25

Densité
Lithologie/Fluide
(g/cm3)
Arenisca 2 644
Calcaire 2 710
Dolomie 2 877
Anhydrite 2 960
Sal 2 040

25
Asquith, G. et Krygowski, D.:Analyse de base des enregistrements de puits, Méthodes AAPG en exploration, Nº 16, 2, L'Association Américaine des Géologues Pétroliers, Tulsa, Oklahoma, pag.
40, 2004.
Problème 5-7. Déterminer le volume de la matrice de l'échantillon de noyau
obtenue d'un réservoir formé par des strates de calcaire. L'échantillon a une
masse de 18,4 grammes.

Solution

m Noyau 18,4 g
VMatrice = = = 6,7897 cm3
ρ Matrice g
2,71
cm3

5.9.2 Expansion de gaz (Loi de Boyle)

Robert Boyle a étudié la relation entre la pression et le volume d'un gaz à

température constante. R. Boyle a déterminé que le produit entre la pression et le
le volume est à peu près constant. Pour un gaz idéal, le produit de la pression
et le volume est exactement constant. La loi de Boyle sous forme d'équation est :

PGazVGazConstante ⇒ P1 V1 = P2 V2

La méthode d'expansion des gaz utilise cette loi pour déterminer le volume de la
matrice des échantillons de noyau. Pour la mesure des pressions et des volumes se
ils utilisent deux cellules comme indiqué dans la figure 5-4. Dans les conditions initiales, il
il y a du gaz dans la cellule A et l'échantillon de roche est positionné dans la cellule B. Ensuite, on
évacuez l'air de la cellule B.

Le conduit qui relie les cellules A et B a un diamètre assez petit de sorte que
que l'on peut négliger le volume de gaz qui se trouve dans ce conduit. Pour
Ainsi, les seuls volumes à prendre en compte sont les volumes des cellules.
lesquelles sont connues. Dans cette condition, le gaz aura une pression initiale (P1)
que se doit à la pression exercée par le gaz piégé dans la cellule A. Dans les
conditions finales on laisse le gaz de la cellule A se dilater en ouvrant la valve. Le gaz
s'étend et remplit la cellule B et les pores interconnectés de l'échantillon (voir figure
5-4). Dans ces conditions, une deuxième pression sera enregistrée (P 2); cette pression sera
inférieur à la pression enregistrée dans les conditions initiales. Alors nous aurons deux
pressions et deux volumes de gaz. Le premier volume de gaz sera égal au volume
de la cellule A, le second volume de gaz sera égal au volume de la cellule A plus
le volume de la cellule B moins le volume de la matrice de l'échantillon.
Représenté sous forme d'équations :

VGaz - 1 =VCellule A

VGaz - 2 =V Cellule B
+V Cellule A
-VMatrice
 Figure 5-4. Étapes de la méthode pour déterminer le volume de la matrice par la méthode d'expansion de gaz.

Introduisant les deux équations dans la loi de Boyle

PV
1 Gaz - 1
= P VGaz
2 -2
⇒ PV
1
Cellule A
=P 2
(V Cellule B
+V Cellule A
− VMatrice )
Déterminer le volume de matrice

− PV
1Cellule A
+P 2
(V Cellule B
+VCellule A )
VMatrice =
P2

Pour appliquer cette méthode, il faut supposer que la température du système reste constante.
constante et qu'il n'y a aucun gain ou perte de gaz, ce dernier signifie que le
nombre de molécules de gaznc'est pareil pendant toute l'expérience.

L'exactitude de cette méthode est généralement élevée, mais elle peut être affectée par le
type de roche et la quantité d'humidité que contient le gaz. Pour réaliser le
l'échantillon d'expérience doit être exempt de composés volatils qui pourraient
contribuer à des erreurs dans les mesures. De plus, le gaz utilisé dans le
l'expérience doit être un gaz sec ou un gaz pur. La contribution des erreurs par parte
L'équipe se limite généralement aux fuites et à l'instabilité de la température. Autre
Un facteur important qui introduit des erreurs dans la mesure est la porosité inefficace.
les pores inefficaces sont isolés des autres pores. Par conséquent, le gaz ne
entre dans ces pores et ce volume est pris comme s'il faisait partie du volume
de la matrice de l'échantillon.

Le gaz couramment utilisé dans les tests est l'hélium parce qu'il se dilate.
rapidement par les pores de l'échantillon. D'autres gaz sont également utilisés, ceux-ci
ils sont sélectionnés en fonction de la façon dont le gaz réagit avec la matrice de la
échantillon.

Problème 5-8. Déterminer la porosité d'un échantillon soumis à une

essai d'expansion de gaz (loi de Boyle). Le gaz utilisé dans l'essai est l'hélium
et est contenu dans un cylindre dont les dimensions sont : diamètre 4 cm et longueur 15
cm. La pression à laquelle le gaz est contenu dans le cylindre est de 54 psia à une
température de 15 ºC. L'échantillon est placé dans un autre compartiment solide
connecté au cylindre par une vanne. Puis, la vanne est ouverte et un enregistrement est effectué.
pression de 52,98 psia et sa température reste constante.

Solution

P1 V1 = P2 V2

⎛ π • 4 2• 15 =⎞ 52,98 psi• V 54 psiπ(• 4• 15 )

54 psi ⎜ ⎟ ⇒ V2 =
⎝4 ⎠
2
52,98 psi

V2 = 192,11 cm3

⎛ π • 4 2• 15 =⎞ 3,61 cm3
V2 = VPoro + VCylindre ⇒ VPoro = 192,11 cm3− ⎜ ⎟
⎝4 ⎠

π π
VTotal = • D 2• h = • 32• 5,1 = 36,05 cm3
4 4

VPoro 3,61 cm3

φ Efficace = = = 0,10 = 10 %
VTotal 36,05 cm3
5.9.3 Méthode de déplacement (désagrégation de l'échantillon)

Cette méthode nécessite le désagrégat des échantillons de roche. Les échantillons sont
désagrégés jusqu'à obtenir des grains de roche lâches. Ces grains fournissent le
volume de la matrice par les méthodes de déplacements décrites
anteriorement comme sont : méthode par gravité et méthode volumétrique.

Cette méthode est considérée comme ayant une plus grande précision par rapport à la méthode
utilisé pour déterminer le volume total car pendant le pesage des grains
Des agrégats de roche, les fluides occupant les pores inefficaces ne sont pas pesés. Cela
Cela signifie que seule la partie solide de l'échantillon est prise en compte.

L'exactitude de cette méthode dépend du processus de désagrégation. Si

durant ce processus, les échantillons sont désagrégés jusqu'à détruire les grains, la
la précision est considérablement réduite.

Problème 5-9. Deux échantillons de roche gréseuse (de la même formation mais
extraits de différents puits) ont été désagrégés et introduits dans l'eau
déplaçant les échantillons A et B de 7,7 et 7,1 cm3respectivement. Avant d'être
désagrégées, les mesures de chaque échantillon ont été prises, le volume total de chaque
3
respectivement. Déterminer la porosité moyenne de la
l'échantillon mesure 9,9 et 8,9 cm
formation de grès à partir des deux échantillons. La couche de grès a un
grosseur de 8 et 9,6 pieds dans les zones d'où les échantillons A et B ont été extraits
respectivement.

Solution

VP VT− VM 9,9− 7,7

φ Un= = = =0,2222
VT VT 9,9

VP VT− VM 8,9− 7,1

φ B= = = =0,2022
VT VT 8,9

∑ (φ • h ) je je
( 22,22je8,0) + ( 20,22je9,6)
φ Prom. = i=1
n
= = 21,13 %
8,0 + 9,6
∑h je
i=1

5.10 Détermination du volume des pores

Le volume de pore des échantillons peut être obtenu par des méthodes préalablement
décrits tels qu'ils sont :
 1.Méthode de déplacement (par gravité)

2.Méthode par expansion de gaz (Loi de Boyle)

5.10.1 Méthode de déplacement

On peut déterminer le volume de pore à la fois par la méthode de gravité et par le

méthode volumétrique.

Méthode par gravité

Dans la déduction de l'équation pour déterminer le volume total de l'échantillon se

trouve l'équation suivante :

WFluide WSam. − WMatrice WSam. − W Seco

VPoro = = =
gjeρ Fluide
gjeρ Fluide
gjeρ Fluide

Cette équation peut être utilisée pour déterminer le volume des pores, en termes de
de masa es :

m Sam. − m Seco
VPoro =
ρ Fluide

Méthode volumétrique

Pour déterminer le volume des pores par cette méthode, certaines doivent être réalisées.
modifications dans l'expérience. La première modification est effectuée dans le
pression appliquée au mercure. Le mercure n'entre pas dans les pores de l'échantillon à
pressions basses. Pour déterminer le volume de pore, la pression est augmentée sur
le mercure. Cette augmentation de pression fait qu'il pénètre dans les pores.

La deuxième modification concerne l'équipement utilisé pour l'expérience.

Pour déterminer le volume des pores, on utilise une cellule connectée directement.
à l'échantillon de noyau (voir figure 5-5). Comme le montre la figure, la cellule B est
solide et n'a d'espace que pour contenir l'échantillon de roche.

Le processus est le même que celui utilisé pour déterminer le volume total de la
muestra. D'abord, on calibre l'équipement et on détermine une ligne de base, puis, on ouvre
la vanne et on applique une pression sur le piston. Les pressions sont élevées, bien au-delà
Au-dessus de la pression atmosphérique, le mercure entre dans le conduit qui communique.
les cellules A et B. Le mercure pénètre dans les pores interconnectés de l'échantillon et
on enregistre la distance parcourue par le piston (la valeur deh). Avec cette valeur, on ...
vous pouvez calculer directement le volume de pore avec l'équation suivante :
 VPoro= π jer 2jeh

Figure 5-5.Détermination du volume de pore par la méthode de déplacement volumétrique.

5.10.2 Méthode d'expansion de gaz

La méthode d'expansion des gaz peut également être utilisée pour déterminer le
volume des pores des échantillons de roche. La procédure est exactement le
même utilisé dans l'expérience pour déterminer le volume de la matrice. Sans
embargo, la configuration de l'équipe est différente ; elle ressemble plus à la configuration
utilisée par la méthode de déplacement volumétrique (figure 5-5).

Pour déterminer le volume de pore par expansion d'un gaz, la cellule B est une
cellule solide qui n'a de l'espace que pour l'échantillon. La cellule A contient le gaz et
quand la vanne est ouverte, le gaz se dilate dans les pores de l'échantillon. La
la pression avant d'ouvrir la vanne est enregistrée et aussi après que le gaz se
s'étendre dans les pores de l'échantillon (voir figure 5-6). Avec ces pressions et le
Le volume de pore peut être déterminé à partir du volume initial connu.
 Figure 5-6. Détermination du volume des pores par la méthode d'expansion de gaz.

Les volumes de gaz dans les conditions initiales et finales sont :

VGaz - 1 =VCellule A

VGaz - 2 =VCellule A
+V Poro

Introduisant les deux équations dans la Loi de Boyle

PV
1 Gaz - 1
= P VGaz
2 -2
⇒ PV
1 Cellule A
=P 2
(V +V
Cellule A Poro
)
Déterminer le volume de pore

PV
1 Cellule A
− PVCellule2 A
VPoro =
P2
Problème 5-10. Détermine la porosité d'un échantillon de carotte de roche calcaire
que s'est passé lors de l'expérience d'expansion de gaz pour mesurer son volume de pore et
volume de matrice. Pendant l'expérience pour déterminer le volume de la
La matrice de l'échantillon a calculé un volume de 21,2 cm3 et a enregistré 100 psia et 64,9
psia comme pressions initiale et finale respectivement. Dans l'expérience pour
Pour déterminer le volume de pore, le même équipement a été utilisé, seul le...
cellule B. Au cours de cette expérience, le même type et la même quantité ont été utilisés
gaz. Dans cette expérience, une pression de 84,9 psia a été enregistrée après avoir laissé
élargir le gaz dans l'échantillon. (La cellule A a un volume de 25,0 cm3)

Solution

PV
1 Cellule A
− PVCellule2 A
VPoro =
P2

100 psiai25,0 cm3− 84,9 psiai 25,0 cm3

VPoro = =3,9662 cm3
84,9 cm3

VPoro 3,9662 cm3

φ= = =0,1576 = 15,76 %
VTotal 3,9662 cm3+21,2 cm3

/ ÊÀi ÛiÊÌ ÃÊ Ì Vi]ÊÛ Ã Ì\Ê

ÜÜÜ° Vi °V ÉÕ V ° Ì