⚡ electrical engineering

Physically-Consistent Fracture Capillary Pressure Type Curves and Relative Permeability: A Unified Aperture-Based Model without Empirical Matching

Cette étude présente un cadre physique unifié pour modéliser la pression capillaire de fracture et la perméabilité relative à l'aide d'un paramètre d'ouverture unique et interprétable, éliminant ainsi le besoin d'ajustements empiriques tout en offrant une alternative robuste et prête pour la simulation aux hypothèses traditionnelles simplistes et aux expériences coûteuses.

Auteurs originaux : Alireza Golsefatan, Seyed Reza Shadizadeh, Hassan Shokrollahzadeh Behbahani

Publié 2026-08-28

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Alireza Golsefatan, Seyed Reza Shadizadeh, Hassan Shokrollahzadeh Behbahani

Article original sous licence CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Profondément sous la surface de la Terre, de vastes réservoirs de pétrole et de gaz résident souvent non pas dans un grès poreux semblable à une éponge humide, mais sont piégés au sein d'un réseau complexe de fissures et de fractures. Ces roches fracturées présentent un défi unique pour les ingénieurs tentant de prédire comment les fluides s'y déplacent. Dans le monde de l'ingénierie des réservoirs, deux forces invisibles dictent l'écoulement du pétrole et de l'eau : la pression capillaire et la perméabilité relative. La pression capillaire est la force qui maintient les fluides en place dans de minuscules espaces, tandis que la perméabilité relative décrit la facilité avec laquelle un fluide peut s'écouler lorsqu'un autre est déjà présent. Pendant des décennies, l'approche standard pour modéliser ces systèmes fracturés s'est appuyée sur une hypothèse pratique mais erronée : que ces fissures sont si larges et uniformes que la pression capillaire est effectivement nulle, et que les fluides glissent les uns sur les autres de manière parfaitement linéaire et prévisible. Cette simplification a facilité les calculs, mais des preuves récentes suggèrent qu'elle est physiquement fausse, entraînant des erreurs significatives dans la prédiction de la quantité de pétrole pouvant être récupérée et de sa vitesse d'écoulement.

Une équipe de chercheurs de l'Université de Technologie du Pétrole en Iran a développé une nouvelle façon unifiée de décrire l'écoulement des fluides dans ces roches fracturées, en abandonnant ces anciens raccourcis. Au lieu de deviner ou d'ajuster des courbes à des données expérimentales avec des nombres arbitraires, ils ont construit un modèle qui part de la forme physique réelle des fissures. Ils ont reconnu qu'aucune fracture n'est parfaitement lisse ou uniforme ; l'écart entre les parois de la roche, appelé ouverture, varie considérablement d'un endroit à l'autre. En traitant cette variation comme un motif statistique spécifique, ils ont créé un cadre qui lie directement la géométrie de la fissure à la manière dont les fluides se comportent à l'intérieur. Leur travail fournit une méthode claire, basée sur la physique, pour calculer à la fois la pression capillaire et la perméabilité relative d'une fracture en utilisant une seule mesure clé : l'écart type de l'ouverture de la fracture. Ce chiffre unique, qui décrit à quel point la largeur de la fissure varie, remplace le besoin de multiples paramètres d'ajustement déroutants qui ont tourmenté les modèles précédents.

Les chercheurs ont commencé par conceptualiser une fracture non pas comme un simple interstice plat, mais comme un réseau où la taille de l'ouverture change constamment le long de sa longueur et de sa largeur. Ils ont déterminé que ces variations suivent une distribution statistique spécifique, similaire à la façon dont les tailles varient dans une grande population. En utilisant cette distribution, ils ont dérivé des relations mathématiques montrant exactement comment l'eau et le pétrole se distribuent à l'intérieur de la fissure. Les ouvertures plus petites ont tendance à retenir le fluide mouillant, comme l'eau, tandis que les ouvertures plus larges permettent au fluide non mouillant, comme le pétrole, de passer. En intégrant ces réalités physiques, ils ont généré un ensemble de « courbes de type ». Ces courbes agissent comme une carte universelle, permettant aux ingénieurs de prédire le comportement d'une fracture simplement en sachant à quel point sa largeur varie. Si la largeur est très constante, l'écoulement se comporte d'une certaine manière ; si la largeur varie considérablement, l'écoulement se comporte différemment. Cette approche élimine le besoin des anciennes courbes en « forme de X » simplistes qui supposaient que les fluides se déplacent en ligne droite indépendamment de la complexité de la fissure.

Pour prouver que leur modèle fonctionnait, l'équipe l'a testé contre un large éventail de données existantes, incluant à la fois des études théoriques et des expériences réelles en laboratoire. Ils ont comparé leurs prédictions aux données de dizaines d'études différentes, couvrant tout, des fractures lisses de laboratoire aux surfaces rocheuses naturelles et rugueuses. Dans les cas où les propriétés physiques de la fracture étaient connues, leur modèle correspondait aux données observées avec une précision remarquable. Dans les cas où les propriétés étaient inconnues, ils ont utilisé une technique d'ajustement pour trouver le meilleur paramètre de variation, qui a ensuite prédit avec succès le comportement des fluides sous différentes conditions. Les résultats étaient cohérents : le nouveau modèle pouvait reproduire des données expérimentales complexes avec une grande précision, surpassant souvent les anciennes méthodes qui reposaient sur de multiples constantes empiriques. Crucialement, le modèle capture également la façon dont la fracture change sous la pression. À mesure que la roche est comprimée par le poids de la terre au-dessus, la fracture se referme et la variation de sa largeur change. Les chercheurs ont montré que leur modèle peut suivre ces changements mécaniques, ajustant les prédictions d'écoulement pour refléter la nouvelle géométrie plus étroite de la fissure.

Les implications de ce travail s'étendent au-delà d'une simple amélioration de calcul. En fournissant un modèle qui repose sur un paramètre unique et physiquement significatif plutôt que sur une collection de chiffres d'ajustement, les chercheurs ont offert un outil qui est à la fois plus simple et plus robuste. Cela signifie que les ingénieurs peuvent désormais intégrer ces descriptions de flux plus réalistes directement dans les simulations informatiques à grande échelle utilisées pour gérer les champs pétroliers. Au lieu de s'appuyer sur l'hypothèse obsolète de canaux ouverts et vides sans résistance interne, ils peuvent désormais simuler la danse complexe des fluides au sein de la géométrie réelle et inégale de la roche. L'étude confirme que la pression capillaire des fractures n'est pas nulle et que la perméabilité relative n'est pas une simple ligne droite ; ce sont des phénomènes complexes, dépendants de la géométrie, qui doivent être pris en compte pour comprendre le véritable potentiel d'un réservoir. Bien que le modèle ait été développé pour des fractures uniques, les auteurs suggèrent que ce cadre pourrait être étendu pour gérer des réseaux entiers de fissures et même des écoulements à trois phases plus complexes, offrant une voie plus claire vers la compréhension de la dynamique cachée de la croûte terrestre fracturée.

Résumé technique : Courbes types de pression capillaire et de perméabilité relative physiquement cohérentes pour les fractures

Énoncé du problème
La modélisation précise de l'écoulement des fluides dans les réservoirs fracturés demeure un défi important en raison de la géométrie complexe des réseaux de fractures et des interactions complexes matrice-fracture. Bien que la pression capillaire ( $P_c$ ) et la perméabilité relative ( $k_r$ ) des fractures soient les relations constitutives fondamentales régissant l'interférence de phases et l'efficacité de déplacement, les pratiques traditionnelles de l'ingénierie des réservoirs reposent souvent sur des hypothèses simplistes. Plus précisément, les simulateurs supposent fréquemment une pression capillaire nulle et un comportement de perméabilité relative linéaire (courbe en X) pour les fractures. Des preuves expérimentales et de modélisation récentes suggèrent que ces hypothèses manquent d'un fondement physique solide, ce qui peut conduire à des erreurs significatives de prédiction de performance. De plus, les modèles existants reposent souvent sur de multiples « paramètres d'ajustement » empiriques qui manquent de signification physique claire ou sont spécifiques à chaque cas, limitant ainsi leur capacité prédictive à travers différentes géométries de fractures. Un fossé critique existe dans le lien direct entre les caractéristiques microstructurales des fractures (spécifiquement la distribution des ouvertures) et ces paramètres d'écoulement.

Méthodologie
Cette étude développe un cadre unifié, fondé sur la physique, pour les modèles de $P_c$ et de $k_r$ des fractures en intégrant les distributions d'ouverture des fractures avec les mécanismes d'écoulement fondamentaux. La méthodologie centrale implique les étapes suivantes :

Caractérisation statistique : Le système de fractures est conceptualisé comme un volume unitaire contenant un réseau entièrement connecté où les fractures individuelles présentent des ouvertures spatialement variables. L'étude utilise une distribution log-normale pour caractériser l'ouverture, la largeur et la longueur des fractures, et une distribution normale pour l'orientation des fractures. L'écart-type de l'ouverture ( $\sigma_b$ ) est identifié comme le paramètre directeur principal.
Dérivation de la saturation de phase : Les saturations mobiles et immobiles pour les phases mouillantes et non-mouillantes sont dérivées en calculant les volumes de phase dans des plages d'ouverture spécifiques définies par des seuils de pression capillaire. Ces plages correspondent à la phase mouillante irréductible, à la phase non-mouillante résiduelle et aux saturations totales de phase.
Développement de courbes types :
- Pression capillaire : La saturation adimensionnelle de la phase mouillante est formulée en fonction de la saturation effective et de la pression capillaire adimensionnelle ( $P_{cDf}$ ). En incorporant la distribution log-normale des ouvertures, des courbes types de pression capillaire de drainage et d'imbibition des fractures sont dérivées en utilisant des fonctions d'erreur (erf). Ces courbes dépendent d'un seul paramètre physique, $\sigma_b$ , et de pressions limites adimensionnelles ( $P_{cDminf}$ , $P_{cDmaxf}$ ).
- Perméabilité relative : Le modèle dérive la perméabilité relative des fractures en intégrant l'équation de Hagen-Poiseuille pour l'écoulement des fluides dans les conduits avec la loi de Darcy pour les milieux poreux. L'intégration tient compte de la distribution des ouvertures et des saturations de phase au cours des processus de drainage et d'imbibition. Des coefficients de tortuosité sont incorporés, avec des exposants déterminés par validation par rapport aux données expérimentales.
Validation et mise en correspondance : Les modèles sont validés par rapport à un large spectre de cadres théoriques (ex. : Pruess et Tsang, Zhou et Wheater, Liu et al.) et d'ensembles de données expérimentales (ex. : Reitsma et Kueper, Huo et al., Chen, Diomampo).
- Lorsque $\sigma_b$ est connu, le modèle génère directement des courbes pour comparaison.
- Lorsque $\sigma_b$ est inconnu, une procédure de mise en correspondance de courbes types est employée. Les données expérimentales de $P_c$ vs saturation sont tracées sur une échelle log-log et mises en correspondance avec les courbes types adimensionnelles développées pour extraire le $\sigma_b$ optimal et les autres paramètres.
Métriques quantitatives : La précision du modèle est évaluée à l'aide de l'erreur quadratique moyenne (MSE) et de l'erreur absolue relative moyenne (AARE).

Contributions clés

Modèle basé sur l'ouverture unifié : L'étude présente un cadre unique qui dérive à la fois la pression capillaire et la perméabilité relative des fractures directement de la distribution d'ouverture des fractures, éliminant le besoin de coefficients d'ajustement empiriques arbitraires.
Paramètre physique unique : Le modèle repose sur l'écart-type de l'ouverture de la fracture ( $\sigma_b$ ) comme seul paramètre d'entrée. Ce paramètre est physiquement interprétable et représente l'hétérogénéité géométrique de la fracture, contrastant avec les modèles traditionnels qui nécessitent plusieurs paramètres d'ajustement spécifiques à chaque cas.
Méthodologie de courbe type : L'introduction de courbes types de pression capillaire de fracture permet une double application : (1) prédire la $P_c$ et la $k_r$ à partir de distributions d'ouverture connues, et (2) caractériser les distributions d'ouverture des fractures à partir de données expérimentales de $P_c$ via la mise en correspondance de courbes types.
Simplicité mathématique et robustesse : Les équations dérivées utilisent des fonctions statistiques standards (fonctions d'erreur) et sont conçues pour une intégration fluide dans les simulateurs de réservoirs commerciaux.

Résultats
Les modèles développés ont démontré une robustesse et une efficacité exceptionnelles à travers divers types de fractures et conditions d'écoulement :

Pression capillaire : Le modèle a réussi à reproduire les modèles théoriques (ex. : Zhou et Wheater) avec un alignement quasi parfait (MSE dans la plage de $10^{-9}$ à $10^{-10}$ ) et a correspondé aux données expérimentales (ex. : Huo et al., Reitsma et Kueper) avec de faibles valeurs d'AARE (typiquement < 20 %). Notamment, le modèle a capturé la tendance physique où l'augmentation du stress effectif conduit à une augmentation du $\sigma_b$ optimisé, reflétant la fermeture des fractures et l'augmentation de l'hétérogénéité.
Perméabilité relative : Le modèle a prédit avec précision les courbes de perméabilité relative pour les processus de drainage et d'imbibition. Il a montré un fort accord avec les limites théoriques (ex. : Brooks et Corey avec $\lambda \to \infty$ ) et les ensembles de données expérimentales (ex. : Chen, Diomampo). Le modèle a réussi à découpler l'interférence de phase des variables expérimentales secondaires, fournissant une validation plus rigoureuse que les tracés traditionnels basés sur la saturation.
Efficacité des paramètres : Dans tous les cas de validation, le modèle a atteint une précision comparable ou supérieure aux modèles complexes existants tout en utilisant uniquement le paramètre unique $\sigma_b$ .

Signification et affirmations
L'article affirme que le cadre proposé offre une alternative fiable et rentable aux expériences de laboratoire coûteuses et chronophages pour déterminer la pression capillaire et la perméabilité relative des fractures. En ancrant les modèles dans la réalité physique des distributions d'ouverture des fractures plutôt que dans l'ajustement de courbes empiriques, l'étude fournit un outil physiquement cohérent pour les ingénieurs de réservoir.

Les auteurs soulignent que la structure mathématique des modèles facilite leur mise en œuvre directe dans les simulateurs de réservoirs commerciaux, permettant des simulations à grande échelle plus réalistes des systèmes fracturés. L'étude conclut que la pression capillaire des fractures n'est ni nulle ni constante, et que la perméabilité relative des fractures s'écarte considérablement des « courbes en X » linéaires ; ces phénomènes doivent être pris en compte à l'aide de modèles basés sur l'ouverture pour capturer précisément l'interférence de phase. Ce travail comble le fossé entre les solutions analytiques complexes et les applications pratiques de l'ingénierie des réservoirs, offrant une méthode standardisée pour caractériser les milieux fracturés sans l'ambiguïté des paramètres d'ajustement arbitraires.

Articles similaires