Seismic analytical reliability analysis of triangular rock wedge slope stability

Resumo Técnico: Análise de Confiabilidade Analítica Sísmica da Estabilidade de Cunha de Rocha Triangular

Definição do Problema
A estabilidade de taludes rochosos em regiões montanhosas é inerentemente estocástica devido às incertezas nos parâmetros geotécnicos, configurações geométricas e carregamento sísmico. Embora existam métodos determinísticos, eles frequentemente falham em quantificar a probabilidade de falha associada a essas variáveis. Pesquisas anteriores abordaram as incertezas na geometria das fraturas ou nas características do material rochoso separadamente, ou utilizaram métodos de simulação numérica como a Simulação de Monte Carlo (MCS) que, embora robustos, podem ser computacionalmente dispendiososos. Uma lacuna significativa permanece na investigação simultânea de incertezas nos parâmetros geométricos, sísmicos e de rocha utilizando uma abordagem analítica que evite os altos custos computacionais das simulações numéricas. Este estudo aborda a necessidade de um método analítico eficiente e preciso para avaliar a confiabilidade de cunhas de rocha triangulares sob condições de deslizamento planar submetidas a forças sísmicas.

Metodologia
O estudo propõe a aplicação do método das Variáveis Aleatórias Conjuntamente Distribuídas (JDRV) para derivar a função de densidade de probabilidade (PDF) do fator de segurança para taludes rochosos. A metodologia envolve os seguintes componentes:

• Estrutura Matemática: O método JDRV utiliza teoremas relativos à transformação de variáveis aleatórias (Teoremas 1–4) para integrar analiticamente as PDFs das variáveis de entrada dentro do modelo de estabilidade do talude. Isso permite a derivação direta da PDF do fator de segurança sem a necessidade de repetidas amostragem aleatória.
• Modelagem Estocástica: Seis parâmetros de entrada são tratados como variáveis aleatórias:
  • Parâmetros Geométricos: Coordenadas da cunha triangular ($y, z$) e o ângulo de superfície de ruptura resultante ($\beta$) e comprimento ($L_{ab}$).
  • Parâmetros da Rocha: Coesão aparente ($c$), ângulo de atrito ($\phi$) e peso específico ($\gamma$).
  • Parâmetro Sísmico: Coeficiente pseudo-estático horizontal ($k_h$).
  • Distribuições: Os parâmetros de rocha e geométricos são modelados usando distribuições normais truncadas, enquanto o coeficiente sísmico é modelado usando uma distribuição exponencialmente truncada.
• Análise Comparativa: Os resultados do JDRV são comparados com outros dois métodos:
  1. Simulação de Monte Carlo (MCS): Utilizada como padrão de referência com 108 simulações para gerar PDFs empíricas e Funções de Distribuição Acumulada (CDFs).
  2. Primeiro Momento de Segunda Ordem (FOSM): Um método aproximado usando expansão de série de Taylor para estimar a média e a variância do fator de segurança.
• Análise de Sensibilidade: Duas abordagens são empregadas para avaliar a influência dos parâmetros:
  1. Aumentar o Coeficiente de Variação (COV) de todos os seis parâmetros estocásticos.
  2. Variar os valores médios de parâmetros específicos (ex: ângulo de atrito, coordenadas geométricas) mantendo outros constantes.

Principais Resultados

• Precisão e Concordância: O estudo demonstra que o método JDRV produz resultados estatisticamente similares ao método FOSM e mostra excelente concordância com o MCS. Especificamente, as CDFs geradas pelo JDRV e pelo MCS exibem um "ajuste perfeito" em todos os três cenários testados (apenas parâmetros de rocha; rocha + sísmico; e parâmetros de rocha + sísmico + geométricos).
• Eficiência Computacional: Um achado primário é a redução significativa no tempo e custo computacional oferecida pelo JDRV em comparação ao MCS. Enquanto o MCS requer milhares de iterações para convergir, o JDRV atinge a convergência com menos iterações devido à sua natureza analítica.
• Desempenho Estatístico: O método JDRV fornece resultados estáveis de média e desvio padrão independentemente do número de simulações, oferecendo maior precisão na estimativa da PDF do fator de segurança em comparação ao ajuste exigido para os outputs do MCS (ex: distribuição epsilon-skew-normal).
• Impacto das Variáveis Estocásticas: Quando todas as variáveis (geométricas, sísmicas e de parâmetros de rocha) são tratadas como estocásticas simultaneamente, a análise revela um índice de confiabilidade menor e uma maior probabilidade de falha em comparação com análises que fixam certas variáveis. Isso ressalta o risco de subestimar a probabilidade de falha ao ignorar incertezas geométricas ou sísmicas.
• Achados de Sensibilidade: A análise de sensibilidade identifica o ângulo de atrito das superfícies de deslizamento como o parâmetro mais influente que afeta a estabilidade do talude rochoso.

Significância e Alegações
O artigo alega que o método JDRV serve como uma ferramenta analítica superior para a análise de confiabilidade de cunhas de rocha sob condições sísmicas. Sua significância reside na capacidade de:

1. Contabilizar simultaneamente as incertezas na geometria, sismicidade e mecânica das rochas dentro de um único framework analítico, uma capacidade frequentemente ausente em estudos anteriores.
2. Reduzir significativamente os encargos computacionais em comparação ao MCS, mantendo alta precisão na determinação das funções de densidade de probabilidade do fator de segurança.
3. Fornecer uma alternativa robusta a métodos aproximados como o FOSM ao capturar a forma completa da PDF, em vez de apenas os dois primeiros momentos.

Os autores enfatizam que tratar todas as variáveis relevantes como parâmetros estocásticos é crítico para uma avaliação de risco realista, pois negligenciar essas incertezas pode levar a uma superestimação da estabilidade do talude. O estudo conclui que o método JDRV, ao combinar relações fundamentais com equações de Jacobiano, oferece uma abordagem confiável, eficiente e precisa para a análise de estabilidade sísmica de cunhas de rocha triangulares.