نظریه پایه و کاربرد توموگرافی لرزهای
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 تکتونیک صفحهای, aktivitas vulkanik, dan gempa bumi sangat dipengaruhi oleh variasi sifat fisik batuan di kedalaman.
1. مفاهیم اساسی امواج لرزهای
امواج لرزهای 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)
برای مثال، امواج ریلی و لاو، که عموماً در ثبت زلزلههای دوردست غالب هستند، پراکندگی دارند و برای ساختار لیتوسفر تا آستنوسفر کمعمق بسیار آموزنده هستند.
سرعت موج لرزهای به پارامترهای الاستیک (مدول) و چگالی بستگی دارد. تغییرات در دما، ترکیب مواد معدنی، فشار، تخلخل، شکستگیها و وجود سیالات یا ذوبهای جزئی بر سرعت تأثیر میگذارند. این مبنای فیزیکی برای نقشهبرداری از تغییرات سرعت است تا سرنخهایی در مورد شرایط زمینشناسی زیرسطحی ارائه دهد.
۲. اصول توموگرافی: مسئلهی رو به جلو و مسئلهی رو به عقب
Tomografi seismik berdiri di atas dua konsep perhitungan:
الف) مسئله رو به جلو (مسئله رو به جلو)
Diberikan model Bumi (misalnya distribusi kecepatan gelombang), kita dapat menghitung prediksi دادههای لرزهای: 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.
از نظر ریاضی، زمان سفر \(T\) را میتوان به صورت زیر نوشت:
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
که در آن \(v(\mathbf{x})\) سرعت موج در موقعیت \(\mathbf{x}\) و \(ds\) عنصر طول مسیر است.
ب) مسئله معکوس
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.
در عمل، وارونگی اغلب به صورت خطی انجام میشود: مدل اولیه کمی آشفته میشود و سپس اختلاف زمان سفر به صورت زیر محاسبه میشود:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
که در آن \(\delta s = \delta(1/v)\) اختلال کندی است. این معادلات در یک دستگاه خطی بزرگ \( \mathbf{d} = \mathbf{Gm} \) مرتب میشوند، سپس با روش حداقل مربعات میرا یا تغییرات آن حل میشوند.
۳. انواع توموگرافی لرزهای
۱) توموگرافی زمان سفر
رایجترین روش از دادههای فاز P و S دریافتی از بسیاری از زلزلهها استفاده میکند. این دادهها، باقیماندههای زمان سفر در برابر یک مدل مرجع (مثلاً IASP91 یا ak135 در مقیاس جهانی) هستند. برای مدلسازی سرعت سهبعدی در پوسته و گوشته بالایی، به ویژه اگر شبکه ایستگاهها متراکم باشد، مناسب است.
۲) توموگرافی موج سطحی
با استفاده از پراکندگی (تغییر در سرعت فاز/گروه با دوره تناوب) امواج ریلی/لاو. این توموگرافی که به ساختارهای کمعمق تا متوسط حساس است، برای نقشهبرداری از ضخامت لیتوسفر، مناطق کمسرعت (LVZ) و تغییرات دما بسیار مفید است.
۳) توموگرافی تلهسیسمیک (توموگرافی تلهسیسمیک)
استفاده از زلزلههای دوردست (تلهلرزه) که امواج آنها از گوشته عبور کرده و سپس توسط یک شبکه محلی ثبت میشوند. مزیت: آنها منابع متعددی دارند و از جهات مختلف میآیند، بنابراین به "روشن کردن" حجم زیرسطحی منطقه مورد مطالعه، به عنوان مثال زیر آتشفشانها یا مناطق فرورانش، کمک میکنند.
4) Tomografi waveform / full waveform inversion (FWI)
استفاده از شکل موج کامل، نه فقط زمان ورود، از نظر تئوری وضوح بالایی را به همراه دارد زیرا به طور کامل از اطلاعات دامنه و فاز استفاده میکند، اما به الزامات محاسباتی بزرگ و یک مدل اولیه خوب برای جلوگیری از تلههای حداقل محلی نیاز دارد.
۴. مراحل کلی پردازش و وارونگی
1. Akuisisi data
مجموعهای از دادههای ثبتشده از لرزهنگارهای دائمی یا موقت. در مقیاس اکتشافی، منبع ممکن است انفجار یا لرزش لرزهای باشد؛ در مقیاس منطقهای-جهانی، منبع اصلی زلزله است.
2. Identifikasi fase dan picking
تعیین زمان رسیدن امواج P، S یا امواج سطحی. کیفیت برداشت، کیفیت مدل را تعیین میکند.
3. Koreksi dan pemodelan awal
اصلاح زمان (انحراف ساعت)، اصلاح ارتفاع و انتخاب مدل مرجع. مدل اولیه میتواند یک بعدی یا سه بعدی ساده باشد.
4. Ray tracing / simulasi gelombang
مسیر پرتوهای مصنوعی یا شکل موجها را برای ساخت ماتریس حساسیت محاسبه کنید.
5. Inversi dan regularisasi
حل سیستم برای بدست آوردن مدل سرعت. منظمسازی بر اساس هدف انتخاب میشود: اینکه آیا ناهنجاریهای شدید برجسته شوند یا روندهای ملایم.
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.
۵. تفسیر مدل سرعت
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.
با این حال، تفسیر باید محتاطانه باشد زیرا سرعت تحت تأثیر عوامل زیادی قرار میگیرد. در حالت ایدهآل، توموگرافی لرزهای با دادههای دیگری مانند گرانش، مگنتوتلوریک (MT)، ژئودزی، زمینشناسی سطحی و سنگشناسی ترکیب میشود.
۶. کاربردهای توموگرافی لرزهای
الف) مطالعه مناطق فرورانش و دینامیک صفحات
توموگرافی جهانی و منطقهای میتواند نقشه صفحه فرورانش گوشته، از جمله هندسه، عمق و تقسیمبندی آن را ترسیم کند. این اطلاعات برای درک منابع زلزلههای بزرگ، مکانیسمهای اتصال صفحات و تکامل تکتونیکی یک منطقه بسیار مهم است.
ب) سیستمهای آتشفشانی و کاهش بلایای طبیعی
در مناطق آتشفشانی، توموگرافی میتواند مناطق کمسرعت مرتبط با محفظههای ماگما، مسیرهای سیال در حال افزایش یا سنگهای تغییر یافته را شناسایی کند. با پایش منظم، تغییرات سرعت (توموگرافی با فاصله زمانی) میتواند شاخصهایی از حرکت سیال/ماگما را که برای هشدار زودهنگام فوران مرتبط است، ارائه دهد.
ج) اکتشاف انرژی و منابع
در مقیاس پوسته کمعمق، توموگرافی در اکتشافات هیدروکربن، زمینگرمایی و معدن استفاده میشود. تغییرات سرعت به نقشهبرداری از سنگشناسی، ساختارهای گسلی، مناطق شکستگی و مخازن کمک میکند. در زمینگرمایی، ادغام توموگرافی با MT اغلب مؤثر است: توموگرافی اطلاعات الاستیک و MT اطلاعات رسانایی سیال را ارائه میدهد.
د) شناسایی گسلهای فعال و خطرات زلزله
توموگرافی محلی میتواند مناطق ضعیف، مناطق آسیبدیده و ناهمگنی اطراف گسلها را آشکار کند. این امر به درک تقسیمبندی گسل، قفلشدگی بالقوه و تغییرات سرعت که بر تقویت موج (اثرات سایت) تأثیر میگذارند، کمک میکند.
ه) ساختار منطقهای پوسته و لیتوسفر
با استفاده از توموگرافی موج سطحی، محققان میتوانند ضخامت لیتوسفر، مرز موهو و ناهنجاریهای آستنوسفر را نقشهبرداری کنند. نتایج، اساس مدلهای ژئودینامیکی، از جمله تشکیل حوضه، کوهزایی و تکامل قارهای را تشکیل میدهند.
۶. محدودیتها و چالشها
توموگرافی لرزهای به شدت به توزیع منابع و ایستگاهها وابسته است. مناطقی که زلزلههای کمی دارند یا ایستگاههای پراکندهای دارند، روشنایی ضعیفی خواهند داشت که منجر به وضوح پایین میشود. علاوه بر این، فرضیات نظریه پرتو میتواند برای ناهمگنی شدید یا فرکانسهای خاص، دقت کمتری داشته باشد. نویز، خطاهای انتخاب و عدم قطعیت مکان مرکز زلزله نیز میتوانند بر نتایج وارونگی تأثیر بگذارند. بنابراین، ارزیابی وضوح و ادغام چند روش برای جلوگیری از تفسیر بیش از حد بسیار مهم است.
بستن
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.
اگر مایل باشید، میتوانم این مقاله را به نسخهای آکادمیکتر (همراه با ارجاعات و کتابشناسی) تبدیل کنم، یا روی یکی از کاربردها (مثلاً توموگرافی برای آتشفشانها یا انرژی زمینگرمایی) تمرکز کنم.
