Théorie de base et application de la tomographie sismique
Tomografi seismik adalah salah satu metode geofisika yang digunakan untuk “memotret” bagian dalam Bumi dengan memanfaatkan gelombang seismik. Seperti halnya CT-scan pada dunia medis yang menyusun citra tiga dimensi tubuh dari perambatan sinar-X, tomografi seismik menyusun model struktur bawah permukaan dari data waktu tempuh, amplitudo, atau bentuk gelombang (waveform) yang direkam oleh jaringan seismometer. Metode ini sangat penting karena interior Bumi tidak dapat diamati secara langsung, sementara dinamika seperti tectonique des plaques, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. Konsep Dasar Ondes sismiques
ondes sismiques adalah gelombang elastik yang merambat melalui medium batuan. Secara umum gelombang seismik dibagi menjadi:
1. Gelombang badan (body waves)
– Gelombang P (primer/kompresi) : merambat paling cepat, dapat melalui padat maupun cair, sensitif terhadap perubahan kecepatan kompresi dan densitas.
– Gelombang S (sekunder/geser) : lebih lambat dari P, tidak dapat merambat dalam fluida, sensitif terhadap kekakuan (shear modulus) batuan.
2. Gelombang permukaan (surface waves)
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
tomographie sismique berdiri di atas dua konsep perhitungan:
a) Problème direct (problème direct)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi données sismiques: waktu tempuh gelombang dari sumber (gempa atau sumber buatan) menuju stasiun. Dalam pendekatan sederhana ray theory, lintasan gelombang diasumsikan sebagai sinar (ray) yang mengikuti prinsip Fermat: gelombang memilih lintasan dengan waktu tempuh minimum.
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
Kebalikannya: kita memiliki data pengamatan (waktu tempuh, residu waktu tempuh, dispersi gelombang permukaan, atau waveform) dan ingin menaksir model kecepatan yang paling cocok. Masalah balik umumnya tidak tunggal (non-unique) dan ill-posed : banyak model dapat sama-sama cocok, dan data memiliki ketidakpastian/noise. Karena itu diperlukan regularisasi, misalnya pemulusan (smoothing), damping, atau batasan geologi.
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) Tomografi waveform / full waveform inversion (Î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. Identifikasi fase dan picking
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. Koreksi dan pemodelan awal
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. Ray tracing / simulasi gelombang
Calculer le trajet des rayons ou des formes d'onde synthétiques pour construire la matrice de sensibilité.
5. Inversi dan regularisasi
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. Evaluasi resolusi dan uji keandalan
Misalnya dengan checkerboard test , spike test, atau analisis covariance/point spread function untuk melihat bagian mana yang benar-benar teresolusi oleh data.
5. Interprétation du modèle de vitesse
Hasil tomografi biasanya berupa peta anomali kecepatan relatif: zona kecepatan tinggi sering diinterpretasikan sebagai batuan lebih dingin, lebih rapat, atau lebih kaku (misalnya slab subduksi). Zona kecepatan rendah dapat menunjukkan temperatur tinggi, batuan teralterasi, zona retakan dengan fluida, atau lelehan parsial—sering terkait gunung api aktif atau astenosfer panas.
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
Tomografi seismik merupakan alat utama untuk menyelidiki struktur internal Bumi dari skala lokal hingga global. Dasar teorinya bertumpu pada hubungan antara sifat elastik batuan dan perambatan gelombang, serta pada penyelesaian masalah balik yang membutuhkan regularisasi dan uji resolusi. Aplikasinya luas: dari pemetaan slab subduksi, sistem magma gunung api, eksplorasi energi, hingga penilaian bahaya gempa. Dengan kemajuan jaringan seismik, komputasi, dan metode inversi (termasuk full waveform inversion), tomografi seismik semakin mampu menghasilkan citra yang lebih tajam dan informatif untuk kebutuhan sains maupun mitigasi bencana.
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).
