Основна теория и приложение на сеизмичната томография
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. Основни понятия Сеизмични вълни
Gelombang seismik 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.
На практика инверсията често се извършва по линеаризиран начин: началният модел се смущава леко и след това разликата във времето за пътуване се изчислява като:
\[
\delta T \approx \int_{\text{ray}} \delta s(\mathbf{x})\, ds
\]
където (Δs = Δ(1/v)) е смущението от бавността. Тези уравнения са подредени в голяма линейна система (d = Gm), след което са решени чрез метода на затихнатите най-малки квадрати или негови варианти.
3. Видове сеизмична томография
1) Томография с времетраене
Най-разпространеният метод използва P- и S-фазови пристигания от много земетресения. Данните са остатъци от времето на пътуване спрямо референтен модел (напр. IASP91 или ak135 в глобален мащаб). Подходящ е за 3D моделиране на скоростта в земната кора и горната мантия, особено ако мрежата от станции е гъста.
2) Повърхностно-вълнова томография
Използвайки дисперсията (промяна във фазовата/групова скорост с периода) на вълните на Рейли/Лав. Чувствителна към плитки до междинни структури, тази томография е много полезна за картографиране на литосферната дебелина, зони с ниска скорост (LVZ) и температурни вариации.
3) Телесеизмична томография (телесеизмична томография)
Използване на отдалечени земетресения (телесеизмични), чиито вълни преминават през мантията и след това се регистрират от локална мрежа. Предимството: те имат множество източници и идват от различни посоки, като по този начин помагат за „осветяване“ на подземния обем на изследваната област, например под вулкани или зони на субдукция.
4) Tomografi waveform / full waveform inversion (FWI)
Използването на цялата форма на вълната, не само на времето на пристигане, теоретично води до висока резолюция, защото използва напълно информацията за амплитудата и фазата, но изисква големи изчислителни изисквания и добър начален модел, за да се избегнат локални минимални капани.
4. Общи етапи на обработка и инверсия
1. Събиране на данни
Събиране на записи от постоянни или временни сеизмометри. В проучвателен мащаб източникът може да бъде експлозия или вибросеизм; в регионален-глобален мащаб основният източник е земетресение.
2. Identifikasi fase dan picking
Определяне на времето на пристигане на P, S или повърхностни вълни. Качеството на пикиране определя качеството на модела.
3. Koreksi dan pemodelan awal
Корекция на времето (отклонение на часовника), корекция на котата и избор на референтен модел. Първоначалният модел може да бъде 1D или прост 3D.
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.
Ако желаете, мога да адаптирам тази статия в по-академична версия (с цитати и библиография) или да се съсредоточа върху едно от приложенията (напр. томография за вулкани или за геотермална енергия).
