A Finite Analytic Method for Two-Dimensional Power-Law Fluid Flow in Heterogeneous Porous Media Using Unstructured Grids

Profondément sous la terre, dans le réseau complexe de roches qui renferment du pétrole, du gaz et de l'eau, les fluides ne se déplacent pas toujours selon des lignes simples et prévisibles. La roche elle-même est rarement uniforme ; c'est une mosaïque de matériaux ayant des capacités très différentes à laisser passer les fluides. Certaines couches sont comme des autoroutes ouvertes, tandis que d'autres sont comme des murs denses. Lorsque les ingénieurs tentent de modéliser comment ces fluides s'écoulent, ils doivent tenir compte de ces changements soudains dans les propriétés de la roche. De plus, les fluides eux-mêmes peuvent être capricieux. Si l'eau s'écoule facilement, de nombreux fluides industriels — tels que les suspensions utilisées dans la fabrication du papier ou les boues de forage utilisées dans les puits de pétrole — se comportent différemment. Ils sont « non-newtoniens », ce qui signifie que leur viscosité et leur résistance à l'écoulement changent en fonction de la force avec laquelle on les pousse. Lorsque ces fluides complexes rencontrent une roche dont le caractère change brusquement, les mathématiques décrivant leur mouvement deviennent incroyablement difficiles, s'effondrant souvent précisément aux points où les propriétés de la roche basculent.

Pendant des décennies, les scientifiques se sont appuyés sur une méthode standard pour calculer ces écoulements, une méthode qui consiste essentiellement à faire la moyenne des propriétés de la roche entre deux points. Cette approche fonctionne assez bien lorsque la roche est relativement constante, mais elle échoue de manière spectaculaire lorsque la roche est très variée. Dans ces cas, la méthode standard sous-estime systématiquement la quantité de fluide qui peut réellement circuler, entraînant des erreurs qui peuvent être massives. Pour résoudre ce problème, des chercheurs de l'Université de science et technologie de Chine ont développé une nouvelle technique de calcul. Au lieu de s'appuyer sur de simples moyennes, cette nouvelle méthode utilise une description mathématique précise de la façon dont le fluide se comporte exactement à la limite où deux types de roches différents se rejoignent. En capturant le comportement réel et net du fluide à ces frontières critiques, la nouvelle méthode offre une image beaucoup plus claire et plus précise de l'écoulement souterrain, même lorsque la roche est extrêmement irrégulière et le fluide complexe.

Les chercheurs se sont concentrés sur un type spécifique de fluide appelé fluide de loi de puissance. Il s'agit de fluides dont la relation entre la force qui les pousse et la vitesse à laquelle ils se déplacent suit une règle non linéaire spécifique. Dans le monde réel, cela couvre un large éventail de substances, des solutions polymères épaisses à certains types de pétrole brut. Le défi surgit lorsque ces fluides circulent à travers des milieux poreux — des roches remplies de minuscules trous — où la perméabilité, ou la facilité d'écoulement, change soudainement d'un endroit à l'autre. Dans de tels environnements, la pression du fluide ne change pas de manière fluide ; au lieu de cela, elle crée un pic aigu ou une « singularité » juste à la frontière entre la roche à écoulement facile et celle à écoulement difficile. Les modèles informatiques traditionnels, qui divisent la roche en une grille de petits triangles, luttent contre ces pics. Ils supposent que la pression change de manière linéaire entre les points, ce qui est physiquement incorrect à ces frontières. En conséquence, les modèles produisent des débits bien trop faibles et nécessitent un nombre énorme de points de grille minuscules pour obtenir ne serait-ce qu'une approximation grossière de la vérité.

Pour corriger cela, l'équipe a développé une Méthode Analytique Finie. L'idée centrale est de cesser de deviner comment le fluide se comporte à la frontière et d'utiliser plutôt une description mathématique exacte de ce comportement. Les chercheurs ont réalisé qu'à proximité du point où les propriétés de la roche changent, la pression du fluide suit un motif courbe spécifique, augmentant ou diminuant d'une manière qui dépend de la nature de la loi de puissance du fluide. Ils ont calculé la forme exacte de cette courbe, appelée fonction propre, pour la partie la plus critique de l'écoulement. Ce calcul impliquait la résolution d'un ensemble complexe d'équations décrivant comment le fluide se déplace dans une zone en forme de coin où les propriétés de la roche diffèrent. En trouvant la solution la plus petite et la plus significative de ces équations, ils ont capturé la manière dominante dont le fluide se comporte précisément au point problématique.

Une fois ce comportement local compris, les chercheurs ont construit un nouveau schéma numérique autour de celui-ci. Au lieu d'utiliser l'ancienne méthode défaillante de moyennage pour relier les points de la grille, ils ont utilisé leur nouvelle solution analytique pour calculer l'écoulement à travers les frontières de chaque cellule triangulaire de leur grille. Cela a permis à l'ordinateur de « savoir » exactement comment le fluide devrait se comporter lorsqu'il se comprime dans le goulot étroit entre les zones de haute et de basse perméabilité. Le résultat est une simulation qui respecte la physique du problème aux points les plus difficiles, plutôt que de lisser les détails par une approximation.

L'équipe a testé cette nouvelle méthode face à plusieurs scénarios exigeants pour voir comment elle se comportait. Ils ont commencé par un motif classique en « damier », où des roches à haute perméabilité et à faible perméabilité étaient disposées en carrés alternés. Dans cette configuration, les propriétés de la roche changeaient d'un facteur de cent. Lorsqu'ils ont comparé leur nouvelle méthode à la méthode traditionnelle, la différence était flagrante. La méthode traditionnelle sous-estimait le débit par des marges énormes, manquant parfois la valeur réelle de près de quatre-vingt-dix pour cent sur une grille grossière. En revanche, la nouvelle méthode restait remarquablement proche de la bonne réponse, même lorsque la grille n'était pas très fine. À mesure qu'ils affinaient la grille, la nouvelle méthode convergeait rapidement vers la bonne réponse, tandis que la méthode traditionnelle peinait à rattraper son retard, même avec des grilles extrêmement fines.

Ils ont également testé la méthode sur un domaine hexagonal régulier et sur un domaine présentant une distribution log-normale aléatoire de la perméabilité, ce qui imite la nature chaotique des véritables réservoirs souterrains. Dans chaque cas, la nouvelle méthode s'est révélée supérieure. Elle gérait les variations extrêmes des propriétés de la roche avec aisance et restait précise, quel que soit le type spécifique de fluide à loi de puissance simulé. L'une des découvertes les plus frappantes fut que la précision de la nouvelle méthode ne dépendait pas fortement des propriétés spécifiques du fluide ou de la roche. Que le fluide soit très épais ou très fluide, ou que le contraste de la roche soit modéré ou extrême, la méthode délivrait des résultats cohérents. Les chercheurs ont constaté qu'un niveau modéré de raffinement de la grille était suffisant pour obtenir une grande précision, les erreurs tombant à moins de huit pour cent même dans les pires scénarios.

Ce travail représente une étape importante dans la façon dont nous simulons l'écoulement des fluides dans des environnements souterrains complexes. En remplaçant une technique de moyennage défaillante par une solution qui respecte la véritable physique de la frontière, les chercheurs ont créé un outil qui est à la fois plus précis et plus efficace. Pour les industries qui dépendent de la compréhension de la manière dont les fluides se déplacent à travers la terre, de la production d'énergie à la gestion environnementale, cela signifie que les simulations peuvent être plus fiables et s'exécuter plus rapidement. La méthode ne nécessite pas la puissance de calcul massive qui était auparavant nécessaire pour obtenir un résultat décent, rendant la modélisation de haute fidélité de systèmes complexes hautement accessibles. Ces résultats suggèrent qu'en observant attentivement les points singuliers où le monde change, nous pouvons construire une meilleure compréhension de l'ensemble.