Grundläggande teori och tillämpning av seismisk tomografi
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 plattektonik, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. Konsep Dasar Seismiska vågor
Seismiska vågor 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)
Till exempel har Rayleigh- och Love-vågor, som generellt sett är dominerande i avlägsna jordbävningsregistreringar, dispersion och är mycket informativa för litosfärens struktur ner till den grunda astenosfären.
Seismiska vågors hastighet beror på elasticitetsparametrar (modul) och densitet. Variationer i temperatur, mineralsammansättning, tryck, porositet, sprickor och närvaron av vätskor eller partiella smältor påverkar hastigheten. Detta är den fysikaliska grunden för att kartlägga hastighetsvariationer för att ge ledtrådar om geologiska förhållanden i underjorden.
2. Principer för tomografi: Framåtproblem och bakåtproblem
Seismisk tomografi berdiri di atas dua konsep perhitungan:
a) Framåtriktat problem (framåtriktat problem)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi seismiska data: 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.
Matematiskt kan restiden \(T\) skrivas som:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
där \(v(\mathbf{x})\) är våghastigheten vid position \(\mathbf{x}\) och \(ds\) är väglängdselementet.
b) Inverst problem
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.
I praktiken utförs inversion ofta linjärt: initialmodellen störs något, och sedan beräknas skillnaden i restid som:
\[
ΔT \approx \int_{ray} \Δs(\mathbf{x})\, ds
\]
där Δs = Δ(1/v) är störningen av långsamheten. Dessa ekvationer är arrangerade i ett stort linjärt system (d = Gm) och löses sedan med den dämpade minstakvadratmetoden eller dess variationer.
3. Typer av seismisk tomografi
1) Resetidstomografi
Den vanligaste metoden använder P- och S-fasdata från många jordbävningar. Data är residualer mot en referensmodell (t.ex. IASP91 eller ak135 på global skala). Lämplig för 3D-hastighetsmodellering i jordskorpan och den övre manteln, särskilt om stationsnätverket är tätt.
2) Ytvågstomografi
Med hjälp av dispersionen (förändring i fas-/grupphastighet med period) av Rayleigh/Love-vågor. Känslig för ytliga till mellanliggande strukturer är denna tomografi mycket användbar för att kartlägga litosfärisk tjocklek, låghastighetszoner (LVZ) och temperaturvariationer.
3) Teleseismisk tomografi (teleseismisk tomografi)
Med hjälp av avlägsna jordbävningar (teleseismiska) vars vågor färdas genom manteln och sedan registreras av ett lokalt nätverk. Fördelen: de har flera källor och kommer från olika håll, vilket hjälper till att "belysa" underjordsvolymen i studieområdet, till exempel under vulkaner eller subduktionszoner.
4) Tomografi waveform / full waveform inversion (FWI)
Att använda hela vågformen, inte bara ankomsttiden, ger teoretiskt sett hög upplösning eftersom den fullt ut utnyttjar amplitud- och fasinformation, men kräver stora beräkningskrav och en bra initialmodell för att undvika lokala minimumfällor.
4. Allmänna steg i bearbetning och inversion
1. Akuisisi-data
Insamling av inspelningar från permanenta eller tillfälliga seismometrar. På prospekteringsskala kan källan vara en explosion eller vibroseism; på regional-global skala är den primära källan en jordbävning.
2. Identifikasi fase dan picking
Bestämning av ankomsttiden för P-, S- eller ytvågor. Kvaliteten på plockningen avgör modellens kvalitet.
3. Koreksi dan pemodelan awal
Tidskorrigering (klockdrift), höjdkorrigering och val av referensmodell. Den initiala modellen kan vara 1D eller enkel 3D.
4. Ray tracing / simulasi gelombang
Beräkna vägen för syntetiska strålar eller vågformer för att konstruera känslighetsmatrisen.
5. Inversi dan regularisasi
Löser systemet för att erhålla en hastighetsmodell. Regularisering väljs baserat på målet: huruvida man vill lyfta fram skarpa avvikelser eller jämna trender.
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. Tolkning av hastighetsmodellen
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.
Tolkningen måste dock vara försiktig eftersom hastigheten påverkas av många faktorer. Helst kombineras seismisk tomografi med andra data såsom gravitation, magnetotellurisk data (MT), geodesi, ytgeologi och petrologi.
6. Seismisk tomografi-tillämpningar
a) Studie av subduktionszoner och plattdynamik
Global och regional tomografi kan kartlägga den mantelsubducerande plattan, inklusive dess geometri, djup och segmentering. Denna information är avgörande för att förstå källorna till stora jordbävningar, plattkopplingsmekanismer och den tektoniska utvecklingen i en region.
b) Vulkaniska system och katastrofbegränsning
I vulkaniska områden kan tomografi identifiera låghastighetszoner associerade med magmakammare, stigande vätskebanor eller förändrade bergarter. Med regelbunden övervakning kan hastighetsförändringar (timelapse-tomografi) ge indikatorer på vätske-/magmarörelse som är relevanta för tidig varning vid utbrott.
c) Energi- och resursutforskning
På grund jordskorpa används tomografi vid prospektering av kolväten, geotermisk energi och gruvdrift. Hastighetsvariationer hjälper till att kartlägga litologi, förkastningsstrukturer, sprickzoner och reservoarer. Inom geotermisk energi är integrationen av tomografi med MT ofta effektiv: tomografi ger elastisk information, MT ger information om fluiders konduktivitet.
d) Karakterisering av aktiva förkastningar och jordbävningsrisker
Lokal tomografi kan avslöja svaga zoner, skadezoner och heterogenitet runt förkastningar. Detta hjälper till att förstå förkastningssegmentering, potentiell låsning och hastighetsvariationer som påverkar vågförstärkning (platseffekter).
e) Regional jordskorpa- och litosfärstruktur
Med hjälp av ytvågstomografi kan forskare kartlägga litosfärens tjocklek, Moho-gränsen och astenosfäriska anomalier. Resultaten utgör grunden för geodynamiska modeller, inklusive bassängbildning, orogenes och kontinental evolution.
7. Begränsningar och utmaningar
Seismisk tomografi är starkt beroende av fördelningen av källor och stationer. Områden med få jordbävningar eller glesa stationer kommer att ha dålig belysning, vilket resulterar i låg upplösning. Dessutom kan antaganden om strålteori vara mindre exakta för stark heterogenitet eller vissa frekvenser. Brus, plockningsfel och osäkerhet kring hypocenterpositionen kan också påverka inversionsresultaten. Därför är upplösningsutvärdering och integration med flera metoder avgörande för att undvika överförtolkning.
Stängning
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.
Om ni önskar kan jag anpassa den här artikeln till en mer akademisk version (med hänvisningar och en bibliografi), eller fokusera på en av tillämpningarna (t.ex. tomografi för vulkaner eller för geotermisk energi).
