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

Resumo Técnico: Deformação Devido a um Par de Falhas Planares Infinitas e um Par de Infinitas-Finitas em um Sólido Linear Padrão Fracionário

Enunciado do Problema
Regiões sismicamente ativas frequentemente contêm múltiplas falhas interagentes onde o deslizamento ao longo de uma falha influencia significativamente a acumulação e a distribuição de tensão em falhas vizinhas. Embora estudos anteriores tenham modelado interações de falhas usando derivadas de ordem inteira dentro de meios viscoelásticos (por exemplo, modelos de Maxwell ou de Sólido Linear Padrão), essas abordagens tradicionais frequentemente negligenciam os efeitos de "memória" inerentes aos materiais geológicos. Este artigo aborda a deformação do sistema litosfera-astenosfera da Terra durante o período aseísmico, modelando dois cenários específicos envolvendo falhas de deslizamento lateral (strike-slip) de fluência (creep) interagentes:

1. Caso-I: A interação entre duas falhas enterradas infinitas.
2. Caso-II: A interação entre uma falha enterrada infinita e uma falha de superfície finita.

O estudo postula que o meio terrestre se comporta como um semiespaço viscoelástico caracterizado por um modelo de Sólido Linear Padrão Fracionário (FSLS), que captura melhor o comportamento reológico dependente do tempo e os efeitos de memória em comparação com modelos de ordem inteira.

Metodologia
Os autores desenvolveram modelos matemáticos teóricos para derivar soluções analíticas para os componentes de deslocamento, tensão e deformação. A metodologia envolve as seguintes etapas principais:

• Formulação Matemática: O meio é definido usando uma equação constitutiva de sólido linear padrão fracionário envolvendo o operador de derivada fracionária de Caputo (${}_{0}D^{\alpha}_{t}$), onde a ordem $\alpha \in (0, 1]$. As relações tensão-deformação incorporam rigidez efetiva ($\mu$) e viscosidade efetiva ($\eta$).
• Equações Governantes: A análise assume deformação quase estática (negligenciando forças de inércia). As equações de movimento de tensão e as condições de contorno são estabelecidas para ambas as geometrias de falha (infinita e finita).
• Técnica de Solução: Os autores empregam a transformada de Laplace para lidar com as equações diferenciais fracionárias dependentes do tempo. As soluções são derivadas usando:
  • O Princípio da Correspondência (Rybicki, 1971).
  • Uma técnica de Função de Green Modificada (Maruyama, 1966) para considerar as condições de descontinuidade através dos planos de falha.
  • A função de Mittag-Leffler ($E_{\alpha}$), que surge naturalmente da transformada inversa de Laplace dos operadores fracionários.
• Cenários Analisados:
  • As soluções são derivadas para o estado antes da fluência da falha (acumulação de tensão tectônica) e após a fluência da falha (deslocamento devido ao deslizamento).
  • Velocidades de fluência ($V_1, V_2$) e funções de fluência dependentes da profundidade são aplicadas para simular o movimento das falhas $F_1$ e $F_2$.
  • Computações numéricas são realizadas utilizando parâmetros derivados de dados geológicos (por exemplo, rigidez $\mu = 3,5 \times 10^{10}$ N/m², viscosidade $\eta = 5 \times 10^{19}$ Pa·s) para gerar resultados gráficos para várias inclinações de falha ($\theta$), ordens fracionárias ($\alpha$) e velocidades de fluência.

Principais Contribuições

1. Modelagem Fracionária de Interação de Falhas: O artigo estende os modelos de interação de falhas existentes ao incorporar o cálculo fracionário, permitindo uma representação mais precisa do comportamento viscoelástico dependente de memória da litosfera.
2. Soluções Analíticas para Geometrias Complexas: O estudo fornece expressões analíticas de forma fechada para os campos de deslocamento, tensão e deformação para duas configurações distintas: um par de falhas enterradas infinitas e um sistema misto de falhas de superfície infinitas e finitas.
3. Quantificação de Efeitos de Memória: A pesquisa demonstra explicitamente como a ordem fracionária $\alpha$ influencia a acumulação e a liberação de tensão, vinculando o parâmetro matemático ao "físico" de memória do meio geológico.

Resultados e Achados Numéricos
A análise numérica, apoiada por representações gráficas, produz diversas observações específicas:

• Magnitudes de Deslocamento: A mudança no deslocamento superficial próximo às falhas é da ordem de $10^{-3}$ m/ano. Isso se alinha bem com dados geofísicos observados (por exemplo, relatórios do USGS de 2–6 mm/ano), enquanto modelos anteriores de ordem inteira previam valores na faixa de $10^{-5}$, mostrando uma discrepância com as observações do mundo real.
• Influência da Geometria da Falha:
  • No Caso-I (duas falhas infinitas), a magnitude do deslocamento próximo à segunda falha é maior do que próxima à primeira.
  • No Caso-II (falhas infinita e finita), a magnitude do deslocamento próximo à falha finita é menor do que próxima à falha infinita.
  • O pico de deslocamento ocorre próximo à traçagem da falha e diminui conforme a distância aumenta, tendendo eventualmente a zero.
• Efeito da Ordem Fracionária ($\alpha$): O estudo constata que a taxa de acumulação e liberação de tensão é altamente dependente de $\alpha$. Valores menores de $\alpha$ (indicando efeitos de memória mais fortes) resultam em magnitudes maiores de acumulação de tensão de cisalhamento e uma relaxação mais lenta. À medida que $\alpha$ aumenta em direção a 1 (aproximando-se do comportamento de ordem inteira), a taxa de liberação de tensão aumenta.
• Força Motriz e Distribuição de Tensão: A força motriz total através das falhas apresenta uma liberação rápida após o início da fluência. A magnitude dessa liberação varia com a inclinação da falha e a velocidade de fluência. Mapas de contorno revelam que a acumulação de tensão é mais significativa perto da superfície livre e das zonas de descontinuidade da falha, enquanto a redução de tensão é proeminente em profundidades correspondentes às localizações das falhas.
• Análise de Deformação: As variações de deformação de cisalhamento ($e_{12}$) encontram-se na faixa de $10^{-4}$, consistentes com dados de eventos do mundo real como o terremoto de Tohoku-oki de 2011.

Significância e Alegações
Os autores alegam que seu modelo fracionário oferece uma representação mais geral e precisa da dinâmica de interação de falhas em comparação com os modelos tradicionais de ordem inteira. Especificamente:

• A abordagem fracionária captura com sucesso os efeitos de memória e a não linearidade dos meios viscoelásticos, que são frequentemente negligenciados em modelos padrão.
• As previsões do modelo para taxas de deslocamento ($10^{-3}$) mostram forte concordância com dados geofísicos observados, sugerindo que o cálculo fracionário é uma ferramenta superior para prever o comportamento sísmico e a deformação do subsolo.
• Os achados indicam que a interação entre falhas influencia o campo de tensão de forma mais extensa do que uma única falha, e que a ordem fracionária $\alpha$ é um parâmetro crítico na determinação da taxa de liberação de tensão.

O artigo conclui que esses modelos teóricos aprimoram a compreensão da deformação do subsolo e seu papel na influência da dinâmica de falhas e na geração de terremotos, fornecendo um arcabouço robusto para analisar períodos aseísmicos em zonas sismicamente ativas.