O cálculo do epicentro de um sismo é um dos procedimentos fundamentais da sismologia prática. A Calculadora de Triangulação de Epicentro de Terremoto é uma ferramenta educacional projetada para estimar a localização do epicentro de um terremoto em um modelo de Terra esférica utilizando observações de três a oito estações de registro.
Ao inserir as coordenadas geográficas e os intervalos de tempo entre as ondas S e P (S−P) para cada estação, além de definir parâmetros de velocidade das ondas e incerteza temporal, o sistema calcula a posição de melhor ajuste. O processo utiliza uma abordagem matemática de mínimos quadrados para minimizar as discrepâncias de distância entre todas as estações, exibindo a qualidade da geometria e os resíduos individuais de cada ponto de observação.
A Física das Ondas Sísmicas e o Intervalo S−P
Durante um terremoto, a liberação repentina de energia gera diferentes tipos de ondas elásticas que se propagam pelo interior da Terra. As duas principais ondas de corpo são:
- Ondas P (Primárias ou Compressionais): São ondas longitudinais que se movem mais rapidamente, comprimindo e dilatando o material rochoso na mesma direção de propagação da onda.
- Ondas S (Secundárias ou de Cisalhamento): São ondas transversais mais lentas, que deslocam o solo perpendicularmente à direção de propagação e não se propagam através de fluidos.
Como as ondas P viajam a uma velocidade superior à das ondas S, a diferença de tempo entre a chegada da onda P e a chegada subsequente da onda S (o intervalo S−P, representado por Δt) aumenta progressivamente à medida que a distância em relação ao foco do terremoto se torna maior.
Sob a premissa de velocidades constantes na crosta, a relação entre a distância da estação ao epicentro (D) e o tempo de percurso obedece à seguinte equação física:
Δt = D(1/Vs − 1/Vp)
Ao isolar a distância, obtemos a fórmula de aproximação utilizada pela ferramenta para determinar o raio de busca a partir de cada estação:
D = Δt ÷ (1/Vs − 1/Vp)
A Geometria da Triangulação e Mínimos Quadrados
Uma única estação sismográfica que registra um intervalo S−P consegue determinar apenas a distância radial até o evento, definindo um círculo de posições possíveis ao seu redor. Duas estações geram dois círculos que se interceptam em, no máximo, dois pontos distintos. Para resolver essa ambiguidade e identificar um único ponto geográfico, é necessária uma terceira estação.
Na sismologia real, as observações contêm pequenas imprecisões de leitura e variações geológicas. Consequentemente, os círculos de distância de três ou mais estações raramente se cruzam em um único ponto perfeito. Para resolver essa inconsistência, a ferramenta aplica uma otimização por mínimos quadrados. O algoritmo busca coordenadas geográficas (latitude e longitude) que minimizem a soma dos quadrados das diferenças entre a distância esférica calculada no modelo e a distância inferida pelo intervalo S−P de cada estação:
minimizar Σ(dᵢ − Dᵢ)² em todas as estações
A distância de superfície (d) entre as coordenadas estimadas do epicentro e cada estação é calculada sobre uma esfera de raio médio constante (R) de 6.371 km, utilizando a relação do ângulo central de haversine:
d = R × ângulo central de haversine, com R = 6371 km
Parâmetros de Entrada e Validação
Para realizar o ajuste do epicentro, o usuário deve fornecer dados consistentes dentro dos limites operacionais do sistema. Caso contrário, mensagens de erro específicas serão exibidas para orientar a correção.
Observações das Estações (Mínimo de 3 e Máximo de 8)
- Nome da estação
- Latitude: Deve estar no intervalo de −90° a 90°. Se violado, exibe:
Estação station: a latitude deve ser de −90° a 90°. - Longitude: Deve estar no intervalo de −180° a 180°. Se violado, exibe:
Estação station: a longitude deve ser de −180° a 180°. - Intervalo S−P: Deve ser maior que 0 e não superior a 2.000 segundos. Se violado, exibe:
Estação station: insira um intervalo S−P maior que 0 e não superior a 2000 segundos.
Parâmetros do Modelo de Propagação
- Velocidade da onda P: Deve ser maior que 0 e não superior a 20 km/s. Deve ser estritamente maior que a velocidade da onda S.
- Velocidade da onda S: Deve ser maior que 0 e não superior a 20 km/s. Deve ser estritamente menor que a velocidade da onda P.
- Incerteza de tempo (1σ): Deve ser maior que 0 e não superior a 60 segundos.
Regras de Consistência e Erros de Interface
- Quantidade de estações: Se o cálculo for iniciado com menos de 3 ou mais de 8 estações, o sistema exibe:
Use entre três e oito estações. - Valores não numéricos: Se um caractere inválido for inserido, exibe:
where: “token” não é um número. - Ordem de velocidade incorreta: Se a velocidade S for maior ou igual à velocidade P, exibe:
A velocidade da onda P deve ser maior que a velocidade da onda S. - Velocidade fora dos limites: Se os limites de velocidade forem violados, exibe:
field deve ser maior que 0 e não superior a 20 km/s. - Incerteza fora dos limites: Se a incerteza temporal violar as regras, exibe:
A incerteza de tempo deve ser maior que 0 e não superior a 60 segundos. - Coordenadas duplicadas: Se duas estações forem posicionadas no mesmo local, exibe:
As estações station e other têm as mesmas coordenadas. Use localizações de estações distintas. - Geometria instável: Se as estações estiverem dispostas em linha reta ou agrupadas de forma muito estreita, impossibilitando a convergência bidimensional, exibe:
Estas estações não limitam uma localização bidimensional estável. Adicione uma estação longe da linha ou do agrupamento atual. - Estouro numérico: Se os cálculos excederem a capacidade computacional, exibe:
Um valor ou resultado intermediário excede o intervalo numérico suportado.
O usuário pode clicar em Carregar exemplo de Tóquio para preencher a ferramenta com dados de teste estruturados ou em Limpar para redefinir todos os campos de entrada.
Resultados e Interpretação de Saída
Após o processamento local, a interface apresenta os resultados estruturados nas seguintes seções:
Solução do Epicentro
Apresenta as coordenadas geográficas calculadas para o Epicentro estimado. A qualidade do ajuste é classificada em três níveis de fidelidade de distância: ajuste de distância próximo, ajuste de distância misto ou ajuste de distância ruim.
A distribuição espacial das estações é avaliada e classificada como geometria de estação forte, geometria de estação moderada ou geometria de estação fraca. Métricas estatísticas adicionais incluem o Resíduo de distância RMS, o Maior resíduo de distância, a Incerteza de posição aproximada 1σ e a Incerteza de distância 1σ por estação.
Esquema da Geometria das Estações
Um diagrama vetorial exibe a disposição espacial relativa das estações e do epicentro calculado. Esta área exibe a legenda obrigatória: Apenas geometria relativa — sem mapa, fronteiras ou dados de terremotos em tempo real..
Ajuste das Estações
Uma tabela detalhada contendo as seguintes colunas para cada estação inserida:
- Estação
- Distância S−P (o raio calculado a partir do tempo de chegada)
- Distância esférica (a distância real do epicentro ajustado até a estação sobre a esfera)
- Resíduo (a diferença matemática entre a distância esférica e a distância S−P)
Fórmula e Substituição
Esta seção detalha a derivação matemática passo a passo para fins didáticos, exibindo as seguintes equações e substituições:
Aproximação de distância: D = Δt ÷ (1/Vs − 1/Vp)Distância de superfície: d = R × ângulo central de haversine, com R = 6371 kmAjuste do epicentro: minimizar Σ(dᵢ − Dᵢ)² em todas as estações- O cálculo do fator de distância:
Fator de distância = vp × vs ÷ (vp − vs) = factor km por segundo de S−P - As derivações individuais de distância para cada estação:
station: interval s × factor km/s = distance km - As coordenadas finais e o resíduo médio:
Melhor ajuste esférico = latitude°, longitude°; resíduo RMS = rms km - A propagação de erro baseada na geometria:
A incerteza de tempo 1σ de timing s resulta em uma incerteza de distância 1σ de distance km; a geometria das estações propaga isso para cerca de position km no eixo horizontal maior.
Método, data e escopo
O painel identifica o método como Ajuste esférico de mínimos quadrados e informa que ele foi revisado em 6 de agosto de 2026. Antes de recalcular, o status pede que as observações sejam editadas. Após o ajuste, ele informa quantas estações foram usadas e, quando os resíduos ou a geometria exigem atenção, exibe um aviso explícito para revisar o resultado.
Limitações do Modelo Didático
Este simulador foi desenvolvido exclusivamente para fins de ensino e demonstração geométrica. Ele adota simplificações físicas importantes que o diferenciam dos sistemas profissionais de monitoramento sísmico:
- Modelo de Velocidade Uniforme: O cálculo assume que as velocidades das ondas P e S são constantes em toda a trajetória. Na Terra real, as velocidades variam significativamente com a profundidade e a composição mineralógica da crosta e do manto.
- Terremoto Raso: O modelo assume que o foco do terremoto (hipocentro) está localizado na superfície da Terra (profundidade zero). Terremotos reais ocorrem em profundidades variáveis, o que exige modelos tridimensionais e curvas de tempo de percurso complexas para determinar a profundidade focal.
- Trajetória Linear dos Raios: O cálculo assume que as ondas viajam em linha reta ao longo da superfície esférica. Na realidade, as ondas sísmicas sofrem refração e reflexão nas interfaces geológicas, gerando caminhos de raio curvos.
Privacidade e Processamento de Dados
A privacidade dos dados inseridos é mantida localmente. Todo o processamento matemático, validação de coordenadas e geração de gráficos ocorrem diretamente no navegador do usuário. Nenhuma coordenada de estação, tempo de chegada ou resultado calculado é enviado para servidores externos ou armazenado fora do dispositivo local.
Perguntas Frequentes (FAQ)
Por que são necessárias três estações?
Uma medição S−P fornece um círculo de distância, não uma direção. Dois círculos podem se encontrar em dois pontos. Uma terceira estação normalmente resolve a escolha, enquanto uma quarta ou mais estações tornam a divergência visível através dos resíduos. Estações distribuídas ao redor da fonte limitam tanto a latitude quanto a longitude melhor do que estações ao longo de uma única linha.
A distância S−P é a distância exata de superfície?
Não. A equação de velocidade constante estima um caminho de propagação comum a partir da diferença do tempo de chegada. Tratar esse raio como distância de superfície esférica assume uma fonte rasa e uma estrutura de velocidade simples. As curvas de tempo de percurso reais dependem da profundidade, da distância e do material atravessado por cada raio.
O que significa um resíduo positivo ou negativo?
O resíduo é a distância esférica ajustada menos a distância inferida do S−P. Um valor positivo significa que o ponto ajustado está mais distante daquela estação do que o seu raio S−P; um valor negativo significa que está mais próximo. Resíduos grandes e mistos mostram que os círculos não compartilham uma interseção limpa sob as velocidades escolhidas.
O que está incluído na incerteza de posição?
Ela propaga o valor único de incerteza de tempo 1σ através do fator de distância S−P e da geometria das estações perto do ponto ajustado. Não inclui um modelo de velocidade incorreto, profundidade focal, fases identificadas incorretamente, desvios de relógio ou a estrutura tridimensional da Terra, portanto não é uma declaração de precisão completa.