Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

Mucho más allá de la superficie de la Tierra, en las rocas porosas que contienen petróleo, gas y calor geotérmico, se libra constantemente una batalla silenciosa entre la piedra sólida y los fluidos atrapados en sus diminutos poros. Cuando los ingenieros intentan extraer estos recursos, a menudo bombean fluidos hacia la roca para separarla, un proceso conocido como fracturación hidráulica. Durante décadas, los científicos han modelado este proceso tratando a la roca como una simple esponja elástica que se estira y vuelve a su forma original. Sin embargo, las rocas reales son más complejas; también pueden doblarse permanentemente, desmoronarse y deslizarse a lo largo de grietas internas, especialmente cuando son comprimidas por altas presiones. Comprender exactamente cómo se rompen estas rocas cuando tanto el esqueleto sólido como el fluido en su interior empujan y tiran es crucial para perforar pozos de forma segura, gestionar el almacenamiento de residuos e incluso predecir eventos geológicos naturales como el agrietamiento de icebergs o el paso del magma a través de la corteza. El desafío radica en crear una descripción matemática que capture la realidad desordenada de la ruptura de la roca sin perderse en cálculos imposibles.

Un equipo de investigadores ha revisado las ecuaciones fundamentales utilizadas para describir este proceso de ruptura, descubriendo que un atajo común utilizado en modelos anteriores conduce a un error significativo. En el mundo de la mecánica de rocas, los científicos utilizan un enfoque de "campo de fase", que trata una grieta no como una línea afilada e invisible, sino como una zona difusa donde el material se debilita gradualmente. Esto permite que las computadoras simulen cómo crece y cambia de forma una grieta. Para hacer esto con precisión, el modelo debe dar cuenta de dos comportamientos distintos: cuando la roca es estirada (tensión) y cuando es comprimida (compresión). En tensión, las diminutas grietas se abren; en compresión, se cierran y las superficies rugosas en su interior se deslizan entre sí, creando fricción. Los investigadores descubrieron que muchos modelos existentes, que se utilizan ampliamente en la industria y la academia, confunden las variables utilizadas para describir la energía del sistema. Específicamente, encontraron que tratar la presión del fluido como una entrada directa en la ecuación de energía, en lugar de la cantidad de fluido en los poros, crea un "salto" o discontinuidad matemática. Este salto significa que el modelo predice que la resistencia de la roca cambia abruptamente y de manera poco realista al cambiar de tracción a compresión, lo que conduce a predicciones inexactas de cómo y dónde se formará una fractura.

Para solucionar esto, los investigadores desarrollaron una nueva formulación consistente que rastrea cuidadosamente la relación entre la deformación de la roca y el contenido de fluido. Construyeron su modelo basándose en una visión microscópica de la roca, imaginándola como una colección de diminutas grietas en forma de moneda llenas de fluido. Al analizar cómo estas microgrietas se abren y cierran bajo diferentes tensiones, derivaron un conjunto de reglas que aseguran que la superficie de resistencia de la roca —el límite matemático que define cuándo se romperá— permanezca suave y continua, independientemente de cómo cambie la presión. Demostraron que, al utilizar un tipo específico de regla de flujo que vincula la deformación permanente de la roca directamente con las fuerzas que actúan sobre ella, e incluir correctamente un término de acoplamiento faltante en la ecuación de energía, el modelo puede describir con precisión la transición de la tensión a la compresión. Esta corrección no es solo un ajuste menor; cambia fundamentalmente la forma en que el modelo calcula la fuerza impulsora que empuja una grieta hacia adelante.

El equipo probó su nuevo modelo contra dos escenarios muy diferentes para ver si se mantenía firme. Primero, simularon un benchmark clásico de fracturación hidráulica donde se inyecta fluido en una roca para crear una grieta. Compararon sus resultados con una solución analítica conocida, que es una respuesta matemática precisa derivada para condiciones ideales. El nuevo modelo coincidió con esta solución mucho más de cerca que las formulaciones mixtas anteriores, prediciendo correctamente la presión necesaria para iniciar la fractura y el ancho de la grieta a medida que crecía. Los modelos anteriores tendían a sobreestimar la presión del fluido necesaria, un fallo que podría conducir a operaciones de perforación ineficientes o inseguras. En la segunda prueba, simularon una muestra de roca siendo comprimida desde dos lados, una condición que imita el entorno profundo bajo tierra donde las rocas están bajo una presión inmensa. Aquí, la roca no solo se rompió; se deformó plásticamente, formando una banda de cizalla donde el material se deslizó sobre sí mismo antes de romperse. El nuevo modelo reprodujo con éxito este comportamiento complejo, capturando tanto las fracturas de cizalla causadas por la compresión mecánica como las fracturas de tensión impulsadas por la inyección de fluido.

Los resultados de estas simulaciones revelan que la deformación plástica de la roca actúa como una barrera significativa para el crecimiento de la fractura. Cuando se permite que la roca se deforme permanentemente, absorbe más energía, lo que significa que la presión del fluido dentro de la grieta debe ser mayor para mantener la fractura abierta y en crecimiento. Esta resistencia adicional cambia todo el perfil de la fractura, haciéndola más estrecha y requiriendo más fuerza para propagarse de lo que sugeriría un modelo puramente elástico. Además, en las pruebas de compresión, el modelo mostró cómo cambia la presión de los poros dentro de la roca a medida que la fractura se forma. A medida que la roca se deforma y crea espacio para que el fluido se mueva, la presión cae localmente, atrayendo más fluido hacia la zona dañada. Esta interacción crea un bucle de retroalimentación donde el flujo de fluido y la falla mecánica de la roca están estrechamente acoplados, una dinámica que el nuevo modelo captura con alta fidelidad.

Al corregir la base matemática de estos modelos, los investigadores han proporcionado una herramienta más confiable para ingenieros y científicos que trabajan con medios porosos. La nueva formulación asegura que la transición entre los diferentes tipos de falla de la roca sea suave y físicamente realista, eliminando los saltos artificiales que plagaban los intentos anteriores. Esta precisión es vital para aplicaciones que van desde la optimización de la producción de petróleo y gas hasta el diseño de un almacenamiento subterráneo seguro para la captura de carbono. El estudio demuestra que incluso en sistemas complejos donde el sólido y el fluido interactúan, un enfoque cuidadoso y consistente de la física subyacente puede revelar la verdadera naturaleza de cómo se rompen los materiales. El modelo ahora está listo para ayudar a predecir el comportamiento de la corteza terrestre con mayor confianza, cerrando la brecha entre la mecánica teórica y la realidad plástica y desordenada del suelo bajo nuestros pies.