Teoria di base e applicazioni della tomografia sismica
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 tettonica delle placche, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. Konsep Dasar Onde sismiche
onde sismiche 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)
Ad esempio, le onde di Rayleigh e di Love, che sono generalmente dominanti nelle registrazioni di terremoti distanti, presentano dispersione e sono molto informative per la struttura della litosfera fino all'astenosfera superficiale.
La velocità delle onde sismiche dipende da parametri elastici (modulo di elasticità) e densità. Variazioni di temperatura, composizione minerale, pressione, porosità, fratture e presenza di fluidi o fusi parziali influenzano la velocità. Questa è la base fisica per la mappatura delle variazioni di velocità, che fornisce indizi sulle condizioni geologiche del sottosuolo.
2. Principi di tomografia: problema in avanti e problema all'indietro
Tomografi seismik berdiri di atas dua konsep perhitungan:
a) Problema in avanti (problema in avanti)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi dati sismici: 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.
Matematicamente il tempo di percorrenza \(T\) può essere scritto come:
\[
T = \int_{\text{raggio}} \frac{ds}{v(\mathbf{x})}
\]
dove \(v(\mathbf{x})\) è la velocità dell'onda nella posizione \(\mathbf{x}\), e \(ds\) è l'elemento di lunghezza del percorso.
b) Problema inverso
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.
In pratica, l'inversione viene spesso eseguita in modo linearizzato: il modello iniziale viene leggermente perturbato e quindi la differenza di tempo di percorrenza viene calcolata come:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
dove \(\delta s = \delta(1/v)\) è la perturbazione di lentezza. Queste equazioni vengono organizzate in un grande sistema lineare \( \mathbf{d} = \mathbf{Gm} \), quindi risolto con il metodo dei minimi quadrati smorzati o sue varianti.
3. Tipi di tomografia sismica
1) Tomografia del tempo di percorrenza
Il metodo più comune utilizza gli arrivi delle fasi P e S di numerosi terremoti. I dati sono residui dei tempi di propagazione rispetto a un modello di riferimento (ad esempio, IASP91 o ak135 su scala globale). Adatto per la modellazione della velocità 3D nella crosta e nel mantello superiore, soprattutto se la rete di stazioni è densa.
2) Tomografia a onde superficiali
Sfruttando la dispersione (variazione della velocità di fase/gruppo in funzione del periodo) delle onde di Rayleigh/Love, questa tomografia, sensibile alle strutture superficiali e intermedie, è molto utile per mappare lo spessore della litosfera, le zone a bassa velocità (LVZ) e le variazioni di temperatura.
3) Tomografia telesismica (tomografia telesismica)
Si utilizzano terremoti distanti (telesismici) le cui onde si propagano attraverso il mantello e vengono poi registrate da una rete locale. Il vantaggio: hanno sorgenti multiple e provengono da diverse direzioni, contribuendo così a "illuminare" il volume sotterraneo dell'area di studio, ad esempio al di sotto di vulcani o zone di subduzione.
4) Tomografi waveform / full waveform inversion (FWI)
L'utilizzo dell'intera forma d'onda, e non solo del tempo di arrivo, consente teoricamente di ottenere un'elevata risoluzione perché sfrutta appieno le informazioni di ampiezza e fase, ma richiede notevoli risorse computazionali e un buon modello iniziale per evitare trappole dovute a minimi locali.
4. Fasi generali di elaborazione e inversione
1. Dati di Akuisisi
Raccolta di registrazioni provenienti da sismometri permanenti o temporanei. Su scala esplorativa, la sorgente può essere un'esplosione o un fenomeno vibrosismico; su scala regionale-globale, la sorgente primaria è un terremoto.
2. Identifikasi fase dan picking
Determinazione del tempo di arrivo delle onde P, S o di superficie. La qualità della selezione determina la qualità del modello.
3. Koreksi dan pemodelan awal
Correzione temporale (deriva dell'orologio), correzione dell'elevazione e selezione del modello di riferimento. Il modello iniziale può essere monodimensionale o tridimensionale semplice.
4. Ray tracing / simulasi gelombang
Calcola il percorso dei raggi o delle forme d'onda sintetiche per costruire la matrice di sensibilità.
5. Inversi dan regularisasi
Risolvere il sistema per ottenere un modello di velocità. La regolarizzazione viene scelta in base all'obiettivo: evidenziare anomalie marcate o smussare le tendenze.
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. Interpretazione del modello di velocità
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.
Tuttavia, l'interpretazione deve essere cauta perché la velocità è influenzata da molti fattori. Idealmente, la tomografia sismica viene combinata con altri dati come quelli gravimetrici, magnetotellurici (MT), geodetici, geologici di superficie e petrologici.
6. Applicazioni della tomografia sismica
a) Studio delle zone di subduzione e della dinamica delle placche
La tomografia globale e regionale permette di mappare la placca del mantello in subduzione, inclusi la sua geometria, la profondità e la segmentazione. Queste informazioni sono fondamentali per comprendere le sorgenti dei grandi terremoti, i meccanismi di accoppiamento delle placche e l'evoluzione tettonica di una regione.
b) Sistemi vulcanici e mitigazione dei disastri
Nelle aree vulcaniche, la tomografia può identificare zone a bassa velocità associate a camere magmatiche, vie di risalita dei fluidi o rocce alterate. Con un monitoraggio regolare, le variazioni di velocità (tomografia time-lapse) possono fornire indicatori del movimento di fluidi/magmi rilevanti per l'allerta precoce delle eruzioni.
c) Esplorazione di energia e risorse
A livello della crosta terrestre superficiale, la tomografia viene utilizzata nell'esplorazione di idrocarburi, geotermia e miniere. Le variazioni di velocità aiutano a mappare la litologia, le strutture di faglia, le zone di frattura e i giacimenti. In ambito geotermico, l'integrazione della tomografia con la tomografia magnetotellurica (MT) è spesso efficace: la tomografia fornisce informazioni sull'elasticità, mentre la MT fornisce informazioni sulla conduttività dei fluidi.
d) Caratterizzazione delle faglie attive e dei rischi sismici
La tomografia locale può rivelare zone di debolezza, zone danneggiate ed eterogeneità intorno alle faglie. Questo aiuta a comprendere la segmentazione delle faglie, il potenziale blocco e le variazioni di velocità che influenzano l'amplificazione delle onde (effetti di sito).
e) Struttura regionale della crosta e della litosfera
Utilizzando la tomografia a onde superficiali, i ricercatori possono mappare lo spessore della litosfera, il confine di Moho e le anomalie astenosferiche. I risultati costituiscono la base per i modelli geodinamici, tra cui la formazione dei bacini, l'orogenesi e l'evoluzione continentale.
7. Limitazioni e sfide
La tomografia sismica dipende fortemente dalla distribuzione delle sorgenti e delle stazioni. Le aree con pochi terremoti o stazioni sparse avranno una scarsa illuminazione, con conseguente bassa risoluzione. Inoltre, le ipotesi della teoria dei raggi possono essere meno accurate in presenza di forte eterogeneità o per determinate frequenze. Anche il rumore, gli errori di rilevamento e l'incertezza nella localizzazione dell'ipocentro possono influenzare i risultati dell'inversione. Pertanto, la valutazione della risoluzione e l'integrazione multimetodo sono cruciali per evitare interpretazioni errate.
Chiusura
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.
Se lo desideri, posso adattare questo articolo in una versione più accademica (con citazioni e bibliografia) oppure concentrarmi su una delle applicazioni (ad esempio, la tomografia per i vulcani o per l'energia geotermica).
