Seismisen tomografian perusteoria ja soveltaminen
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 laattatektoniikka, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. Konsep Dasar Seismiset aallot
Seismiset aallot 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)
Esimerkiksi Rayleigh- ja Love-aallot, jotka ovat yleensä hallitsevia kaukaisten maanjäristysten tallenteissa, ovat dispersioltaan erilaisia ja ne ovat erittäin informatiivisia litosfäärin rakenteesta aina matalaan astenosfääriin asti.
Seismisen aallon nopeus riippuu kimmoparametreista (moduulista) ja tiheydestä. Lämpötilan, mineraalikoostumuksen, paineen, huokoisuuden, murtumien sekä nesteiden tai osittaisten sulamien läsnäolon vaihtelut vaikuttavat nopeuteen. Tämä on fysikaalinen perusta nopeusvaihteluiden kartoittamiselle, jotta saadaan vihjeitä maanalaisista geologisista olosuhteista.
2. Tomografian periaatteet: eteenpäin suuntautuva ongelma ja taaksepäin suuntautuva ongelma
Seisminen tomografia berdiri di atas dua konsep perhitungan:
a) Eteenpäin suuntautuva ongelma (eteenpäin suuntautuva ongelma)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi seismiset tiedot: 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.
Matemaattisesti matka-aika \(T\) voidaan kirjoittaa muodossa:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
jossa v(x) on aallon nopeus kohdassa x ja ds on reitin pituuselementti.
b) Käänteinen ongelma
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.
Käytännössä inversio suoritetaan usein linearisoidulla tavalla: alkuperäistä mallia muutetaan hieman ja sitten matka-aikaero lasketaan seuraavasti:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
jossa \(\delta s = \delta(1/v)\) on hitaushäiriö. Nämä yhtälöt järjestetään suureksi lineaariseksi yhtälöryhmäksi \( \mathbf{d} = \mathbf{Gm} \) ja ratkaistaan sitten vaimennettua pienimmän neliösumman menetelmää tai sen variaatioita käyttäen.
3. Seismisen tomografian tyypit
1) Matka-ajan tomografia
Yleisin menetelmä käyttää P- ja S-vaiheen saapumisia useista maanjäristyksistä. Tiedot ovat matka-ajan jäännöksiä referenssimallia (esim. IASP91 tai ak135 globaalilla tasolla) vasten. Sopii 3D-nopeusmallinnukseen kuoressa ja ylävaipassa, erityisesti jos asemaverkosto on tiheä.
2) Pinta-aaltotomografia
Käyttää Rayleigh/Love-aaltojen dispersiota (vaiheen/ryhmän nopeuden muutos periodin funktiona). Tämä tomografia on herkkä matalille ja keskikokoisille rakenteille, ja se on erittäin hyödyllinen litosfäärin paksuuden, matalan nopeuden vyöhykkeiden (LVZ) ja lämpötilavaihteluiden kartoittamiseen.
3) Teleseisminen tomografia (teleseisminen tomografia)
Käyttämällä kaukaisia maanjäristyksiä (teleseismisiä), joiden aallot kulkevat vaipan läpi ja tallentuvat sitten paikalliseen verkkoon. Etu: niillä on useita lähteitä ja ne tulevat eri suunnista, mikä auttaa "valaisemaan" tutkimusalueen maanalaista tilavuutta, esimerkiksi tulivuorten tai subduktiovyöhykkeiden alla.
4) Tomografi waveform / full waveform inversion (FWI)
Koko aaltomuodon käyttäminen, ei pelkästään saapumisaikaa, tuottaa teoriassa korkean resoluution, koska se hyödyntää täysin amplitudi- ja vaiheinformaatiota, mutta vaatii suuria laskentavaatimuksia ja hyvän alkumallin paikallisten minimiansojen välttämiseksi.
4. Käsittelyn ja inversion yleiset vaiheet
1. Akuisisi data
Pysyvien tai väliaikaisten seismometrien tallenteiden kerääminen. Tutkimusmittakaavassa lähde voi olla räjähdys tai vibroseismi; alueellisella tai globaalilla tasolla ensisijainen lähde on maanjäristys.
2. Identifikasi fase dan picking
P-, S- tai pinta-aaltojen saapumisajan määrittäminen. Poiminnan laatu määrää mallin laadun.
3. Koreksi dan pemodelan awal
Aikakorjaus (kellon siirtymä), korkeuskorjaus ja referenssimallin valinta. Alkuperäinen malli voi olla 1D- tai yksinkertainen 3D-malli.
4. Ray tracing / simulasi gelombang
Laske synteettisten säteiden tai aaltomuotojen polku herkkyysmatriisin muodostamiseksi.
5. Inversi dan regularisasi
Järjestelmän ratkaiseminen nopeusmallin saamiseksi. Regularisointi valitaan tavoitteen perusteella: halutaanko korostaa teräviä poikkeamia vai tasoittaa trendejä.
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. Nopeusmallin tulkinta
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.
Tulkinnassa on kuitenkin oltava varovainen, koska nopeuteen vaikuttavat monet tekijät. Ihannetapauksessa seisminen tomografia yhdistetään muihin tietoihin, kuten painovoimaan, magnetotelluuriseen (MT), geodesiaan, pintageologiaan ja petrologiaan.
6. Seismisen tomografian sovellukset
a) Subduktiovyöhykkeiden ja laattadynamiikan tutkimus
Globaali ja alueellinen tomografia voivat kartoittaa vaipan subduktiolaatan, mukaan lukien sen geometrian, syvyyden ja segmentoinnin. Nämä tiedot ovat ratkaisevan tärkeitä suurten maanjäristysten lähteiden, laattojen kytkentämekanismien ja alueen tektonisen kehityksen ymmärtämiseksi.
b) Tulivuorenpurkausjärjestelmät ja katastrofien lieventäminen
Tulivuorialueilla tomografialla voidaan tunnistaa magmakammioihin, nouseviin nestevirtoihin tai muuttuneisiin kiviin liittyviä matalan nopeuden vyöhykkeitä. Säännöllisellä seurannalla nopeusmuutokset (timelapse-tomografia) voivat antaa indikaattoreita nesteen/magman liikkeestä, jotka ovat olennaisia purkauksen varhaisen varoittamisen kannalta.
c) Energian ja luonnonvarojen etsintä
Matalan maankuoren mittakaavassa tomografiaa käytetään hiilivetyjen, geotermisen energian ja kaivosten etsinnässä. Nopeusvaihtelut auttavat kartoittamaan litologiaa, siirrosrakenteita, murtumavyöhykkeitä ja tekoaltaita. Geotermisessa energiassa tomografian ja MT:n integrointi on usein tehokasta: tomografia tarjoaa tietoa elastisuudesta, MT tarjoaa tietoa nesteen johtavuudesta.
d) Aktiivisten siirrosten ja maanjäristysvaarojen karakterisointi
Paikallinen tomografia voi paljastaa heikkoja alueita, vaurioalueita ja heterogeenisyyttä siirrosten ympärillä. Tämä auttaa ymmärtämään siirrosten segmentoitumista, mahdollista lukittumista ja nopeusvaihteluita, jotka vaikuttavat aaltojen vahvistumiseen (paikkavaikutukset).
e) Alueellinen maankuoren ja litosfäärin rakenne
Pinta-aaltotomografian avulla tutkijat voivat kartoittaa litosfäärin paksuutta, Moho-rajaa ja astenosfäärin poikkeavuuksia. Tulokset muodostavat perustan geodynaamisille malleille, mukaan lukien altaiden muodostuminen, orogeneesi ja mannerten kehitys.
7. Rajoitukset ja haasteet
Seisminen tomografia on erittäin riippuvainen lähteiden ja asemien jakautumisesta. Alueilla, joilla on vähän maanjäristyksiä tai harvat asemat, valaistus on huono, mikä johtaa alhaiseen resoluutioon. Lisäksi sädeteorian oletukset voivat olla epätarkempia voimakkaan heterogeenisyyden tai tiettyjen taajuuksien tapauksessa. Kohina, poimintavirheet ja hypokeskuksen sijainnin epävarmuus voivat myös vaikuttaa inversiotuloksiin. Siksi resoluution arviointi ja monimenetelmäinen integrointi ovat ratkaisevan tärkeitä ylitulkinnan välttämiseksi.
Sulkeminen
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.
Voin halutessasi muokata tämän artikkelin akateemisempaan versioon (viittauksineen ja lähdeluetteloineen) tai keskittyä johonkin sovelluksista (esim. tulivuorten tai geotermisen energian tomografia).
