Grundlæggende teori og anvendelse af 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 pladetektonik, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. Konsep Dasar Seismiske bølger
Seismiske bølger 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)
For eksempel har Rayleigh- og Love-bølger, som generelt er dominerende i fjerne jordskælvsoptagelser, spredning og er meget informative for litosfærens struktur ned til den lave asthenosfære.
Seismiske bølgers hastighed afhænger af elasticitetsparametre (modul) og densitet. Variationer i temperatur, mineralsammensætning, tryk, porøsitet, sprækker og tilstedeværelsen af væsker eller delvise smelter vil påvirke hastigheden. Dette er det fysiske grundlag for kortlægning af hastighedsvariationer for at give spor om geologiske forhold i undergrunden.
2. Tomografiens principper: Fremadrettet problem og bagudrettet problem
Seismisk tomografi berdiri di atas dua konsep perhitungan:
a) Fremadrettet problem (fremadrettet problem)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi seismiske 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.
Matematisk kan rejsetiden \(T\) skrives som:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
hvor \(v(\mathbf{x})\) er bølgehastigheden ved position \(\mathbf{x}\), og \(ds\) er vejlængdeelementet.
b) Omvendt 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 praksis udføres inversion ofte lineært: den oprindelige model perturberes en smule, og derefter beregnes rejsetidsforskellen som:
\[
ΔT \approx \int_{ray} \Δs(\mathbf{x})\, ds
\]
hvor Δs = Δ(1/v) er langsomhedsforstyrrelsen. Disse ligninger er arrangeret i et stort lineært system (d = Gm) og derefter løst ved hjælp af den dæmpede mindste kvadraters metode eller dens variationer.
3. Typer af seismisk tomografi
1) Rejsetidstomografi
Den mest almindelige metode bruger P- og S-fasedata fra mange jordskælv. Dataene er rejsetidsresidualer i forhold til en referencemodel (f.eks. IASP91 eller ak135 på global skala). Velegnet til 3D-hastighedsmodellering i jordskorpen og den øvre kappe, især hvis stationsnetværket er tæt.
2) Overfladebølgetomografi
Brug af dispersionen (ændring i fase/gruppehastighed med periode) af Rayleigh/Love-bølger. Denne tomografi er følsom over for overfladiske til mellemliggende strukturer og er meget nyttig til kortlægning af litosfærisk tykkelse, lavhastighedszoner (LVZ'er) og temperaturvariationer.
3) Teleseismisk tomografi (teleseismisk tomografi)
Ved hjælp af fjerne jordskælv (teleseismiske), hvis bølger bevæger sig gennem kappen og derefter registreres af et lokalt netværk. Fordelen: de har flere kilder og kommer fra forskellige retninger, hvilket hjælper med at "oplyse" undergrundsvolumenet i undersøgelsesområdet, for eksempel under vulkaner eller subduktionszoner.
4) Tomografi waveform / full waveform inversion (FWI)
Brug af den komplette bølgeform, ikke kun ankomsttidspunktet, giver teoretisk set høj opløsning, fordi den fuldt ud udnytter amplitude- og faseinformation, men kræver store beregningskrav og en god initialmodel for at undgå lokale minimumsfælder.
4. Generelle stadier af bearbejdning og inversion
1. Akuisisi-data
Indsamling af optagelser fra permanente eller midlertidige seismometre. På en udforskningsskala kan kilden være en eksplosion eller vibroseisme; på en regional-global skala er den primære kilde et jordskælv.
2. Identifikasi fase dan picking
Bestemmelse af ankomsttidspunktet for P-, S- eller overfladebølger. Kvaliteten af udvælgelsen bestemmer modellens kvalitet.
3. Koreksi dan pemodelan awal
Tidskorrektion (urdrift), elevationskorrektion og valg af referencemodel. Den oprindelige model kan være 1D eller simpel 3D.
4. Ray tracing / simulasi gelombang
Beregn banen for syntetiske stråler eller bølgeformer for at konstruere følsomhedsmatricen.
5. Inversi dan regularisasi
Løsning af systemet for at opnå en hastighedsmodel. Regularisering vælges baseret på målet: om der skal fremhæves skarpe anomalier eller udjævnede tendenser.
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. Fortolkning af hastighedsmodellen
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.
Fortolkning skal dog være forsigtig, da hastigheden påvirkes af mange faktorer. Ideelt set kombineres seismisk tomografi med andre data såsom tyngdekraft, magnetotellurisk (MT), geodæsi, overfladegeologi og petrologi.
6. Seismisk tomografi-applikationer
a) Undersøgelse af subduktionszoner og pladedynamik
Global og regional tomografi kan kortlægge den kappe-subducerende plade, herunder dens geometri, dybde og segmentering. Denne information er afgørende for at forstå kilderne til store jordskælv, pladekoblingsmekanismer og den tektoniske udvikling i en region.
b) Vulkaniske systemer og katastrofeafbødning
I vulkanske områder kan tomografi identificere lavhastighedszoner forbundet med magmakamre, stigende væskebaner eller ændrede bjergarter. Med regelmæssig overvågning kan hastighedsændringer (time-lapse-tomografi) give indikatorer for væske-/magmabevægelse, der er relevante for tidlig varsling af udbrud.
c) Energi- og ressourceefterforskning
På den overfladiske jordskorpeskala anvendes tomografi i efterforskning af kulbrinter, geotermisk energi og minedrift. Hastighedsvariationer hjælper med at kortlægge litologi, forkastningsstrukturer, sprækkezoner og reservoirer. I geotermisk energi er integrationen af tomografi med MT ofte effektiv: tomografi giver elastisk information, MT giver information om væskeledningsevne.
d) Karakterisering af aktive forkastninger og jordskælvsfarer
Lokal tomografi kan afsløre svage zoner, skadeszoner og heterogenitet omkring forkastninger. Dette hjælper med at forstå forkastningssegmentering, potentiel låsning og hastighedsvariationer, der påvirker bølgeforstærkning (stedseffekter).
e) Regional skorpe- og litosfærisk struktur
Ved hjælp af overfladebølgetomografi kan forskere kortlægge litosfærisk tykkelse, Moho-grænsen og astenosfæriske anomalier. Resultaterne danner grundlag for geodynamiske modeller, herunder bassindannelse, orogenese og kontinental evolution.
7. Begrænsninger og udfordringer
Seismisk tomografi er meget afhængig af fordelingen af kilder og stationer. Områder med få jordskælv eller sparsomme stationer vil have dårlig belysning, hvilket resulterer i lav opløsning. Desuden kan antagelser i stråleteorien være mindre nøjagtige for stærk heterogenitet eller bestemte frekvenser. Støj, plukkefejl og usikkerhed omkring hypocenterplacering kan også påvirke inversionsresultater. Derfor er opløsningsevaluering og integration af flere metoder afgørende for at undgå overfortolkning.
Lukker
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.
Hvis du ønsker det, kan jeg tilpasse denne artikel til en mere akademisk version (med citater og en bibliografi) eller fokusere på en af anvendelserne (f.eks. tomografi til vulkaner eller geotermisk energi).
