Seisminės tomografijos pagrindinė teorija ir taikymas
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 tektonik lempeng, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. Pagrindinės sąvokos Seisminės bangos
Seisminės bangos 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)
Pavyzdžiui, Rayleigh ir Love bangos, kurios paprastai dominuoja tolimųjų žemės drebėjimų įrašuose, turi dispersiją ir yra labai informatyvios apie litosferos struktūrą iki pat seklios astenosferos.
Seisminių bangų greitis priklauso nuo elastingumo parametrų (modulio) ir tankio. Temperatūros, mineralinės sudėties, slėgio, poringumo, įtrūkimų ir skysčių ar dalinių lydinių buvimo pokyčiai turės įtakos greičiui. Tai yra fizikinis pagrindas greičio pokyčių žemėlapiams sudaryti, siekiant gauti užuominų apie požemines geologines sąlygas.
2. Tomografijos principai: tiesioginė ir galinė sritys
Tomografi seismik berdiri di atas dua konsep perhitungan:
a) Problema į priekį (problema į priekį)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi seisminiai duomenys: 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.
Matematiškai kelionės laikas \(T\) gali būti užrašytas taip:
\[
T = ∫_{\text{ray}} ∫frac{ds}{v(\mathbf{x})}
\]
kur \(v(\mathbf{x})\) yra bangos greitis \(\mathbf{x}\) pozicijoje, o \(ds\) yra kelio ilgio elementas.
b) Atvirkštinė problema
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.
Praktiškai inversija dažnai atliekama tiesiškai: pradinis modelis šiek tiek perturbuojamas, o tada kelionės laiko skirtumas apskaičiuojamas taip:
\[
∫ΔT apytiksliai ∫Δs(x)₀, ds
\]
kur \(\delta s = \delta(1/v)\) yra lėtumo perturbacija. Šios lygtys yra išdėstytos į didelę tiesinę sistemą \( \mathbf{d} = \mathbf{Gm} \) ir tada sprendžiamos slopinamu mažiausių kvadratų metodu arba jo variacijomis.
3. Seisminės tomografijos tipai
1) Kelionės laiko tomografija
Dažniausiai naudojamas daugelio žemės drebėjimų P ir S fazių atvykimas. Duomenys yra sklidimo laiko liekanos, palyginti su etaloniniu modeliu (pvz., IASP91 arba ak135 pasauliniu mastu). Tinka 3D greičio modeliavimui plutoje ir viršutinėje mantijoje, ypač jei stočių tinklas yra tankus.
2) Paviršinių bangų tomografija
Naudojama Rayleigh/Love bangų dispersija (fazės/grupės greičio pokytis priklausomai nuo periodo). Ši tomografija yra jautri seklioms ir tarpinėms struktūroms, todėl labai naudinga kartografuojant litosferos storį, mažo greičio zonas (LVZ) ir temperatūros pokyčius.
3) Teleseisminė tomografija (teleseisminė tomografija)
Naudojant tolimus žemės drebėjimus (teleseisminius), kurių bangos sklinda per mantiją ir vėliau įrašomos vietinio tinklo. Privalumas: jie turi kelis šaltinius ir sklinda iš įvairių krypčių, taip padėdami „apšviesti“ tyrimo srities požeminį tūrį, pavyzdžiui, po ugnikalniais ar subdukcijos zonomis.
4) Tomografi waveform / full waveform inversion (FWI)
Naudojant visą bangos formą, o ne tik atvykimo laiką, teoriškai gaunama didelė skiriamoji geba, nes visiškai išnaudojama amplitudės ir fazės informacija, tačiau tam reikalingi dideli skaičiavimo reikalavimai ir geras pradinis modelis, kad būtų išvengta vietinių minimalių spąstų.
4. Bendrieji apdorojimo ir inversijos etapai
1. Duomenų rinkimas
Įrašų iš nuolatinių arba laikinų seismometrų rinkimas. Žvalgybos mastu šaltinis gali būti sprogimas arba vibroseizmas; regioniniu-pasauliniu mastu pagrindinis šaltinis yra žemės drebėjimas.
2. Identifikasi fase dan picking
P, S arba paviršiaus bangų atvykimo laiko nustatymas. Modelio kokybę lemia signalo rinkimo kokybė.
3. Koreksi dan pemodelan awal
Laiko korekcija (laikrodžio poslinkis), aukščio korekcija ir etaloninio modelio pasirinkimas. Pradinis modelis gali būti 1D arba paprastas 3D.
4. Ray tracing / simulasi gelombang
Apskaičiuokite sintetinių spindulių arba bangų formų kelią, kad sukurtumėte jautrumo matricą.
5. Inversi dan regularisasi
Sistemos sprendimas greičio modeliui gauti. Reguliarizavimas parenkamas atsižvelgiant į tikslą: ar išryškinti ryškias anomalijas, ar išlyginti tendencijas.
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. Greičio modelio interpretavimas
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.
Vis dėlto interpretuojant reikia būti atsargiam, nes greičiui įtakos turi daug veiksnių. Idealiu atveju seisminė tomografija derinama su kitais duomenimis, tokiais kaip gravitacija, magnetotellūra (MT), geodezija, paviršiaus geologija ir petrologija.
6. Seisminės tomografijos taikymas
a) Subdukcijos zonų ir plokščių dinamikos tyrimas
Visuotinė ir regioninė tomografija gali sudaryti mantiją subdukuojančios plokštės žemėlapį, įskaitant jos geometriją, gylį ir segmentaciją. Ši informacija yra labai svarbi norint suprasti didelių žemės drebėjimų šaltinius, plokščių sąveikos mechanizmus ir regiono tektoninę evoliuciją.
b) Vulkaninės sistemos ir nelaimių padarinių švelninimas
Vulkaninėse vietovėse tomografija gali nustatyti mažo greičio zonas, susijusias su magmos kameromis, kylančių skysčių keliais arba pakitusiomis uolienomis. Reguliariai stebint greičio pokyčius (laipsniško vaizdo tomografija), galima gauti skysčių / magmos judėjimo rodiklius, svarbius ankstyvam išsiveržimo įspėjimui.
c) Energijos ir išteklių žvalgyba
Sekliame plutos mastu tomografija naudojama angliavandenilių, geoterminėje ir kasybos žvalgyboje. Greičio pokyčiai padeda sudaryti litologijos, lūžių struktūrų, lūžių zonų ir rezervuarų žemėlapius. Geoterminėje energijoje tomografijos integravimas su MT dažnai yra efektyvus: tomografija suteikia informaciją apie elastingumą, MT – informaciją apie skysčių laidumą.
d) Aktyviųjų lūžių ir žemės drebėjimų pavojų apibūdinimas
Vietinė tomografija gali atskleisti silpnąsias zonas, pažeidimo zonas ir nevienalytiškumą aplink lūžius. Tai padeda suprasti lūžių segmentaciją, galimą užsiblokavimą ir greičio pokyčius, kurie turi įtakos bangų stiprinimui (vietos poveikiui).
e) Regioninė plutos ir litosferos struktūra
Naudodami paviršinių bangų tomografiją, tyrėjai gali sudaryti litosferos storio, Moho ribos ir astenosferos anomalijų žemėlapius. Gauti rezultatai sudaro geodinaminių modelių, įskaitant baseinų formavimąsi, orogenezę ir žemynų evoliuciją, pagrindą.
7. Apribojimai ir iššūkiai
Seisminė tomografija labai priklauso nuo šaltinių ir stočių pasiskirstymo. Vietovėse, kuriose mažai žemės drebėjimų arba stočių yra nedaug, apšvietimas bus prastas, todėl skiriamoji geba bus maža. Be to, spindulių teorijos prielaidos gali būti mažiau tikslios esant dideliam heterogeniškumui arba tam tikriems dažniams. Triukšmas, parinkimo klaidos ir hipocentro vietos neapibrėžtumas taip pat gali turėti įtakos inversijos rezultatams. Todėl skiriamosios gebos įvertinimas ir kelių metodų integravimas yra labai svarbūs siekiant išvengti perteklinės interpretacijos.
Uždarymas
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.
Jei pageidaujate, galiu pritaikyti šį straipsnį į akademinę versiją (su citatomis ir bibliografija) arba sutelkti dėmesį į vieną iš taikymo sričių (pvz., ugnikalnių ar geoterminės energijos tomografija).
