Residual Saturation under Pressure-Controlled Drainage

Résumé Technique : Saturation Résiduelle sous Drainage à Pression Contrôlée

Énoncé du Problème
L'article traite d'une question fondamentale non résolue en physique des milieux poreux : le comportement de la saturation résiduelle lors d'un drainage à pression contrôlée dans la limite de taille infinie. Alors que les modèles standards de percolation par invasion (IP) prédisent que la saturation de la phase envahissante tend vers zéro à mesure que la taille du système tend vers l'infini (laissant une fraction finie de fluide défendant piégée uniquement en raison de la nature fractale du cluster de percée), il reste incertain si un processus à pression contrôlée — qui permet à la phase envahissante d'accéder à tous les pores connectés à l'entrée à une pression capillaire donnée — résulte en une saturation résiduelle finie de la phase défendante dans la limite thermodynamique. Cette distinction est cruciale pour comprendre les phénomènes à l'échelle du réservoir, tels que l'injection de l'eau, où la question se pose de savoir si la saturation résiduelle en huile reste finie lorsque la région envahie approche les dimensions du réservoir.

Méthodologie
Les auteurs formulent le drainage à pression contrôlée comme une Percolation de Liens avec Piégeage (BPT - Bond Percolation with Trapping) sur un graphe de réseau de pores.

• Représentation par Graphe : Le milieu poreux est modélisé comme un graphe non orienté où les nœuds représentent des régions de pores et les arêtes représentent des cols avec des rayons d'entrée assignés (seuils capillaires).
• Définition du Processus : Contra-irement à l'IP standard, qui progresse chemin par chemin via le chemin localement le plus accessible, la BPT opère sous une contrainte de pression globale. À mesure que la pression capillaire augmente, le « graphe actif » (comprenant les cols dont les seuils d'entrée sont inférieurs à la pression actuelle) s'étend.
• Mécanisme de Piégeage : L'invasion est définie par la connectivité à l'entrée sur le graphe actif. Le piégeage se produit lorsqu'une partie de la phase défendante perd sa connectivité avec l'ensemble de sortie. L'algorithme met à jour le système en : (1) augmentant la pression pour activer de nouveaux cols ; (2) identifiant la région envahie (pores actifs connectés à l'entrée) ; (3) identifiant la région défendue (pores non envahis connectés à la sortie) ; et (4) classant tous les pores restants déconnectés comme « piégés ».
• Simulations : L'étude compare la BPT à la percolation par invasion (IP) standard sur divers réseaux 2D (carré, triangulaire, Voronoi) et réseaux 3D (diamant, cubique simple avec des nombres de coordination variables). Les simulations utilisent des rayons de cols distribués selon une loi log-normale et des moyennes d'ensemble pour analyser les effets de taille finie.

Contributions Clés et Résultats

1. Saturation Résiduelle Finie dans la Limite Infinie :
   La conclusion principale est que le drainage à pression contrôlée conduit à un état envahi qui remplit l'espace, où la saturation résiduelle de la phase défendante tend vers une limite finie et non nulle ($0 < S_\infty < 1$) lorsque $L \to \infty$. Cela contraste nettement avec l'IP standard, où la saturation envahie tend vers zéro à la percée ($S_{IP} \to 0$) car le cluster envahisseur est fractal et contourne de larges régions du fluide défendant.

2. Mise à l'Échelle de Taille Finie en Deux Dimensions :
   En 2D, l'écart de la saturation résiduelle par rapport à sa limite thermodynamique suit une loi de mise à l'échelle de taille finie :
   $$S_r(L) \approx S_\infty - \left(\frac{L_0}{L}\right)^\delta$$
   où l'exposant $\delta \approx 0,25$. Cet exposant est universel pour différents types de réseaux (carré, triangulaire, Voronoi) et est indépendant des détails microscopiques. La convergence lente ($\Delta S_r \sim L^{-0,25}$) implique que l'obtention de volumes élémentaires représentatifs en 2D nécessite des tailles de système extrêmement grandes (par exemple, doubler la taille du système ne réduit le biais de taille finie que d'un facteur $\approx 1/16$).

3. Effets de Dimensionnalité et de Nombre de Coordination en Trois Dimensions :
   En 3D, les corrections de taille finie décroissent nettement plus vite, avec un exposant $\delta \approx 0,75$. Par conséquent, le volume représentatif requis pour approcher la limite de taille infinie est beaucoup plus petit en 3D qu'en 2D. De plus, bien que la saturation résiduelle asymptotique reste finie en 3D, sa valeur dépend du nombre de coordination du réseau ; une connectivité plus élevée offre davantage de voies d'invasion alternatives, réduisant la fraction piégée, sans toutefois l'éliminer totalement.

4. Distinction avec la Percolation par Invasion :
   L'article clarifie que la différence entre la BPT et l'IP n'est pas simplement une question de forme de courbe ou de dynamique de rupture (burst), mais une différence fondamentale d'états asymptotiques. L'IP produit un cluster envahisseur fractal lors de la percée, tandis que la BPT produit une région envahie qui remplit l'espace avec une fraction piégée finie coexistant dans la limite infinie.

Signification
Cet article établit un lien théorique direct entre la théorie de la percolation et les relations pression-saturation en démontant que la connectivité et le piégeage sont les facteurs déterminants de la saturation résiduelle sous drainage à pression contrôlée. En montrant qu'une fraction finie de l'espace poreux reste piégée même dans la limite de système infini, ce travail étend le modèle standard de percolation par invasion au-delà de l'état de percée. Cela fournit un cadre rigoureux pour interpréter les lois constitutives macroscopiques et suggère que le problème du « volume élémentaire représentatif » est très sensible à la dimensionnalité, étant bien plus restrictif en 2D qu'en 3D. Les résultats impliquent que, dans les scénarios à pression contrôlée, une fraction finie du fluide défendant est irrécupérable quel que soit la taille du système, une conclusion ayant des implications directes pour la redistribution des fluides et les processus de récupération dans le sous-sol.