Elastic wave propagation in fractured media with spring-type and frictional contact deformation laws

El suelo bajo nuestros pies rara vez es un bloque sólido e ininterrumpido. En muchos lugares, desde profundos reservorios de petróleo hasta la corteza misma de la Tierra, la roca está atravesada por fracturas: grietas planas y delgadas que pueden extenderse por kilómetros. Cuando la energía viaja a través de este paisaje fracturado, ya sea una onda sísmica de un terremoto o un pulso de sonido utilizado para inspeccionar un material, las grietas cambian la historia. Dispersan la energía, debilitan la señal y, a veces, incluso provocan que las caras de la roca se froten entre sí o se deslicen. Comprender exactamente cómo estas ondas interactúan con las grietas es vital para los geólogos que intentan mapear lo que hay bajo tierra y para los ingenieros que garantizan la seguridad de las estructuras. Sin embargo, la física de esta interacción es notoriamente difícil de precisar. Una simple grieta no es solo un hueco; es un límite complejo donde las caras de la roca pueden presionarse entre sí, deslizarse con fricción o abrirse y cerrarse bruscamente, todo mientras la onda pasa a través de ella en una fracción de segundo.

Durante décadas, los científicos han dependido de modelos simplificados para predecir estos comportamientos, tratando a menudo las grietas como si estuvieran hechas de resortes invisibles que se estiran y comprimen. Aunque son útiles, estos modelos fallan cuando las ondas son lo suficientemente fuertes como para obligar a las caras de la roca a tocarse, deslizarse o separarse por completo. No pueden capturar el movimiento de adherencia-deslizamiento (stick-slip) donde la fricción mantiene la roca en su lugar hasta que la fuerza es demasiado grande, ni pueden describir con precisión el endurecimiento no lineal que ocurre cuando una grieta se cierra por compresión. Un nuevo estudio realizado por investigadores de la Universidad de Bergen, el Centro de Investigación de Noruega y la Universidad Técnica de Múnich aborda estas lagunas mediante la construcción de un marco computacional integral que simula ondas elásticas moviéndose a través de roca fracturada con un realismo sin precedentes. En lugar de depender de suposiciones simplificadas, el equipo desarrolló un método que puede manejar cuatro niveles distintos de complejidad, que van desde el comportamiento básico tipo resorte hasta la mecánica de contacto total, donde se calculan la fricción y la apertura y cierre físico de las grietas en tiempo real.

Los investigadores construyeron un entorno digital unificado donde la roca y las grietas se tratan como un sistema único e interconectado. En sus simulaciones, la roca se representa como un volumen tridimensional, mientras que las fracturas se modelan como superficies bidimensionales incrustadas en ella. Esto permite que la computadora rastree cómo una onda que viaja a través de la roca golpea una grieta, cómo se deforma la grieta y cómo esa deformación envía ondulaciones de vuelta a la roca. El equipo probó cuatro formas específicas en las que las grietas podrían comportarse. Los dos primeros modelos utilizaron leyes de tipo resorte, donde la fuerza a través de la grieta es proporcional a cuánto se mueve la roca, con una versión que permite una relación lineal simple y otra que utiliza una regla no lineal más compleja que tiene en cuenta el hecho de que las grijas son más difíciles de comprimir a medida que se cierran. Los otros dos modelos introdujeron la realidad de la fricción y el contacto. En estos escenarios, las caras de las grietas pueden tocarse y deslizarse entre sí, gobernadas por leyes de fricción similares a las que impiden que los neumáticos de un coche patinen, mientras que también se aplica la regla de que las caras de la roca no pueden atravesarse entre sí.

Para asegurar que su nuevo marco fuera preciso, el equipo realizó una serie de pruebas rigurosas. Simularon ondas golpeando una sola grieta y compararon sus resultados computacionales con soluciones matemáticas conocidas para casos simples. Encontraron que su método podía predecir el comportamiento de la onda con alta precisión, coincidiendo con las expectativas teóricas tanto para la velocidad de la onda como para la cantidad de energía que pasaba o rebotaba. Luego llevaron el sistema más allá, simulando ondas lo suficientemente fuertes como para causar que las grietas se abrieran, se pegaran y se deslizaran. En estos escenarios más complejos, donde las caras de las grietas se frotaban entre sí con un coeficiente de fricción de 1.0, las simulaciones mostraron que el método seguía siendo estable y producía resultados consistentes. Los investigadores verificaron que su enfoque podía manejar los cambios repentinos en la física que ocurren cuando una grieta pasa de estar adherida a deslizarse, una transición que suele confundir a los modelos más simples.

Cuando el equipo comparó los cuatro modelos diferentes uno al lado del otro, las diferencias en los patrones de onda resultantes fueron sorprendentes. Las simulaciones revelaron que la elección del modelo importa significemente. Cuando se permitía que una grieta se cerrara y se endureciera según la ley no lineal, transmitía mucha más energía de la onda que una grieta modelada con un simple resorte lineal. Esto se debe a que el modelo no lineal tiene en cuenta el hecho de que, a medida que la grieta se cierra, las caras de la roca se presionan más fuertemente, haciendo que la grieta sea efectivamente más rígida y más transparente para la onda. Además, la inclusión de la fricción cambió la forma de los frentes de onda. En los modelos friccionales, la onda que emergía al otro lado de la grieta era más amplia y distorsionada que en las versiones sin fricción. Los investigadores observaron que la fricción causaba que las caras de las grietas se mantuvieran adheridas durante períodos más largos, retrasando la liberación de energía y alterando la forma en que la onda se dispersaba.

El estudio no se detuvo en las grietas individuales. Los investigadores aplicaron su marco a un dominio tridimensional que contiene seis fracturas que se intersectan, creando una pequeña red de grietas con diferentes propiedades. En este entorno complejo, asignaron diferentes límites de cierre máximo a diferentes conjuntos de fracturas, lo que significa que algunas grietas podían cerrarse más herméticamente que otras. La simulación mostró que el campo de ondas se volvió altamente heterogéneo. Las fracturas que eran más rígidas y podían cerrarse más herméticamente reflejaban más energía y transmitían menos, creando sombras distintas en el patrón de la onda. Los investigadores pudieron rastrear exactamente cuánto se cerraba cada grieta y cuánto se deslizaba, revelando que incluso pequeñas diferencias en las propiedades físicas de las grietas conducían a respuestas de deformación notablemente diferentes. Por ejemplo, una grieta que se permitía cerrar más herméticamente exhibía más comportamiento de "adherencia", manteniendo sus caras unidas durante más tiempo antes de deslizarse, lo que a su vez afectaba la onda que pasaba a través de ella.

Las implicaciones de estos hallazgos se extienden más allá de los resultados inmediatos. Al demostrar que un único marco computacional puede manejar desde un simple comportamiento de resorte hasta un complejo contacto friccional, los investigadores han proporcionado una herramienta que puede usarse para interpretar datos sísmicos del mundo real con mayor fidelidad. La capacidad de simular cómo las ondas se dispersan y atenúan en presencia de grietas friccionales, de apertura y de cierre significa que los geocientíficos pueden distinguir mejor entre diferentes tipos de formaciones rocosas y fracturas llenas de fluidos. El estudio confirma que ignorar los aspectos no lineales y friccionales de la mecánica de fractura puede conducir a errores significativos en la predicción de cómo se mueve la energía a través del subsuelo. Si bien el trabajo actual se centra en interacciones puramente mecánicas, el marco está construido sobre una base que es compatible con modelos de flujo de fluidos, lo que sugiere un camino hacia simulaciones aún más realistas donde el agua o el petróleo que se mueven a través de las grietas interactúan con las ondas que pasan. El trabajo es una demostración de que, al aceptar la plena complejidad de cómo se comportan las grietas de la roca, podemos construir una imagen más clara del mundo oculto bajo nuestros pies.