Théorie de base et application de la tomographie sismique
La tomographie sismique est une méthode géophysique permettant de « photographier » l'intérieur de la Terre grâce aux ondes sismiques. À l'instar d'un scanner en médecine, qui crée une image tridimensionnelle du corps à partir de rayons X, la tomographie sismique modélise les structures du sous-sol à partir des données de temps de propagation, d'amplitude et de forme d'onde enregistrées par un réseau de sismomètres. Cette méthode est essentielle car l'intérieur de la Terre est inaccessible à l'observation directe, alors que des phénomènes dynamiques tels que la tectonique des plaques , l'activité volcanique et les séismes sont fortement influencés par les variations des propriétés physiques des roches en profondeur.
1. Concept de base des ondes sismiques
Les ondes sismiques sont des ondes élastiques qui se propagent dans la roche. On distingue généralement deux types d'ondes sismiques :
1. Ondes de volume
– Ondes P (primaires/de compression) : se propagent le plus rapidement, peuvent traverser les solides et les liquides, sont sensibles aux changements de vitesse de compression et de densité.
– Ondes S (secondaires/de cisaillement) : plus lentes que les ondes P, ne peuvent pas se propager dans les fluides, sensibles à la rigidité (module de cisaillement) de la roche.
2. Ondes de surface
Par exemple, les ondes de Rayleigh et de Love, qui sont généralement dominantes dans les enregistrements de séismes lointains, présentent une dispersion et sont très informatives sur la structure de la lithosphère jusqu'à l'asthénosphère superficielle.
La vitesse des ondes sismiques dépend des paramètres élastiques (module d'élasticité) et de la densité. Les variations de température, de composition minérale, de pression, de porosité, de fractures et la présence de fluides ou de phases de fusion partielle influent sur cette vitesse. C'est sur ce principe physique que repose la cartographie des variations de vitesse, permettant d'obtenir des informations sur les conditions géologiques du sous-sol.
2. Principes de la tomographie : problème direct et problème inverse
La tomographie sismique repose sur deux concepts de calcul :
a) Problème direct (problème direct)
À partir d'un modèle de la Terre (par exemple, la distribution de la vitesse des ondes), nous pouvons effectuer des prédictions. données sismiquesLe temps de propagation d'une onde depuis sa source (séisme ou source artificielle) jusqu'à la station est appelé temps de parcours. Dans une approche simplifiée par la théorie des rayons, le trajet de l'onde est assimilé à un rayon qui suit le principe de Fermat : l'onde emprunte le chemin de temps de parcours minimal.
Mathématiquement, le temps de trajet \(T\) peut s'écrire comme suit :
\[
T = \int_{\text{rayon}} \frac{ds}{v(\mathbf{x})}
\]
où \(v(\mathbf{x})\) est la vitesse de l'onde à la position \(\mathbf{x}\), et \(ds\) est l'élément de longueur de trajet.
b) Problème inverse
L'inverse est également vrai : nous disposons de données d'observation (temps de propagation, résidus de temps de propagation, dispersions des ondes de surface ou formes d'onde) et souhaitons estimer le modèle de vitesse le plus adapté. Le problème inverse est généralement non unique et mal posé : plusieurs modèles peuvent convenir, et les données comportent des incertitudes ou du bruit. Par conséquent, une régularisation est nécessaire, telle que le lissage, l'amortissement ou l'intégration de contraintes géologiques.
En pratique, l'inversion est souvent réalisée de manière linéarisée : le modèle initial est légèrement perturbé, puis la différence de temps de trajet est calculée comme suit :
\[
\delta T \approx \int_{\text{rayon}} \delta s(\mathbf{x})\, ds
\]
où \(\delta s = \delta(1/v)\) représente la perturbation de lenteur. Ces équations sont organisées en un grand système linéaire \( \mathbf{d} = \mathbf{Gm} \), puis résolues par la méthode des moindres carrés amortis ou ses variantes.
3. Types de tomographie sismique
1) Tomographie du temps de trajet
La méthode la plus courante utilise les arrivées des phases P et S de nombreux séismes. Les données correspondent aux résidus de temps de propagation par rapport à un modèle de référence (par exemple, IASP91 ou ak135 à l'échelle mondiale). Cette méthode est particulièrement adaptée à la modélisation 3D de la vitesse dans la croûte et le manteau supérieur, notamment lorsque le réseau de stations est dense.
2) Tomographie par ondes de surface
Utilisant la dispersion (variation de la vitesse de phase/de groupe en fonction de la période) des ondes de Rayleigh/Love, cette tomographie, sensible aux structures superficielles à intermédiaires, est très utile pour cartographier l'épaisseur de la lithosphère, les zones de faible vitesse (ZFV) et les variations de température.
3) Tomographie télésismique (tomographie télésismique)
On utilise les séismes lointains (télésismiques) dont les ondes se propagent à travers le manteau et sont ensuite enregistrées par un réseau local. L'avantage : ils proviennent de sources multiples et de directions variées, permettant ainsi de mieux visualiser le sous-sol de la zone étudiée, par exemple sous les volcans ou les zones de subduction.
4) Tomographie de forme d'onde / inversion de forme d'onde complète (Îles-du-Prince-Édouard)
L'utilisation de la forme d'onde complète, et non seulement du temps d'arrivée, permet théoriquement d'obtenir une haute résolution car elle exploite pleinement les informations d'amplitude et de phase, mais elle nécessite d'importantes ressources de calcul et un bon modèle initial pour éviter les pièges des minima locaux.
4. Étapes générales du traitement et de l'inversion
1. Données Akuisisi
Collecte d'enregistrements provenant de sismomètres permanents ou temporaires. À l'échelle de l'exploration, la source peut être une explosion ou un vibrosismique ; à l'échelle régionale ou mondiale, la source principale est un séisme.
2. Identification et sélection des phases
Détermination du temps d'arrivée des ondes P, S ou de surface. La qualité du modèle dépend de la précision de la détection.
3. Correction et modélisation initiales
Correction temporelle (dérive d'horloge), correction d'altitude et sélection du modèle de référence. Le modèle initial peut être unidimensionnel ou 3D simple.
4. Tracé de rayons / simulation d'ondes
Calculer le trajet des rayons ou des formes d'onde synthétiques pour construire la matrice de sensibilité.
5. Inversion et régularisation
La résolution du système permet d'obtenir un modèle de vitesse. Le choix de la régularisation dépend de l'objectif : mettre en évidence les anomalies marquées ou les tendances progressives.
6. Évaluation de la résolution et tests de fiabilité
Par exemple, avec le test en damier, le test de pic ou l'analyse de la fonction de covariance/d'étalement du point pour voir quelles parties sont réellement résolues par les données.
5. Interprétation du modèle de vitesse
La tomographie produit généralement une carte des anomalies de vitesse relative : les zones de vitesse élevée sont souvent interprétées comme des roches plus froides, plus denses ou plus rigides (par exemple, des plaques en subduction). Les zones de faible vitesse peuvent indiquer des températures élevées, des roches altérées, des zones de fractures remplies de fluides ou une fusion partielle, souvent associées à des volcans actifs ou à une asthénosphère chaude.
Cependant, l'interprétation doit être prudente car la vitesse est influencée par de nombreux facteurs. Idéalement, la tomographie sismique est combinée à d'autres données telles que la gravimétrie, la magnétotellurie (MT), la géodésie, la géologie de surface et la pétrologie.
6. Applications de la tomographie sismique
a) Étude des zones de subduction et de la dynamique des plaques
La tomographie globale et régionale permet de cartographier la plaque mantellique en subduction, notamment sa géométrie, sa profondeur et sa segmentation. Ces informations sont essentielles pour comprendre les sources des grands séismes, les mécanismes de couplage des plaques et l'évolution tectonique d'une région.
b) Systèmes volcaniques et atténuation des catastrophes
Dans les zones volcaniques, la tomographie permet d'identifier les zones de faible vitesse associées aux chambres magmatiques, aux remontées de fluides ou aux roches altérées. Un suivi régulier des variations de vitesse (tomographie en temps réel) peut fournir des indicateurs de mouvements de fluides/magma, essentiels pour l'alerte précoce aux éruptions.
c) Exploration énergétique et des ressources
À l'échelle de la croûte superficielle, la tomographie est utilisée dans l'exploration des hydrocarbures, de la géothermie et des mines. Les variations de vitesse permettent de cartographier la lithologie, les structures de failles, les zones de fractures et les réservoirs. En géothermie, l'intégration de la tomographie et de la magnétotellurique (MT) est souvent efficace : la tomographie fournit des informations sur l'élasticité, tandis que la MT renseigne sur la conductivité des fluides.
d) Caractérisation des failles actives et des risques sismiques
La tomographie locale permet de révéler les zones de faiblesse, les zones endommagées et l'hétérogénéité autour des failles. Ceci contribue à la compréhension de la segmentation des failles, des blocages potentiels et des variations de vitesse qui influencent l'amplification des ondes (effets de site).
e) Structure crustale et lithosphérique régionale
Grâce à la tomographie par ondes de surface, les chercheurs peuvent cartographier l'épaisseur de la lithosphère, la limite de Moho et les anomalies de l'asthénosphère. Ces résultats constituent la base de modèles géodynamiques, notamment pour la formation des bassins, l'orogenèse et l'évolution continentale.
7. Limites et défis
La tomographie sismique est fortement dépendante de la distribution des sources et des stations. Les zones à faible activité sismique ou à faible densité de stations présentent un éclairage insuffisant, ce qui entraîne une faible résolution. De plus, les hypothèses de la théorie des rayons peuvent être moins précises en cas de forte hétérogénéité ou pour certaines fréquences. Le bruit, les erreurs de pointage et l'incertitude quant à la localisation de l'hypocentre peuvent également influencer les résultats de l'inversion. Par conséquent, l'évaluation de la résolution et l'intégration de plusieurs méthodes sont essentielles pour éviter toute surinterprétation.
Clôture
La tomographie sismique est un outil essentiel pour étudier la structure interne de la Terre, de l'échelle locale à l'échelle globale. Son fondement théorique repose sur la relation entre les propriétés élastiques des roches et la propagation des ondes, ainsi que sur la résolution de problèmes inverses nécessitant une régularisation et des tests de résolution. Ses applications sont nombreuses : de la cartographie des plaques subductées et des systèmes magmatiques volcaniques à l'exploration énergétique et à l'évaluation des risques sismiques. Grâce aux progrès réalisés dans les réseaux sismiques, le calcul et les méthodes d'inversion (notamment l'inversion de forme d'onde complète ), la tomographie sismique est de plus en plus capable de produire des images plus nettes et plus informatives, à des fins scientifiques et de prévention des catastrophes.
Si vous le souhaitez, je peux adapter cet article en une version plus académique (avec des citations et une bibliographie), ou me concentrer sur l'une des applications (par exemple la tomographie pour les volcans ou pour l'énergie géothermique).
