Le fonctionnement de la triangulation d'épicentre
Le calculateur de triangulation d'épicentre de séisme est un outil pédagogique gratuit conçu pour estimer l'épicentre d'un séisme sur un modèle de Terre sphérique en utilisant les observations de trois à huit stations d'enregistrement. En saisissant les coordonnées géographiques et les intervalles d'arrivée des ondes S−P pour chaque station, ainsi que les vitesses des ondes et l'incertitude temporelle définies par l'utilisateur, l'outil calcule la position de l'épicentre offrant le meilleur ajustement.
Pour y parvenir, le calculateur applique un traitement par moindres carrés afin de minimiser les écarts de distance sur l'ensemble des stations. Les résultats affichent l'épicentre estimé, la qualité de la géométrie des stations et les résidus individuels de chaque station pour évaluer la cohérence des observations.
Pour faciliter la prise en main, l'utilisateur peut cliquer sur Charger l'exemple de Tokyo afin de remplir automatiquement les champs avec des données d'exemple, ou sur Effacer pour réinitialiser l'ensemble des entrées.
Paramètres d'entrée et de configuration
Le calcul repose sur deux ensembles de paramètres saisis par l'utilisateur:
1. Observations des stations (de 3 à 8 stations)
Pour chaque station, les champs suivants doivent être renseignés:
- Nom de la station (identifiant textuel).
- Latitude: Doit être comprise entre −90° et 90°.
- Longitude: Doit être comprise entre −180° et 180°.
- Intervalle S−P: Doit être supérieur à 0 et d'au plus 2000 secondes.
2. Modèle de propagation
Ces paramètres définissent les propriétés physiques du milieu traversé par les ondes:
- Vitesse de l'onde P: Doit être supérieure à 0, d'au plus 20 km/s, et strictement supérieure à la vitesse de l'onde S.
- Vitesse de l'onde S: Doit être supérieure à 0, d'au plus 20 km/s, et strictement inférieure à la vitesse de l'onde P.
- Incertitude de temps (1σ): Doit être supérieure à 0 et d'au plus 60 secondes.
Formules mathématiques et modélisation
Le calculateur utilise un modèle simplifié où la Terre est assimilée à une sphère parfaite de rayon constant R = 6371 km. Les calculs s'effectuent selon les étapes et formules suivantes, présentées dans la section Formule et substitution de l'outil:
Approximation de la distance
Les ondes P (primaires/compression) voyagent plus rapidement que les ondes S (secondaires/cisaillement). L'écart temporel entre leur arrivée croît proportionnellement à la distance de la station par rapport à la source. Sous l'hypothèse de vitesses constantes, la distance de propagation D est calculée à partir de la différence de temps d'arrivée Δ t par la formule:
Approximation de la distance: D = Δt ÷ (1/Vs − 1/Vp)
Le facteur de distance est ainsi dérivé de la manière suivante:
Facteur de distance = vp × vs ÷ (vp − vs) = factor km par seconde d'intervalle S−P
Pour chaque station, la distance individuelle est obtenue par:
station: interval s × factor km/s = distance km
Distance de surface et ajustement
La distance de surface d entre les coordonnées d'une station et l'épicentre sur la sphère est calculée via la formule de la trigonométrie sphérique:
Distance de surface: d = R × angle central haversine, avec R = 6371 km
L'ajustement de l'épicentre consiste à trouver les coordonnées géographiques qui minimisent la somme des carrés des écarts entre la distance de surface théorique (dᵢ) et la distance déduite de l'intervalle S−P (Dᵢ) pour toutes les stations:
Ajustement de l'épicentre: minimiser Σ(dᵢ − Dᵢ)² sur toutes les stations
Le résultat final affiche les coordonnées optimales et le résidu quadratique moyen (RMS):
Meilleur ajustement sphérique = latitude°, longitude°; résidu RMS = rms km
Analyse des résultats et de la géométrie
L'outil fournit plusieurs indicateurs pour évaluer la fiabilité de la solution géométrique:
- Rayon calculé: La distance théorique estimée pour chaque station à partir de son intervalle S−P.
- Ajustement de distance: Qualifié selon la précision globale en ajustement de distance étroit, ajustement de distance mixte ou ajustement de distance médiocre.
- Géométrie des stations: Évaluée en géométrie des stations forte, géométrie des stations modérée ou géométrie des stations faible. Une bonne répartition spatiale des stations autour de l'événement est essentielle pour contraindre efficacement la position.
- Résidus de distance: Le tableau d'ajustement affiche pour chaque station la Distance S−P calculée, la Distance sphérique finale et le Résidu (différence entre la distance sphérique ajustée et la distance S−P). Un résidu positif indique que l'épicentre ajusté est plus éloigné de la station que son rayon S−P théorique; un résidu négatif indique qu'il est plus proche.
- Incertitude de position: L'outil calcule la propagation de l'incertitude temporelle saisie par l'utilisateur:
Une incertitude de temps 1σ de timing s donne une incertitude de distance 1σ de distance km; la géométrie des stations propage cela à environ position km sur le grand axe horizontal.
Un schéma visuel simplifié est généré sous l'intitulé Schéma de la géométrie des stations, accompagné de la mention: « Géométrie relative uniquement — sans carte, frontières ni données de séisme en temps réel. »
Règles de validation et messages d'erreur
Pour garantir la stabilité des calculs numériques, l'outil applique des règles de validation strictes et affiche des messages d'erreur spécifiques en cas de saisie incorrecte:
| Règle / Cas limite | Message d'erreur affiché |
|---|---|
| Nombre de stations hors limites (< 3 ou > 8) | Utilisez entre trois et huit stations. |
| Saisie non numérique | where: « token » n'est pas un nombre. |
| Latitude hors limites | Station station: la latitude doit être comprise entre −90° et 90°. |
| Longitude hors limites | Station station: la longitude doit être comprise entre −180° et 180°. |
| Intervalle S−P hors limites | Station station: saisissez un intervalle S−P supérieur à 0 et d'au plus 2000 secondes. |
| Vitesse d'onde hors limites | field doit être supérieure à 0 et d'au plus 20 km/s. |
| Ordre des vitesses incorrect | La vitesse de l'onde P doit être supérieure à la vitesse de l'onde S. |
| Incertitude de temps hors limites | L'incertitude de temps doit être supérieure à 0 et d'au plus 60 secondes. |
| Coordonnées de stations identiques | Les stations station et other ont les mêmes coordonnées. Utilisez des emplacements de stations distincts. |
| Alignement ou regroupement instable | Ces stations ne permettent pas de contraindre une position bidimensionnelle stable. Ajoutez une station éloignée de la ligne ou du groupe actuel. |
| Dépassement des capacités de calcul | Une valeur ou un résultat intermédiaire dépasse la plage de nombres prise en charge. |
Confidentialité et traitement des données
Toutes les coordonnées de stations, les temps d'intervalle S−P et les résultats de calcul restent sur cet appareil — rien n'est téléversé. Le traitement s'effectue localement dans le navigateur Web de l'utilisateur.
Foire aux questions (FAQ)
La distance S−P est-elle la distance de surface exacte?
Non. L'équation à vitesse constante estime un trajet de propagation commun à partir de la différence de temps d'arrivée. Traiter ce rayon comme une distance de surface sphérique suppose une source peu profonde et une structure de vitesse simple. Les courbes de temps de trajet réelles dépendent de la profondeur, de la distance et des matériaux traversés par chaque rayon.
Pourquoi trois stations sont-elles nécessaires?
Une seule mesure S−P donne un cercle de distance, pas une direction. Deux cercles peuvent se croiser en deux points. Une troisième station permet normalement de lever l'ambiguïté, tandis qu'une quatrième station ou plus rend l'écart visible à travers les résidus. Des stations réparties autour de la source contraignent mieux la latitude et la longitude que des stations alignées.
Que signifie un résidu positif ou négatif?
Le résidu est la distance sphérique ajustée moins la distance déduite de l'intervalle S−P. Une valeur positive signifie que le point ajusté est plus éloigné de cette station que son rayon S−P; une valeur négative signifie qu'il est plus proche. Des résidus mixtes élevés indiquent que les cercles ne partagent pas une intersection nette avec les vitesses choisies.
Qu'est-ce qui est inclus dans l'incertitude de position?
Elle propage la valeur unique de l'incertitude de temps 1σ à travers le facteur de distance S−P et la géométrie des stations près du point ajusté. Elle n'inclut pas les erreurs de modèle de vitesse, la profondeur focale, les erreurs d'identification des phases, les biais d'horloge des stations ou la structure tridimensionnelle de la Terre, et ne constitue donc pas une mesure de précision complète.