בסיס תיאורטי ויישומים של טומוגרפיה סייסמית

תיאוריה בסיסית ויישום של טומוגרפיה סייסמית

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)
לדוגמה, גלי ריילי ולאב, שבדרך כלל דומיננטיים בצילומי רעידות אדמה רחוקות, הם בעלי פיזור, והם אינפורמטיביים מאוד לגבי מבנה הליתוספירה עד לאסטנוספרה הרדודה.

מהירות גלים סייסמיים תלויה בפרמטרים אלסטיים (מודול) ובצפיפות. שינויים בטמפרטורה, בהרכב המינרלים, בלחץ, בנקבוביות, בסדקים ובנוכחות נוזלים או התכות חלקיות ישפיעו על המהירות. זהו הבסיס הפיזי למיפוי שינויים במהירות כדי לספק רמזים לגבי תנאים גיאולוגיים תת-קרקעיים.

2. עקרונות הטומוגרפיה: בעיה קדימה ובעיה אחורה

טומוגרפיה סייסמית 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.

בפועל, היפוך מבוצע לעתים קרובות באופן ליניארי: המודל הראשוני מופרע מעט, ולאחר מכן הפרש זמן הנסיעה מחושב כך:
\[
ΔT \approx \int_{ray} \deltas(x)\, ds
\]
כאשר Δs = Δ(1/v) היא הפרעת האיטיות. משוואות אלו מסודרות למערכת לינארית גדולה (d = Gm), ולאחר מכן נפתרות באמצעות שיטת הריבועים הפחותים המרוסנים או וריאציות שלה.

3. סוגי טומוגרפיה סייסמית

1) טומוגרפיה בזמן נסיעה
השיטה הנפוצה ביותר משתמשת בנתוני פאזה P ו-S מרעידות אדמה רבות. הנתונים הם שאריות זמן נסיעה כנגד מודל ייחוס (למשל, IASP91 או ak135 בקנה מידה עולמי). מתאים למידול מהירות תלת-ממדי בקרום ובמעטפת העליונה, במיוחד אם רשת התחנות צפופה.

2) טומוגרפיה של גלי שטח
שימוש בפיזור (שינוי במהירות פאזה/קבוצה עם מחזור) של גלי ריילי/לאב. טומוגרפיה זו, הרגישה למבנים רדודים עד בינוניים, שימושית מאוד למיפוי עובי ליתוספרי, אזורים בעלי מהירות נמוכה (LVZs) ושינויי טמפרטורה.

3) טומוגרפיה טלסייסמית (טומוגרפיה טלסייסמית)
שימוש ברעידות אדמה רחוקות (טלזיסמיות) שגלן עובר דרך המעטפת ואז מוקלט על ידי רשת מקומית. היתרון: יש להן מקורות מרובים והן מגיעות מכיוונים שונים, ובכך מסייעות "להאיר" את נפח תת-הקרקע של אזור המחקר, למשל מתחת להרי געש או אזורי סדוקציה.

4) Tomografi waveform / full waveform inversion (FWI)
שימוש בצורת הגל המלאה, לא רק בזמן ההגעה, מניב תיאורטית רזולוציה גבוהה משום שהוא מנצל באופן מלא את מידע האמפליטודה והפאזה, אך דורש דרישות חישוביות גדולות ומודל התחלתי טוב כדי להימנע ממלכודות מינימום מקומיות.

4. שלבים כלליים של עיבוד והיפוך

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.

5. פרשנות מודל המהירות

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), גיאודזיה, גיאולוגיה של פני השטח ופטרולוגיה.

6. יישומי טומוגרפיה סייסמית

א) לימוד אזורי סחיטה ודינמיקת לוחות
טומוגרפיה גלובלית ואזורית יכולה למפות את הלוח החותך את המעטפת, כולל הגיאומטריה שלו, עומקו וחלוקתו לסגמנטציה. מידע זה חיוני להבנת מקורות רעידות אדמה גדולות, מנגנוני צימוד לוחות והאבולוציה הטקטונית של אזור.

ב) מערכות געשיות והפחתת אסונות
באזורים געשיים, טומוגרפיה יכולה לזהות אזורים בעלי מהירות נמוכה הקשורים לתאי מאגמה, מסלולי נוזלים עולים או סלעים שעברו שינוי. בעזרת ניטור קבוע, שינויי מהירות (טומוגרפיה בהילוך מהיר) יכולים לספק אינדיקטורים לתנועת נוזלים/מגמה הרלוונטיים להתרעה מוקדמת על התפרצות.

ג) חיפושי אנרגיה ומשאבים
בקנה מידה רדוד של קרום, טומוגרפיה משמשת בחיפושי פחמימנים, גיאותרמיים וכרייה. שינויי מהירות מסייעים במיפוי ליתולוגיה, מבני שבר, אזורי שברים ומאגרים. בגיאותרמי, שילוב של טומוגרפיה עם מטאורולוגיה יעיל לעיתים קרובות: טומוגרפיה מספקת מידע אלסטי, מטאורולוגיה מספקת מידע על מוליכות נוזלים.

ד) אפיון שברים פעילים וסכנות רעידות אדמה
טומוגרפיה מקומית יכולה לחשוף אזורים חלשים, אזורי נזק והטרוגניות סביב שברים. זה עוזר להבין את פילוח השברים, פוטנציאל הנעילה ושינויי המהירות המשפיעים על הגברת הגל (השפעות אתר).

ה) מבנה קרום וליתוספרי אזורי
באמצעות טומוגרפיה של גלי שטח, חוקרים יכולים למפות את עובי הליטוספירה, את גבול המוהו ואת אנומליות האסתנוספירה. התוצאות מהוות את הבסיס למודלים גיאודינמיים, כולל היווצרות אגן, אורוגנזה ואבולוציה יבשתית.

7. מגבלות ואתגרים

טומוגרפיה סייסמית תלויה במידה רבה בפיזור המקורות והתחנות. אזורים עם מעט רעידות אדמה או תחנות דלילות יהיו בעלי תאורה ירודה, וכתוצאה מכך רזולוציה נמוכה. יתר על כן, הנחות תורת הקרניים יכולות להיות פחות מדויקות עבור הטרוגניות חזקה או תדרים מסוימים. רעש, שגיאות בחירה ואי ודאות במיקום היפומרכז יכולים גם הם להשפיע על תוצאות היפוך. לכן, הערכת רזולוציה ואינטגרציה רב-שיטות הן קריטיות כדי למנוע פרשנות יתר.

סְגִירָה

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.

אם תרצו, אוכל לעבד מאמר זה לגרסה אקדמית יותר (עם ציטוטים וביבליוגרפיה), או להתמקד באחד היישומים (למשל טומוגרפיה להרי געש או לאנרגיה גיאותרמית).