COMMISSION DES COMMUNAUTES EUROPEENNES

énergie

Etude des phénomènes de transfert

de chaleur et de masse qui accompagnent
la réinjection dans un gisement
àhauteenthalpie

1979 EUR 6627 FR

 COMMISSION DES COMMUNAUTES EUROPEENNES

énergie

Etude des phénomènes de transfert

de chaleur et de masse qui accompagnent
la réinjection dans un gisement
à haute enthalpie

A. Bories et A. Podaire

Rapport établi par le

Groupe d'études sur les
milieux poreux, Toulouse (France)
Contrat n° 249-77

1979 EUR 6627 FR

 NOTATIONS

i = l(liquide), v(vapeur),
ε porosité,
p masse volumique de la phase i,
i
saturation en phase i,
vecteur vitesse de la phase i,
terme source de phase i,
débit de production en masse de phase i,
Κ perméabilité intrinsèque,
Κ perméabilité relative de la phase i,
ri
viscosité cinématique de la phase i,
pression,
g vecteur accélération de la pesanteur
P. masse volumique de la roche,
enthalpie de la roche,
h. enthalpie de la phase i (valeur prise sur la courbe de change­
ι ment de phase),
: vecteur flux de chaleur (dû à la diffusion thermique),
+ ρ S : densité du mélange diphasique,
g g

P1 S h + ρ S h
"—*—*■ : enthalpie du mélange diphasique,
p
l Sl + p
g Sg

Ρ : pression capillaire,
c
λ* : conductivité thermique,
Τ : température,
C : capacité calorifique de la roche,
pr
ι : énergie interne,
υ : température de saturation,
: pression de saturation,
C : capacité calorifique massique de la phase i,
¿£ : chaleur latente de vaporisation,
Ρ : pression,
Χ,Υ,Ζ : dimensions spatiales,
Y : titre massique de vapeur,
m : débit masse en liquide du fluide,
m : débit masse en vapeur du fluide,
m : débit massique total du fluide,
t : dimension temps,
o : coefficient d'expansion volumique de la vapeur,
φ : flux de chaleur provenant de 1'extérieur,
μ. : viscosité dynamique de la phase i.
INTRODUCTION

D'une manière g é n é r a l e , l ' é n e r g i e s t o c k é e d a n s un g i s e m e n t géother-

mqiue p e u t ê t r e d é c o m p o s é e en d e u x p a r t i e s :

- l ' é n e r g i e contenue dans l e f l u i d e q u i c o n s t i t u e l a fraction

d ' é n e r g i e r é c u p é r a b l e l o r s de l ' e x p l o i t a t i o n ,

- l ' é n e r g i e s t o c k é e dans l a m a t r i c e s o l i d e e t dont s e u l e une f a i -

b l e p a r t i e p e u t ê t r e r é c u p é r é e à l a s u i t e des t r a n s f e r t s de c h a l e u r e t
de masse q u i se m a n i f e s t e n t au cours du p r o c e s s u s de changement de phase.

Lorsque l e s gisements s o n t de type à vapeur dominante, ce d e r n i e r

type d ' é n e r g i e n ' e s t p r a t i q u e m e n t pas r é c u p é r é au cours de l a p r o d u c -
t i o n , compte t e n u notamment des f a i b l e s v a l e u r s de l a c h a l e u r volumique
de l a vapeur d ' e a u .

Sans e n t r e r dans l e s d é t a i l s d ' a n a l y s e q u i p e r m e t t r a i e n t d'affiner

c e s c o n c l u s i o n s , on c o n ç o i t donc q u ' i l p u i s s e ê t r e i n t é r e s s a n t d'amélio-
r e r , non seulement l e t a u x de r é c u p é r a t i o n , mais également l a l o n g é v i -
t é de l ' e x p l o i t a t i o n en p r é l e v a n t une p l u s grande q u a n t i t é de c h a l e u r
sur la roche.

C ' e s t dans ce b u t q u ' a é t é récemment proposée l a mise en oeuvre

de p r o c é d u r e s de r é i n j e c t i o n de f l u i d e r é c u p é r é e t condensé au s e i n
du gisement.

Si au p l a n des p r i n c i p e s , une t e l l e p r o c é d u r e p e u t sembler s a t i s -

f a i s a n t e , i l c o n v i e n t de s o u l i g n e r t o u t e f o i s q u ' i l n ' e x i s t e , à l ' h e u r e
a c t u e l l e , aucun r é s u l t a t p e r m e t t a n t de p r é c i s e r l e s c o n d i t i o n s de s a
mise en o e u v r e .

Le programme de r e c h e r c h e q u i a é t é e n t r e p r i s avec l e s o u t i e n de
l a Communauté Européenne, dont l a première phase a é t é développée au
c o u r s de l ' e x e r c i c e 1978/1979, s ' i n s c r i t dans l e cadre de ces p r é o c c u -
p a t i o n s . L ' é t u d e e f f e c t u é e concerne e s s e n t i e l l e m e n t :

- une é t u d e b i b l i o g r a p h i q u e s u r l e s mécanismes physiques de base

concernant l e s t r a n s f e r t s couplés de c h a l e u r e t de masse avec ou sans
changement de phase en m i l i e u poreux e t l e u r m o d é l i s a t i o n mathématique,
 « l a mise en évidence des paramètres c a r a c t é r i s t i q u e s dont dépen-
dent l e s phénomènes,
- l a d e s c r i p t i o n d'un p r o j e t d ' i n s t a l l a t i o n expérimentale.
MODELISATION MATHEMATIQUE DES TRANSFERTS DE CHALEUR ET DE MASSE
EN MILIEU POREUX

ETUDE PARAMETRIQUE

1.1. EQUATIONS DU PHENOMENE

Le modèle mathématique que nous envisageons de mettre en oeuvre pour

l'interprétation des expériences est inspiré de celui proposé par LUIKOV [_4J
pour l'étude des phénomènes de transfert de chaleur et de masse en milieu
poreux. Ce modèle ayant déjà fait l'objet d'une analyse approfondie [lj ,
nous ne présenterons ici que les résultats essentiels qui peuvent être adop­
tés, compte tenu de la nature particulière de nos expériences.

Dans la suite du texte, les grandeurs, afférentes à la vapeur, au

liquide et au solide, seront respectivement indicées i = l , 2 , 3 o u v , 1, r.

I.1. 1. Présentation globale des équations

I.1.1.1. Equation_de_continuité

Elle s'écrit en toute généralité pour un milieu poreux indéformable

ε *- (p. S.) = ­ V. (p. v\) + I.

dt

p. V. représente le flux convectif de i et sera explicité ultérieurement,

I. le terme source, exprimé en masse par unité de volume et de temps,
traduit la création de l'une ou l'autre phase, par suite des phéno­
mènes de vaporisation condensation.
 Σ I, = O I. ­ ­ I
i l v
i=l

Σ S,
i ­ 1 S,
1 + Sν ­ 1
i=l

1.1.1.2. Conservation de l'enthalgie

En admettant que l'équilibre local entre le fluide sous ces diffé­

rentes phases est réalisé, l'ensemble squelette solide plus fluide peut
être assimilé à un milieu continu fictif équivalent. Avec cette hypothèse,
qui est celle retenue en première approche de cette étude, l'équation de
conservation de 1'enthalpie s'écrit ;

2 2
— [d­ε) P r h + e Σ ρ S h ] = ­ V . [j + Σ h ρ V ] + —
q
3t i=l i=l dt

h. étant 1'enthalpie massique de la phase i,

h., p. , V. le flux de chaleur convectif volumique,

3 le flux de chaleur dû à la diffusion thermique.

Il est important de rappelée que cette formulation suppose négligea­

bles les effets de dissipation visqueuse.

Compte tenu des relations :

"2 2 3h.
ε Σ — (ρ s h.) ­ Σ e p. s — + h ε — (p. S.)
i­1 it i=l 3t 3t

et :

a +
e — (ρ S ) ­ ­ V . (ρ V ) + I
at
le bilan d'enthalpie peut encore s'écrire :

2 3h.
— ((1 - ε) ρ _h ) + Σ ρ ε S — i - V . [λ* 7 τ]
3t i=l at

2
-> dP
- Σ (p.V. . V h . + h . I.) + —
x x l x x
i=l dt

1.1.1.3. Equation_de_mouvement

Les débits massiques en phase vapeur et liquide seront explicités

au moyen de l'équation classique de Darcy en utilisant le concept de per­
méabilité relative K . :
rx
K K
rl r **■ ι
Ρ i V. = v. [V Ρ. - ρ ± g ]
1
les pressions de l'eau liquide et de la vapeur étant reliées par l'expres­
sion :

P, = Ρ ­ Ρ
1 ν cap

où Ρ est la pression capillaire,

cap

1.1.1.4. Equations d'état

Pour fermer complètement le système d'équations, il convient

d'y adjoindre les relations qui explicitent les variations des propriétés
des fluides et de la matrice poreuse quant à son aptitude à permettre
1'écoulement.

. Variation de la masse volumique des fluides : p. = f(P.,T)

. Comportement de la vapeur,
 . Lois d'évolution des viscosités : u.(P,T)

. Lois d'évolution des perméabilités relatives : K . (S,)

ri 1

1.1.1.5. Remarques_­_HYpothèses_simglificatrices

Les travaux effectués récemment dans le Groupe sur les phénomènes

de thermomigration, ont permis de constater que, d'une manière générale,
les transferts de type diffusionnel associés à l'existence des effets
capillaires, sont, d'une part, relativement faibles, d'autre part, ne se
manifestent que très lentement. Si l'on prend pour temps de référence la
durée moyenne des expériences effectuées par CADY [2] deux à trois heu­
res, il est donc vraisemblable que pour des conditions de fonctionnement
d'aussi courte durée, les effets capillaires ne joueront qu'un rôle se­
condaire par rapport aux mécanismes d'écoulement dus à l'expansion de la
phase vapeur. Dans les conditions de nos expériences, les termes qui leur
sont associés pourront être négligés en première approximation. Les flux
massiques seront alors explicités par les relations suivantes :

Κ Κ ,
­*■ rl
Ρλ V = — VP
ν
ι
κκ
·*■ rv
ρ ν = — VP
ν ν
V
ν
l ' é c o u l e m e n t du f l u i d e s ' e f f e c t u a n t dans l e p l a n horizontal.

1.2. MOD ELE MATHEMATIQUE

1.2.1. Equations

Compte tenu des hypothèses choisies, le modèle mathématique que

nous proposons d'utiliser, sera le suivant :
 2
3 ­
continuité :ε — (p. S.) = ­ V . ( p . V.) + I . avec Σ I . = 0 I, = ­ I
.. i l i i ι ­ ι 1 ι ν

3 t . 1=1
Κ ΚJ
·> ri
mouvement : ρ . V. = ­ VP
ι ι
V
i

2 3h
3 ­
énergie : (l­ε) p h + Σ ρ. ε S. (—­) = V . (λ V Τ)
r r X X
3t i=l 3t

2 2

­ Σ (ρ .V.) . V h . ­ Σ h. I + —
• 1 1 1 1 · 1 1 i J t
x=l 1=1 dt

­ état : p. = ρ.(Ρ,Τ) μ. = y.(Ρ,Τ) Κ , = Κ .(S.)

ι ι ι ι ri η 1

que l ' o n peut é c r i r e , en u t i l i s a n t une combinaison de l ' é q u a t i o n de con­

t i n u i t é avec l ' é q u a t i o n du mouvement e t l e s expressions de :

h = C Τ
r r

ÜjL­eS. Í Ü + $ vh, = c , 2L+SÌ ­ α τ > SL

dt ι Λ3t. i ι ι dt ρ. dt

et Σ h. I . = (h, ­ h ) I ,
. . ι ι 1 v i
ι=1

KK KK
ε — [ p L s + (i ­ s x ) p v J = V. [ ( — £ ! · + —ΞΣ) ν p j
3t ν, ν
1 ν
KK C KK c
* 3T „ ,,*„„,., rl 1 rv v, „„ „„ <J? ^ . dP
(pO Λ — = ν . (λΛ VT) + (—£=­=· + — 2 L i j vp . VT + ^C ï
3t ν, ν dt
 10

" l'expansion volumique du liquide est supposée très faible par rapport
à celle de la phase vapeur pour laquelle : a y_ —
~ Τ

­ (hv ­ hL) = X„ chaleur latente de vaporisation.

(PC)* « (pC)r (1 ­ ε) + ε [(ρ0)χ S;L + (pC)v (1 ­ S ^ ]

l'équation d'évolution de la pression étant exprimée par la relation carac­

téristique Ρ (Τ) de la courbe de saturation.

1.2.2. Grandeurs adimensionnelles

En posant :

* Je le
Χ , Y , Ζ = X/L, Y/L, Z/L, L longueur de référence

t* = t/t , P* = P/P , T* = Τ/Τ , (λ*)* = λ*/λ*

o o o o

(pC) /(pC)o = (pC ) , pL = P1/P0 , Pv = Pv/P0 '

v, = v , / v , ν = ν /v
1 l o v ν o

avec Ρ et Τ : conditions thermodynamiques initiales,

it it
λ , (pC) : conductivité thermique e t chaleur volumique du milieu poreux
saturé, à Ρ , Τ ,
o o
ν : viscosité cinématique de l'eau à Ρ , Τ ,
ρ : masse volumique de l'eau à Ρ , Τ ,
o o o
λ

* o
a
«= 7
q
<p<
 11

et, en introduisant les grandeurs réduites dans les équations précédentes,

on obtient (I étant homogène à une variation de masse par unité de temps)

e ^_ [(p* S,, + „ - β]? pj]- ^ V * . [,itt -& ,* p*]

3t ρ L v V, V
o o l v

( pc*)* ¿4 =^-l2 v
* · « λ *>* v* T *> + KP
° k 2 ί^Γ1 + " ^
3t L v (pC) L ν V
ο ο ν ν
3
7 * ρ * . y* τ * +
a?
'ο*- * + ο__ dP_
Τ (pC) (pC) Τ dt
ο ο ο ο

m e t t a n t a i n s i en é v i d e n c e l e s g r a n d e u r s a d i m e n s i o n n e l l e s :

P t K a* t K P t C K P t C
e o o g o o o i o o v
' 2 ' 2 ' * 2 ' * 2
ρ L V L ν (pC) L v (pC) L''
O O o o o o

Ρ *χ ρ
*o *^­ o
T (pC)* (PC)*T
o o o o

La constante de temps t étant indéterminée, si l'on pose :

L2
a
q

i l vient :

P K
o
KP
o o
P C ,
i
KP
o o
p C
v
L ρ
o
•¿
¿ζ, Ρ
o
ρ ν a* ρ ν a* (pC)* ρ v a* (pC)* T (pC)* (pC)*T
o o q o o q o o o q o o o o o

avec Κ Ρ /ρ ν = a cdiffusivité massique en phase liquide,

o o o m
 12

On peut encore é c r i r e :

Ρ
o
, a/a* , to C ) / ( ρ Ο * , to C ) / ( p C ) * , L 3 ρ^ Χ / Τ (pC)* ,
m q o l o o v o o o o
(PC)* τ
o o
d'où l e s grandeurs adimensionnelles :

e . V a' . K C l)/(P C, : . <po c.,,«*,* . ¿ ρ.*/*.**« - V

mq o l o o v o o o o ,pCj* T
o o

Le fluide saturant étant donné, transfert de chaleur et de masse dépen­

dront des paramètres suivants :

­ porosité e ,
­ nature du milieu poreux,
2
­ granulometrie (a étant proportionnel ã d diamètre des grains),
­ saturation initiale,
­ conditions thermodynamiques initiales,
­ conditions d'exploitation.
 13

ETUDE CRITIQUE DE DIFFERENTES APPROCHES

DE RESOLUTION NUMERIQUE

Ce bref exposé critique a pour but de dégager les caractéris-

tiques communes et les originalités particulières des différents modè-
les utilisés pour décrire les transferts de chaleur et de masse en milieu
poreux, puis d'en souligner les lacunes éventuelles afin de préciser le
cadre de notre étude.

Pour plus de complétude, nous présenterons les équations régis-

sant les phénomènes physiques, et leurs variantes, en insistant sur les
hypothèses s implificatrices et les relations conduisant à la formulation
définitive.

Nous soulignerons en conclusion les éléments des travaux anté-

rieurs qui peuvent, dès à présent, être utilisés dans le contexte de notre
étude.
 14

I. PRESENTATION DES EQUATIONS

1.1. Préliminaire

Nous traiterons le cas d'un réservoir soumis, en certains points

du domaine, à des processus d'injection ou/et de production de fluide.

Les puits de production/injection sont localisés en des points

ou lieux géométriques discrets du réservoir. Les valeurs des termes cor-
respondants de débits massique et thermique ne sont pas réparties avec une
densité continue, ou constante, à l'intérieur de ce réservoir, mais impo-
sées en certains points ou lieux géométriques de celui-ci.

En mécanique des milieux continus, les équations de continuité

et de conservation de l'énergie sont établies par bilans respectifs des
transferts de masse et de chaleur sur un volume élémentaire, de dimensions
microscopiques...

A tout volume élémentaire, considéré comme point physique ou maté-

riel ducontinu, correspond bijectivement un point mathématique, limite de
ce volume [l]

Supposons un puits de production ponctuel.

Le volume élémentaire sur lequel on procède au bilan massique ou

thermique, et le point physique associé au point mathématique où est loca-
lisé le puits, ont, d'après ce qui précède, des échelles identiques : la
prise en compte directe des termes de production s'avère alors impossible
en raisonnant sur des volumes microscopiques. Nous ne pourrons faire in-
tervenir ces termes, qu'une fois effectué un changement dimensionnel, nous
plaçant à une échelle macroscopique du continu, infiniment grande par rap-
port à la taille d'une particule.

Nous incluerons ensuite les quantités notifiant les divers prélè-

vements ou apports de fluide, dans la formulation de l'équilibre global au
niveau d'un bloc macroscopique.
 15

Les fonctions spatiales q., inhérentes à ces quantités, présen­

teront une ou plusieurs discontinuités (nullité dans tout l'espace sauf
aux lieux géométriques d'implantation des puits).

1.2. Equation de continuité

Pour chaque phase i du liquide, nous écrirons

— (ερ S ) = ­ V . (p V ) + I ­ q
3t

avec : Σ I = 0

1.3. Equation du mouvement

Nous supposerons l'équation classique de Darcy, applicable, soit

Κ Κ
Piv. EL [ y p - p g]
v
i

1.4. Equation de conservation de 1'enthalpie

En définissant un terme moyen de conduction/dispersion pour

toutes les phases, ce qui revient à assimiler l'ensemble squelette
solide + fluide à un milieu continu fictif équivalent, nous sommes con­
duits à :

— Γ (l­ε) p h + ε Σ p.S.h.1 = ­ V.LΓ J + Σ p.V.h.H Σ q. h. + —

.w L r r îii­l q . ι ι i­I _ ^i ι
dt i i. dt
 16

Cette formulation suppose négligeables les effets de dissipation

visqueuse.

1.5. Equation d'état

ρ (Ρ,Τ) i = liquide ou vapeur

En sommant les équations de continuité concernant chaque phase,

puis en substituant dans la nouvelle équation la relation donnée par la
loi de Darcy, nous obtenons :

Κ , Κ ρ,Κ , ρ Κ
3_ [ερ] . ν . [ κ (­LL+ ­L2.) yp ­ κ (­LÜ+ ­SJS, J ] ­ Σ q
at ν, ν ν, ν ι
l g l g

En remplaçant dans l'équation de l'énergie les vitesses par les

valeurs données par la loi de Darcy, nous écrirons :

2­ [ (ΐ-ε) p h + eph 1 = - 7.J + V. Γ κ (­LL­Ì. + ­ÏSL2, VP

*­ r r J α
r r q
at ν. νΛ
i g

ρ, Κ ,h, ρ κ h _ ,_

νΊ v i dt
1 g

Nous regrouperons sous les termes "pondérés" ρ et h les expres­

sions :

p = S
l pl + S
g pg

S. p, h, + S p h
h = J- ι ι g g g
S p + S rp
i i g g
 17

Hypothèses simplificatrices

a) Le milieu considéré ne comprenant qu'un seul composant fluide, la somme

des saturations de chacune des phases égale 1. Soit :

S, + S = 1
i g

b) Les pressions du liquide et de la vapeur sont liées par :

P = p ­ p
1 v cap

P étant la pression générée par les effets capillaires,

cap

Les études faites relativement à l'importance des effets capil­

laires montrent que ceux­ci sont faibles et ne se manifestent que très
lentement. En conséquence, nous négligerons les termes qui leur sont as­
sociés. D'où :

P, = P = P
1 v

ç_) Nous supposerons le mouvement du fluide dans le milieu poreux suffisam­

ment lent et les surfaces de contact des différentes phases suffisamment
grandes, pour que l'équilibre local entre phases puisse être considéré com­
me établi instantanément.

Cela permet d'introduire un flux moyen de chaleur, dû à la dif­

fusion thermique, sous la forme :

­ * ■ *

J = ­ λ VT
q

D'après SOMERTON, la conductivité thermique λ est une fonction

de la température, de la porosité, de la saturation en liquide. Seule
l'influence de la saturation, dont le rôle est largement prépondérant,
sera prise en considération.
 10

d) Nous supposerons l ' e n t h a l p i e de l a r o c h e f o n c t i o n de l a t e m p é r a t u r e

h = C . T
r
pr

Choix des inconnues du problème

Nous imagerons notre discussion en raisonnant sur le diagramme

pression­enthalpie (p, h ) . Nous exprimerons 1'enthalpie h par :

h* = h S + h, S,
g g 1 1

(A une pression p donnée, correspondent les valeurs de h et h prises sur

la courbe d'équilibre).

La courbe de changement de phase délimite sur ce diagramme trois

zones distinctes (pour une valeur de la pression inférieure à celle du
point critique) :

» une zone monophasique d'eau sous pression,

­ une zone diphasique eau­vapeur,

■ une zone monophasique de vapeur surchauffée.

Considérons les relations intervenant entre les quatre grandeurs

caractéristiques : pression (ρ) , saturation en liquide (S ) , température
in
(Τ), enthalpie (h ) dans chacune des régions citées précédemment.
­ En milieu monophasique (saturation en liquide égale à 0 ou 1) ,
la connaissance de la pression et de la température détermine une valeur
unique de 1'enthalpie h , soit :

Sh = 0 ; 1

h* = f (ρ,Τ)
 19

­ En milieu diphasique, pression et température sont liées par

la courbe de saturation.

A une pression et une saturation en liquide données, correspond

une enthalpie h unique. Soit :

P = P(T)

h* = fipfSj

Quel que soit le nombre de phases du milieu, deux des quatre

grandeurs introduites suffisent à décrire l'état thermodynamique du sys­
tème.

Pour traduire le changement de phase, nous déterminerons un cou­

ple d'inconnues indépendantes, à comportement strictement monotone, l'une
par rapport à l'autre (assurant la bijection de la relation entre ces
d eux inconnues) dans les trois régions précitées.

Nous excluerons donc de notre choix des variables telles que la

saturation en liquide (constante en zone monophasique) et la température
(liée à la pression dans la zone diphasique).

Seules restent alors la pression p et 1'enthalpie h .

Rien ne contrarie le choix de la pression comme première inconnue

du modèle.
χ
Par contre, 1'enthalpie h n'apparaît pas explicitement dans l'é­
quation de conservation de l'énergie.
La grandeur enthalpie intervient dans cette équation sous la
forme pondérée :

p
h l Sl hl + p
g Sg h g
p
l Sl + p
g Sg
 20

Nous prendrons l'expression de h comme deuxième inconnue du

problème. Puis nous relierons numériquement la valeur de h trouvée à une
valeur correspondante (et unique) de h , afin de définir un critère de
changement de phase.

La connaissance quantitative du couple (p,h) déterminera en­

tièrement les caractéristiques thermodynamiques du milieu.

Remarques

a) TORONYI [ 5] et KRUGER [ 4 J traitent le problème à pression et sa­

turation inconnues· Ce faisant, et à cause du choix de la saturation
comme deuxième inconnue, ils limitent leurs possibilités à la résolu­
tion du cas diphasique.

b) B ROWNELL & Co [θ] travaillent quant à eux à densité du fluide et

énergie interne inconnues. Ces variables étant en correspondance biuni­
voque avec la pression p et 1'enthalpie h :

P = PL SL + p S = Px(p,h) ε^ρ,ΐι) + p (p,h) S (p,h) = p(p,h)

£
U = h* - = f(p,h) — E
P(P/h)

les deux couples d'inconnues seront équivalents, le choix du

couple (p,h) semblant toutefois plus explicite.

Relations

Nous constatons une quasi-unanimité des auteurs, dont les

approches traitent des trois régions de la courbe pression-enthalpie,
sur la forme des relations entre paramètres et inconnues. Ces relations
 21

peuvent, par contre, différer quantiativement (problèmes de lissage...).

Emanerons les expressions sans expliciter numériquement :

a) Porosité

En développant en séri^ de Taylor autour d'une porosité initiale

ε = ε (ρ)

b) Densités

Dans la zone d'eau sous pression :

P = p
l i(p'h)

Dans la zone à vapeur dominante :

p s = p s (p,h)

Dans la zone diphasique :

ρ = ρ χ 3 χ + p v S v = p(p,h)

c) Saturations

Dans les zones monophasiques, les saturations en liquide sont

triviales.

Dans la zone diphasique :

Sx = S ^ h )

d) Enthalpies

Dans la zone diphasique, les enthalpies associées à chaque phase

prendront leurs valeurs sur la courbe de changement de phase :

h = h (p)
g g

hx = V ? )
 22

e) Température

Elle sera fonction :

- de la pression et de 1'enthalpie dans les régions monophasiques

Τ = T(p,h)

- de la pression dans la région diphasique :

T = T(p) = Ts(pJ
s

f) Viscosités

Dans les régions monophasiques

y = μ
i ι(ρ'^

Dans la région diphasique ;

μ = μ
1 ι(ρ)

g) Perméabilités

La littérature utilise généralement les équations de COREY

pour représenter les perméabilités relatives en fonction de la saturation
en liquide. Ces expressions analytiques dérivent d'un modèle physique
idéal pour un processus de déplacement par drainage où le phénomène de
vaporisation domine celui de condensation.

Soit Κ . = Κ . (S.) = Κ .(p,h)

ri ri 1 ri v

h) Conductivité thermicjue

Elle s'exprimera en fonction de la saturation.

Soit : λ* = λ* (SL) = λ* (ρ,η)

 23

1.6. Modèle retenu

La démarche suivie conduit au système final de deux équations diffé

rentielles aux dérives partielles couplées :

et. (p,h) ­2­ + α (p,h) — = V . Γ et, (p,h) Vp + et. (p,h) g 1 ­ Σ q.

1 ¿ J 4 i
3t 3t i

β (p,h) ­2­ + β (p,h) — = V . Γ β (p,h) Vp + β (p,h) g + β (p,h) Vh] ­ Σq,h.

3 t ^ a t i 1 1
dP
+ —
dt
avec :

α
ι = P — + εA
dp 3p h

a
2
= ε A
3h p

Ρ κ
a
3
=κ ("*♦ g r 9)
V
l v
g
3. = (Ph ­ p C T) ^ + th (^­) + (l­ε) p c Ä
p p
dp 3p h 3p h

β_ = e h (^­) + ερ + (l­ε) p c Ä
r pr
3h p 3h ρ

Κ ,h, K h * ,3T,
β.
3
= κ ( ­ ^ ­ ^ + ­ ^ ) + λ <Γ*
3 p h
ν, ν
ι g
ρ Κ ,h, ρ Κ h
34 = _ κ ( JL£LJ, + g rgg}
ν. ν
ι g
ß5 = χ* Ä
3h ρ
 24

Conditions aux limites

Le problème consiste à fixer des conditions aux limites réalistes,

possédant une souplesse d'adaptabilité suffisante, compte tenu du choix
du modèle physique, et une certaine simplicité mathématique.

Deux options s'offrent à notre examen.

Dans le premier cas, nous pouvons considérer les frontières du

domaine comme des éléments actifs du modèle. Les limites interviennent
alors pleinement dans la dynamique du problème. Elles sont perméables
aux transferts et permettent d'injecter ou de produire du fluide. Les con­
ditions aux limites consistent généralement en des valeurs de masse pro­
duite et de pression (ou de gradient de pression) données (KRUGER [4j ) .

Traitons rapidement le cas de la condition de masse produite.

Elle s'exprime par :

2 2
Κ ,ρ, Κ ρ , Κ ,ρ, Κ p
Κ [ ( _τΐΐ! + _τώ, 3E. (_LL_L+^SUL) S .£] . f(t)
ν, ν 3n ν, v
1 g 1 g

soit : ct(p,h) - 2 - - ß(p,h) g . η = f(t)

3n

Cette condition n'est pas indépendante du système d'équations dif­

férentielles, simulé. En effet, son comportement est régi par les valeurs
de certains quantités en des points internes au réservoir. Nous tiendrons
compte de cette expression en tant qu'équation supplémentaire, ajoutée au
système principal et résolue simultanément à ce dernier.

Dans le deuxième cas, les limites du réservoir n'influencent pas

les phénomènes étudiés, c'est à dire que les conditions composées sur
ces surfaces n'introduisent pas de quantités agissantes, mais servent seu­
lement à "fermer le problème".
 25

Les conditions, aux limites traduiront alors les nullités des

flux de masse et de chaleur à travers les frontières du modèle.

La condition de flux de masse nul correspond à un réservoir non

réalimenté par le milieu environnant.

La condition de flux de chaleur nul signifie que les quantités de

chaleur gagnées ou perdues par le réservoir sont petites en comparaison
de celles générées par le transfert de chaleur convectif forcé, lors de
la production.

Ces relations s'explicitent par :

ct(p,h) ^E ­ ß(p,h) g . n = 0 |
{L)
3n

2à­ 0 I
(Σ)
3n '

La condition portant sur le flux de masse dépend du système prin­

cipal. Nous devrons donc, comme pour toute condition à caractère non li­
néaire, la joindre à ce dernier et tenter une résolution globale.

Le deuxième groupe des conditions aux limites ne présente pas

d'avantages marquants sur le premier. La seule hypothèse fortement simpli­
ficatrice consiste à supprimer les effets gravitationnels (soit en les négli­
geant, soit en travaillant dans un plan orthogonal au vecteur accélération
de la pesanteur). Nous nous ramenons alors à un problème de NEUMANN classique
(MERCER [β] ) . Mais n'est­ce­pas là esquiver un aspect fort intéressant de
l'étude ?

Dans une première approche, nous traiterons le problème en suppo­

sant les parois du réservoir imperméables aux transferts. Nous localiserons
alors les apports ou pertes de masse et de chaleur en des points particu­
liers du réservoir.
 26

Conditions initiales

Dans le cas de milieux diphasiques, nous définirons des distri-

butions de pression et de saturation (conduisant de fait à des valeurs
initiales de p (pression) et h (enthalpie "pondérée")).

Dans le cas de milieux monophasiques, nous fixerons des distri-

butions de pression et de température, déterminant elles aussi des va-
leurs initiales de p et h.

II. TRAVAUX ANTERIEURS

A partir du système d'équations (I) , deux types d'approches

peuvent être envisagés. Uneraodélisation"locale"dans laquelle les va-
riables indépendantes sont paramétrées en fonctions des variables d'es-
pace et de temps, et qui constitue la représentation la plus détaillée.
Une modélisation "globale" prenant en compte la dépendance des paramètres
par rapport au temps seul, la plus largement utilisée jusqu'à présent.

Nous tenterons, en nous appuyant sur certains exemples éprouvés,

de dégager quelques traits généraux de ces deux sortes de modèles.

II. 1. Rapide critique des modèles globaux [_2 J [3j

D'une manière générale, ces modèles ne permettent de comprendre

ni le fonctionnement des réservoirs géothermiques, ni les mécanismes
mis en jeu au cours de leur exploitation.

Ils permettent uniquement d'accéder à. une estimation des réserves

énergétiques disponibles et de la productivité des gisements en fonction
de la porosité et de la composition des fluides.
 27

Le modèle [2] permet en outre, d'envisager deux sortes de con­

ditions :

­ des phases liquide et vapeur séparées avec abaissement du ni­

veau liquide en cours de production,

■■ un mélange liquide­vapeur,

la récupération ayant lieu en phase vapeur dans les deux cas.

Il s'avère nécessaire de supposer les conditions thermodynamiques

uniformes dans le gisement (hypothèse fondamentale).

L'on est ainsi conduit à estimer certains paramètres, soit à par­

tir de mesures, soit en étudiant de manière fondamentale certains phéno­
mènes (conditions d'équilibre liquide­vapeur en milieu poreux, transferts
de chaleur solide­fluide).

Cependant, plusieurs ensembles de valeurs, sans pour autant être

exacts, peuvent restituer une histoire de la production : on peut mettre
en doute le caractère représentatif des données de champ.

II.2. Modélisation locale

Cette présentation numérique distingue deux ensembles de varia­

bles :

­ un ensemble de doublets (ou triplets) caractéristiques de la

position d'un point dans l'espace mathématique(bi ou tridimensionnel),

­ un ensemble de singlets de variables de temps.

Les éléments de chacun de ces ensembles décrivent entièrement

les pavés dans lesquels ils sont définis.

La base du traitement numérique consistera en une discrétisa­

tion de ces ensembles, c'est à dire en une approximation du problème
réel dans un espace discret.
 28

La dimension temps sera traitée par une méthode implicite de dif-

férences finies. Nous optimiserons stabilité et convergence du modèle en
adaptant un paramètre de pondération, semblable à celui de la méthode
de CRANK NICHOLSON.

Deux possibilités se présentent pour la prise en compte des di-

mensions spatiales :

a) une formulation variationnelle du problème postulant la déter-

mination d'une fonctionnelle ; le minimum de cette fonctionnelle donne
la solution du système proposé.

Nous pouvons partager l'espace en éléments finis et écrire le

système discret correspondant.

La difficulté principale réside en l'introduction de conditions

aux limites non linéaires dans la fonctionnelle d'un système déjà peu
agréable.

b) un calcul discret par un schéma aux différences finies.

Dans ce cas, la souplesse géométrique du modèle est restreinte

par les conditions strictes de mailiage.

L'écriture mathématique correspondant aux deux options réduit

le système à une équation matricielle non linéaire de la forme :

A[X] .X=B[XJ

Cette équation sera résolue par une méthode itérative de type

NEWTON-RAPHSON.

Critère de stabilité

Les auteurs consultés se réfèrent, pour la stabilité du modèle, à

un nombre empirique définissant un pas de temps maximal (Throughput ratio
ou R^,).
 29

L'évaluation explicite des perméabilités relatives en fonction

de la saturation cause d'importantes instabilités en milieu multiphasi-
que.

En utilisant des approximations du type de celles proposées par

COREY, une erreur minime par rapport à la valeur exacte des perméabilités
relatives se répercute de façon peu contrôlable sur le reste du systè-
me et peut grandement affecter la stabilité de la solution. Le IL·, per-
met d'analyser ce phénomène.

L'expression explicite de chaque terme de débit par rapport à

chacune des variables dépendantes du problème résulte d'une stabilité
conditionnelle. Cette stabilité conditionnelle conduit à une expression
donnant un pas de temps maximal. L'utilisation d'un pas de temps supérieur
à ce maximum produirait une divergence dans les calculs.

Explicitons la quantité R,.

Pour une phase i :

Q. At
Vi" V
—
' pSi

soit le rapport du volume total de phase i passant dans le bloc pendant un

pas de temps At au volume de phase i dans le bloc.

Pratiquement, la stabilité d'une formulation d'un écoulement

multiphasique est notifiée par le rapport plus grossier.

soit le rapport du débit total en volume passant à travers le bloc pen-

dant le pas de temps At, au volume des pores dans le bloc.

Il est toutefois possible de vérifier la stabilité du modèle

local en constatant que sa dégénérescence doit s'accorder avec la solution
donnée par un modèle global.
 30

II.3. Problèmes particuliers traités

a) TORONYI_ [_s]

TORONYI simule le cas d'un réservoir hydrothermique à vapeur

dominante par une technique de différences finies.

Le puits de production s'étend sur un bloc entier du maillage,

situé à l'intérieur du réservoir.

Le débit de chaleur cédée par le réservoir dans le puits s'ex­

prime par :

q = q k + q. h
^H g g 1 1

et le débit de masse perdue par :

q ­ q, + q g

Le taux de production de vapeur est déterminé grâce au rapport

q
σ = —­*­ avec σ = σ (p,h)
g q g g v

Les valeurs de h et h (prises sur la courbe de changement

de phase pour une pression p donnée) étant calculées dans le simulateur
le débit q étant fixé au puits de production, nous connaîtrons par suite
les quantités σ , puis q et enfin q .

Il suffit donc de fixer le débit massique de production en un

point interne au réservoir pour exhiber la totalité des caractéristiques
inhérentes au puits.

Ce problème est repris dans l'étude entreprise par MERCER.

 31

b) MERCER [β]

MERCER tente aussi la simulation d'un réservoir à eau chaude,

en définissant le débit ponctuel à la production sous la forme :

q' = - Σ q( x k 'y k ) i <x - x

k)(y - yk)
k=l

avec : q décharge massique du réservoir au puits ponctuel k localisé

au point (Xk/Yk)»
δ fonction de DIRAC
1 nombre de puits.

Ces modèles ne semblent pas être exploités par rapport aux pos­
sibilités qu'ils offrent. Dans son compte­rendu, MERCER insiste princi­
palement sur l'aspect numérique du problème.

c) COATS [7]
Le travail de COATS présente quelques originalités quant au trai­
tement de certaines variables intervenant dans la formulation du problème.
Mais point dans la mise en équations, ni dans la méthode d'approxima­
tion numérique. Il simule le comportement d'un réservoir géothermique
sous des régimes établis ou semi­établis. Peu d'allusions sont faites
concernant les cas spécifiquement transitoires.

COATS innove en traitant simultanément puits et réservoir tandis

que la littérature antérieure [ 5j découple ces deux entités (le débit de
production (fixé) du réservoir devenant le débit d'entrée du puits). Cette
inclusion du puits dans le volume constituant le réservoir permet la ré­
partition du débit total désiré parmi les différentes couches de péné­
tration du tubage.

Par suite, COATS introduit un critère de "délivrance" imposant

une valeur minimale du gradient de pression créé entre le réservoir et
une ouverture de fonctionnement du puits ; au­dessous de cette valeur,
le puits ne produira pas.
 32

Le traitement numérique d'un puits multicouches requiert des

simplifications semi-implicites ; pour éviter des calculs trop impor-
tants, les dérivées partielles et variables intervenant dans l'expres-
sion du débit de production seront évaluées au pas de temps précédent
l'itération en cours. Ce mode de calcul est classique pour les problèmes
de régime non permanent.

L'écoulement dans le puits est supposé s'effectuer sans glissement

(ce qui se traduit par l'égalité de la saturation et de l'écoulement frac-
tionnel en volume), en régime turbulent pleinement développé. La gravité
se manifestant seulement dans les écoulements verticaux et non dans les
écoulements horizontaux, il est nécessaire d'introduire deux ensembles
de perméabilités relatives.

Grâce à ce modèle, COATS étudie les phénomènes physiques suivants t

1_° E_ffet^ d_'un_flash_de_ vaEeH.r_sH.r_lã. E?k^YJFi\Pi-e-ßlyi- EHL'-S.

Les équations définissant les débits ne tiennent pas compte de

l'expansion volumique, à volume massique constant, de la phase vapeur.
Pour pallier à ce fait, est défini un facteur de débit f, fonction de
la saturation en gaz et de la pression en un point de la circonférence du
puits d'une part, de la pression en un point infiniment voisin du réser-
voir, d'autre part.

2? C_ompa_raison_ entre modèle_numérique_et_ résolution analytique

Les tests du modèle par un problème monodimensionnel s'avérant
généralement peu probants, un problème radial sert d'essai préliminaire
à la stabilité du simulateur et à l'appréciation de l'erreur de tronca-
ture sur la dimension temps. L'auteur expose un mode de maillage radial
à partir de la spécification d'un rayon imaginaire à l'intérieur du
puits et d'un rayon extérieur définissant la limite du domaine.

3_° Etude_ d_'un_prob_lème_ b_id_imen_si_onnel_à_un_ P_uit:L

Sont plus particulièrement observées les influences sur la récupé-
ration des positions respectives de l'ouverture du puits de production
et de la base du tubage.
 33

4_° £imulat_io_n_d_|_un_ £éseryoir_ à_ jnatriçes_f racturées_

L'écoulement dans le système à blocs matriciels + fissures est
modélisé en considérant une formation sans fracture avec une grande perméa-
bilité effective reflétant la conductivité du système.

Toutefois, les formations étudiées sont très faiblement fractu-

rées (1 % du volume total est occupé par les fractures), ce qui enrobe
d'évidences les conclusions quant à la nature des transferts dans les
fissures (prépondérance des forces de gravité...).

5° Extraction_de_ chaleur, à_ partir_de_ _la_roche_chaude_ et_sèche

On effectue une injection d'eau froide dans un réservoir à matri-
ces + fractures chauffé et sec. COATS examine l'influence de la longueur
de la fracture, et de la localisation de l'injection sur les grandeurs
caractéristiques de production (température de sortie, fraction énergétique
récupérée...) .

III. REMARQUES GENERALES - CONCLUSION

D'une façon générale, ces modèles numériques ne considèrent

les phénomènes physiques que pour examiner la manière dont ils inter-
viennent sur la production.

L'esprit des travaux étudiés est soumis à cet objectif : prédo-

minance de l'exploitation avec recensement puis analyse des différents
facteurs physiques ou techniques influençants.

Nous prendrons le contrepied de cette optique. Le but avoué de

notre futur travail est la caractérisation des transferts de chaleur et
de masse entre les différentes phases du milieu considéré.

Nous nous attacherons en particulier, à mettre en évidence la

progression du front*de vaporisation. Nous utiliserons les quantités
 34

privilégiées par les auteurs consultés comme "pseudo" paramètres

du système, pouvant réagir de concert avec l'évolution des phénomènes,
et non comme finalités à notre étude.

Nous tenterons d'organiser notre travail de la manière suivante ¡

a) modélisation d'un réservoir horizontal (forces de gravité né-

gligées) avec puits de production ponctuel,

b) modélisation d'un réservoir en "élévation". Introduction de

puits de production multicouches,

c) à partir des expériences précédentes, extensions possibles :

- tridimensionnel,
- réalimentations diverses,
 35

PROJET D'INSTALLATION EXPERIMENTALE

La formulation théorique du problème physique ayant été établie, il

convient de préciser maintenant les conditions expérimentales qui seront
réalisées pour simuler les phénomènes que nous envisageons d'étudier. Pour
des raisons de facilité , la simulation sera conduite sur une installation
de configuration géométrique simple, dans laquelle transferts de chaleur et
de masse seront unidirectionnels.

Au point de vue pratique, cette configuration peut être assimilée à une

portion du domaine au sein duquel se développe la double phase par suite de
la décompression.

/M/1Î
.·,. ·„· , .! . <· ,1. ;]

'X^r.'^'::^.
•c^.^­';,«­ i,..··.·, i .'"«-* ·.'·*' 7·ν·ν

:l:
 36

Dans une telle approche, les conditions seront supposées telles que,
sur la frontière latérale, les flux thermique et massique seront nuls, qu'à
l'amont Ρ, T seront imposées ou bien les densités des flux seront nulles et
qu'à l'aval, la densité de flux thermique sera nulle et la pression imposée.

II.1. DESCRIPTION DE L'INSTALLATION EXPERIMENTALE

Elle comprendra (figure 1) :

­ une cellule cylindrique contenant le milieu poreux,

­ un dispositif d'alimentation en eau,

­ les dispositifs de mesure des grandeurs physiques essentielles :

saturations (absorption de rayons γ) , pressions, températures (capteurs) ,
débits massiques (condenseur),

·■ les appareils destinés à l'enregistrement de ces grandeurs.

II. 1.1. Cellule d'essais^

La cellule expérimentale sera conçue pour fonctionner dans des condi­

tions thermodynamiques moyennes :

P ^ 10 bars T ^ 180°C
o — o —

aisément réalisables en laboratoire, ne nécessitant pas, compte tenu des

volumes mis en jeu, de tests du comportement mécanique exigés par le Service
des Mines et permettant de concilier partiellement certains éléments contra­
dictoires tels que : bonne tenue en pression (nécessité de parois épaisses)
et faibles perturbations thermiques dues à l'enveloppe en contact avec le
milieu poreux (parois peu conductrices).
 37

Afin de limiter les effets thermiques parasites que pourrait entraîner

la mise en place d'embouts à l'extrémité aval, nous conviendrons de réali­
ser une cellule symétrique, condition qui permettra de satisfaire à la nul­
lité des flux dans le plan de symétrie χ = 0 du modèle et d'assurer une meil­
leure simulation du cas réel.

Compte tenu de la conception générale choisie, la cellule expérimentale

dont une représentation schématique est donnée sur la Figure 3, se composera
d'un tube de 200 cm de longueur, 5 cm de diamètre, contenant le milieu po­
reux. Pour limiter les effets thermiques radiaux perturbateurs, dus à la
conduction longitudinale à travers l'enveloppe, celle­ci sera constituée par
un tube en résine époxy de 1 cm d'épaisseur. La tenue mécanique de ce tube
sera renforcée au moyen d'un tube métallique extérieur en durai de 5 mm d'é­
paisseur.

Le milieu poreux mis en place dans la cellule sera limité longitudina­

lement par des embouts métalliques pourvus de pastilles poreuses (figure 2).
Les différences de dilatation entre le tube de durai et l'enveloppe de résine
seront compensées au moyen de ressorts tarés.

^extérieur

Pour limiter les fuites thermiques vers l'extérieur ainsi que pour
assurer une plus grande rapidité aux opérations de réchauffage du massif
avant expérience, nous munirons le tube métallique destiné à renforcer l'en­
veloppe de résine de résistances chauffantes dont l'alimentation électrique
sera réalisée par l'intermédiaire d'une chaîne de régulation.

Cette régulation s'effectuera à deux niveaux : soit par rapport à la

température ce consigne correspondant à la condition initiale, soit en cours
d'expérience par rapport à la température régnant sur l'axe du milieu poreux.
Dans ce cas, le rôle de la régulation consistera à minimiser l'écart de
température entre milieu poreux et enveloppe extérieure. Le réglage de cet écart
peut être réduit au maximum à +_ 1,25°C . Le tube de durai sera isolé
de l'environnement par des coquilles de laine de verre.
 33

Etant donné qu'en cours d'expérience la régulation ne s'effectuera que

dans une section du modèle (non prise en compte d'éventuels transferts longi-
tudinaux provenant des extrémités) et que, par ailleurs, même avec (D = 0
il sera difficile d'évaluer la contribution de l'enveloppe sur le transfert
de chaleur , il s'avérera nécessaire d'apprécier ces phénomènes.

II.1.1.2. Alimentation_en eau et 'nise_en_gression_du_massif

L'alimentation en eau du milieu poreux ainsi que sa mise en pression

seront réalisées au moyen d'un réservoir étanche à l'intérieur duquel une
enveloppe déformable remplie de gaz sera maintenue à pression constante, à
l'aide d'un dispositif bouteille-détendeur. L'eau, éventuellement préchauf-
fée, sera injectée dans le milieu poreux à travers les embouts métalliques
portés à la température de fonctionnement T par circulation d'un fluide
thermostaté (figure 1).

II.1.1.3. Mise_en_groduction

Elle sera réalisée par ouverture d'une vanne micrométrique placée à

l'extrémité aval du condenseur X = 0. Le pilotage de la production s'effec-
tuera par laminage du condensât à travers la vanne, ce qui permettra d'éviter
l'apparition de phénomènes d'instabilité.

L'orifice de sortie des fluides, dont le diamètre est de 5 mm, sera

pourvu d'une capsule poreuse maintenant le milieu poreux en place et loca-
lisé sur une génératrice du plan horizontal diamétral à égale distance (X = 0)
des extrémités de la cellule (figure 3)·
 39

I I . 1 . 2 . Moyens de mesures

■^lil.----. ' -ËBurË_É2f_gEaD^su£s_E!}ysi3y55_

T i_Pi_§
L'évolution des conditions thermodynamiques T(X,t) et P(X,t) sera
appréciée au moyen de 24 thermocouples et de 12 capteurs de pression disposés
le long d'une génératrice de la cellule (figure 4). Les indications délivrées
par les sondes de température et les capteurs de pression seront suivies, soit
sur enregistreur potentiométrique, soit sur centrale de mesure avec sortie
numérique.

Les profils de saturation S(X,t) permettant de suivre la manière dont

évolue le changement de phase au sein du milieu poreux, pourront être déter­
minés par une méthode désormais courante fondée sur l'absorption d'un rayon­
nement γ.

Le choix de la source et son activité ont été effectués de telle sorte,

que compte tenu des coefficients d'absorption des différents matériaux tra­
versés par le rayonnement (enveloppes et milieu poreux, soit sec, soit com­
plètement saturé), le contraste reste optimal dans une gamme de variation de la
porosité comprise entre 15 et 38 % et afin que les mesures restent signifi­
catives pour des temps de comptages faibles(^ 10 x) par rapport à la durée
moyenne d'une expérience (2 heures).

La prise en considération de ces contraintes nous a conduit à choisir

une source d'Americium 241 de 1 Curie [3j

La source et le détecteur seront montés sur un chariot mobile se dépla­

çant au­dessus du modèle.

De conception classique, la chaîne de comptage associée au détecteur

est schématiquement présentée sur la figure 5.

11^1.2^2. Mesure des_débits_liquide­vap_eur

A la sortie de la cuve réservoir, le débit d'injection dans le modèle

pourra être mesuré à l'aide d'un débitmètre à flotteur.
 40

Le débit global de sortie (eau + vapeur) sera déterminé par pesée

après passage du mélange diphasique à travers un condenseur. Nous utiliserons
d'ailleurs le même dispositif pour mesurer et suivre la loi d'évolution
de la quantité de chaleur produite au cours de la détente (enthalpie des
fluides produits).

iîiiiii^· 2ÉÏË£5}iDa£i22_55_llËD£î}alBiË_Ê55_ÎlliideË_B£°duit5 :
calorimètre ­
condenseur

La méthode que nous avons retenue pour déterminer le titre du mélange

sortant du milieu poreux est une méthode calorimétrique fondée sur la con­
densation et le refroidissement des fluides produits. Une méthode fondée sur
la vaporisation du mélange aurait pu également être envisagée. Elle n'a pas
été rentenue pour l'instant, car elle présente de plus grandes difficultés
de mise en oeuvre.

Principe de la mesure

Le principe du dispositif de mesure est représenté sur le schéma ci­

dessous. Il s'agit d'un échangeur pouvant fonctionner à co ou contre­courant.
Le mélange d'eau et de vapeur issu du milieu poreux traverse le condenseur
dont le refroidissement est assuré par une circulation d'eau. Les températu­
res d'entrée et de sortie, ainsi que les débits masse du mélange et du fluide
réfrigérant étant mesurés, si l'on suppose que les deux phases sont animées
de la même vitesse dans le circuit de récupération, il est. possible de dé­
terminer le titre du mélange y = m /m + m = m /m et par suite les rapports
v e ν v
des débits.

t 'f!

Π
f2
 41

Si h désigne les enthalpies en phase liquide, X. la chaleur latente

de vaporisation à pression constante, 1 et 2 respectivement à l'entrée et
à la sortie de l'appareil et m le débit masse de l'eau de refroidissement,
la seule équation du bilan thermique suffit à déterminer Y.

Vhl ­ V +
V ^ + h
l ­ h2) = m (h
f f2 ­ hfl>

soit :

m(h ­ h 2 + Y¿£) = m f (h ­h )

d'où :

m
1 f
Y= [ (h
Zp ­ f2 ­ V ­hl + h 2 3

Pour fixer les idées, choisissons T = 150°C, T = 20°C, m = 0,2 et

0,4 g/s de débit masse eau + vapeur d'eau et m = 10 g/s. En se plaçant dans
le cas contre­courant avec T = T = 20°C donc h„ = h ., l'expression ci­
dessus devient :

mf
Y = [ ­ (h.. ­ h2) ­ h. + h. ]
t m

et

h f 2 ­ h2 = — (h1 ­ h 2 + γ.^:
m
f
Dans notre c a s h ­ h = 547 J/g »C = 211!
¿C 2115 J/g

h ­ h 2 = (m/10) (547 + 2115 Y)

 42

y = 0 y = 1

m m/10 (hf2­h2)/4,18 TX °c (hf2­h2)/4,18 T

f2 C f2

0,2 0,02 2,62 22,62 12,74 32,74

0,4 0,04 5,23 25,23 25,47 45,47

Calcul d'erreur sur y

Si l'on appelle Ah f l'erreur sur h ­ h et Ah l'erreur sur h~ ­ h ,

l'erreur sur Y s'écrit :

Y = - ~
1 m

[ Ah + —
f
Ahf + (hf2 ­ hfl)
^f ΛΜ
(—­ + — +2At)
""f -£.
oC m m m m.t eL

où M représente la masse d'eau produite pendant t servant à déterminer m,

et M la masse d'eau de refroidissement mesurée pendant le même intervalle
de temps pour calculer m f . Etant donné que pendant nos expériences m f sera
maintenu constant, la mesure de ce débit peut être extrêmement précise et
ΔΜ /m = 0.

Si avec les données précédentes, on adopte t = 60 s avec At = 1/2 s ¡

,ΔΜ _ , ,. ΔΜ, ,,
(— + 2 At) m /m.t = (1 + — ) /6 m
m m

Si : T ­ T sont mesurés à 0,5°C près, donc également h ­ h . .

T ­ T sont mesurés à 2°C donc également h. ­ h„

­3 vJ?
et T à 2°C près, soit une erreur relative de 3.10 sur <*­ on obtient
pour m = 0,2, avec m /m = 10/0,2 = 50 ; ΔΜ/m = 0,5/0,2 (en supposant l'er­
reur de mesure sur la masse condensée de 0,5 g ) , soit :

(1 + ΔΜ/πι)/6 m = (1 + 2,5)/l,2 = 2,92

\
 43

les résultats suivants :

m
Y yAÍ¿ h h
f2" fl
m
f
­ A h f 2,92 Ah [ 1 ¿Il ΔΥ
(h h
f2 ­ fl}

0 0 2,62 25 7,65 2 34,65 0,016 0,016

1 0,003 12,74 25 37,6 2 64,6 0,03 0,03

résultats qui montrent que sans prendre de dispositions draconiennes, un tel

dispositif permettra d'accéder à des mesures satisfaisantes.
 VW-p^-φ
ι r'i
Ju,
ii
WMzzm ^////i///777i\^\\V\\\W
» ï ·:: ·.:-·.·:·· · . ¡
Ov vw, V é y ,
p\WvY\\^f
1 |:

p-V?, W77T/Z7ρ ~ · - ι χ . : . . . . Ι - : ι ■ ' . ..i........L -.J p^.:....).,--..^-^^;

®4^
7
: cellule
: thermocouples et
capteurs ¿e pression
: eriouts de préchauffaço
(«) Pression
Θ Température
: calorimètre condenseur
L
et :·;' : source et détecteur
rcäurii c¿5 saturations A
Θ
cuve ce r.iso en pression
t:.er:-o£ti-t Po
r­sçulûticn do t.4np<£ra­
t u r e ce Cuve c'è m/se
¿ovólos:;«
e t :; : e n r e g i s t r e u r s
©
pression e t te­.pératuro en pression

SCHEMA DE PRINCIPE DE L'INSTALLATION

C5
 45
Λ'ΓΛ 11
Fig.2
£01
1
Π e 'Ο s orí ' c í ^
/, Kli Ρ η s î i 11 o poreuse
\ /
R'

m '7 \ Γ\ ' * · · · ·« · · . * · · · · . · ··*'.·

Μ \ :J ι:. - ι*·'-'. •.·ΐ·.·Ίν\^Λ»ι\ ' / V
\'■ / / / /
z■7^
CCT
Milieu-poreux

Thermocouple

χ-

Condenseur
Vanne micrométrique
Capsule poreuse
de sortie
en acier inox
 POSITION DES CAPTEURS

IMA : thermocouples milieu poreux

TIA : thermocouples cellule
TDA : thermocouples enveloppe durai

fenêtre de scrutation position d'as colliers chauffants

'par rayons X
DS5

ç
TMA12I TMA11 ITMA10 TMASO
t u £ C Ι Β ï
û ι - Γ «Û
TMAQ8 'ΤΓ.1Α07 TMAOSj TMA05 ΤΓν*Α04 ΤΓΛΑ03
TMA02 TiYiAOI J J
ΓίΊ i ! i 2 u S40 ! CÕ0
-i
780
5—ι
700^
h "
620 _>: Ü_,—
40 400 ' 330
1
•co 22C
140 ! GO Γ
poreux

TIAS Tl Ä7 TIA6 TIJAS TIA* TPA3 TIA2 TI Al

isolant £20 6G0 440 200
1
— ι —
!-
r
Dur TDA7 TDA5 Τ D A3 Τ D Al
550 220
 43
NG. 5

Alimentation
stabilisée collimateurs
BT 7 000 F
HT 7 127

PréampU.
1
cristal Haï source Am
Amplification
Mise en forme
Enregistreur indication
le to me tre linéaire du taux de
(intégrateur) comptage

Sélection
d'amplitude
(seuil et fenêtre)

Echelle double avec Imprimante

horloge (comptage)

Programmateur
d'amplitude

Schema synoptique de la chaîne de comptage

 49

CONCLUSION

Les bases e s s e n t i e l l e s à la poursuite d'une étude sur les phéno-

mènes de changement de phase avec ou sans réinjection permettant d ' a p -
p r é c i e r le caractère r e p r é s e n t a t i f des modèles mathématiques de d e s c r i p -
tion des phénomènes, nous p a r a i s s e n t , à l'heure a c t u e l l e , acquises.

En vue d'une poursuite du programme dans l a perspective d'un pro-

longement du c o n t r a t , nous tenons à souligner qu'au cours de ces der-
n i e r s mois, un c e r t a i n nombre de r é a l i s a t i o n s complémentaires ont été
effectuées. Elles concernent principalement le lancement de l'exécution
de l ' i n s t a l l a t i o n d ' e s s a i , l ' é l a b o r a t i o n de programmes numériques con-
cernant les l o i s d'évolution des p r o p r i é t é s physiques du fluide en
fonction des conditions thermodynamiques e t le démarrage d'un premier
programme de simulation numérique du système d'équations des phénomènes,
dans le cas u n i d i r e c t i o n n e l .
 50

REFERENCES BIBLIOGRAPHIQUES

P. KRUGER
Geothermal Energy Development
Gov. Rep. Announce (1977)

D.B. LOMBARD
D e v e l o p m e n t and u t i l i z a t i o n o f g e o p r e s s u r e d g e o t h e r m a l resource
in t h e Couth C e n t r a l United S t a t e s
G e o l o g y o f A l t e r n a t e E n e r g y R e s o u r c e s (1977)

HOUSTON GEOL. SOC.

S e l e c t e d Geopressured Energy B i b l i o g r a p h y
Geology of A l t e r n a t e Energy R e s o u r c e s p p . 337-342 (1977)

SMITH M . , POTTER R . , BROWN D . , AAMODT T .

I n d u c t i o n and growth of f r a c t u r e s i n h o t rock
S t a n f o r d U n i v e r s i t y P r e s s p p . 251-268 (1973)

R.D. Mc FARLAND, H.D. MURPHY

E x t r a c t i n g E n e r g y from H y d r a u l i c a l l y f r a c t u r e d g e o t h e r m a l
reservoirs
11th I n t e r s o c . Energy C o n v e r s i o n Eng. Conf. p p . 828-835 (1977)

A . C . GRINGARTEN, H . J . RAMEY
E f f e c t o f h i g h volume v e r t i c a l f r a c t u r e s on g e o t h e r m a l steam
well behavior
2nd J . N . D e v e l o p & Use o f G e o t h e r m a l R e s o u r c e Sump. p . 1759-176
(1976)

R. CELATI, P . SQUARCI, G.C. STEFEANI, L . TOFFI

A n a l y s i s of water l e v e l s and r e s e r v o i r p r e s s u r e measurements
in geothermal w e l l s
2nd U.N. D e v e l o p & Use o f G e o t h e r m a l R e s o u r c e s Symp. p p . 1 5 8 3 -
1590 (1976)

S.K. GARG, T . R . BLAKE, D.H. BROWNELL, A . H . NAYFEIH, J . W . PRITCHETT

Simulation of f l u i d - r o c k i n t e r a c t i o n s in a geothermal b a s i n
Gov. R e p . Announce (1977)

A. HUNSBEDT, P . KRUGER, AL. LONDON

Laboratory s t u d i e s of f l u i d p r o d u c t i o n from a r t i f i c i a l l y frac-
t u r e d geothermal r e s e r v o i r s
SPE 6 5 5 3 (1977)
 51

K.H. CHEN
Measurement of water content of porous media under geothermal
system conditions
Gov. Rep. Announce (1977)

A. HUNSBEDT, P. KRUGER, A.L. LONGON

A laboratory model of stimulated geothermal reservoirs
Gov. Rep. Announce (1977)

B.C. Mc CABE, E. ZAJAC

High temperature geothermal energy system
Magma Energy Inc. (1976)

W.F. FRANZ, H.V. HESS

R e c o v e r y from g e o t h e r m a l reservoirs
T e x a c o I n c . (1977)

W.E. BRIGHAM, P . KRUGER

H e a t and mass t r a n s f e r i n p o r o u s r o c k fragments
S t a n f o r d U n i v . R e p o r t D e c . 1975

H.K. CHEN
Measurement of water content of porous media under geothermal
system conditions
Tulsa N° 236742 (1976)

C R . WILSON, P.A. WITHERSPOON

Steady state flux in rigid networks of fractures
Waters Resources (1974)

STANFORD GEOTHERMAL PROGRAM

Stanford University, Stanford, California

R.M. TORONYI
Two phase, two dimensional simulation of geothermal reservoir
and the wellbore system
SPE 5521 (1975)

JW. MERCER, C R . FAUST

Simulation of water and vapor dominated hydrothermal reservoirs
SPE 5520 (1975)

J.W. MERCER
Finite element approach to the modelling of hydrothermal systems
Ph. D. Thesis Illinois (1973)

C R . FAUST, J.W. MERCER

Mathematical modelling of geothermal' systems
2nd U.N. Symposium on Development and Use of Geothermal Resources
(1975)
 s¿

T.J. LASSETER, P . A . WITHERSPOON, M . J . LIPPMANN

The n u m e r i c a l s i m u l a t i o n o f h e a t a n d mass t r a n s f e r i n m u l t i ­
d i m e n s i o n a l two p h a s e g e o t h e r m a l r e s e r v o i r s
2nd U.N. Symposium on D e v e l o p m e n t a n d Use o f G e o t h e r m a l R e s o u r c e s
(1975)

C . L . Mc MICHAEL, G.W. THOMAS

R e s e r v o i r S i m u l a t i o n by G a l e r k i n ' s Method
S o c . P e t . E n g . J . (1973)

N . D . SHUTLER
Numerical three phase model of the two dimensional steam flood
process
SPE 2798 (1970)

K.H. COATS, W.D. GEORGE, CHIEH CHU, Β.E. MARCUM

Three dimensional simulation to steamflooding
SPE 4500

P. CHENG, L. TECKCHANDANI
Numerical solutions for transient heating and withdrawal of
fluid in a liquid dominated geothermal reservoir
Gov. Rep. Announce (1977)

R.N. HORNE, M.J. O' SULLIVAN

Numerical modeling of the Wairakei geothermal reservoir
SPE 6356 (1977)

C R . FAUST, J.W. MERCER

Finite difference model of two dimensional, single and two
phase heat transport in a porous medium, vertion I
U.S. Geol. Surv. Open File Remp. (1977)

OF. ISOKRARI
Numerical simulation of U.S. gulf coast geothermal geo pressured
reservoirs
Diss. Abstr. Int. Sect. Β (1977)

K.H. COATS
Geothermal reservoir Modelling
52 Annu. SPE of AIME Fall Tech. Conf. SPE 6892 (1977)

WHITING­RAMEY
Application of material energy balances to geothermal steam pro­
duction
Journal of Petroleum Technology (1969)

E. BRIGHAM, B. MORROW
P/Z behavior for geothermal steam reservoirs
SPE Paper N° 4899 (1974)
 53

P. CHENG
Simularity solutions for mixed convection from horizontal
impermeable surfaces in saturated porous media
Gov. Rep. Announce (1977)

P. CHENG, W.C CHAN

Similarity solutions for convection of groundwater adjacent
to horizontal impermeable surfaces with axisymmetric tempera-
ture distribution
Gov. Rep. Announce (1977)

P. CHENG
Combined free and forced boundary layer flows about inclined
surfaces in a porous medium
Gov. Rep. Announce (1977)

P. CHENG, I.D. CHANG

Buoyancy induced flows in a saturated porous medium adjacent
to impermeable horizontal surfaces
Gov. Rep. Announce (1977)

C R . FAUST, J.W. MERCER

Theoretical analysis of fluid flow and energy transport in
hydrothermal systems
W.S. Geo!. Survey Open File Rep. (1977)

A.F. MOENCH
Simulation of steam transport in vapor dominated geothermal
reservoirs
U.S. Geol. Survey Open File Rep. 76-607 (1976)

D.R. KASSOY
Some remarks on t h e continuum d e r i v a t i o n of t h e e q u a t i o n s for
flow i n t h e r m a l l y a c t i v e porous media
Nat. S c i . Found. & Base Western Reserve - Univ. Heat & Mass
T r a n s f e r in Porous Media Workshop P r o c . p p . 88-102 (1976)

A.V. LUIKOV
Heat and Mass transfer in capillary porous bodies
Pergamon Press (1966)

J. BEAR
Dynamics of Fluids in porous media
American Elsevier
Environmental Science Series

P. KRUGER, H.J. RAMEY,Jr

Stimulation and reservoir engineering of geothermal resources
Stanford Geothermal Program - Stanford Univ. - Stanford
Californie
 54

D.M. BROWNELL Jr, S.K. GARG, J.W. PRITCHETT

Computer Simulation of geothermal reservoirs
Ventura, California, 45th Annual California Regional Meeting
of the S.P.E. ­ 1975 ­ Paper SPE 5381

S. BORIES, K. GHALLEΒ
T r a n s f e r t s de c h a l e u r e t de masse dans l e s m i l i e u x poreux
non s a t u r é s
1972 ­ G.E. I.F.P.­I.M.F ­ ENSEEIHT ­ TOULOUSE

G. CADY
Model Studies of geothermal fluid production
Ph. D. ­ Stanford University ­ 1969

GALAUP
C o n t r i b u t i o n à l ' é t u d e des méthodes de mesure en écoulement
d i p h a s i q u e ­ A p p l i c a t i o n à l ' a n a l y s e s t a t i s t i q u e des é c o u l e ­
ments à b u l l e s
Thèse de D o c t e u r ­ I n g é n i e u r ­ Grenoble ­ 1975
CDNA06627FRC