⚡ electrical engineering

Sensitivity Analysis of Shield Tunnel-Induced Ground Settlement to Stratigraphic Spatial Variability: From Random Field Simulation to Causal Decoupling

Este artículo presenta un marco integrado que combina la simulación de campo aleatorio, experimentos numéricos controlados y regresión simbólica para cuantificar cómo la variabilidad espacial geológica en el granito meteorizado afecta los asentamientos inducidos por túneles de escudo, identificando el módulo de deformación de los estratos excavados y de sobrecarga como los factores dominantes y derivando un modelo predictivo físicamente interpretable.

Autores originales: Chenlong Yue, Jian Zhang, Tugen Feng, Haibo Wang, Bin Zheng

Publicado 2026-08-06

📖 5 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Chenlong Yue, Jian Zhang, Tugen Feng, Haibo Wang, Bin Zheng

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

Imagina el suelo bajo nuestros pies no como un bloque de roca sólido e inalterable, sino como un pastel de capas gigante y desordenado. A veces el glaseado es grueso, otras veces es delgado; a veces el bizcocho es suave y esponjoso, otras veces es duro como un ladrillo. Esta es la realidad de la geología, especialmente cuando intentamos excavar túneles para metros subterráneos a través de él. Durante décadas, los ingenieros han intentado predecir cuánto se hundirá el terreno (asentamiento) cuando una gigantesca máquina tuneladora (TBM) se abre paso a través de él. Solían adivinar basándose en fórmulas simples o trataban al suelo como si fuera perfectamente uniforme, como si fuera un solo sabor de helado. Pero ahora sabemos que el suelo está lleno de sorpresas. Tiene "variabilidad espacial", que es una forma elegante de decir que el suelo cambia sus propiedades de un lugar a otro, incluso en distancias cortas. Si excavas en un punto, el terreno puede ser rígido; te mueves unos pocos metros y puede ser blando. Esta imprevisibilidad es un gran dolor de cabeza para los ingenieros porque, si calculan mal, el suelo podría hundirse demasiado, agrietando las aceras y los edificios de arriba. Así que, la gran pregunta es: ¿cómo predecimos este hundimiento cuando el suelo es tan desordenado y cambiante?

Este artículo se sumerge en esa realidad desordenada, mirando específicamente los túneles de escudo (aquellos que usan un cabezal de corte giratorio gigante) que atraviesan granito meteorizado en Shenzhen, China. Los investigadores querían averiguar exactamente cómo afecta la "desordenada" naturaleza del terreno al hundimiento y, lo que es más importante, qué parte específica del terreno importa más. No se limitaron a adivinar; construyeron un sofisticado gemelo digital del túnel utilizando un potente programa informático llamado PLAXIS 3D. Simularon el terreno como un "campo aleatorio", que es como lanzar dados para crear miles de versiones diferentes y realistas de las capas subterráneas, cada una con resistencias y rigideces ligeramente distintas. Descubrieron que, con solo cambiar la disposición natural de las rocas y el suelo, la cantidad que el terreno se hundió podía variar aproximadamente 3,38 mm. Eso puede parecer insignificante, pero en el mundo de la tunelería, es una diferencia enorme que podría significar la diferencia entre un túnel seguro y un edificio agrietado.

La parte más emocionante de su descubrimiento es cómo descubrieron por qué se hunde el terreno. Trataron el suelo como un equipo de dos partes: el "estrato excavado" (la capa que el túnel está cortando realmente) y la "sobrecarga" (las capas de tierra y roca situadas encima del túnel). Realizaron 37 experimentos computacionales diferentes, cambiando una cosa a la vez, como un científico ajustando una receta. Descubrieron que la rigidez (qué tan duro o blando es) de la capa que el túnel está cortando es la jefa absoluta. Si esa capa es más rígida, el terreno se hunde menos. La rigidez de las capas superiores también importa, pero es lo segundo más importante. Sorprendentemente, descubrieron que otros factores, como qué tan "pegajoso" es el suelo (cohesión) o el ángulo en el que se desliza (ángulo de fricción), apenas marcaron la diferencia. Incluso el número exacto de capas o cómo estaban apiladas no importaba mucho, siempre que la rigidez general fuera la misma. Es como darse cuenta de que, cuando estás empujando una caja pesada, el peso de la caja importa mucho más que el color de la alfombra sobre la que está apoyada.

Para que esto sea útil para los ingenieros del mundo real, el equipo utilizó un truco matemático ingenioso llamado "regresión simbólica". Piensa en esto como un robot que observa todos sus datos informáticos e intenta escribir una frase sencilla y legible (una ecuación) que explique los resultados, en lugar de dar simplemente una respuesta de "caja negra". Llegaron a una fórmula que utiliza solo tres números —la rigidez de la capa del túnel, la rigidez de las capas superiores y la pegajosidad de la capa del túnel— para predecir cuánto se hundirá el terreno. Esta fórmula fue sorprendentemente precisa, acertando la predicción con menos de medio milímetro de diferencia en promedio. También utilizaron una herramienta llamada SHAP para verificar su trabajo, lo que confirmó que la rigidez de la capa del túnel es, de hecho, la superestrella del espectáculo.

El artículo esencialmente argumenta contra la vieja idea de que necesitamos conocer cada pequeño detalle del terreno para hacer una buena predicción. En su lugar, sugieren que si nos enfocamos en obtener correctamente la "rigidez" de las capas adecuadas, podemos hacer mejores predicciones sin perdernos en los detalles. Excluyeron explícitamente la idea de que el ángulo de fricción o el número exacto de capas sean actores principales en este tipo específico de terreno de granito. Aunque sus resultados se basan en simulaciones por computadora y no en excavaciones físicas en el mundo real, la consistencia entre sus experimentos controlados y su nueva fórmula les otorga una alta confianza en sus hallazgos. No pretenden haber resuelto todos los problemas de tunelería del universo, pero han proporcionado una forma mucho más clara, simple y precisa de entender cómo la variabilidad natural del terreno afecta la seguridad de los túneles, ayudando a los ingenieros a construir metros más seguros en el futuro.

Resumen Técnico: Análisis de Sensibilidad del Asentamiento del Terreno Inducido por Túneles de Escudo ante la Variabilidad Espacial Estratigráfica

Planteamiento del Problema
La excavación de túneles con escudo en entornos geológicos complejos, particularmente en formaciones de granito meteorizado caracterizadas por interfaces estratigráficas irregulares y contrastes de rigidez pronunciados, introduce una incertidumbre significativa en las predicciones del asentamiento del terreno. Los análisis deterministas tradicionales suelen asumir propiedades de suelo homogéneas o dependen de valores medios, fallando al capturar la naturaleza dual de la incertidumbre geológica: la variabilidad discreta de las secuencias estratigráficas y la variabilidad espacial continua de los parámetros geotécnicos dentro de los estratos. Además, los análisis de sensibilidad existentes frecuentemente tratan todo el perfil del suelo como un sistema unificado, pasando por alto la distinción mecánica fundamental entre el estrato excavado (la fuente de la deformación) y la sobrecarga (la vía de transmisión de la deformación). Este estudio aborda la necesidad de cuantificar, explicar y predecir la incertidumbre del asentamiento derivada de estas variabilidades espaciales, específicamente en el contexto del Ferrocarril Interurbano de Shenzhen.

Metodología
Los autores desarrollaron un marco integrado que combina la simulación de campos aleatorios, experimentos numéricos controlados y regresión simbólica:

Modelado Numérico y Validación: Se estableció un modelo de elementos finitos 3D utilizando PLAXIS 3D, simulando un túnel de escudo con equilibrado de presión de tierras (EPB) de 8.0 m de diámetro y una profundidad de enterramiento promedio de 20 m. El modelo utilizó el modelo constitutivo de Mohr–Coulomb para el suelo y modelos elásticos lineales para los componentes estructurales. Se validó frente a datos de monitoreo de campo y el método empírico de Peck, demostrando una alta precisión en la reproducción de las cubetas de asentamiento superficial en granito meteorizado.
Marco de Acoplamiento de Dual-Incertidumbre: Para simular la variabilidad espacial geológica, el estudio empleó un enfoque de dos módulos:
- Simulación Estratigráfica: Se utilizaron Cadenas de Markov Acopladas (CMC) para generar 11 realizaciones geológicas independientes, capturando la incertidumbre discreta de las secuencias estratigráficas (por ejemplo, interfaces de meteorización fluctuantes).
- Asignación de Parámetros: Se aplicaron Campos Aleatorios Condicionales (CRF) con funciones de autocorrelación exponencial dentro de cada realización para asignar valores espacialmente variables de módulo de deformación ( $E$ ), cohesión ( $c$ ) y ángulo de fricción ( $\phi$ ), capturando la variabilidad continua de los parámetros.
Experimentos Paramétricos Controlados: Se llevó a cabo una serie de 37 experimentos numéricos de un solo parámetro a la vez (OVAT) para desacoplar los efectos de los parámetros geológicos. El estudio varió independientemente $E$ , $c$ y $\phi$ tanto para el estrato excavado como para la sobrecarga, así como el número de capas, para aislar sus contribuciones causales al asentamiento superficial máximo ( $S_{max}$ ).
Regresión Simbólica e Interpretación: El conjunto de datos se expandió utilizando el Muestreo por Hipercubo Latino (LHS) para cubrir el espacio de parámetros de manera más uniforme. Posteriormente, se empleó la regresión simbólica para derivar una ecuación de predicción explícita e interpretable. Finalmente, se aplicó el análisis SHAP (Explicaciones de Valores de Shapley Aditivos) al modelo resultante para verificar la clasificación de sensibilidad y asegurar la interpretabilidad física.

Resultos Clave

Impacto de la Variabilidad Espacial: Las 11 realizaciones de campos aleatorios produjeron asentamientos superficiales máximos que oscilaron entre $-11.25$ mm y $-14.63$ mm, una variación de $3.38$ mm. Esto confirma que la variabilidad espacial geológica por sí sola puede causar una incertidumbre de asentamiento prácticamente significativa, la cual los modelos deterministas basados en valores medios pueden subestimar.
Clasificación de Sensibilidad: Los experimentos controlados establecieron una jerarquía de sensibilidad clara: $E_{exc}$ $E_{e x c}$ (módulo de deformación del estrato excavado) $> E_{upper}$ $> E_{u pp er}$ (módulo de deformación de la sobrecarga) $> c_{exc}$ $> c_{e x c}$ (cohesión del estrato excavado).
- Las variaciones en $E_{exc}$ (35–75 MPa) causaron un cambio de asentamiento de $3.72$ mm.
- Las variaciones en $E_{upper}$ (5–30 MPa) causaron un cambio de $2.57$ mm.
- Las variaciones en $c_{exc}$ (25–48 kPa) causaron un cambio de $1.09$ mm.
- El ángulo de fricción ( $\phi$ ) y la complejidad de la estratificación tuvieron efectos insignificantes, con variaciones de asentamiento inferiores a $0.18$ mm.
Roles Mecánicos: El estudio identificó que el estrato excavado actúa como la fuente de deformación, donde tanto la rigidez como, en menor medida, la resistencia gobiernan la generación de la deformación. Por el contrario, la sobrecarga actúa principalmente como una vía de transmisión de la deformación, donde la rigidez es el único factor determinante; los parámetros de resistencia ( $c$ y $\phi$ ) en la sobrecarga no mostraron influencia medible en el asentamiento superficial.
Modelo Predictivo: La ecuación de regresión simbólica derivada, que incorpora $E_{exc}$ , $E_{upper}$ y $c_{exc}$ , alcanzó un $R^2$ de $0.806$, un MAE de $0.265$ mm y un RMSE de $0.454$ mm. El análisis SHAP en este modelo reprodujo la clasificación de sensibilidad ( $E_{exc} \gg E_{upper} > c_{exc}$ ), confirmando su consistencia física.

Significancia y Reivindicaciones
El artículo afirma avanzar en el análisis del asentamiento de túneles de escudo, pasando de una simple caracterización de la incertidumbre a un enfoque de desacoplamiento causal basado en mecanismos. Sus principales contribuciones incluyen:

Innovación Metodológica: El desarrollo de un marco de dual-incertidumbre que modela simultáneamente la incertidumbre de la secuencia estratigráfica (vía CMC) y la variabilidad de los parámetros dentro del estrato (vía CRF), ofreciendo una representación más realista de las formaciones de granito meteorizado que los modelos estocásticos deterministas o de fuente única.
Perspectiva Mecanicista: El desacoplamiento explícito del estrato excavado y la sobrecarga revela que ambos desempeñan roles mecánicos fundamentalmente diferentes. Esto desafía los análisis de sensibilidad global convencionales que tratan el perfil del suelo como una unidad única, proporcionando una base más físicamente significativa para la evaluación de parámetros.
Utilidad Práctica: El estudio demuestra que el módulo de deformación es el parámetro dominante que controla el asentamiento, significativamente más que la cohesión o el ángulo de fricción. Propone una ecuación de regresión simbólica explícita e interpretable que sirve como una alternativa transparente a los modelos de aprendizaje automático de "caja negra". Esta herramienta facilita la predicción rápida de asentamientos preliminares, el cribado de parámetros geológicos y la evaluación de riesgos sin sacrificar la interpretabilidad física.

Los autores posicionan este trabajo como un paso hacia el modelado predictivo informado por la física en la ingeniería geotécnica, cerrando la brecha entre las complejas simulaciones estocásticas y la toma de decisiones prácticas en ingeniería.

Más como este