Calculadora de Triangulação de Epicentro de Sismo

O Funcionamento da Triangulação de Epicentros

A localização do epicentro de um sismo é um dos procedimentos fundamentais da sismologia prática. A "Earthquake Epicenter Triangulation Calculator" é uma ferramenta pedagógica online concebida para estimar a localização do epicentro de um sismo num modelo de Terra esférica, recorrendo a observações registadas por um conjunto de três a oito estações sismográficas.

Ao introduzir as coordenadas geográficas de cada estação e o respetivo intervalo de tempo entre a chegada das ondas S e P (intervalo S−P), o utilizador consegue calcular a posição que melhor se ajusta aos dados observados. O cálculo baseia-se num tratamento matemático por mínimos quadrados que minimiza as discrepâncias de distância em todas as estações sismográficas. A ferramenta apresenta a estimativa do epicentro, a força da geometria da rede de estações e os resíduos individuais de cada ponto de observação, revelando o nível de consistência dos dados introduzidos.

O processamento de todos os dados introduzidos — incluindo as coordenadas das estações, os tempos S−P e os resultados calculados — ocorre localmente no dispositivo do utilizador. Nenhuma informação é enviada para servidores externos, garantindo que os dados permanecem estritamente no seu navegador.

---

Parâmetros de Entrada e Validação de Dados

Para efetuar a modelação do epicentro, o utilizador deve configurar os parâmetros do modelo de propagação e fornecer as observações de cada estação sismográfica. O sistema permite a utilização de um mínimo de 3 estações e um máximo de 8 estações. Se tentar realizar o cálculo fora destes limites, a ferramenta apresenta a mensagem de erro: Utilize entre três e oito estações.

1. Parâmetros do Modelo de Propagação

Estes parâmetros definem as propriedades físicas do meio de propagação uniforme e a incerteza associada às leituras temporais:

• Velocidade da onda P (fieldPVelocity): Deve ser superior a 0 e não superior a 20 km/s. Deve também ser estritamente superior à velocidade da onda S. Caso contrário, é exibido o erro: A velocidade da onda P deve ser superior à velocidade da onda S. ou Vₚ/Vₛ deve ser superior a 0 e não superior a 20 km/s. (onde Vₚ/Vₛ corresponde ao nome do campo em falta).
• Velocidade da onda S (fieldSVelocity): Deve ser superior a 0 e não superior a 20 km/s, e inferior à velocidade da onda P.
• Incerteza temporal (1σ) (fieldTimingSigma): Representa o desvio padrão na leitura dos tempos de chegada. Deve situar-se entre 0 e 60 segundos. Valores fora deste intervalo geram o erro: A incerteza temporal deve ser superior a 0 e não superior a 60 segundos..

2. Observações das Estações

Para cada estação sismográfica (identificada por um nome ou etiqueta), devem ser introduzidos os seguintes dados:

• Latitude: Deve situar-se no intervalo de −90° a 90°. Se violado, surge o erro: Estação 1: a latitude deve situar-se entre −90° e 90°. (onde 1 é o nome da estação).
• Longitude: Deve situar-se no intervalo de −180° a 180°. Se violado, surge o erro: Estação 1: a longitude deve situar-se entre −180° e 180°. (onde 1 é o nome da estação).
• Intervalo S−P: O tempo decorrido entre a chegada da onda primária e da onda secundária. Deve ser superior a 0 e não superior a 2000 segundos. O incumprimento gera o erro: Estação 1: insira um intervalo S−P superior a 0 e não superior a 2000 segundos. (onde 1 é o nome da estação).

Se for introduzido um valor não numérico em qualquer campo, o sistema indica: input: “abc” não é um número. (onde input identifica o local do erro e abc o texto inválido introduzido). No caso de introduzir coordenadas idênticas para estações diferentes, o sistema impede o cálculo com o erro: As estações 1 e 2 têm as mesmas coordenadas. Utilize localizações de estações distintas. (onde 1 e 2 representam os nomes das estações duplicadas). Se as estações estiverem dispostas de forma linear ou num agrupamento demasiado estreito, impossibilitando a convergência matemática, é exibido o aviso: Estas estações não delimitam uma localização bidimensional estável. Adicione uma estação afastada da linha ou do agrupamento atual.. Erros de processamento numérico extremo resultam na mensagem: Um valor ou resultado intermédio excede o intervalo numérico suportado..

O utilizador pode clicar em Carregar exemplo de Tóquio para preencher automaticamente a ferramenta com dados de teste, ou em Limpar para repor todos os campos.

A Física e a Matemática do Ajuste do Epicentro

As ondas P (compressionais) viajam a uma velocidade superior à das ondas S (de corte ou cisalhamento). Consequentemente, à medida que a frente de onda se afasta do foco do sismo, o intervalo temporal entre a chegada de ambas as ondas aumenta proporcionalmente à distância percorrida.

Sob a premissa de velocidades constantes adotada por este modelo pedagógico, a diferença de tempo de chegada (Δ t) satisfaz a relação física:

Δ t = D (1 / Vₛ - 1 / Vₚ)

Ao isolar a distância (D), obtemos a fórmula de aproximação utilizada pela ferramenta para calcular o raio de distância de cada estação ao epicentro:

Aproximação da distância: D = Δ t ÷ (1 / Vₛ - 1 / Vₚ)

Fórmula e substituição:
- Distância de superfície: d = R × ângulo central de haversine, com R = 6371 km
- Ajuste do epicentro: minimizar Σ(dᵢ − Dᵢ)² em todas as estações

O cálculo do raio esférico assume a Terra como uma esfera perfeita com um raio constante (R) de 6371 km. Como as observações reais contêm erros de leitura e o modelo de velocidade simplificado introduz desvios, as circunferências de distância calculadas para cada estação raramente se cruzam num único ponto exato. Para resolver esta inconsistência, a ferramenta aplica um ajuste por mínimos quadrados, cujo objetivo é encontrar as coordenadas geográficas que minimizam a soma dos quadrados das diferenças entre a distância esférica calculada sobre a superfície (dᵢ) e a distância teórica obtida pelo intervalo S−P (Dᵢ) para todas as estações ativas:

Minimizar Σ (dᵢ - Dᵢ)²

---

Interpretação dos Resultados e Resíduos

Após o processamento, a ferramenta apresenta os resultados estruturados em secções claras:

Ajuste das Estações

Apresenta uma tabela detalhada com as seguintes colunas para cada estação:

• Estação (columnStation)
• Distância S−P (columnObserved): O raio calculado a partir do intervalo de tempo.
• Distância esférica (columnFitted): A distância real medida sobre a superfície da esfera entre a estação e o epicentro estimado pelo ajuste.
• Resíduo (columnResidual): Calculado como a distância esférica ajustada menos a distância inferida pelo intervalo S−P (dᵢ - Dᵢ). Um resíduo positivo indica que o epicentro calculado se encontra mais afastado da estação do que o raio S−P sugeria; um resíduo negativo indica que está mais próximo.

Qualidade da Solução

A ferramenta classifica a qualidade do ajuste obtido através de indicadores qualitativos e quantitativos:

• Qualidade do ajuste de distância: Classificado como ajuste de distância próximo, ajuste de distância misto ou ajuste de distância fraco.
• Qualidade da geometria das estações: Classificado como geometria de estações forte, geometria de estações moderada ou geometria de estações fraca. Uma boa distribuição angular das estações em redor do epicentro garante uma geometria forte, enquanto estações concentradas num único quadrante fragilizam a solução espacial.
• Resíduo de distância RMS e Maior resíduo de distância.
• Incerteza de posição 1σ aproximada e Incerteza de distância 1σ por estação.

Abaixo dos resultados, é apresentado um Esquema da geometria das estações que ilustra graficamente a disposição relativa das estações e do epicentro calculado, acompanhado da legenda: Apenas geometria relativa — sem mapa, fronteiras ou dados de sismos em tempo real..

---

Limitações do Modelo Pedagógico

Este simulador foi desenvolvido para fins didáticos e de demonstração em sala de aula. Os cálculos assumem um sismo de origem superficial (foco na superfície) e velocidades de propagação de onda constantes e homogéneas em toda a crosta.

Na sismologia profissional, a determinação de catálogos sísmicos oficiais utiliza modelos de velocidade tridimensionais e estratificados, que consideram a profundidade focal do sismo (hipocentro), a curvatura das trajetórias dos raios sísmicos através das diferentes camadas da Terra, desvios de tempo das estações e incertezas associadas à identificação de fases adicionais. A incerteza de posição calculada nesta ferramenta propaga estritamente o valor de incerteza temporal (1σ) introduzido pelo utilizador, ignorando erros sistemáticos do modelo de velocidade ou desvios instrumentais.

---

Perguntas Frequentes

Porque são necessárias três estações?
Uma única medição S−P fornece apenas um raio de distância, definindo uma circunferência de possíveis localizações em redor da estação. Duas estações geram duas circunferências que se cruzam em dois pontos distintos. Uma terceira estação é necessária para eliminar a ambiguidade e identificar o ponto de interseção único. A adição de quatro ou mais estações permite aplicar o método dos mínimos quadrados, tornando visíveis os resíduos e aumentando a precisão do ajuste face a erros de leitura.

A distância S−P é a distância de superfície exata?
Não. A equação de velocidade constante estima uma trajetória de propagação comum a partir da diferença no tempo de chegada. Tratar esse raio como distância de superfície esférica assume uma fonte superficial e uma estrutura de velocidades 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 a partir do intervalo S−P. Um valor positivo significa que o ponto ajustado está mais afastado dessa estação do que o seu raio S−P; um valor negativo significa que está mais próximo. Resíduos grandes e mistos indicam que as circunferências não partilham uma interseção limpa sob as velocidades escolhidas.

O que está incluído na incerteza de posição?
Esta propaga o valor único de incerteza temporal 1σ através do factor de distância S−P e da geometria das estações perto do ponto ajustado. Não inclui erros de modelo de velocidade incorreto, profundidade focal, fases mal identificadas, desvios de relógio ou a estrutura tridimensional da Terra, pelo que não constitui uma garantia de precisão completa.