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

En lo profundo del subsuelo, en la compleja red de rocas que contienen petróleo, gas y agua, los fluidos no siempre se mueven en líneas simples y predecibles. La roca en sí rara vez es uniforme; es un mosaico de materiales con capacidades muy distintas para permitir el paso de fluidos. Algunas capas son como autopistas abiertas, mientras que otras son como muros densos. Cuando los ingenieros intentan modelar cómo fluyen estos fluidos, deben tener en cuenta estos cambios repentinos en las propiedades de la roca. Además, los fluidos mismos pueden ser complicados. Mientras que el agua fluye fácilmente, muchos fluidos industriales —como las pastas utilizadas en la fabricación de papel o los lodos de perforación utilizados en pozos petroleros— se comportan de manera diferente. Son "no newtonianos", lo que significa que su espesor y su resistencia al flujo cambian dependiendo de qué tan fuerte se les presione. Cuando estos fluidos complejos encuentran una roca que cambia su carácter abruptamente, las matemáticas que describen su movimiento se vuelven increíblemente difíciles, fallando a menudo precisamente en los puntos donde las propiedades de la roca cambian.

Durante décadas, los científicos han dependido de una forma estándar de calcular estos flujos, un método que esencialmente promedia las propiedades de la roca entre dos puntos. Este enfoque funciona lo suficientemente bien cuando la roca es relativamente consistente, pero falla drásticamente cuando la roca es altamente variada. En esos casos, el método estándar subestima sistemáticamente cuánto fluido puede realmente moverse, lo que genera errores que pueden ser masivos. Para resolver esto, investigadores de la Universidad de Ciencia y Tecnología de China han desarrollado una nueva técnica computacional. En lugar de depender de promedios simples, este nuevo método utiliza una descripción matemática precisa de cómo se comporta el fluido justo en el borde donde se encuentran dos tipos de roca diferentes. Al capturar el comportamiento real y nítido del fluido en estas fronteras críticas, el nuevo método proporciona una imagen mucho más clara y precisa del flujo subterráneo, incluso cuando la roca es extremadamente irregular y el fluido es complejo.

Los investigadores se centraron en un tipo específico de fluido conocido como fluido de ley de potencia. Estos son fluidos donde la relación entre la fuerza que los empuja y la velocidad con la que se mueven sigue una regla no lineal específica. En el mundo real, esto abarca una amplia gama de sustancias, desde soluciones poliméricas espesas hasta ciertos tipos de petróleo crudo. El desafío surge cuando estos fluidos fluyen a través de medios porosos —rocas llenas de diminutos agujeros— donde la permeabilidad, o la facilidad de flujo, cambia repentinamente de un lugar a otro. En tales entornos, la presión del fluido no cambia suavemente; en cambio, crea un pico agudo o una "singularidad" justo en el límite entre la roca de flujo fácil y la de flujo difícil. Los modelos computacionales tradicionales, que dividen la roca en una cuadrícula de pequeños triángulos, tienen dificultades con estos picos. Asumen que la presión cambia en línea recta entre puntos, lo cual es físicamente incorrecto en estas fronteras. Como resultado, los modelos producen tasas de flujo demasiado bajas y requieren un número enorme de puntos de cuadrícula diminutos para obtener incluso una aproximación rudimentaria de la verdad.

Para solucionar esto, el equipo desarrolló un Método Analítico Finito. La idea central es dejar de adivinar cómo se comporta el fluido en la frontera y, en su lugar, utilizar una descripción matemática exacta de ese comportamiento. Los investigadores se dieron cuenta de que, cerca del punto donde las propiedades de la roca cambian, la presión del fluido sigue un patrón curvo específico, aumentando o disminuyendo de una manera que depende de la naturaleza de la ley de potencia del fluido. Calcularon la forma exacta de esta curva, conocida como función propia (eigenfunction), para la parte más crítica del flujo. Este cálculo implicó la resolución de un conjunto complejo de ecuaciones que describen cómo se mueve el fluido en un área con forma de cuña donde las propiedades de la roca difieren. Al encontrar la solución más pequeña y significativa de estas ecuaciones, capturaron la forma dominante en que el fluido se comporta justo en el punto problemático.

Una vez comprendido este comportamiento local, los investigadores construyeron un nuevo esquema numérico alrededor de él. En lugar de utilizar el antiguo y defectuoso método de promediado para conectar los puntos de la cuadrícula, utilizaron su nueva solución analítica para calcular el flujo a través de las fronteras de cada celda triangular de su cuadrícula. Esto permitió que la computadora "supiera" exactamente cómo debería comportarse el fluido al comprimirse a través del estrecho cuello entre las zonas de alta y baja permeabilidad. El resultado es una simulación que respeta la física del problema en los puntos más difíciles, en lugar de suavizarlos mediante una aproximación.

El equipo probó este nuevo método contra varios escenarios desafiantes para ver cómo funcionaba. Comenzaron con un patrón clásico de "tablero de ajedrez", donde rocas de alta permeabilidad y baja permeabilidad se disponían en cuadrados alternos. En esta configuración, las propiedades de la roca cambiaban por un factor de cien. Cuando compararon su nuevo método con el tradicional, la diferencia fue evidente. El método tradicional subestimaba el flujo por márgenes enormes, perdiendo a veces casi el noventa por ciento del valor real en una cuadrícula gruesa. En contraste, el nuevo método se mantuvo notablemente cerca de la respuesta correcta, incluso cuando la cuadrícula no era muy fina. A medida que hacían la cuadrícula más fina, el nuevo método convergía a la respuesta correcta muy rápidamente, mientras que el método tradicional luchaba por alcanzarla, incluso con cuadrículas extremadamente finas.

También probaron el método en un dominio hexagonal regular y en un dominio con una distribución log-normal aleatoria de permeabilidad, que imita la naturaleza caótica de los reservorios subterráneos reales. En cada caso, el nuevo método resultó superior. Manejó las variaciones extremas en las propiedades de la roca con facilidad y se mantuvo preciso independientemente del tipo específico de fluido de ley de potencia que se simulara. Uno de los hallazgos más sorprendentes fue que la precisión del nuevo método no dependía fuertemente de las propiedades específicas del fluido o de la roca. Ya fuera que el fluido fuera muy espeso o muy fluido, o que el contraste de la roca fuera leve o extremo, el método entregaba resultados consistentes. Los investigadores descubrieron que un nivel moderado de refinamiento de la cuadrícula era suficiente para lograr una alta precisión, con errores que caían a menos del ocho por ciento incluso en los peores escenarios.

Este trabajo representa un avance significativo en la forma en que simulamos el flujo de fluidos en entornos subterráneos complejos. Al reemplazar una técnica de promediado defectuosa con una solución que respeta la verdadera física de la frontera, los investigadores han creado una herramienta que es tanto más precisa como más eficiente. Para las industrias que dependen de comprender cómo se mueven los fluidos a través de la tierra, desde la producción de energía hasta la gestión ambiental, esto significa que las simulaciones pueden ser más confiables y ejecutarse más rápido. El método no requiere la enorme potencia computacional que se necesitaba anteriormente para obtener una respuesta decente, haciendo que el modelado de alta fidelidad de sistemas heterogéneos complejos sea mucho más accesible. Los hallazgos sugieren que, al observar de cerca los puntos singulares donde el mundo cambia, podemos construir una mejor comprensión del todo.