地震層析成像的基本理論與應用
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. 斷層掃描原理：正向問題與反向問題
Tomografi seismik berdiri di atas dua konsep perhitungan:
a) 正向問題（正向問題）
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\) 是路徑長度元素。
b) 逆問題
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} \)，然後用阻尼最小二乘法或其變體求解。
3. 地震層析成像的類型
1）旅行時間斷層掃描
最常用的方法是利用多次地震的P波和S波到時資料。這些數據是相對於參考模型（例如，全球範圍內的IASP91或ak135）的走時殘差。此方法適用於地殼和上地函的三維速度建模，尤其適用於台網密集的情況。
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. 地震層析成像應用
a) 俯衝帶與板塊動力學研究
全球和區域層析成像技術可以繪製地函俯衝板塊的分佈圖，包括其幾何形狀、深度和分段情況。這些資訊對於理解大地震的震源、板塊耦合機制、區域構造演化至關重要。
b) 火山系統與災害減緩
在火山地區，層析成像技術可以識別與岩漿房、上升流體通道或蝕變岩相關的低速區域。透過定期監測，速度變化（延時層析成像）可以提供與火山爆發早期預警相關的流體/岩漿運動指標。
c) 能源和資源勘探
在淺層地殼尺度上，層析成像技術應用於油氣、地熱和礦產勘探。速度變化有助於繪製岩性、斷層構造、裂縫帶和儲層圖。在地熱勘探中，層析成像與大地電磁法（MT）的結合通常非常有效：層析成像提供彈性訊息，而大地電磁法提供流體電導率資訊。
d) 活動斷層和地震災害的特徵
局部層析成像可以揭示斷層周圍的薄弱帶、損傷帶和非均質性。這有助於了解斷層分段、潛在的閉鎖以及影響波放大（場地效應）的速度變化。
e) 區域地殼與岩石圈結構
利用面波層析成像技術，研究人員可以繪製岩石圈厚度、莫霍面邊界和軟流圈異常圖。這些結果構成了地球動力學模型的基礎，包括盆地形成、造山運動和大陸演化。
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.
如果您願意，我可以將這篇文章改編成更學術的版本（附引文和參考書目），或重點介紹其中一個應用（例如火山斷層掃描或地熱能斷層掃描）。
