Baz teyorik ak aplikasyon tomografi sismik

Teyori debaz ak aplikasyon tomografi sismik

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 plak, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.

1. Konsèp debaz yo Vag sismik

Vag sismik 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)
Pa egzanp, ond Rayleigh ak Love, ki jeneralman dominan nan anrejistreman tranblemanntè ki lwen, gen dispèsyon, epi yo bay anpil enfòmasyon sou estrikti litosfè a jis nan astenosfè fon an.

Vitès ond sismik la depann de paramèt elastik yo (modil) ak dansite. Varyasyon nan tanperati, konpozisyon mineral, presyon, porosit, frakti, ak prezans likid oswa fonn pasyèl ap afekte vitès la. Sa a se baz fizik pou trase varyasyon vitès yo pou bay endikasyon sou kondisyon jeyolojik anba tè.

2. Prensip Tomografi: Pwoblèm Devan ak Pwoblèm Dèyè

Tomografi sismik berdiri di atas dua konsep perhitungan:

a) Pwoblèm pou pi devan (pwoblèm pou pi devan)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi done sismik: 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.

Matematikman, nou ka ekri tan vwayaj la \(T\) kòm:
\[
T = \int_{\text{reyon}} \frac{ds}{v(\mathbf{x})}
\]
kote \(v(\mathbf{x})\) se vitès vag la nan pozisyon \(\mathbf{x}\), epi \(ds\) se eleman longè chemen an.

b) Pwoblèm envès
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.

An pratik, souvan yo fè envèsyon an nan yon fason linearize: yo deranje modèl inisyal la yon ti kras, epi answit yo kalkile diferans tan vwayaj la kòm:
\[
ΔT \approx \int_{\text{reyon}} Δs(\mathbf{x})\, ds
\]
kote \(\delta s = \delta(1/v)\) se perturbasyon ralentisman an. Ekwasyon sa yo ranje nan yon gwo sistèm lineyè \( \mathbf{d} = \mathbf{Gm} \), answit rezoud pa metòd kare minimòm amortisman an oswa varyasyon li yo.

3. Kalite Tomografi Sismik

1) Tomografi tan vwayaj
Metòd ki pi komen an itilize arive faz P ak S ki soti nan anpil tranblemanntè. Done yo se rezidyèl tan vwayaj kont yon modèl referans (pa egzanp, IASP91 oswa ak135 sou yon echèl mondyal). Li apwopriye pou modèl vitès 3D nan kwout la ak manto siperyè a, sitou si rezo estasyon an dans.

2) Tomografi vag sifas
Lè l sèvi avèk dispèsyon (chanjman nan vitès faz/gwoup avèk peryòd) ond Rayleigh/Love yo. Tomografi sa a sansib a estrikti ki soti nan fon rive nan entèmedyè, li trè itil pou trase epesè litosferik, zòn ki gen vitès ki ba (LVZ), ak varyasyon tanperati.

3) Tomografi teleseismik (tomografi teleseismik)
Yo itilize tranblemanntè ki lwen (teleseismik) ki gen vag ki vwayaje nan manto a epi ki anrejistre pa yon rezo lokal. Avantaj la: yo gen plizyè sous epi yo soti nan divès direksyon, kidonk ede "eklere" volim anba tè zòn etid la, pa egzanp anba vòlkan oswa zòn subdiksyon.

4) Tomografi waveform / full waveform inversion (FWI)
Lè w itilize fòm ond lan nèt, pa sèlman lè rive a, teyorikman ou jwenn yon rezolisyon ki wo paske li itilize enfòmasyon anplitid ak faz yo nèt, men li mande anpil kondisyon kalkil ak yon bon modèl inisyal pou evite pyèj minimòm lokal yo.

4. Etap Jeneral Pwosesis ak Envèsyon

1. Akuisisi data
Koleksyon anrejistreman ki soti nan sismomèt pèmanan oswa tanporè. Sou yon echèl eksplorasyon, sous la kapab yon eksplozyon oswa yon vibroseism; sou yon echèl rejyonal-mondyal, sous prensipal la se yon tranblemanntè.

2. Identifikasi fase dan picking
Detèmine tan rive vag P, S, oubyen vag sifas yo. Kalite seleksyon an detèmine kalite modèl la.

3. Koreksi dan pemodelan awal
Koreksyon tan (derivasyon revèy), koreksyon elevasyon, ak seleksyon modèl referans. Modèl inisyal la ka 1D oswa 3D senp.

4. Ray tracing / simulasi gelombang
Kalkile trajektwa reyon sentetik yo oubyen fòm ond yo pou konstwi matris sansiblite a.

5. Inversi dan regularisasi
Rezoud sistèm nan pou jwenn yon modèl vitès. Yo chwazi regilarizasyon an selon objektif la: si wi ou non pou mete aksan sou anomali byen file oswa tandans lis.

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. Entèpretasyon Modèl Vitès la

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.

Sepandan, entèpretasyon an dwe fèt avèk prekosyon paske vitès la enfliyanse pa anpil faktè. Idealman, tomografi sismik la konbine avèk lòt done tankou gravite, mayetotelurik (MT), jeodezi, jewoloji sifas, ak petwòli.

6. Aplikasyon Tomografi Sismik

a) Etid sou zòn subdiksyon ak dinamik plak yo
Tomografi mondyal ak rejyonal ka trase yon kat jeyografik plak subduktif manto a, ki gen ladan jeyometri li, pwofondè li ak segmentasyon li. Enfòmasyon sa a enpòtan pou konprann sous gwo tranblemanntè yo, mekanis kouplaj plak yo, ak evolisyon tektonik yon rejyon.

b) Sistèm vòlkanik ak rediksyon dezas
Nan zòn vòlkanik yo, tomografi ka idantifye zòn ki gen vitès ki ba ki asosye avèk chanm magma, chemen likid k ap monte, oswa wòch ki modifye. Avèk siveyans regilye, chanjman vitès yo (tomografi time-lapse) ka bay endikatè mouvman likid/magma ki enpòtan pou avètisman bonè sou eripsyon.

c) Eksplorasyon enèji ak resous
Nan echèl kwout tè a ki pa twò fon, yo itilize tomografi nan eksplorasyon idrokarbur, jewotèmal, ak min. Varyasyon vitès ede trase litoloji, estrikti fay, zòn frakti, ak rezèvwa. Nan jewotèmal, entegrasyon tomografi ak MT souvan efikas: tomografi bay enfòmasyon elastik, MT bay enfòmasyon sou konduktivite likid.

d) Karakterizasyon fay aktif ak danje tranblemanntè
Tomografi lokal la ka revele zòn fèb, zòn domaj, ak eterogeneite otou fay yo. Sa ede konprann segmentasyon fay yo, potansyèl bloke, ak varyasyon vitès ki enfliyanse anplifikasyon vag (efè sit).

e) Estrikti rejyonal kwout la ak litosfè a
Lè yo itilize tomografi vag sifas, chèchè yo ka trase yon kat jeyografik sou epesè litosferik la, limit Moho a, ak anomali astenosferik yo. Rezilta yo fòme baz pou modèl jeodinamik yo, tankou fòmasyon basen, orojènèz, ak evolisyon kontinantal.

7. Limit ak Difikilte

Tomografi sismik depann anpil de distribisyon sous ak estasyon yo. Zòn ki gen kèk tranblemanntè oswa estasyon ki ra yo pral gen move ekleraj, sa ki lakòz yon rezolisyon ki ba. Anplis de sa, sipozisyon teyori reyon yo ka mwens egzak pou gwo eterogeneite oswa sèten frekans. Bri, erè nan seleksyon, ak ensètitid nan pozisyon iposant lan kapab enfliyanse rezilta envèsyon yo tou. Se poutèt sa, evalyasyon rezolisyon ak entegrasyon plizyè metòd enpòtan pou evite twòp entèpretasyon.

Penutup

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.

Si ou vle, mwen ka adapte atik sa a nan yon vèsyon pi akademik (avèk sitasyon ak yon bibliyografi), oubyen konsantre sou youn nan aplikasyon yo (pa egzanp, tomografi pou vòlkan oubyen pou enèji jeotèmal).

Kite yon kòmantè