🔬 physics

Elastic self-averaging does not ensure representative failure behavior in field-constrained basalt fracture networks

Este estudio demuestra que, si bien las propiedades elásticas de las redes de fractura de basalto pueden estabilizarse en ciertas escalas, el volumen elemental representativo (VER) para el comportamiento de falla es significativamente mayor y distinto debido a la influencia dominante del declive de la resistencia sobre los cambios de rigidez, lo que requiere un enfoque multicriterio para definir los tamaños de modelos representativos.

Autores originales: Ming Li, Xiyong Chen

Publicado 2026-10-07

📖 7 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Ming Li, Xiyong Chen

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

En lo profundo del subsuelo, donde se excavan enormes cavernas para centrales hidroeléctricas, la roca rara vez es un bloque sólido e ininterrumpido. En su lugar, es un complejo rompecabezas de fragmentos de piedra mantenidos unidos por una red de grietas, fallas y juntas. Los ingenieros necesitan saber cómo se comportará esta roca fracturada cuando se excave en ella o se le apliquen cargas pesadas. Para realizar estas predicciones, dependen de un concepto llamado "volumen elemental representativo". En términos sencillos, este es el tamaño mínimo de una muestra de roca necesario para representar con precisión el comportamiento de todo el macizo rocoso. Si una muestra es demasiado pequeña, podría omitir una grieta crítica y dar una falsa sensación de seguridad; si es demasiado grande, los cálculos se vuelven imposiblemente costosos. Durante décadas, la suposición estándar ha sido que una vez que una muestra de roca es lo suficientemente grande como para mostrar una rigidez promedio estable, es lo suficientemente grande para predecir todo lo demás, incluyendo cómo y cuándo fallará la roca.

Un nuevo estudio desafía esta creencia largamente sostenida. Investigadores de la Universidad de Ciencias y Tecnología de Sichuan investigaron si una muestra de roca que parece estable en términos de su rigidez es realmente estable en términos de su resistencia y capacidad para absorber energía antes de romperse. Se centraron en el basalto, una roca volcánica común que se encuentra en las cavernas subterráneas del suroeste de China. Al crear cientos de modelos computacionales de roca fracturada a diferentes tamaños, que variaban desde cubos de cinco metros hasta enormes bloques de cincuenta metros, rastrearon el camino desde una roca elástica y estable hasta una que se desmorona y falla. Su trabajo revela que la roca puede parecer perfectamente estable y representativa en cierto tamaño, pero aun así puede ocultar una tendencia peligrosa a fallar mucho antes y de forma más violenta de lo que sugeriría una muestra más pequeña.

Los investigadores comenzaron construyendo una biblioteca digital de la estructura interna de la roca. Utilizaron datos del mundo real de un túnel de exploración para crear dos tipos de modelos de roca fracturada: uno que representa la roca "débilmente meteorizada", que ha estado expuesta al aire y al agua y contiene más grietas, y otro que representa la roca "fresca", que es más profunda e intacta. Generaron doscientas cuarenta redes tridimensionales únicas de fracturas para cada tipo de roca, variando el tamaño del cubo virtual de cinco a cincuenta metros. Primero, probaron cómo se deformaba la roca bajo presión sin llegar a romperse. Descubrieron que la rigidez promedio de la roca, y el número de grietas por volumen, se volvieron estadísticamente estables una vez que el cubo alcanzó los diez metros de tamaño. En este punto, el número promedio de grietas y el área superficial de las fracturas dejaron de cambiar significamente a medida que la muestra crecía.

Sin embargo, la estabilidad en el número de grietas no significa que la muestra esté lista para el diseño de ingeniería. Los investigadores notaron una restricción geométrica crucial: las grietas más grandes en sus modelos tenían treinta metros de largo. Un cubo de diez metros es simplemente demasiado pequeño para contener una gricción de treinta metros sin cortarla en los bordes. Por lo tanto, argumentaron que un modelo debe tener al menos treinta metros de ancho para ser una muestra "admisible en los límites", lo que significa que puede contener el rango completo de tamaños de fractura presentes en la roca real sin una truncación artificial. Esto estableció un nuevo y mayor umbral para lo que cuenta como un tamaño de muestra válido, elevando el límite de los diez metros a los treinta metros.

Con las reglas geométricas establecidas, el equipo pasó a la parte más crítica del experimento: simular cómo se rompe realmente la roca. Tomaron un modelo cuidadosamente calibrado para cada tipo y tamaño de roca, lo seccionaron en secciones bidimensionales y los sometieron a una compresión intensa, imitando el estrés de una excavación subterránea. Observaron todo el proceso, desde el empuje inicial hasta el momento en que la roca alcanzó su resistencia máxima y comenzó a desmoronarse. Los resultados fueron sorprendentes. Mientras que la rigidez de la roca cambió muy poco a medida que los modelos crecían de cinco a cincuenta metros, la resistencia de la roca disminuyó drásticamente. El esfuerzo máximo que la roca podía soportar disminuyó casi un treinta por ciento a medida que el tamaño aumentaba. Aún más sorprendente fue la energía que la roca podía absorber antes de fallar. Esta energía, que representa el trabajo que la roca puede realizar antes de romperse, cayó más de un cincuenta por ciento en los modelos más grandes en comparación con los más pequeños.

El estudio explica esta discrepancia observando cómo falla la roca. En los modelos más pequeños, las grietas eran menos numerosas y más aisladas. A medida que los modelos crecían, contenían fracturas más largas y más oportunidades para que las pequeñas grietas se conectaran entre sí. Los investigadores observaron que en los modelos más grandes, las grietas locales se conectarían a través de toda la muestra, formando un camino continuo que causaría el colapso de la roca mucho antes. Este proceso, conocido como localización de la fractura, significó que la roca alcanzó su punto de ruptura a un nivel de esfuerzo mucho menor, a pesar de que su rigidez inicial permaneciera igual. Los investigadores calcularon que la pérdida de la capacidad de energía se debió casi por completo a esta caída en la resistencia, más que a un cambio en la rigidad de la roca. Encontraron que la relación entre la resistencia y la energía seguía una regla matemática precisa donde una pequeña caída en la resistencia conduce a una caída masiva en la capacidad de energía, un hecho que fue confirmado por sus simulaciones.

El equipo también examinó la dirección en la que se cargaba la roca. Incluso en los modelos más grandes de cincuenta metros, la roca no se comportaba de la misma manera en todas las direcciones. Dependiendo de la orientación de las grietas respecto a la presión, la roca podía ser significativamente más fuerte o más débil. La resistencia promedio pudo haberse estabilizado alrededor de los cuarenta metros, pero la diferencia entre las direcciones más fuertes y las más débiles seguía siendo sustancial. Esto significa que un modelo que parece "representativo" en promedio aún podría ser peligrosamente engañoso si la estructura de ingeniería real resulta alinearse con la dirección más débil. El estudio demostró que la estructura interna de la roca seguía siendo anisotrópica, o dependiente de la dirección, incluso en las escalas más grandes probadas.

En última instancia, la investigación sugiere que no existe un único "número mágico" para el tamaño de una muestra de roca que funcione para todos los propósitos de ingeniería. Una muestra de diez metros puede ser suficiente para contar las grietas, y una muestra de treinta metros puede ser necesaria para incluir el rango completo de tamaños de fractura, pero ninguna es suficiente para predecir cómo fallará la roca bajo carga. Crucialmente, el estudio encontró que incluso en la escala más grande de cincuenta metros probada, los modelos no lograron resolver la convergencia de la absorción de energía o el comportamiento de las direcciones de carga más conservadoras (más débiles). Para comprender la absorción de energía y el riesgo de colapso repentino, los ingenieros pueden necesitar modelos incluso más grandes que cincuenta metros, o deben aceptar que sus predicciones llevarán un margen de incertidumbre significativo. El estudio concluye que la estabilidad de la rigidez de una roca es un indicador deficiente de su estabilidad ante el fallo. Una roca puede parecer sólida y predecible en una escala pequeña, pero puede ocultar una red compleja de conexiones que causarán que se haga añicos inesperadamente cuando la escala sea lo suficientemente grande como para permitir que esas conexiones se formen. Esta distinción es vital para la seguridad de las estructuras subterráneas, recordando a los ingenieros que el comportamiento de un macizo rocoso no se trata solo de cuántas grietas tiene, sino de cómo están dispuestas esas grietas y cómo interactúan cuando se aplica la presión.

Resumen Técnico: Auto-promediado Elástico y Falla Representativa en Redes de Fractura de Basalto

Planteamiento del Problema
El concepto de un Volumen Elemental Representativo (REV, por sus siglas en inglés) es central en la mecánica de rocas, aunque a menudo se trata como una constante de material única. En macizos rocosos fracturados, la representatividad geométrica (por ejemplo, la densidad de fractura) no implica necesariamente representatividad mecánica (por ejemplo, resistencia, absorción de energía o localización de la falla). Un problema crítico no resuelto es si la estabilización de una propiedad elástica media asegura que la trayectoria completa de carga —incluyendo la resistencia máxima, la energía pre-pico y la organización de la falla post-pico— también haya convergido. Estudios previos han reducido la representatividad únicamente al módulo elástico o a la resistencia máxima, pasando potencialmente por alto el desacoplamiento del escalamiento de la rigidez, la resistencia y la energía. Además, no está claro cómo se comportan las redes de fractura con restricciones de campo, particularmente aquellas con fracturas largas comparables al tamaño del dominio, con respecto al auto-promediado y la admisibilidad de los límites.

Metodología
El estudio utiliza un marco numérico jerárquico que combina la generación de Redes de Fractura Discreta (DFN) con simulaciones del Método de Elementos Finitos-Discretos (FDEM).

Contexto Geológico: Los modelos se basan en basalto fresco y débilmente meteorizado de una caverna de energía hidroeléctrica subterránea en el suroeste de China. Las familias de fracturas se calibraron mediante mapeo de campo y escaneo láser 3D, resultando en dos poblaciones con distintas intensidades de fractura ( $P_{32}$ ) pero con una longitud máxima de fractura compartida de 30 m.
Escala y Conjunto: Se analizaron seis tamaños de cubo (5, 10, 20, 30, 40 y 50 m). Para el análisis geométrico y elástico, se generaron 20 realizaciones independientes de DFN para cada tamaño y tipo de roca (240 realizaciones en total). Para el análisis de falla mecánica, un DNF calibrado en 3D por cada tipo de roca y tamaño se seccionó en tres planos ortogonales y se cargó en dos direcciones, produciendo 72 respuestas completas de esfuerzo-deformación.
Análisis Elástico: Se utilizó el software DFNLab para calcular la complacencia direccional de límite inferior (asumiendo juntas de deslizamiento libre) para todas las 240 realizaciones con el fin de evaluar el auto-promediado elástico y la disminución de la varianza.
Análisis de Falla: Las simulaciones FDEM emplearon un Modelo de Zona Cohesiva Extrínseco (ECZM) para simular el agrietamiento inducido por esfuerzo, permitiendo que las fracturas preexistentes interactúen con nuevas grietas de tracción y corte. Las métricas incluyeron resistencia máxima ( $\sigma_p$ ), módulo secante ( $E_{50}$ ), energía pre-pico, relación de grietas y dimensión fractal.
Marco Analítico: El estudio introduce un factor de forma de energía adimensional ( $\eta$ ) para descomponer el escalamiento de la energía en contribuciones de resistencia, rigidez y forma de la curva. Además, define una jerarquía de criterios de representatividad: (1) estacionariedad estadística de la intensidad de fractura, (2) admisibilidad de los límites relativa al tamaño máximo de fractura, (3) convergencia de la media elástica y la varianza, (4) convergencia de la resistencia máxima, y (5) convergencia de la energía y la localización.

Resultos Clave

Convergencia Desacoplada: La intensidad geométrica ( $P_{30}$ y $P_{32}$ ) alcanzó la estacionariedad estadística a los 10 m (CV < 10%). Sin embargo, debido a la longitud máxima de fractura de 30 m, los dominios de 10 m sufrieron un sesgo de truncamiento de frontera. La primera escala de admisibilidad de frontera fue de 30 m.
Auto-promediado Elástico: El módulo elástico medio de límite inferior se estabilizó entre 40 y 50 m (cambios < 0.7%). Sin embargo, el coeficiente de variación (CV) continuó disminuyendo significativamente (de ~7.5% a ~2.2% para la roca débilmente meteorizada), y el exponente de auto-promediado ( $\beta$ ) fue de 0.43–0.47, muy por debajo del valor ideal de 1.5 para 3D. Esto indica que, aunque la media se estabiliza, la incertidumbre decae lentamente y el sistema no ha alcanzado un régimen asintótico real dentro de las escalas probadas.
Escalamiento de Resistencia y Energía: Mientras que el módulo secante promedio direccional ( $E_{50}$ ) fue casi insensible a la escala (cambiando solo un 3.6–4.0% de 5 a 50 m), la resistencia máxima disminuyó entre un 26.7% y 30.0%, y la densidad de energía pre-pico disminuyó entre un 47.3% y 55.3%.
Descomposición de la Energía: El factor de forma de energía $\eta$ se mantuvo casi constante (~1.21). La validación cruzada confirmó que la tendencia de la energía está gobernada principalmente por la disminución cuadrática de la resistencia de pico ( $\sigma_p^2$ ), no por cambios en la rigidez. La identidad exacta $\alpha_W = 2\alpha_\sigma - \alpha_E + \alpha_\eta$ se mantuvo, siendo la resistencia la contribución dominante al escalamiento de la energía.
Localización de la Falla: A medida que la escala aumentó, la relación de grietas y la dimensión fractal aumentaron significamente (la relación de grietas subió entre un 167–234%). El modo de falla cambió de un daño de punta aislado a trayectorias macroscópicas complejas que abarcan fronteras, involucrando una coalescencia significativa de modo mixto. A los 50 m, la anisotropía direccional ( $A_c \approx 1.43–1.46$ ) persistió, y la meseta de resistencia media no garantizó la convergencia para direcciones de diseño específicas.

Significado y Reivindicaciones
El artículo argumenta que un único valor escalar de REV es insuficiente para los macizos rocosos fracturados. En su lugar, propone un marco de representatividad dependiente de la etapa donde diferentes propiedades convergen a diferentes escalas:

10 m: Suficiente para la estacionariedad estadística de la intensidad de fractura ( $P_{30}/P_{32}$ ).
30 m: La escala mínima para la geometría de admisibilidad de frontera (sin truncamiento de la cola de fractura de 30 m).
~40 m: Suficiente para la estabilización de la media elástica de límite inferior y la resistencia máxima promedio direccional.
>50 m (No resuelto): La absorción de energía, la localización post-pico y las respuestas direccionales conservadoras no han convergido hacia el límite de 50 m.

El estudio concluye que el auto-promediado elástico no garantiza un comportamiento de falla representativo. Un tamaño de modelo que parece representativo para la rigidez o la resistencia media puede subestimar significativamente la pérdida de energía y sobreestimar la capacidad de deformación debido a la dependencia cuadrática de la energía respecto a la resistencia. Los autores enfatizan que la selección del modelo de ingeniería debe ser específica para la respuesta: los dominios más pequeños pueden ser suficientes para evaluaciones elásticas, pero se requieren dominios más grandes o envolventes de incertidumbre explícitas para el análisis de margen de deformación, liberación de energía y trayectoria de falla. Los hallazgos resaltan que la estacionariedad geométrica y la admisibilidad de la frontera son condiciones necesarias pero no suficientes para la homogeneización mecánica en redes de fractura con restricciones de campo.

Más como este