⚡ electrical engineering

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

Este estudio presenta un marco unificado basado en la física para modelar la presión capilar de fractura y la permeabilidad relativa utilizando un único parámetro de apertura interpretable, eliminando la necesidad de ajustes empíricos al tiempo que ofrece una alternativa robusta y lista para simuladores frente a los supuestos tradicionalmente simplificados y los costosos experimentos.

Autores originales: Alireza Golsefatan, Seyed Reza Shadizadeh, Hassan Shokrollahzadeh Behbahani

Publicado 2026-08-28

📖 5 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Alireza Golsefatan, Seyed Reza Shadizadeh, Hassan Shokrollahzadeh Behbahani

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Esta es una explicación generada por IA del artículo a continuación. No ha sido escrita ni avalada por los autores. Para mayor precisión técnica, consulte el artículo original. Leer descargo de responsabilidad completo

Bajo la superficie de la Tierra, vastos depósitos de petróleo y gas a menudo residen no en arenisca porosa como una esponja mojada, sino atrapados dentro de una compleja red de grietas y fisuras. Estas rocas fracturadas presentan un desafío único para los ingenieros que intentan predecir cómo se mueven los fluidos a través de ellas. En el mundo de la ingeniería de yacimientos, dos fuerzas invisibles dictan el flujo de petróleo y agua: la presión capilar y la permeabilidad relativa. La presión capilar es la fuerza que mantiene los fluidos en su lugar dentro de espacios diminutos, mientras que la permeabilidad relativa describe con qué facilidad fluye un fluido cuando otro ya está presente. Durante décadas, el enfoque estándar para modelar estos sistemas fracturados se ha basado en un supuesto conveniente pero erróneo: que estas grietas son tan anchas y uniformes que la presión capilar es efectivamente cero, y que los fluidos se deslizan uno al lado del otro de una manera perfectamente lineal y predecible. Esta simplificación ha facilitado los cálculos, pero la evidencia reciente sugiere que es físicamente incorrecta, lo que conduce a errores significativos al predecir cuánto petróleo se puede recuperar y con qué rapidez fluirá.

Un equipo de investigadores de la Universidad Tecnológica del Petróleo de Irán ha desarrollado una nueva forma unificada de describir el flujo de fluidos en estas rocas fracturadas que abandona estos viejos atajos. En lugar de adivinar o ajustar curvas a datos experimentales con números arbitrarios, construyeron un modelo que parte de la forma física real de las grietas. Reconocieron que ninguna fractura es perfectamente lisa o uniforme; el espacio entre las paredes de la roca, conocido como apertura, varía drásticamente de un lugar a otro. Al tratar esta variación como un patrón estadístico específico, crearon un marco que vincula la geometría de la grieta directamente con el comportamiento de los fluidos en su interior. Su trabajo proporciona un método claro, basado en la física, para calcular tanto la presión capilar como la permeabilidad relativa de una fractura utilizando solo una medida clave: la desviación estándar de la apertura de la fractura. Este único número, que describe cuánto varía el ancho de la grieta, reemplaza la necesidad de múltiples y confusos parámetros de ajuste que han plagado los modelos anteriores.

Los investigadores comenzaron conceptualizando una fractura no como un simple espacio plano, sino como una red donde el tamaño de la abertura cambia constantemente a lo largo de su longitud y anchura. Determinaron que estas variaciones siguen una distribución estadística específica, similar a cómo varían las estaturas en una gran población. Utilizando esta distribución, derivaron relaciones matemáticas que muestran exactamente cómo el agua y el petróleo se distribuyen dentro de la grieta. Las aberturas más pequeñas tienden a retener el fluido mojante, como el agua, mientras que las aberturas más grandes permiten el paso del fluido no mojante, como el petróleo. Al integrar estas realidades físicas, generaron un conjunto de "curvas tipo". Estas curvas actúan como un mapa universal, permitiendo a los ingenieros predecir el comportamiento de una fractura simplemente conociendo cuánto varía su anchura. Si la anchura es muy constante, el flujo se comporta de una manera; si la anchura varía significamente, el flujo se comporta de otra. Este enfoque elimina la necesidad de las antiguas y simplificadas curvas en "forma de X" que asumían que los fluidos se mueven en línea recta independientemente de la complejidad de la grieta.

Para demostrar que su modelo funcionaba, el equipo lo probó contra una amplia gama de datos existentes, incluyendo tanto estudios teóricos como experimentos de laboratorio del mundo real. Compararon sus predicciones con datos de docenas de estudios diferentes, que cubrían desde fracturas de laboratorio lisas hasta superficies de roca rugosas y naturales. En los casos donde las propiedades físicas de la fractura eran conocidas, su modelo coincidió con los datos observados con una precisión notable. En los casos donde las propiedades eran desconocidas, utilizaron una técnica de ajuste para encontrar el mejor parámetro de variación, el cual predijo con éxito el comportamiento del fluido bajo diferentes condiciones. Los resultados fueron consistentes: el nuevo modelo podía replicar datos experimentales complejos con alta precisión, superando a menudo los métodos antiguos que dependían de múltiples constantes empíricas. Crucialmente, el modelo también captura cómo la fractura cambia bajo presión. A medida que la roca es comprimida por el peso de la tierra arriba, la fractura se cierra y la variación en su anchura cambia. Los investigadores demostraron que su modelo puede rastrear estos cambios mecánicos, ajustando las predicciones de flujo para reflejar la nueva y más estrecha geometría de la grieta.

Las implicaciones de este trabajo se extienden más allá de una simple mejora en el cálculo. Al proporcionar un modelo que depende de un único parámetro físicamente significativo en lugar de una colección de números de ajuste, los investigadores han ofrecido una herramienta que es tanto más simple como más robusta. Esto significa que los ingenieros ahora pueden integrar estas descripciones de flujo más realistas directamente en las simulaciones computacionales a gran escala utilizadas para gestionar campos petroleros. En lugar de depender de la hipótesis obsoleta de que las fracturas son canales abiertos y vacíos sin resistencia interna, ahora pueden simular la intrincada danza de los fluidos dentro de la geometría real e irregular de la roca. El estudio confirma que la presión capilar de la fractura no es cero y que la permeabilidad relativa no es una simple línea recta; son fenómenos complejos, dependientes de la geometría, que deben tenerse en cuenta para comprender el verdadero potencial de un yacimiento. Aunque el modelo fue desarrollado para fracturas individuales, los autores sugieren que este marco podría expandirse para manejar redes completas de grietas e incluso flujos de tres fases más complejos, ofreciendo un camino más claro para comprender la dinámica oculta de la corteza fracturada de la Tierra.

Resumen Técnico: Curvas de Tipo de Presión Capilar y Permeabilidad Relativa de Fracturas Físicamente Consistentes

Planteamiento del Problema
El modelado preciso del flujo de fluidos en yacimientos fracturados sigue siendo un desafío significativo debido a la compleja geometría de las redes de fracturas y las intrincadas interacciones matriz-fractura. Si bien la presión capilar de fractura ( $P_c$ ) y la permeabilidad relativa ( $k_r$ ) son las relaciones constitutivas fundamentales que gobiernan la interferencia de fases y la eficiencia de desplazamiento, las prácticas tradicionales de ingeniería de yacimientos suelen basarse en supuestos excesivamente simplificados. Específicamente, los simuladores frecuentemente asumen una presión capilar nula y un comportamiento de permeabilidad relativa lineal (curva en X) para las fracturas. La evidencia experimental y de modelado reciente sugiere que estos supuestos carecen de una base física sólida, lo que puede conducir a errores significativos en la predicción del desempeño. Además, los modelos existentes suelen depender de múltiples "parámetros de ajuste" empíricos que carecen de un significado físico claro o que son específicos para cada caso, lo que limita su capacidad predictiva en diferentes geometrías de fractura. Existe una brecha crítica en el vínculo directo entre las características microestructurales de la fractura (específicamente la distribución de la apertura) y estos parámetros de flujo.

Metodología
Este estudio desarrolla un marco unificado y basado en la física para los modelos de $P_c$ y $k_r$ de fractura, integrando las distribuciones de apertura de las fracturas con la mecánica de flujo fundamental. La metodología central consiste en los siguientes pasos:

Caracterización Estadística: El sistema de fractura se conceptualiza como un volumen unitario que contiene una red totalmente conectada donde las fracturas individuales exhiben aperturas espacialmente variables. El estudio emplea una distribución lognormal para caracterizar la apertura, el ancho y la longitud de la fractura, y una distribución normal para la orientación de la fractura. La desviación estándar de la apertura ( $\sigma_b$ ) se identifica como el parámetro gobernante principal.
Derivación de la Saturación de Fase: Las saturaciones móviles e inmóviles para las fases húmeda y no húmeda se derivan calculando los volúmenes de fase dentro de rangos de apertura específicos definidos por umbrales de presión capilar. Estos rangos corresponden a la fase húmeda irreducible, la fase no húmeda residual y las saturaciones totales de fase.
Desarrollo de Curvas de Tipo:
- Presión Capilar: La saturación adimensional de la fase húmeda se formula como una función de la saturación efectiva y la presión capilar adimensional ( $P_{cDf}$ ). Al incorporar la distribución de apertura lognormal, se derivan las curvas de tipo de Presión Capilar de Drenaje e Imbibición utilizando funciones de error (erf). Estas curvas dependen de un único parámetro físico, $\sigma_b$ , y de presiones límite adimensionales ( $P_{cDminf}$ , $P_{cDmaxf}$ ).
- Permeabilidad Relativa: El modelo deriva la permeabilidad relativa de la fractura integrando la ecuación de Hagen-Poiseuille para el flujo de fluidos en conductos con la ley de Darcy para medios porosos. La integración tiene en cuenta la distribución de la apertura y las saturaciones de fase a través de los procesos de drenaje e imbibición. Se incorporan coeficientes de tortuosidad, con exponentes determinados mediante la validación contra datos experimentales.
Validación y Ajuste: Los modelos se validan contra un amplio espectro de marcos teóricos (por ejemplo, Pruess y Tsang, Zhou y Wheater, Liu et al.) y conjuntos de datos experimentales (por ejemplo, Reitsma y Kueper, Huo et al., Chen, Diomampo).
- Cuando $\sigma_b$ es conocido, el modelo genera curvas directamente para comparación.
- Cuando $\sigma_b$ es desconocido, se emplea un procedimiento de ajuste de curvas de tipo. Los datos experimentales de $P_c$ vs. saturación se grafican en una escala log-log y se ajustan contra las curvas de tipo adimensionales desarrolladas para extraer el $\sigma_b$ óptimo y otros parámetros.
Métricas Cuantitativas: La precisión del modelo se evalúa mediante el Error Cuadrático Medio (MSE) y el Error Absoluto Relativo Promedio (AARE).

Contribuciones Clave

Modelo Unificado Basado en la Apertura: El estudio presenta un marco único que deriva tanto la presión capilar como la permeabilidad relativa de la fractura directamente de la distribución de la apertura de la fractura, eliminando la necesidad de coeficientes de ajuste empíricos arbitrarios.
Parámetro Físico Único: El modelo depende de la desviación estándar de la apertura de la fractura ( $\sigma_b$ ) como el único parámetro de entrada. Este parámetro es físicamente interpretable y representa la heterogeneidad geométrica de la fractura, contrastando con los modelos tradicionales que requieren múltiples parámetros de ajuste específicos para cada caso.
Metodología de Curvas de Tipo: La introducción de curvas de tipo de presión capilar de fractura permite una aplicación dual: (1) predecir $P_c$ y $k_r$ a partir de distribuciones de apertura conocidas, y (2) caracterizar las distribuciones de apertura de las fracturas a partir de datos experimentales de $P_c$ mediante el ajuste de curvas de tipo.
Simplicidad Matemática y Robustez: Las ecuaciones derivadas utilizan funciones estadísticas estándar (funciones de error) y están diseñadas para una integración fluida en simuladores de yacimientos comerciales.

Resultados
Los modelos desarrollados demostraron una robustez y eficacia excepcionales a través de diversos tipos de fractura y condiciones de flujo:

Presión Capilar: El modelo replicó con éxito los modelos teóricos (por ejemplo, Zhou y Wheater) con una alineación casi perfecta (MSE en el rango de $10^{-9}$ a $10^{-10}$ ) y coincidió con los datos experimentales (por ejemplo, Huo et al., Reitsma y Kueper) con valores bajos de AARE (típicamente < 20%). Notablemente, el modelo capturó la tendencia física donde el aumento del esfuerzo efectivo conduce a un aumento en la $\sigma_b$ optimizada, reflejando el cierre de la fractura y el aumento de la heterogeneidad.
Permeabilidad Relativa: El modelo predijo con precisión las curvas de permeabilidad relativa tanto para los procesos de drenaje como de imbibición. Mostró una fuerte concordancia con los límites teóricos (por ejemplo, Brooks y Corey con $\lambda \to \infty$ ) y con conjuntos de datos experimentales (por ejemplo, Chen, Diomampo). El modelo logró desacoplar la interferencia de fases de las variables experimentales secundarias, proporcionando una validación más rigurosa que las gráficas tradicionales basadas en la saturación.
Eficiencia de Parámetros: En todos los casos de validación, el modelo alcanzó una precisión comparable o superior a la de los modelos complejos existentes utilizando únicamente el parámetro único $\sigma_b$ .

Significancia y Reivindicaciones
El artículo afirma que el marco propuesto ofrece una alternativa confiable y rentable a los costosos y laboriosos experimentos de laboratorio para determinar la presión capilar y la permeabilidad relativa de las fracturas. Al fundamentar los modelos en la realidad física de las distribuciones de apertura de las fracturas en lugar de en el ajuste de curvas empíricas, el estudio proporciona una herramienta físicamente consistente para los ingenieros de yacimientos.

Los autores enfatizan que la estructura matemática de los modelos facilita su implementación directa en simuladores de yacimientos comerciales, permitiendo simulaciones a gran escala más realistas de sistemas fracturados. El estudio concluye que la presión capilar de la fractura no es ni cero ni constante, y que la permeabilidad relativa de la fractura se desvía significativamente de las "curvas en X" lineales; estos fenómenos deben contabilizarse utilizando modelos basados en la apertura para capturar con precisión la interferencia de fases. El trabajo cierra la brecha entre las soluciones analíticas complejas y las aplicaciones prácticas de la ingeniería de yacimientos, ofreciendo un método estandarizado para caracterizar medios fracturados sin la ambigüedad de los parámetros de ajuste arbitrarios.

Más como este