Deformation Due to a Pair of Infinite and a Pair of infinite-finite planar faults in Fractional Standard Linear Solid

Resumen Técnico: Deformación debida a un par de fallas planas infinitas y un par de fallas planas infinitas-finitas en un Sólido Lineal Estándar Fraccionario

Planteamiento del Problema
Las regiones sísmicamente activas contienen frecuentemente múltiples fallas interactuantes donde el deslizamiento a lo largo de una falla influye significativamente en la acumulación y distribución de esfuerzos en las fallas vecinas. Mientras que los estudios previos han modelado las interacciones de fallas utilizando derivadas de orden entero dentro de medios viscoelásticos (por ejemplo, modelos de Maxwell o de Sólido Lineal Estándar), estos enfoques tradicionales a menudo pasan por alto los efectos de "memoria" inherentes a los materiales geológicos. Este artículo aborda la deformación del sistema litosfera-astenosfera de la Tierra durante el periodo aseísmico mediante el modelado de dos escenarios específicos que involucran fallas de desgarre (strike-slip) de movimiento lento (creep) que interactúan:

1. Caso-I: La interacción entre dos fallas enterradas infinitas.
2. Caso-II: La interacción entre una falla enterrada infinita y una falla de superficie finita.

El estudio postula que el medio terrestre se comporta como un semiespacio viscoelástico caracterizado por un modelo de Sólido Lineal Estándar Fraccionario (FSLS), el cual captura mejor el comportamiento reológico dependiente del tiempo y los efectos de memoria en comparación con los modelos de orden entero.

Metodología
Los autores desarrollaron modelos matemáticos teóricos para derivar soluciones analíticas para los componentes de desplazamiento, esfuerzo y deformación. La metodología involucra las siguientes etapas clave:

• Formulación Matemática: El medio se define utilizando una ecuación constitutiva de sólido lineal estándar fraccionario que emplea el operador de derivada fraccionaria de Caputo (${}_{0}D^{\alpha}_{t}$), donde el orden $\alpha \in (0, 1]$. Las relaciones de esfuerzo-deformación incorporan rigidez efectiva ($\mu$) y viscosidad efectiva ($\eta$).
• Ecuaciones Gobernantes: El análisis asume una deformación cuasiestática (negligenciando las fuerzas inerciales). Se establecen las ecuaciones de movimiento de esfuerzo y las condiciones de contorno tanto para geometrías de falla infinitas como finitas.
• Técnica de Solución: Los autores emplean la transformada de Laplace para manejar las ecuaciones diferenciales fraccionarias dependientes del tiempo. Las soluciones se derivan utilizando:
  • El Principio de Correspondencia (Rybicki, 1971).
  • Una técnica de Función de Green Modificada (Maruyama, 1966) para dar cuenta de las condiciones de dislocación a través de los planos de falla.
  • La función de Mittag-Leffler ($E_{\alpha}$), que surge naturalmente de la transformada inversa de Laplace de los operadores fraccionarios.
• Escenarios Analizados:
  • Se derivan soluciones tanto para el estado antes del creep de la falla (acumulación de esfuerzo tectónico) como después del creep de la falla (dislocación debido al deslizamiento).
  • Se aplican velocidades de creep ($V_1, V_2$) y funciones de creep dependientes de la profundidad para simular el movimiento de las fallas $F_1$ y $F_2$.
  • Se realizan computaciones numéricas utilizando parámetros derivados de datos geológicos (por ejemplo, rigidez $\mu = 3.5 \times 10^{10}$ N/m², viscosidad $\eta = 5 \times 10^{19}$ Pa·s) para generar resultados gráficos para diversos ángulos de inclinación de falla ($\theta$), órdenes fraccionarios ($\alpha$) y velocidades de creep.

Contribuciones Clave

1. Modelado Fraccionario de la Interacción de Fallas: El artículo extiende los modelos de interacción de fallas existentes al incorporar el cálculo fraccionario, permitiendo una representación más precisa del comportamiento viscoelástico dependiente de la memoria de la litosfera.
2. Soluciones Analíticas para Geometrías Complejas: El estudio proporciona expresiones analíticas de forma cerrada para los campos de desplazamiento, esfuerzo y deformación para dos configuraciones distintas: un par de fallas enterradas infinitas y un sistema mixto de una falla infinita y una falla de superficie finita.
3. Cuantificación de los Efectos de Memoria: La investigación demuestra explícitamente cómo el orden fraccionario $\alpha$ influye en la acumulación y liberación de esfuerzos, vinculando el parámetro matemático con la "memoria" física del medio geológico.

Resultos y Hallazgos Numéricos
El análisis numérico, respaldado por representaciones gráficas, arroja varias observaciones específicas:

• Magnitudes de Desplazamiento: El cambio en el desplazamiento superficial cerca de las fallas es del orden de $10^{-3}$ m/año. Esto concuerda bien con los datos geofísicos observados (por ejemplo, los informes del USGS de 2–6 mm/año), mientras que los modelos previos de orden entero predecían valores en el rango de $10^{-5}$, mostrando una discrepancia con las observaciones del mundo real.
• Influencia de la Geometría de la Falla:
  • En el Caso-I (dos fallas infinitas), la magnitud del desplazamiento cerca de la segunda falla es mayor que cerca de la primera.
  • En el Caso-II (falla infinita y finita), la magnitud del desplazamiento cerca de la falla finita es menor que cerca de la falla infinita.
  • El desplazamiento alcanza picos cerca de la traza de la falla y disminuye a medida que aumenta la distancia, tendiendo eventualmente a cero.
• Efecto del Orden Fraccionario ($\alpha$): El estudio encuentra que la tasa de acumulación y liberación de esfuerzo depende altamente de $\alpha$. Valores más bajos de $\alpha$ (que indican efectos de memoria más fuertes) resultan en magnitudes más altas de acumulación de esfuerzo cortante y una relajación más lenta. A medida que $\alpha$ aumenta hacia 1 (aproximándose al comportamiento entero), la tasa de liberación de esfuerzo aumenta.
• Fuerza Impulsora y Distribución de Esfuerzos: La fuerza impulsora total a través de las fallas muestra una liberación rápida tras el inicio del creep. La magnitud de esta liberación varía con la inclinación de la falla y la velocidad de creep. Los mapas de contorno revelan que la acumulación de esfuerzo es más significativa cerca de la superficie libre y las zonas de dislocación de la falla, mientras que la reducción de esfuerzo es prominente en las profundidades correspondientes a las ubicaciones de las fallas.
• Análisis de Deformación: Se encuentran variaciones de la deformación cortante ($e_{12}$) en el rango de $10^{-4}$, consistentes con datos de eventos reales como el terremoto de Tohoku-oki de 2011.

Significancia y Reivindicaciones
Los autores afirman que su modelo fraccionario ofrece una representación más general y precisa de la dinámica de interacción de fallas en comparación con los modelos tradicionales de orden entero. Específicamente:

• El enfoque fraccionario logra capturar con éxito los efectos de memoria y la no linealidad de los medios viscoelásticos, que a menudo se pasan por alto en los modelos estándar.
• Las predicciones del modelo respecto a las tasas de desplazamiento ($10^{-3}$) muestran una fuerte concordancia con los datos geofísicos observados, lo que sugiere que el cálculo fraccionario es una herramienta superior para predecir el comportamiento sísmico y la deformación del subsuelo.
• Los hallazgos indican que la interacción entre fallas influye en el campo de esfuerzos de manera más extensa que una sola falla, y que el orden fraccionario $\alpha$ es un parámetro crítico en la determinación de la tasa de liberación de esfuerzo.

El artículo concluye que estos modelos teóricos mejoran la comprensión de la deformación del subsuelo y su papel al influir en la dinámica de las fallas y la generación de terremotos, proporcionando un marco robusto para analizar los periodos aseísmicos en zonas sísmicamente activas.