Modèles mathématiques et numériques avancés pour la simulation du polymère dans les réservoirs pétroliers
| Auteur / Autrice : | Guissel Lagnol Dongmo Nguepi |
| Direction : | Quang Huy Tran, Christophe Berthon, Benjamin Braconnier, Christophe Preux |
| Type : | Thèse de doctorat |
| Discipline(s) : | Mathématiques et Informatique |
| Date : | Soutenance le 13/12/2021 |
| Etablissement(s) : | université Paris-Saclay |
| Ecole(s) doctorale(s) : | École doctorale Sciences et technologies de l'information et de la communication (Orsay, Essonne ; 2015-....) |
| Partenaire(s) de recherche : | Laboratoire : Institut français du pétrole Énergies nouvelles (Rueil-Malmaison, Hauts-de-Seine) |
| Référent : CentraleSupélec (2015-....) | |
| graduate school : Université Paris-Saclay. Graduate School Informatique et sciences du numérique (2020-….) | |
| Jury : | Président / Présidente : Pauline Lafitte-Godillon |
| Examinateurs / Examinatrices : Claire Chainais-Hillairet, Paul Vigneaux, Bijan Mohammadi, Rodolphe Turpault | |
| Rapporteurs / Rapporteuses : Claire Chainais-Hillairet, Paul Vigneaux | |
| DOI : | 10.70675/9e8a52a3zf233z4cebz828cz4e55822ec5b2 |
Résumé
Une technique efficace pour accroître la production d’un champ pétrolier consiste à y injecter un mélange d’eau et de polymère. La viscosité du polymère réduit en effet la mobilité de l’eau, qui pousse alors mieux l’huile, d’où un taux d’extraction plus élevé. La simulation numérique d’un tel procédé de récupération d’hydrocarbures revêt donc d’une importance capitale. Or, malgré des décennies de recherche, la modélisation des écoulements avec polymère en milieu poreux et sa résolution numérique demeurent un sujet difficile. D’une part, les modèles habituellement employés par les ingénieurs de réservoir présentent, au mieux, des singularités de type résonance qui les rend faiblement hyperboliques. Ce défaut donne lieu à certaines complicationsD’une part, les modèles habituellement employés par les ingénieurs de réservoir présentent, au mieux, des singularités de type résonance qui les rend faiblement hyperboliques. Ce défaut donne lieu à certaines complications mais reste acceptable. Au pire, quand on veut incorporer l’effet du volume de pore inaccessible (IPV), les modèles deviennent non hyperboliques, ce qui aggrave les instabilités numériques susceptibles d’apparaître.D’autre part, les schémas numériques classiques ne conduisent pas à des résultats satisfaisants. Sans IPV, la diffusion excessive autour de l’onde de contact fait perdre les informations pertinentes. Avec IPV, l’existence des valeurs propres complexes crée des instabilités exponentielles au niveau continu qu’il faut traiter au niveau discret sous peine d’arrêt prématuré du code.L’objectif de cette thèse est de remédier à ces difficultés. Au niveau des modèles, nous analysons plusieurs lois d’IPV et établissons une équivalence entre deux d’entre elles. Nous proposons de surcroît des conditions suffisantes raisonnables sur la loi d’IPV en vue de l’hyperbolicité faible du système d’écoulement. Au niveau des schémas pour le problème sans IPV, nous préconisons une correction afin d’améliorer la précision des discontinuités de contact. Pour le problème avec IPV,nous élaborons une méthode de relaxation qui garantit la stabilité des calculs quelle que soit la loi IPV.