Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Ecoulement des fluides en milieux poreux et rhéologie

1 Physique des roches

1.1 INTRODUCTION

L’étude des roches couvre, selon les auteurs, un domaine plus ou moins vaste. Il s’agit, dans ce
manuel, de la description du milieu poreux (grains, pores et fissures) et de la quantification de
ses principales propriétés physiques. Une roche peut être définie comme un assemblage de
minéraux ayant hérité de liaisons plus ou moins fortes au cours de son histoire géologique. La
description d'une roche se fait sur un échantillon observé à la loupe et éventuellement au
microscope optique polarisant ou électronique à balayage. On décrit la texture, les minéraux
présents, la taille des grains, les fissures et les altérations.

1.2 LES PRINCIPALES CATÉGORIES DE ROCHES

Le nom d'une roche est donné en fonction de sa composition minéralogique, de sa texture et de

son mode de formation. Ce dernier conduit à distinguer trois grandes familles de roches :
sédimentaires, métamorphiques et magmatiques.

TABLEAU 1 : PRINCIPALES ROCHES MAGMATIQUES

Les roches magmatiques, résultent de la solidification de magmas. Ceux qui parviennent

directement à la surface forment les roches volcaniques ou effusives, les laves, dans lesquelles
on n'observe à l'oeil nu que très peu de cristaux (en raison d’un refroidissement rapide). Au
contraire la solidification en profondeur donne les roches plutoniques dont la lenteur du
refroidissement permet la croissance des cristaux, bien visibles donc à l’œil nu. Les minéraux
sont essentiellement des silicates : quartz, feldspaths alcalins et plagioclases, amphiboles, micas
etc..

Les classifications font intervenir la texture de la roche (arrangement et taille des cristaux) et la
composition minéralogique. Associées aux roches métamorphiques, les roches magmatiques,
dont les plus répandues sont les granites, forment l’essentiel de la croûte terrestre continentale,
bien qu’elles soient souvent cachées sous les roches sédimentaires.

1
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Les roches sédimentaires sont formées à la surface de la terre sur le sol ou au fond de l'eau. Elles
résultent de la désagrégation des roches préexistantes. Leur dépôt en couches initialement
proches de l'horizontale résulte :

- de l'action des agents d'érosion et de transport (eaux, crues, glaciers, vent) qui déposent
finalement des roches détritiques ;

- de l'activité des êtres vivants (roches organogènes) ;

- de phénomènes physico-chimiques (roches hydrochimiques et salines formées essentiellement

par précipitation).

Ces actions se combinent souvent et beaucoup de roches sédimentaires sont en fait d'origine
mixte). Ces sédiments se consolident et se cimentent lorsqu'ils sont enfouis en profondeur
(l’augmentation de la température et de la pression favorisant des réactions entre les minéraux et
les fluides interstitiels). Les roches sédimentaires couvrent 75 % de la surface des continents et
la quasi-totalité des fonds océaniques, mais leur épaisseur est limitée.
Il y a de nombreuses classifications des roches sédimentaires fondées sur la composition
chimique, la granularité des constituants ou le mode de formation qui sont présentées dans les
ouvrages de géologie: grès, calcaires, roches argileuses, évaporites, les roches carbonées
(pétroles, charbons et lignites).

1.3 LE MILIEU POREUX

1.3.1 DÉFINITION ET MESURE

Dans les roches la matière minérale forme un squelette solide qui ne remplit pas tout l’espace, et
dont le complément est appelé vide. La proportion de vide est appelée porosité (du grec Ποροσ,
passage, ou détroit, entre les îles). La forme des vides, leur taille, leur répartition, leurs liaisons
ou au contraire l’isolement de certains, influent sur le comportement mécanique et sur les
propriétés de couplage.

[Link]ésentation d’un milieu poreux

2
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

La porosité ɸ est, par définition, le rapport du volume des vides Vv au volume total Vt :

1.4 PRINCIPALES MÉTHODES DE MESURE DE LA POROSITÉ

Volume total

1. Mesure directe des dimensions de l'éprouvette

2. Poussée d'Archimède dans le mercure, qui ne mouille pas la roche et ne pénètre pas dans les
pores sans une forte pression

Volume des vides

3. Absorption d'un fluide mouillant par saturation sous vide

Volume du solide

4. Mesure au pycnomètre de la masse volumique du solide, après broyage

5. Compressibilité des gaz parfaits, en réalisant une détente isotherme en reliant une enceinte V1
(contenant l'échantillon) de volume V1 à la pression P1 à une enceinte vide de volume V2
(pycnomètre à hélium). Après détente, la pression d'équilibre P2 permet de calculer le volume
d’une masse connue de poudre ou de roche.
La porosité est de l’ordre du centième pour certains marbres et quartzites, du dixième pour
beaucoup de roches sédimentaires, elle peut atteindre 0,5 pour certaines craies et tufs. La forme
des vides est représentée en première approximation par un coefficient de forme, rapport de la
plus petite dimension à la plus grande. Ceci permet de distinguer les vides de type pore dont le
coefficient de forme est entre 10-1 et 1) et les vides de type fissure, dont le coefficient de forme
est très faible, entre 10-2 et 10-4, et dont l’épaisseur est négligeable.
La part des pores dans le volume des vides est prépondérante, mais par contre les fissures
contribuent pour l’essentiel à la surface spécifique.
L'espace poreux peut être étudié de manière directe par observation au microscope optique ou
électronique, éventuellement après remplissage par un produit colorant.
Les techniques d’analyse d'images permettent ensuite une estimation quantitative de la porosité.

3
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure 2. Porosité en fonction de la nature de la roche et la profondeur

Calcul de la porosité

Un échantillon de roche propre et sec pesant 425 g a été saturé à 100% avec de la saumure
de gravité spécifique (y) = 1,07. Le nouveau poids est de 453 g. L'échantillon de base mesure 12
cm de long et 4 cm de diamètre. Calculer la porosité de la roche échantillon.

Solution

Le volume global de l’échantillon est

Le volume des pores est

4
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Au fur et à mesure que la profondeur augmente, la porosité démunie, cela est dû à l’effet
conjugué de la pression et de la température lors de la 1iere formation géologique des
roches

Figure 3. Relation moyenne entre la porosité et la profondeur de l'enfouissement pour les

schistes et les températures et les profondeurs pour la genèse du pétrole et du gaz

1.4.1 SATURATION DES FLUIDES

La porosité d'une roche réservoir est très importante car elle permet de mesurer la
capacité de cette roche à stocker des fluides (pétrole, gaz et eau).
Il est tout aussi important de déterminer dans quelle mesure les pores de la roche sont
remplis de fluides spécifiques. Cette propriété est appelée saturation en fluide et est
exprimée en tant que fraction, ou pourcentage, du volume total de pores occupé par
l'huile, le gaz ou l'eau. Ainsi, par exemple, la saturation en huile So est égale à:

5
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Des expressions similaires peuvent être écrites pour le gaz et l'eau. Il est évident que:
Idéalement, en raison de la différence de densité de fluide, un réservoir de pétrole est
formé de telle sorte que, du haut vers le bas du lit de sable, il y aura du gaz, de l'huile et
de l'eau. L'eau de fossile, cependant, est presque toujours trouvée à travers le réservoir
de pétrole. L'eau de fossile est l'eau de mer piégée dans les espaces poreux des
sédiments pendant leur déposition et lithification, bien avant que le pétrole ne migre
dans la roche réservoir.
Outre la densité, la mouillabilité et la tension interfaciale se combinent pour modifier la
manière dont les trois fluides sont distribués dans le réservoir.

La quantité d'eau fossile présente dans l'espace poreux varie de 100% en dessous de la
zone d'huile à théoriquement nulle à des hauteurs au-dessus du niveau d'eau libre.
Cependant, dans les cas pratiques, une teneur presque constante en eau fossile (Si,)
existe au-dessus de la zone de transition. L'amplitude de Si et la hauteur de la zone de
transition dépendent de la taille et la texture des pores. Les valeurs élevées de Siw sont
indicatives de petites tailles de pores. La zone de transition correspond à la zone de
saturation en eau variable. Les puits forés dans cette zone produiront des hydrocarbures
et de l'eau, et les puits forés au-dessus de cette zone, c'est-à-dire dans la zone de
saturation en eau, ne produiront que des hydrocarbures.

UTILISATION QUANTITATIVE DE LA POROSITÉ

L'une des méthodes les plus simples de calcul de la teneur en huile du réservoir est
appelée la méthode volumétrique. L'expression mathématique pour l'huile-in-place
initiale (N en bbl) par cette méthode est:

6
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

L'équation 1.4 donne le volume d'huile contenu dans la roche poreuse aux conditions de
pression et de température du réservoir. Cependant, l'huile de surface ou "réservoir de
stockage" finalement vendue par le producteur est différente du volume de liquide qui
était sous terre.
La différence est due aux changements dans les propriétés de l'huile à mesure que la
pression diminue de la haute pression et température souterraine à la pression et
température de surface. Cette réduction de p et T provoque la sortie de certains
composants volatils (évaporation), ce qui provoque un rétrécissement du volume de
liquide. Cette réduction de volume est exprimée par le facteur de volume de formation
d'huile, Boi. Ainsi, l'huile de réservoir de stockage initialement en place est:

Où Boi est en barils de gisement par baril de réservoir de stockage ou bbl/STB.

Dans cette équation, Soi est remplacé par (1 - Si,), où Si, est la saturation irréductible ou
eau fossile. Ceci implique qu'aucun gaz libre n'est présent dans l'espace poreux. Étant
donné qu'aucun réservoir de pétrole n'est homogène, il faut faire la moyenne des facteurs
As, h, ɸ et Siw. La constante 7.758 devient 10.000 si As, h sont exprimés en hectares (1
hectare = 10.000 m2) et m, respectivement, et N en m3.

EXEMPLE

Calculer le pétrole initial in place (N) d'un réservoir d'huile si A = 1,600 acres,
H = 32 ft, ɸ = 22%, Siw = 20% et Boi = 1,23 bbl / STB.

SOLUTION

En utilisant l'équation 1.5, nous avons:

N = 7,758 (1,600) (32) (0,22) (1- 0,20) / 1,23 = 56,8 x l06 STB

Une expression semblable à l'équation 1.5 peut être dérivée pour estimer le gaz initial in
situ. Dans ce cas, il est commode d'exprimer le volume de gaz en pieds cubes. Dans des
conditions standard, c'est-à-dire Psc = 14,7 psia et Tsc = 60 ° F, le gaz initial en place
dans un réservoir volumétrique est donné par:

Où Bgi, le facteur de volume initial de formation de gaz en ft3 / SCF, est calculé à partir
de:

7
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Le facteur initial de déviation de gaz (également appelé compressibilité), Zi, est calculé
à la pression initiale pi du réservoir de gaz. Ce facteur explique la différence entre les
volumes de gaz réel et idéal.
La température T du réservoir est en degré Rankin (R °).

EXEMPLE

Un réservoir de gaz volumétrique présente les caractéristiques suivantes:

A = 1,320 acres T = 200°F

Pi = 4,000 psia h = 45ft
Zi = 0.916 ɸ = 0.175
Siw = 0.23

SOLUTION

Le facteur de volume initial de formation de gaz est:

Le gaz initial en place est:

LA PERMÉABILITÉ

En plus d'être poreux, une roche réservoir doit avoir la capacité de permettre aux fluides
pétroliers de s'écouler à travers ses pores interconnectés. La capacité de la roche à
conduire des fluides est appelée perméabilité. Cela indique que les roches non poreuses
n'ont pas de perméabilité. La perméabilité d'une roche dépend de sa porosité effective,
par conséquent, elle est affectée par la taille du grain de la roche, la forme du grain, la
distribution granulométrique, l'empaquetage du grain et le degré de consolidation et de
cimentation.
Le type d'argile ou de matériau de ciment entre les grains de sable affecte également la
perméabilité, surtout là où l'eau fraîche est présente. Certaines argiles, en particulier les
smectiques (bentonites) et les montmorillonitiques, gonflent dans de l'eau douce et ont
tendance à bloquer partiellement ou complètement les espaces poreux.

8
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure [Link]éabilités d’un fluide à travers un matériau poreux

Figure5. perméabilité inter-particulaire et de fracture

L'ingénieur français Henry Darcy a développé une équation de flux de fluide qui est
devenue depuis un des outils mathématiques standard de l'ingénieur pétrolier.
Cette équation est exprimée sous forme différentielle comme suit

Où :

u = fluid velocity, cm/s.

q = flow rate cm3/s.
k = permeability of the porous rock, Darcy (0.986923 µm2).
Ac = cross-sectional area of the rock, cm2.
µ = viscosity of the fluid, centipoises (cP).
l = length of the rock sample, cm.
𝑑𝑝
= pressure gradient in the direction of the flow, atm/cm.
𝑑𝑙

9
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure6. Représentation schématique de la loi de Darcy

Un Darcy est relativement une haute perméabilité. La perméabilité de la plupart des

roches réservoir pétrolier est inférieure à un Darcy. Ainsi, une plus petite unité de
perméabilité, le millidarcy (mD), est largement utilisée dans l'industrie pétrolière et
gazière. En unités SI, le micromètre carré (um2) est utilisé à la place du m2.

La perméabilité, k, dans l'équation 1.8 est appelée perméabilité "absolue" si la roche est
saturée à 100% avec un seul fluide (ou phase), tel que de l'huile, du gaz ou de l'eau. En
présence de plus d'un fluide, la perméabilité est appelée perméabilité «effective» (ko, kg,
ou kw étant respectivement une perméabilité à l'huile, au gaz ou à l'eau). Les fluides du
réservoir interagissent les uns avec les autres pendant leur mouvement à travers les
canaux poreux de la roche, par conséquent la somme des perméabilités efficaces de
toutes les phases sera toujours inférieure à la perméabilité absolue.
En présence de plus d'un fluide dans la roche, le rapport de la perméabilité effective de
toute phase à la perméabilité absolue de la roche est connu comme la perméabilité
"relative" de cette phase. Par exemple, la perméabilité relative de l'huile, du gaz et de
l'eau serait kro = ko / k, krg = kg / k, krw = kw / k, respectivement.

CLASSIFICATION DE LA PERMÉABILITÉ

Les réservoirs pétroliers peuvent avoir une perméabilité primaire, qui est également
connue sous le nom de perméabilité de la matrice, et de perméabilité secondaire.
La perméabilité de la matrice est née au moment du dépôt et de la lithification
(durcissement) des roches sédimentaires. La perméabilité secondaire résulte de
l'altération de la matrice de roches par compactage, cémentation, fracturation, et la
solution.

10
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure7. Effets de cimentation et compaction sur la perméabilité

Alors que le compactage et la cimentation réduisent généralement la perméabilité,

comme le montre la figure 8, la fracturation et la solution tendent à l'augmenter. Dans
certaines roches réservoir, en particulier les carbonates à faible porosité, la perméabilité
secondaire fournit le conduit d'écoulement principal pour la migration de fluide.

Porosité = 36% Porosité = 36%

KH = 1 000 mD KH = 100 mD
KV = 600 mD KV = 25 mD

Figure8. Grains de sable sans Figure9. Grains de sable sans

argile Matériaux de argile Matériaux de
cimentation cimentation
Ces deux figures montrent les effets du matériau de cimentation de l'argile sur la
porosité et la perméabilité

11
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

FACTEURS AFFECTANT LA MAGNITUDE DE PERMÉABILITÉ

La perméabilité des roches de réservoir pétrolier peut varier de 0,1 à 1,000 millidarcies
ou plus, comme le montre le tableau 1.9. La qualité d'un réservoir telle que déterminée
par la perméabilité, en mD, peut être jugée comme suit: pauvre si k <1, juste si 1 < k
<10, modérée si 10 < k < 50, bonne si 50 < k < 250 et très bonne si k> 250 mD. Dans les
champs de l'est du Texas, la perméabilité peut atteindre 4,600 mD. Les réservoirs ayant
une perméabilité inférieure à 1 mD sont considérés comme «imperméables». De telles
valeurs de faible perméabilité se retrouvent généralement dans les matrices calcaires et
dans les réservoirs de gaz à faible perméabilité de l'ouest des États-Unis. Les techniques
de stimulation telles que la fracturation hydraulique et l'acidification augmentent la
perméabilité de ces roches et permettent l'exploitation de tels réservoirs à faible
perméabilité, considérée comme non économique, il y a seulement 50 ans, les roches
ayant une perméabilité de 50 mD ou moins étaient considérées comme serrées.

Tableau 2. Perméabilité et Porosité de quelques réservoirs d’huile aux USA

Les facteurs affectant l'amplitude de la perméabilité dans les sédiments sont:

(a) Forme et taille des grains de sable: Si la roche est composée de gros grains plats
uniformément disposés avec la plus grande dimension horizontale, comme illustré à la
figure 10, sa perméabilité horizontale (Kh) sera très élevée alors que la perméabilité
verticale (Kv) sera moyenne à grande. Si la roche est composée principalement de gros
grains arrondis, sa perméabilité sera considérablement élevée et de même magnitude
dans les deux directions,
Comme le montre la figure 11. La perméabilité des roches réservoir est généralement
plus faible, surtout dans le sens vertical, si les grains de sable sont petits et de forme
irrégulière la figure 12.

12
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure 10. effets des gros grains plats sur la Figure 11. effets de gros grains arrondis sur la
perméabilité perméabilité

Figure 12. effets des petits grains irréguliers sur la perméabilité

La plupart des réservoirs de pétrole entrent dans cette catégorie. Les réservoirs à
perméabilité directionnelle sont appelés anisotropes. L'anisotropie affecte grandement
les caractéristiques d'écoulement du fluide de la roche. La différence de perméabilité
mesurée parallèlement et verticalement au plan de stratifcation est la conséquence de
l'origine du sédiment, parce que les grains se déposent dans l'eau avec leurs côtés les
plus longs et les plus plats dans une position horizontale. Le compactage subséquent du
sédiment augmente l'ordre des grains de sable de sorte qu'ils se trouvent généralement
dans la même direction

(b) Laminage: Les minéraux lamellaires tels que la muscovite et les laminés de schiste,
agissent comme des barrières à la perméabilité verticale. Dans ce cas, le rapport KH / Kv
est généralement compris entre 1,5 et 3 et peut dépasser 10 pour certaines roches
réservoir. Parfois, cependant, Kv est plus élevé que KH en raison de fractures ou de
jonctions verticales et de canaux de solutions verticales. Les joints ne font obstacle à la
perméabilité horizontale que s'ils sont remplis d'argile ou d'autres minéraux.

13
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

L'importance des minéraux argileux comme déterminant de la perméabilité est souvent

liée non seulement à leur abondance, mais aussi à leur minéralogie et à la composition
des fluides poreux. Si les minéraux argileux qui recouvrent les surfaces des grains se
dilatent et / ou se détachent en raison des changements de la chimie des fluides poreux
ou de l'invasion du filtrat de boue, la perméabilité sera considérablement réduite.

(c) Cimentation: La figure 8 montre que la perméabilité et la porosité des roches

sédimentaires sont influencées par l'étendue de la cimentation et l'emplacement du
matériau de cimentation dans l'espace poreux.

(d) Fractures et solutions: Dans les roches de grès, la fracturation n'est pas une cause
importante de la perméabilité secondaire, sauf lorsque les grès sont intercalés avec des
schistes, des calcaires et des dolomies. Dans les carbonates, la solution de minéraux par
percolation des eaux acides de surfaces et souterraines lorsqu'elles passent le long des
pores primaires, les fissures, les fractures et les plans de stratification, augmentent la
perméabilité de la roche du réservoir. Les perméabilités horizontale et verticale sont
égales dans de nombreux réservoirs de carbonate.

CORRELATION DE KOZENY
Kozeny a dérivé une des corrélations les plus fondamentales et les plus populaires
exprimant la perméabilité en fonction de la porosité et de la surface spécifique.
Considérons un échantillon de roche poreuse de section transversale A et de longueur L
comme étant constitué d'un nombre n de tubes capillaires droits en parallèle, les espaces
entre les tubes étant scellés par un matériau de cimentation. Si les tubes capillaires sont
tous du même rayon r (cm) et longueur L (cm), le débit q (cm3/s) à travers ce faisceau de
tubes, selon l'équation de Poiseuille, est:

Où la perte de pression Δp sur la longueur L est exprimée en dynes / cm2.

14
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure 13. Relation typique perméabilité-porosité pour diverses roches

La loi de Darcy peut aussi approximer le flux de fluides à travers ces n Capillaires:

10

Où AC est la surface totale de la section transversale, y compris les zones cimentées, de ce

faisceau de tubes capillaires.

Egalisons les équations 9 et 10 et résolvons k donnes:

11

Par définition, la porosité est :

12

En remplaçant AC = nΠr2 / Φ de l'équation 12 dans l'équation 11, on obtient une relation plus
simple entre perméabilité et porosité pour les pores de même taille et de rayons égaux à r:

13

Où k est en cm2 (1cm2 = 1,013 x 108 darcys) ou en µm2 (1 mD = 9,871 x 10-4 µm2) et Φ est une
fraction.

Soit SvP la surface interne par unité de volume de pore, où la surface As pour n tubes capillaires
est n (2ΠrL) et le volume poreux Vp est n (Πr2L)

14

Soit Svgr la surface spécifique d'un matériau poreux ou la surface totale exposée dans l'espace
poreux par unité de volume de grain. Pour un faisceau de tubes capillaires, la surface totale
exposée, At, est équivalente à la surface interne As; Et le volume du grain, Vgr. Est égal à :
Ac L (1-Φ). Ainsi

15
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

15

La combinaison des équations 14 et 15 donne:

16

L'équation 15 peut être écrite comme:

17

En remplaçant SvP par l'équation 17, on obtient:

18

Après détermination de la surface spécifique par unité de volume de pore, SvP, à partir de
données capillaires ou d'analyse d'image pétrographique (PIA), l'équation 16 est utilisée pour
obtenir SVgr.

EXEMPLE

Un échantillon carotté d'une formation de grès uniforme a une perméabilité de 480 mD et une
porosité de 0,17. Estimer:
(a) le rayon moyen de la gorge des pores de la carotte;
(b) surfaces spécifiques SvP et Svgr

SOLUTION

(a) En supposant que les canaux d'écoulement dans l'échantillon de la carotte peuvent être
représentés par un faisceau de tubes capillaires, le rayon de gorge de pore peut être estimé à
partir de l'équation 13. Tout d'abord, la perméabilité est convertie de mD en µm2:

La résolution de l'équation 13 pour r donne:

16
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

(b) La surface spécifique par unité de volume de pore est donnée par l'équation 14:

La surface spécifique par unité de volume de grain peut être estimée à l'aide de l'équation 16:

Toutes les équations ci-dessus utilisées pour établir la relation entre la perméabilité et la porosité
(équation 18) sont basées sur l'hypothèse que la roche poreuse peut être représentée par un faisceau
de tubes capillaires rectilignes. Cependant, la longueur de trajet moyenne qu'une particule fluide
doit parcourir est en fait supérieure à la longueur L de l'échantillon de carotte. Le départ d'un
milieu poreux étant constitué par un faisceau de tubes capillaires droits peut être mesuré par le
coefficient de tortuosité ɩ qui est exprimé comme:

19

Où La est le chemin d'écoulement réel et L la longueur de la carotte. L'équation 19 est préférée ici
parce que dans la plupart des expériences de laboratoire, le produit du facteur de résistivité de
formation (F) et de la porosité est lié au rapport La / L par la corrélation suivante:

17
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

20

L'exposant C est la constante de corrélation, qui va de 1,7 à 2. On notera que cette plage est assez
proche de celle du facteur de cimentation m.

Ainsi, pour un faisceau de tubes capillaires tortueux, la loi de Poiseuille devient:

21

En combinant l'équation 21 avec l'équation 10 et en utilisant la même approche que ci-dessus, on

peut montrer que les équations 13, 17 et 18 deviennent respectivement:

22

23

24

Wyllie et Spangler ont suggéré que le facteur 2 soit remplacé par un paramètre plus général, à
savoir le facteur de forme des pores Kps. Carman a indiqué que le produit Kps ɩ peut être
approximé de 5 pour la plupart des matériaux poreux. L'équation 24 pour les roches poreuses peut
alors être écrite comme suit:

25

L'équation 25 est la forme la plus populaire de l'équation de Kozeny, même si, dans la roche
poreuse actuelle, Kpsɩ est variable et beaucoup plus grande que 5.

EXEMPLE

Un sac de sable de grains fins uniformes à une porosité effective de 0,[Link] diamètre moyen des
grains est d'environ 1/8 mm. Calculez la perméabilité de ce sac de sable non consolidé.

18
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

SOLUTION

La surface spécifique des grains peut être estimée, en supposant que les grains sont sphériques,
comme suit:

Où dgr est en cm. 1/8 mm = 0,0125 cm.

Maintenant, en utilisant l'équation de Kozeny, 24, la perméabilité peut être estimée:

L'équation de Carman-Kozeny peut également être utilisée pour estimer la perméabilité:

Le changement de la constante de 2 à 5 donne une variation de 40% de la valeur de k.

LE CONCEPT DES UNITES DE DEBIT

THÉORIE MATHÉMATIQUE DES UNITÉS DE DÉBIT

Selon Tiab, une unité d'écoulement hydraulique est un corps continu à travers un volume de
réservoir spécifique qui possède pratiquement des propriétés pétrophysiques et fluides cohérentes,
qui caractérisent uniquement sa communication statique et dynamique avec le puits de forage.
Tiab, Tiab et al., Et Amaefule et al. Ont développé une technique pour identifier et caractériser une

19
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

formation ayant des caractéristiques hydrauliques similaires, ou des unités d'écoulement, basées
sur les mesures microscopiques d'échantillons de roche.
Cette technique est basée sur une équation de Kozeny-Carman modifiée et sur le concept de rayon
hydraulique moyen. La forme générale de l'équation 25 est:

26

où:

K = perméabilité, µm2,
Φe = porosité effective,
Svgr = surface spécifique par unité de volume de grain,
ɩ = tortuosité du chemin d'écoulement, et
KT = Kpsɩ = facteur de zonage efficace.

L'équation 26 peut s'écrire:

27

Où ΦR est:

28

Le paramètre KT, appelé ici facteur de zonage efficace au niveau des pores, est fonction de la taille
et de la forme des pores, de la taille et de la forme des grains, de la distribution des pores et des
grains, de la tortuosité, de la cimentation et du type de système de pores, e.g. Intergranulaire,
intercristalline, vacuolaire ou fracturée.
Le paramètre KT pour une formation de grès homogène peut être estimé à partir de:

29

L'indice de lithologie J1 est déterminé à partir des données de pression capillaire.

Les données expérimentales montrent que le tracé de la fonction J de Leverett, J (S*W), par rapport
à la saturation d'eau normalisée S*W sur un graphique log-log donne une droite selon l'équation
suivante:

La lithologie est la nature des roches formant un objet, ensemble, ou couche géologique.

20
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

30

Où J1 est l'intersection de la droite (extrapolée) à S*W = 1, comme le montre la figure 14. La

saturation en eau normalisée est définie comme suit:

31

L'indice de répartition de la taille des pores λ est la pente de la ligne. L'indice de lithologie J 1 va de
0,44 pour un grain sphérique non consolidé à 0,20 pour une formation de grès propre consolidée
avec une distribution homogène des pores. Le tableau 3.2 présente les valeurs typiques de J 1 et KT
pour différentes formations. Des valeurs élevées de J 1 se retrouvent généralement dans des
réservoirs à perméabilité élevée, tandis que de faibles valeurs de J 1 correspondent à des réservoirs à
faible perméabilité.

Figure 14. La détermination de KT à partir d'un diagramme de J par rapport à S*W

21
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Tableau 3. Valeurs typiques de J1 et KT pour plusieurs formations

AIRE DE SURFACE SPÉCIFIQUE

La surface spécifique peut être estimée par au moins trois techniques:

Le procédé d'adsorption de gaz, l'analyse d'image pétrographique (PIA) et la résonance magnétique
nucléaire (RMN). La méthode de base pour mesurer la surface à partir de la technique d'adsorption
de gaz consiste à déterminer la quantité de gaz inerte, typiquement de l'azote, de l'argon ou du
krypton, nécessaire pour former une couche épaisse d'une molécule sur la surface d'un échantillon
à température cryogénique. La zone de l'échantillon est alors calculée en utilisant la zone connue,
d'autres considérations, devant être occupée par chaque molécule de gaz dans ces conditions. Le
procédé d'adsorption de gaz est largement utilisé dans la détermination de la surface spécifique des
matériaux poreux. Il doit toutefois être limité aux milieux poreux qui n'ont pas de grandes surfaces
spécifiques et où les grains de la matrice sont singulièrement lisses et réguliers, c'est-à-dire
sphéricité > 0,7 et rondeur > 0,5 comme le montre la figure 15.

Un diagramme log-log de Svgr, (cm -1) par rapport au diamètre moyen du grain dgr, (cm) a donné la
corrélation suivante:

32

Le numérateur est en fait le facteur de forme du grain Kgs, comme indiqué dans l'équation 33. Le
diamètre moyen des grains peut être obtenu à partir de plusieurs méthodes: analyse par tamisage,
PIA, compacteur et micromètre. Cette corrélation est applicable aux grains ayant une sphéricité ≥
10,7 et de rondeur ≥ 30,5. La forme générale de l'équation 32 est la suivante:

33

22
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure 15. Rondeur et sphéricité de formes différentes

La surface spécifique des pores peut être déterminée à partir de :

34

Où Lp et Ap sont respectivement le périmètre des pores et la section transversale des pores. En

utilisant la PIA, un facteur de forme de pores planaire fps peut être déterminé comme suit:

35

Le facteur fps, illustré au tableau 4, indique le périmètre Lp si la fonction planaire ou 2D était un

cercle. La plage de fps est 3,75 (sphéricité de 0,5 ou moins et rondeur de 0,3 ou moins) à 5,84
(sphéricité de 0,5 ou moins et rondeur supérieure à 0,5). Il n'y a pas de relation pratique entre le
facteur de forme de pores fps, 2D et le facteur de forme de pores Kps, 3D. Pour un pore sphérique
idéal fps = 1 et Kps = 6.

23
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Tableau 4. facteur de forme des pores 2D

La méthode de résonance magnétique nucléaire ou RMN semble être actuellement la technique la

plus précise pour estimer la surface spécifique. Dans ce cas, les surfaces spécifiques Svgr et Spv
sont obtenues à partir de:

36

37

où:
ANMR = surface spécifique RMN du matériau sec, m2 / g
Pm = densité grain-matrice, g / cm3
SVgr = surface spécifique par unité de volume de grain, m2 / cm3

24
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure 16. la saturation en eau irréductible en fonction de la zone de surface (SPV) dans la
formation de grès

Les valeurs de Spv et de Svgr obtenues à partir de la RMN sont généralement supérieures aux
valeurs obtenues par la PIA ou la technique d'adsorption de gaz. Plusieurs études ont montré que la
surface spécifique, mesurée avec l'une de ces trois méthodes, est liée à la saturation irréductible de
l'eau ou simplement à la saturation en eau par une relation de la forme générale:

38

Où a et b sont des constantes de corrélation. Zemanek a étudié les réservoirs de grés à faible
résistivité et a trouvé que les surfaces mesurées par la technique de RMN étaient quantitativement
en accord avec les saturations d'eau irréductible à partir des données de la courbe de pression
capillaire.
Les figures 15 et 16 démontrent une bonne corrélation entre Spv et Sw et Swi, obtenant
respectivement:

39

40

Où Sw, et Sw, sont exprimés en pour cent et Spv en m2 / cm3. Les valeurs de Spv utilisées dans
les figures 15 et 16 ont été obtenues à partir de l'équation 37.

25
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

APPLICATIONS DE LA LOI DE DARCY

Ce chapitre décrit les caractéristiques du flux de fluides à travers des matériaux géologiques
poreux. Les pores, ou conduits d'écoulement, sont des capillaires complexes et interconnectés et
des canaux de tailles variables. Le flux de fluides compressibles et incompressibles à travers des
roches poreuses est décrit par la loi de Darcy et ses dérivés. Le cas le plus simple d'écoulement de
fluide à travers un milieu poreux est l'écoulement linéaire d'un fluide monophasé sous un gradient
de pression constant, connu sous le nom d'écoulement linéaire en régime stationnaire. Lorsque
deux fluides sont présents dans un milieu poreux, l'écoulement en régime permanent se produit
sous un gradient de pression constant seulement lorsque les saturations en fluides restent
constantes. Si les saturations changent par rapport au temps (par exemple, si la saturation en eau
augmente alors que la saturation en huile diminue), le flux de fluides se caractérise par un
écoulement à l'état instable.

Pour exprimer la quantité de fluide qui s'écoule à travers un système de roche poreuse de
géométrie et de dimensions spécifiées, comme celle montrée à la figure 17. Il est nécessaire
d'intégrer la loi de Darcy sur les limites du système poreux. Cette loi, sous sa forme différentielle
simple, est:

41

où:
V = vitesse apparente d'écoulement du fluide, cm / s.
K = perméabilité de la roche poreuse.
P = viscosité du fluide en écoulement, centipoise.
𝑑𝑝
= gradient de pression dans la direction d'écoulement, atm / cm.
𝑑𝑥
X = distance dans le sens d'écoulement, toujours positive, cm.

Figure 17. Système d'écoulement linéaire typique dans un échantillon.

Cette relation empirique unidimensionnelle a été développée par l'ingénieur français Henry Darcy
en 1856 alors qu'il étudiait le flux d'eau à travers des filtres à sable pour la purification de l'eau.
(1) La loi de Darcy pourrait être étendue à des fluides autres que l'eau,
(2) la constante de proportionnalité est en fait le rapport de mobilité Up, et
(3) La loi de Darcy est indépendante de la direction du flux dans champ gravitationnel de la Terre.
Le gradient dp/dx est la force d'entraînement et peut être dû à des gradients de pression de fluide et
/ ou à des gradients hydrauliques (gravitationnels).

26
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Dans les réservoirs de pétrole avec un grand bouchon de gaz en expansion et des caractéristiques
importantes de drainage par gravité, cependant, les gradients gravitationnels sont importants et
doivent être pris en compte lors de l'analyse de la performance du réservoir.

FLUX LINEAIRE DE FLUIDES INCOMPRESSIBLES

Les hypothèses suivantes sont nécessaires au développement des équations de base décrivant le
flux linéaire de fluides incompressibles ou légèrement compressibles à travers des milieux poreux:
(1) il existe des conditions d'écoulement stationnaire;
(2) la roche poreuse est saturée à 100% avec le fluide s'écoulant; Cependant, une phase fixe et
immobile peut être présente et est souvent, comme c'est le cas pour l'écoulement d'huile dans une
roche poreuse contenant une saturation d'eau irréductible, ou dans le cas d'un écoulement d'huile
avec une phase gazeuse immobile de saturation inférieure à la saturation de gaz critique;
(3) la viscosité du fluide s'écoulant est constante;
(4) les conditions isothermes prévalent;
(5) la roche poreuse est homogène et isotrope;
(6) la porosité et la perméabilité sont constantes, c'est-à-dire indépendantes de la pression;
(7) l'écoulement est des effets de turbulence laminaire, c'est-à-dire négligeables; et
(8) les forces de gravité sont négligeables.
Avec ces restrictions à l'esprit, la vitesse apparente est :

42

Où q est le débit volumétrique (cm3 / sec) et A est la section transversale perpendiculaire à la

direction d'écoulement (cm2). La vitesse réelle (va) est déterminée en divisant la vitesse apparente
(v) par la porosité de la roche (Φ). Si une saturation d'eau irréductible, Siw, est présente, la vitesse
réelle dans un réservoir humide est:

43

La combinaison des équations 42 et 43 donne:

44

En séparant et intégrant les variables entre les limites 0 et L, et P1 et P2, on obtient l'expression
suivante pour le débit volumétrique:

45

L'équation 45 est l'équation de débit linéaire classique utilisée dans les calculs de débit de fluide.
Cette expression est écrite dans les unités fondamentales qui définissent l'unité de Darcy. La
transformation en unités de champs de pétrole couramment utilisées, telle que q = bbl/D, A = ft2,
p = psi, L = ft et k = mD donne:

27
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

46

Dans les unités SI, on obtient un débit de 1 m3 / s pour un fluide s'écoulant à travers un milieu
poreux avec une perméabilité de 1 µm2, une surface de section transversale de 1 m2 et une viscosité
du fluide de 1 Pa x s sous un gradient de pression de 1012 Pa / m.

EXEMPLE

Un échantillon de noyau cylindrique de 10 cm de long a été soumis à un essai d'écoulement

linéaire en laboratoire sous une pression différentielle de 3,4 atm en utilisant un fluide de viscosité
2,5 cp. Le diamètre du noyau est de 4 cm. On obtient un débit de 0,35 cm3 / s. Calculer la
perméabilité de cet échantillon de noyau.

SOLUTION

La figure 17 est une représentation schématique de l'échantillon de base. En utilisant l'équation 45,
la perméabilité de l'échantillon de noyau est:

Pour utiliser l'équation 46, il faut d'abord convertir les données en unités de champs pétrolifères.

La perméabilité de cet échantillon en utilisant l'équation 46 est:

28
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

En raison des nombreux systèmes unitaires utilisés par l'industrie, il est très important que les
ingénieurs pétroliers puissent convertir des unités d'un système à un autre.
Pour estimer la pression en un point quelconque d'un système d'écoulement linéaire, l'équation 44
est intégrée entre les limites de 0 et de x, et p1 et p, respectivement, ce qui donne:

47

D'après l'équation 45, on obtient:

48

En divisant l'équation 47 par l'équation 48 et en résolvant pour la pression variable, p, on obtient:

49

Cette équation indique que le comportement sous pression d'un système d'écoulement linéaire
pendant l'écoulement stationnaire est une droite en fonction de la distance.

FLUX LINEAIRE DE GAZ

Considérons le même système d'écoulement linéaire de la figure 17, sauf que le fluide s'écoulant
est maintenant du gaz naturel. Cependant, comme le gaz se dilate à mesure que la pression
diminue, le gradient de pression augmente vers l'extrémité aval et, par conséquent, le débit q n'est
pas constant mais est une fonction de p. En supposant que la loi de Boyle est valide (facteur de
déviation de gaz z = 1) et un débit massique constant, c'est-à-dire que pq est constante, on peut
écrire:

50

Où les indices indiquent le point de mesure, q le débit moyen et p la pression moyenne. La

combinaison de cette relation avec la loi de Darcy, c'est-à-dire l'équation 45, donne:

51

Où µg est la viscosité du gaz dans les unités CP. Séparer les variables et intégrant entre p1et p2, et 0
et L donne:

29
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

52

53

L'expression du débit moyen qui suit peut être obtenue en combinant les équations 50 et 53:

54a

Si l'on suppose que la pression moyenne p est égale à (p1 + p2) / 2, l'équation 54a se réduit à:

54b

54b équation est identique à l'équation 45, ce qui donne le débit volumétrique des fluides
incompressibles. Par conséquent, la loi pour l'écoulement linéaire de gaz idéal est la même que
pour un liquide, tant que le débit de gaz est exprimé en fonction de la pression arithmétique.
Pour inclure l'effet des variations du facteur d'écart de gaz, z, des conditions normales de pression,
psc et température, p, à la pression moyenne, et la température, T,

55

En combinant les équations 54 et 55, et en résolvant pour qsc:

56

Où qsc est en cm3 / sec. Dans la mesure où p = (p1 + p2) / 2, l'équation 56 devient:

57

La conversion des unités de Darcy en unités de champ pratique et en supposant que zsc = 1
à psc = 14,7 psia et Tsc = 60° F ou 520 °R donne:

58

30
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

où:
qsc = débit volumétrique aux conditions standard, SCF / D.
µg = viscosité des gaz, cP
K = perméabilité de la roche réservoir, mD.
A = section transversale, ft2.
T = température moyenne du réservoir de gaz,°R
z = facteur moyen d'écart de gaz à T et p sans dimension.
L = longueur du corps de sable, ft.
Δp2 = p12 – p22; Psia2.

Si le débit moyen est exprimé en ft3 / jour à la pression moyenne P et la température moyenne T et
d'autres variables sont exprimées en unités de champs pétrolifères, l'équation 54 devient:

59

L'équation suivante est utile pour déterminer le débit volumétrique de sortie q2 à la pression p2, qui
est généralement la pression atmosphérique dans une expérience de laboratoire:

60

Où q2 est en cm3 / s. Si des unités de champs pétroliers pratiques sont utilisées dans l'équation 60,
où q2 est exprimée en ft3 / jour:

61

EXEMPLE

Un tuyau horizontal de 2 po de diamètre intérieur et 12 po de long est rempli d'un sable de 24% de
porosité. Ce sable a une saturation d'eau irréductible de 28% et une perméabilité au gaz de 245
mD. La viscosité du gaz est de 0,015 cP.
(a)Quelle est la vitesse réelle du gaz (en cm / sec) sous une pression différentielle de 100 psi?
(b)Quel est le débit moyen du gaz en ft3 / D et cm3 / sec?

SOLUTION

(a) La vitesse réelle va du gaz s'écoulant peut être calculée à partir de l'équation 7.3 où ɸ = 0.24,
SWi = 0.28 et la vitesse apparente v peuvent être obtenus à partir de la loi de Darcy. Dans la mesure
où k = 0.245 Darcy, µg = 0,015 cP, L = 12 po x 2.54 cm / po. = 30.48 cm et Δp = (100 psi) / (14,7
psi / atm) = 6,80 atm, la vitesse apparente est égale à:

Et la vitesse réelle est:

31
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

(b) Le débit volumétrique moyen de gaz à travers ce sac en ft3 / D est obtenu à partir de l'équation
7.19, où k = 254 mD, L = 1 ft, et A = Π(1/12)2 = 0.0128 ft2:

En supposant un taux de débit constant, q, une équation de distribution de pression le long d'un
corps de sable linéaire peut également être dérivée en combinant les lois de Boyle et Darcy et en
intégrant p1 et p, et 0 et x. En remplaçant L par x dans l'équation 52 et intégrant:

62

D'après l'équation 53, on peut obtenir:

63

En divisant l'équation 62 par l'équation 63 et en résolvant la pression variable p, on obtient:

64

LA LOI DE DARCY ET DE POISEUILLE

La loi de Darcy pour le flux linéaire des fluides incompressibles dans les roches poreuses et
perméables et l'équation de Poiseuille pour le débit capillaire liquide sont assez semblables. La
forme générale de la loi de Poiseuille pour l'écoulement d’un liquide visqueux à travers des tubes
capillaires est:

65

où:
R = rayon du tube capillaire, cm.
L = longueur du tube capillaire, cm.
µ = viscosité du fluide en circulation, poise

32
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Δp = chute de pression, dynes / cm2 (= 1,0133 x l06 atm). Si les canaux conducteurs de fluide dans
un milieu poreux peuvent être représentés par un faisceau de tubes capillaires parallèles de
diamètres différents, alors le débit à travers ce système est:

66

où:
nj = nombre de tubes de rayon rj.
N = nombre de groupes de tubes de rayons différents.

Cette expression peut être réécrite comme suit:

67

Où C est le coefficient d'écoulement Il est évident que l'équation 66 est

semblable à l'équation 67 (loi de Darcy) où le coefficient C est équivalent à la perméabilité. Ainsi:

68

Où A est la superficie totale de la section transversale, comme illustré à la figure 18.

Figure 18. Le système d'écoulement de Poiseuille pour capillaires droites.

En remplaçant A = ΠR2 (R le rayon de la section transversale totale) dans l'équation 68, on obtient:

33
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

69

FLUX LINEAIRE A TRAVERS LES FRACTURES ET LES CANAUX

Les réservoirs d'huile avec fracture-porosité matricielle contiennent également des canaux de
solution. La matrice (porosité inter-granulaire) est habituellement de faible perméabilité et contient
la majeure partie de l'huile (96% -99%). Considérant que ces fractures et canaux de solution ne
peuvent contenir un volume important d'huile, généralement inférieur à 4% de l'huile totale dans
un réservoir, ils sont très importants pour la réalisation de taux de production économique.

Figure 19a. Roche carbonatée Figure 19b. Roche carbonatée

montrant une porosité: A, vacuoles; B, montrant une porosité provenant
canaux communs; C, canaux plans; D, de fracturation et fissuration
canal de solution

La porosité de fracture est commune dans de nombreuses roches sédimentaires et est formée par
une défaillance structurelle de la roche sous des charges causées par diverses formes de
diastrophisme, telles que le pliage et la faille.

FLUX TRAVERSANT LES FRACTURES

La signification des fractures comme porteurs de fluide peut être évaluée en considérant une seule
fracture s'étendant à une certaine distance dans le corps de la roche et s'ouvrant dans le puits de
forage, comme le montre la figure 20.
Rappelant l'équation d'hydrodynamique classique pour l'écoulement à travers des fentes
d'autorisations fines et de largeur unitaire telles que rapportées par Croft et Kotyakhov

70

où:
h = hauteur (ou épaisseur) de la fracture, cm.
wf = largeur de la fracture, cm.

34
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Figure 20. Modèle linéaire pour le débit de rupture.

L = longueur de la fracture, cm.

P = viscosité du fluide, Poise.
Δp = chute de pression (p1 - p2), dynes / cm2.

La vitesse réelle du fluide traversant la fracture est ainsi:

71

En supposant que la porosité de la fracture est l'unité et que la saturation d'eau connée (eau de
formation) dans la fracture est nulle, la vitesse réelle (selon la loi de Darcy où Δp est exprimé en
dynes / cm2, k dans Darcy, p dans Poise et L en cm) est:

72

En combinant les équations 71 et 72et en résolvant la perméabilité de la fracture kf (où wf = cm et

kf = Darcy):

73

Les fractures sont classées comme ouvertes (espace ouvert visible), fermées (non visible espace
ouvert dans la section mince), partiellement, ou complètement rempli.
De nombreux réservoirs carbonatés présentent des fractures avec un certain degré de remplissage,
qui peuvent être constitués de cristaux de calcite, de dolomie, de pyrite, de gypse, etc. L'équation
73 n'est donc valable que pour le cas où la fracture est totalement ouverte et propre de toute
particule de remplissage, c'est-à-dire, ɸf = 1.

35
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

La perméabilité réelle de la fracture peut être déterminée à partir de l'équation de la vitesse réelle
du fluide s'écoulant à travers la fracture:

74

Où ɸf est la porosité fractionnaire de la fracture et Swcf est la saturation en eau de formation dans la
fracture. Par définition:

75

La vitesse apparente de l'équation 74 est:

76

Où la vitesse réelle est exprimée par le Taux actuel du débit d'écoulement de fluide à travers la
rupture divisée par la surface de la fracture, ou:

77

Et le débit q est exprimé par l'équation 70. En remplaçant Va et q dans l'équation 76, on obtient:

78

Assimilant cette expression par la loi de Darcy (équation 72) et en résolvant pour la perméabilité
effective de la rupture (en Darcy), on obtient:

79

EXEMPLE

Un bloc cubique d'une roche carbonatée avec un système de porosité intercristallin-intergranulaire

a une porosité de matrice de 19%. La perméabilité de la matrice est de 1 mD. Calculer:
(a) la perméabilité de la fracture si chaque pied carré contient une fracture dans le sens de
l'écoulement du fluide, et
(b) le taux débit à travers la fracture et le système matrice- rupture.
La largeur de la fracture est de 2,5 x 103 po, la viscosité du fluide s'écoulant est de 1,5 cP et Δp à
travers ce bloc est de 10 psia.

SOLUTION

(a) La perméabilité d'une fracture est estimée à partir de l'équation 79, où :

36
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Il est évident à partir de cette valeur de perméabilité extrêmement élevée que les fractures
contribuent sensiblement à la récupération du pétrole à partir de formations serrées qui autrement
seraient non commerciales. Cette contribution est en fait encore plus élevée car un pied carré de
roche carbonatée est généralement susceptible de contenir plus d'une fracture.
(b) Le débit à travers la fracture peut seulement être estimé à partir de la Loi Darcy (équation 46),
où

Le débit à travers la matrice est également obtenu à partir de l'équation 46, où la perméabilité de la
matrice est 1mD et

Ainsi :

Le débit total à travers le bloc est:

L'importance de la fracture sur la productivité des réservoirs peut être mieux appréciée en termes
de pourcentage de contribution par rapport au débit total, qui pour ce cas est de :

0,533 / 0,54 = 98,6%.

Le volume d'huile contenu dans les fractures et la matrice est:

80

Où Vom est le volume d'huile contenue dans la matrice et Vof est le volume d'huile contenu dans les
fractures, qui peut être estimé à partir de l'équation suivante:

81

82

où:

37
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Vo = huile-sur-place, m3.
As = surface de production de la formation, m2
ɸm = porosité fractionnaire de la matrice seulement.
h = épaisseur moyenne de formation, m.
ɸf = porosité fractionnaire des fractures seulement.
Swm = saturation en eau dans la matrice.
Swf = saturation en eau dans les fractures.

Le volume récupérable d'huile est:

83

où:
Em = facteur de récupération de la matrice, fraction.
Ef = facteur de récupération des fractures, fraction.

Si la perméabilité de la matrice est négligeable, c'est-à-dire inférieure à 0,1 mD,

La perméabilité moyenne du système d'écoulement matrice-fractura est obtenu de:

84

où:
Kf = perméabilité à la rupture
Km = perméabilité de la matrice
A = surface totale de la section transversale
nf = nombre de fractures par unité de surface
wf = largeur de fracture,
h = hauteur de rupture.

La perméabilité moyenne du réservoir de carbonate de l'exemple ci-dessus peut être estimée à

partir de l'équation 7.45:
où:

Dans la mesure où:

Par conséquent:

38
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

ECOULEMENT A TRAVERS LES CANAUX DE SOLUTION

Craft et Hawkins et Aguilera ont combinés la loi de Poiseuille pour l'écoulement visqueux des
liquides à travers des tubes capillaires avec la loi de Darcy pour l'état d'équilibre des écoulements
linéaires des fluides incompressibles pour estimer la perméabilité des canaux de solution (Figure
21). Le débit volumétrique actuel du fluide s'écoulant à travers un tube capillaire de rayon r, et la
longueur L est:

Figure 21. Système de canal-matrice.

85

A partir de la loi de Darcy, en supposant que la porosité du canal, ɸC, et la saturation irréductible
de l'eau, Siwc, sont égales à l'unité et à zéro, respectivement, le débit est égal à:

86

En comparant les équations 85 et 86, et en résolvant la perméabilité du rendement des canaux de

solution:
87

Où:

kC = perméabilité du canal de solution, Darcy. rC = rayon du canal tubulaire, cm.

Le développement de la porosité dans certains réservoirs carbonatés est dû au lessivage des roches
carbonatées par la minéralisation des eaux. Les précipités de cette eau circulante peuvent remplir
les pores et les canaux existants avec une variété de fines particules (sel, anhydrite et gypse),

39
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

rendant la porosité du canal de solution inférieure à l'unité. De plus, la saturation en eau de ces
canaux, qui se sont formé en raison de la circulation de l'eau n’est pas nulle. Ainsi, la surface réelle
ouverte à l'écoulement est:

88

Et l'équation 87 devient:

89

où:
ɸC = Porosité du canal de solution.
Siwc = saturation irréductible de l'eau dans le canal.

La perméabilité moyenne d'un système d'écoulement de canal-matrice peut être calculée à partir de
l'équation:

90

où:
KC = perméabilité des canaux, Darcy,
Km = perméabilité de la matrice, Darcy,
A = surface de section transversale, cm2,
n, = nombre de canaux par unité de surface, et
rC, = rayon du canal de solution, cm.

EXEMPLE

Un échantillon cubique d'une formation calcaire a une perméabilité matricielle de 1 mD et contient

5 canaux-solution par ft2. Le rayon de chaque canal est de 0,05 cm. Calculer:
(a) la perméabilité canal-solution prenant pour hypothèse une porosité-cavité de 3% et une
saturation irréductible de l'eau dans ces canaux égale à 18%; et
(b) la perméabilité moyenne de cette roche.

SOLUTION

(a) La perméabilité du canal-solution peut être obtenue à partir de équation 89:

En utilisant l'équation 87, c'est-à-dire en supposant que ɸC = 1 et Siwc = 0, la perméabilité d'un

canal est de 31,500 Darcy, ce qui est plus de 40 fois la valeur de kC, obtenue à partir de l'équation
7.50 et donc irréaliste.
(B) La perméabilité moyenne de ce bloc contenant 5 canaux est estimée à partir de l'équation 90,
où :

40
 Cours : Ecoulement des fluides en milieux poreux et rhéologie – Hamdi Abdelkader

Cet exemple illustre l'importance d'estimer la saturation irréductible de l'eau et la porosité des
canaux de solution et des fractures. Ces paramètres jouent un rôle important dans la détermination
du pétrole en place dans les pores et les fractures des vaisseaux et leur ignorance peut entraîner une
surestimation de la capacité de production des puits des réservoirs de carbonate.

41