Calcul de vitesse de propagation des ondes P et S
Calculez rapidement les vitesses sismiques à partir d’une distance et de temps d’arrivée, ou à partir des propriétés mécaniques du milieu. Cet outil est conçu pour l’analyse pédagogique, géophysique et géotechnique.
Résultats
Renseignez vos données puis cliquez sur Calculer pour obtenir les vitesses des ondes P et S, le rapport Vp/Vs et une visualisation graphique.
Guide expert du calcul de vitesse de propagation des ondes P et S
Le calcul de la vitesse de propagation des ondes P et S est un élément central de la sismologie, de la géophysique appliquée, de la géotechnique et même du contrôle non destructif des matériaux. Dans le sous-sol, les ondes sismiques transmettent des informations très précieuses sur la rigidité, la densité, la fracturation, la saturation en fluide et la nature des couches traversées. Une estimation correcte des vitesses permet non seulement de localiser les séismes, mais aussi d’interpréter l’architecture de la croûte, d’évaluer des sites d’infrastructure et d’améliorer les modèles de réponse sismique locale.
Les ondes P, dites ondes primaires, sont des ondes de compression. Elles provoquent des variations de volume dans le milieu et se déplacent généralement plus vite que les ondes S. Les ondes S, dites ondes secondaires, sont des ondes de cisaillement. Elles ne se propagent pas dans les fluides car ceux-ci ne supportent pas le cisaillement à long terme. Cette différence physique explique pourquoi la vitesse P est toujours supérieure à la vitesse S dans les solides, et pourquoi l’absence de propagation des ondes S dans les liquides est une signature importante en géophysique interne.
1. Formules fondamentales pour le calcul des vitesses sismiques
Deux grands cadres de calcul sont utilisés en pratique :
- Approche temps de trajet : si l’on connaît la distance parcourue et le temps d’arrivée d’une onde, la vitesse se calcule directement par le rapport distance sur temps.
- Approche mécanique du milieu : si l’on connaît les modules élastiques et la densité, on peut estimer la vitesse intrinsèque des ondes.
Dans le cas le plus direct :
- Mesurer ou estimer la distance entre la source et le capteur.
- Identifier le temps d’arrivée de l’onde P.
- Identifier le temps d’arrivée de l’onde S.
- Calculer Vp = d / tP et Vs = d / tS.
Dans le cas mécanique, les formules sont les suivantes :
- Vp = √((K + 4G/3) / ρ)
- Vs = √(G / ρ)
Où K représente le module de compressibilité, G le module de cisaillement et ρ la densité massique. Les unités doivent être cohérentes. Si K et G sont en pascals et ρ en kg/m³, alors la vitesse obtenue sera en m/s. En pratique, les géophysiciens travaillent souvent avec des valeurs en GPa pour les roches et convertissent ensuite en pascals avant le calcul.
2. Pourquoi les ondes P vont plus vite que les ondes S
La vitesse des ondes dépend de la résistance du milieu à la déformation et de sa densité. Une onde P exploite la rigidité volumique et la rigidité au cisaillement. Elle bénéficie donc à la fois du module K et du module G. À l’inverse, l’onde S ne dépend que du module de cisaillement G. Par conséquent, sauf cas pathologique ou erreur de mesure, la relation Vp > Vs est toujours vérifiée dans un solide isotrope élastique.
Le rapport Vp/Vs est particulièrement utile. Il renseigne sur la nature des matériaux, leur porosité, leur teneur en fluides et parfois leur degré d’altération. Dans les roches cristallines compactes, ce rapport est souvent proche de 1,7 à 1,9. Dans des formations plus saturées, argileuses ou altérées, il peut augmenter.
3. Exemples typiques de vitesses selon les matériaux
Les valeurs suivantes sont des ordres de grandeur fréquemment utilisés dans la littérature académique et professionnelle. Elles peuvent varier selon la pression, la température, la porosité, la fracturation et la saturation en eau.
| Matériau | Vp typique | Vs typique | Vp/Vs approximatif |
|---|---|---|---|
| Eau | ≈ 1480 m/s | 0 m/s | Non défini |
| Sédiments meubles saturés | 1500 à 2200 m/s | 100 à 500 m/s | 3,0 à 8,0 |
| Grès consolidé | 2500 à 4500 m/s | 1400 à 2500 m/s | 1,6 à 1,9 |
| Calcaire | 4000 à 6500 m/s | 2200 à 3400 m/s | 1,7 à 1,9 |
| Granite | 5000 à 6200 m/s | 3000 à 3600 m/s | 1,7 à 1,8 |
| Basalte | 5500 à 7000 m/s | 3200 à 4000 m/s | 1,7 à 1,9 |
On voit immédiatement l’intérêt opérationnel de ces vitesses. Une vitesse S très faible près de la surface peut révéler des sols meubles, potentiellement plus sensibles à l’amplification sismique. À l’inverse, des vitesses élevées suggèrent souvent un substratum compétent ou une roche compacte.
4. Exemple de calcul simple avec les temps d’arrivée
Supposons qu’une source sismique soit à 50 km d’une station. L’onde P arrive après 8,3 s et l’onde S après 14,3 s. Les vitesses sont alors :
- Vp = 50 000 / 8,3 ≈ 6024 m/s
- Vs = 50 000 / 14,3 ≈ 3497 m/s
- Vp/Vs ≈ 1,72
Ces résultats sont cohérents avec une roche consolidée de type cristallin ou volcanique dense. Une fois les vitesses obtenues, on peut estimer la nature du milieu, comparer les valeurs à des bases de référence et vérifier la qualité de la mesure. Si le rapport Vp/Vs est inhabituellement élevé, il faut envisager une saturation en fluide, une erreur de picking ou un fort contraste lithologique.
5. Exemple de calcul à partir des propriétés élastiques
Considérons un matériau avec K = 36 GPa, G = 45 GPa et ρ = 2700 kg/m³. Après conversion en pascals :
- K = 36 × 109 Pa
- G = 45 × 109 Pa
- ρ = 2700 kg/m³
On obtient alors :
- Vs = √(45 × 109 / 2700) ≈ 4082 m/s
- Vp = √((36 × 109 + 4 × 45 × 109 / 3) / 2700) ≈ 6086 m/s
Ces valeurs sont très plausibles pour une roche dure et compacte. Cette méthode est particulièrement utile en mécanique des roches, en acoustique des matériaux et en modélisation géophysique quand on dispose de paramètres issus d’essais de laboratoire.
6. Données de référence utilisées en géophysique et en ingénierie
Le tableau ci-dessous rassemble quelques statistiques et repères opérationnels largement cités dans les domaines sismique et géotechnique.
| Indicateur | Valeur ou plage courante | Interprétation pratique |
|---|---|---|
| Vp dans l’eau à température ambiante | ≈ 1480 m/s | Propagation compressive uniquement, pas d’onde S |
| Vs30 site classe B | 760 à 1500 m/s | Roche ou sol très dense selon les codes de construction |
| Vs30 site classe C | 360 à 760 m/s | Sol dense à très dense ou roche tendre |
| Vs30 site classe D | 180 à 360 m/s | Sol rigide mais potentiellement amplificateur |
| Rapport Vp/Vs roches compactes | 1,6 à 1,9 | Valeur attendue en milieu consolidé |
| Rapport Vp/Vs matériaux très saturés | > 2,0 dans certains cas | Peut signaler fluides, altération ou faible rigidité au cisaillement |
7. Application au zonage sismique et au génie civil
En ingénierie sismique, la vitesse des ondes S est souvent encore plus stratégique que celle des ondes P. La raison est simple : les mouvements de cisaillement gouvernent une partie importante de la réponse dynamique des structures. Le paramètre Vs30, qui représente la vitesse moyenne des ondes S dans les 30 premiers mètres du sol, est couramment utilisé pour classer les sites et estimer leur comportement lors d’un séisme. Des valeurs faibles de Vs30 sont généralement associées à une amplification plus marquée des secousses.
Dans les études de site, les ingénieurs croisent souvent plusieurs méthodes : réfraction sismique, MASW, downhole, crosshole, essais pressiométriques et données géologiques. Le calcul des vitesses P et S ne doit donc jamais être interprété isolément. Il doit être replacé dans un cadre géologique cohérent, avec contrôle des incertitudes, de la topographie, des effets de couche et des erreurs instrumentales.
8. Principales sources d’erreur dans le calcul des vitesses
- Distance mal estimée : la trajectoire réelle de l’onde peut être plus longue que la distance horizontale mesurée.
- Picking imprécis : identifier exactement le premier arrivage P ou S n’est pas toujours trivial sur un sismogramme bruité.
- Anisotropie : certaines roches ne propagent pas les ondes à la même vitesse dans toutes les directions.
- Hétérogénéité du milieu : un calcul distance sur temps donne souvent une vitesse moyenne, pas une vitesse locale.
- Unités incohérentes : une confusion entre km et m, ou entre GPa et Pa, produit des erreurs majeures.
9. Bonnes pratiques pour obtenir un calcul fiable
- Vérifier systématiquement les unités d’entrée avant de calculer.
- Comparer les résultats avec des plages réalistes pour le matériau étudié.
- Analyser le rapport Vp/Vs en plus des vitesses absolues.
- Utiliser plusieurs capteurs ou plusieurs trajets quand c’est possible.
- Documenter la méthode de mesure, la géométrie et les hypothèses.
10. Comment interpréter le rapport Vp/Vs
Le rapport Vp/Vs est un excellent indicateur synthétique. Un rapport faible à modéré, de l’ordre de 1,6 à 1,9, est fréquent pour des roches compactes. Un rapport plus élevé peut être observé dans des matériaux plus poreux, altérés, fissurés ou riches en fluides. En sismologie du réservoir, ce rapport joue un rôle important pour distinguer les changements de saturation. En géotechnique, il peut aider à repérer des couches meubles ou mal consolidées.
Cela dit, ce ratio ne doit pas être surinterprété. Deux milieux très différents peuvent présenter des rapports proches si plusieurs paramètres se compensent. Il faut donc le combiner avec la densité, la géologie, les logs de forage, les essais de laboratoire et d’autres mesures géophysiques.
11. Sources institutionnelles et académiques à consulter
Pour approfondir, il est recommandé de consulter des sources reconnues. Voici quelques références extérieures utiles :
- USGS.gov pour des explications sur les ondes sismiques, les séismes et les méthodes d’analyse.
- earthquake.usgs.gov pour les ressources pédagogiques et techniques sur la propagation des ondes sismiques.
- serc.carleton.edu pour des contenus universitaires sur la sismologie, les matériaux terrestres et l’interprétation des vitesses.
12. Conclusion
Le calcul de vitesse de propagation des ondes P et S constitue une base indispensable pour comprendre le comportement mécanique du sous-sol et interpréter correctement les signaux sismiques. Qu’il soit réalisé à partir de temps d’arrivée ou de propriétés élastiques, il fournit des informations robustes sur la nature du milieu. L’onde P révèle la réponse du matériau à la compression, tandis que l’onde S met en évidence sa rigidité au cisaillement. Ensemble, ces deux vitesses, ainsi que leur rapport, offrent une lecture très riche des caractéristiques géologiques et géotechniques.
Pour une utilisation professionnelle, il faut toujours associer les résultats numériques à une validation terrain, à un contrôle des unités, à des comparaisons avec des valeurs de référence et à une interprétation géologique rigoureuse. L’outil ci-dessus permet un calcul rapide et visuel, mais la qualité des conclusions dépend avant tout de la qualité des données d’entrée et du contexte scientifique dans lequel elles sont analysées.