Analyse asymptotique du problème de Riemann pour les écoulements compositionnels polyphasiques en milieux poreux et applications aux réservoirs souterrains, Asymptotic analyse of Riemann problem for multiphase compositional flow in porous media with application to subterranean reservoirs

AVERTISSEMENT



Ce document est le fruit d’un long travail approuvé par le jury de
soutenance et mis à disposition de l’ensemble de la communauté
universitaire élargie.
Il est soumis à la propriété intellectuelle de l’auteur au même titre que sa
version papier. Ceci implique une obligation de citation et de
référencement lors de l’utilisation de ce document.
D’autre part, toute contrefaçon, plagiat, reproduction illicite entraîne une
poursuite pénale.

Contact SCD INPL: mailto:scdinpl@inpl-nancy.fr




LIENS




Code de la propriété intellectuelle. Articles L 122.4 e la propriété intellectuelle. Articles L 335.2 – L 335.10
http://www.cfcopies.com/V2/leg/leg_droi.php
http://www.culture.gouv.fr/culture/infos-pratiques/droits/protection.htm
INSTITUT NATIONAL POLYTHECHNIQUE DE LORRAINE
École Nationale Supérieure d’Électricité et de Mécanique
Laboratoire d’Énergétique et de Mécanique Théorique et Appliquée
École Doctorale EMMA
THÈSE
Présentée en vue de l’obtention du grade de
DOCTEUR DE L’I.N.P.L.
Spécialité :
Mécanique-Énergétique
Par :
Anahita ABADPOUR
Analyse asymptotique du problème de Riemann pour les
écoulements compositionnels polyphasiques en milieux poreux et
applications aux réservoirs souterrains
Soutenue publiquement le 4 Décembre 2008
Composition du jury :
Rapporteur M. Franck PLOURABOUÉ
Rapporteur M. Stéphane ZALESKI
Examinateur M. Christian MOYNE
Invité M. Peppino TERPOLILLI
Invité Mme. Irina PANFILOVA
Ditrecteur de thése M. Mikhail PANFILOVAvant Propos
Je tiens à remercier tout particulièrement Monsieur Mikhail PANFILOV, Professeur à
l’INPL, pour le travail productif et stimulant que nous avons mené ensemble.
Je remercie l’entreprise Total pour la prise en charge de son financement ainsi que pour
l’intérêt qu’ils ont porté pour les progrès des élèves dans les pays qu’ils travaillent.
Je tiens à remercier Monsieur C. MOYNE, Directeur du Laboratoire d’Énergétique et de
Mécanique Théorique et Appliquée (LEMTA-ENSEM-INPL) de me avoir accueillie dans
son laboratoire.
Je remercie également M. Franck PLOURABOUÉ, Directeur de rechertche CNRS à l’ins-
titut Mécanique des Fluides, Toulouse, M. Stéphane ZALESKI, Professeur à l’université
Paris-6, M. Christian MOYNE, Directeur de recherche CNRS à l’institute National Poly-
technique deLorraine, MmeIrinaPANFILOVA Doctoràl’Institut NationalPolytechnique
deLorraineetM.Peppino TERPOLILLI Doctoraucentre scientifique deTotalàPaupour
avoir accepté de juger mon travail en participant au jury.
Je remercie tout mes amis du Laboratoire pour leur accueil et leur soutien durant ces
années de travail.
Je tiens particulièrement à remercier mon conjoint, Vahid, pour me supporter avec tous les
hauts et les bas de ces années et merci aussi à mes parents, Shahrbanoo et Ghasem, pour
me décerner l’enthousiasme de l’apprentissage.Qualités et spécialités des membres du jury
Franck PLOURABOUÉ Directeur de rechertche CNRS à l’institut Mécanique
Rapporteur des Fluides, Toulouse. Group d’Étude sur les Milieux
Poreux.
Spécialité : Milieux fluides et réactifs
Stéphane ZALESKI Professeur,L’universitéPierreetMarieCURIE(Paris-
Rapporteur 6),Institut d’Alembert. Tour 55, case 162, UPMC
75252 PARIS Cedex 05, France.
Spécialité : Mécanique des fluides
Christian MOYNE Directeur de recherche CNRS, l’institute National Po-
Examinateur lytechnique de Lorraine, ENSEM, LEMTA, BP 160,
2 avenue de la Foret de Haye, 54502, Vandoeuvre lès
Nancy, France. Spécialité : Mécanique et énergétique,
milieux poreux
Peppino TERPOLILLI Doctor des Mathématiques, Ingénieur au centre scien-
Invité tifique de Total, Avenue Larribau 64018 Pau cedex,
France.
Spécialité : Spécialité : Mathématiques appliquées
Irina PANFILOVA Doctor, L’institute National Polytechnique de Lor-
Invité raine, ENSEM, LEMTA, BP 160, 2 avenue de la Foret
de Haye, 54502, Vandoeuvre lès Nancy, France.
Spécialité : Mécanique des fluides et milieux poreux
Mikhail PANFILOV Professeur, L’institute National Polytechnique de Lor-
Directeur de thèse raine, ENSEM, LEMTA, BP 160, 2 avenue de la Foret
de Haye, 54502, Vandoeuvre lès Nancy, France.
Spécialité : Génie Mécanique, Génie civilRésumé
l’écoulement et le transport de deux ou plusieurs phases dans les milieux poreux sont
observés dans la récupération assistée du pétrole (EOR : Enhanced Oil Recovery) dont
l’objectif est déplacer efficacement les hydrocarbures par l’injection de fluides qui solubi-
lisent et / ou mobilisent la phase oléique, en la poussant vers les puits de production.
Dans la récupération du pétrole, de l’eau ou du gaz sont injectés dans le réservoir en créant
les fronts miscibles ou partiellement miscibles entre les fluides.
Aucoursdel’injection, lescomposants dufluide injecté se dissolvent partiellement oucom-
plètement dans la phase initiale. Ce processus, qui détermine un échange de masse entre
les phases, conduit à des changements souhaitables dans les propriétés de l’huile telle que
la une réduction de sa viscosité et / ou l’évaporation partielle ou complete.
En réalité, la situation est souvent complexée par la présence d’un troisième fluide, autre
que le fluide injecté et le fluide déplacé. Il s’agit de l’eau résiduelle au venue de l’aquifère
adjacent lors de l’exploitation du réservoir elle joue un rôle important dans le processus de
déplacement. En particulier, une partie du gaz injecté dissolvant dans l’eau peut réduire
l’efficacité de EOR. Un autre exemple de présence de l’eau est WAG : l’injection artificielle
des tranches d’eau alternées avec le gaz. Cette technique permet, d’une part, augmenter
l’efficacité de EOR ( à couse de l’augmentation moyenne de la viscosité du fluide dépla-
çant),maisaussiderésoudreleproblème destockagedesgazacides(CO )danslescouches2
géologiques.
Desconsidérableseffortssontétéappliquéspourdévelopper lathéorieanalytique del’écou-
lement simultanée de deux phases partiellement miscibles dans les milieux poreux sont
7connues dans la littérature scientifique. La modélisation d’un écoulement unidimension-
nel des systèmes tri-phasiques entièrement immiscibles (pétrole/eau/gaz) a été également
intensivement développée. En même temps, trèspeuderésultatssontdisponible sur l’écou-
lement partiellement miscible en trois phases.
L’hydrodynamique de ces processus est fortement couplée avec la thermodynamique du
système multiphasique multicompositionnel, ce qui rend difficile l’analyse mathématique
des ces système et leur modélisation numérique.
Dans la pratique, tousles systèmes multiphasiques sont multicomposants et contiennentN
composants chimiques tels que les hydrocarbures à partir de CH jusqu’à C H et plus4 10 22
lourds même, ainsi que H O, CO , N et parfois le sulfure.2 2 2
Le modèle mathématique complet d’écoulement multiphasique multi-compositionnel
consiste aux équations de bilan de la quantité du mouvement pour chaque phase et aux
équations de bilan de la masse pour chaque composant chimique présent dans le système.
Cemodèledoitêtrecomplétéauxrelationsthermodynamiquesdécrivantl’équilibreetl’état
des phases.
L’hypothèse de l’équilibre local est traditionnellement appliquée dans le cadre de cette
théorie, étant basée sur le fait que la surface de contact entre les phases en milieux poreux
est très grande.
Par conséquent, l’écoulement d’un tel mélange est décrit par N équations différentielles
nonlinéaires de conservation de masse pour chaque composant etN+2 relations transcen-
dantales thermodynamiques dans le cas d’écoulement diphasique, et 2N +3 relations dans
le cas d’écoulement triphasique. L’analyse qualitative efficace d’un tel modèle est pratique-
ment impossible.
Si on considère que seule la partie thermodynamique, on constate que, la théorie ther-
modynamique de ces systèmes est bien développée dans les ouvrages fondamentaux. Sa
description mathématique, qui est souvent appelée le modèle compositionnel, consiste aux
équations d’état des phases, équations d’état rhéologiques des phases et équations d’équi-
libre reliant les potentiels chimiques de chaque composant dans chaque phase, les pres-
sions et les t