地震層析成像的基本理論與應用
地震層析成像是一種利用地震波「拍攝」地球內部結構的地球物理方法。它類似醫學中的CT掃描，利用X射線產生人體三維影像，地震層析成像則利用地震儀網路記錄的傳播時間、振幅或波形資料來模擬地下結構。這種方法至關重要，因為地球內部結構無法直接觀測，而板塊構造、火山活動和地震等地質動力學又深受深層岩石物理性質變化的影響。
1. 地震波的基本概念
地震波是彈性波，它能在岩石中傳播。地震波通常分為：
1. 體波
– P 波（縱波/壓縮波）：傳播速度最快，可以穿過固體和液體，對壓縮速度和密度的變化很敏感。
– S 波（次波/剪切波）：速度比 P 波慢，不能在流體中傳播，對岩石的剛度（剪切模量）敏感。
2. 表面波
例如，瑞利波和勒夫波通常在遠震記錄中占主導地位，它們具有頻散特性，對於了解岩石圈直至淺層軟流圈的結構非常有幫助。
地震波速度取決於彈性參數（模量）和密度。溫度、礦物成分、壓力、孔隙度、裂縫以及流體或部分熔融物的存在都會影響速度。這是繪製速度變化圖以提供有關地下地質條件線索的物理基礎。
2. 斷層掃描原理：正向問題與反向問題
地震層析成像技術基於兩個計算概念：
a) 正向問題（正向問題）
給定一個地球模型（例如，波速分佈），我們可以計算地震資料的預測值：波從震源（地震或人工震源）到台站的傳播時間。在簡單的射線理論方法中，假設波的傳播路徑是一條遵循費馬原理的射線：波會選擇傳播時間最短的路徑。
從數學角度來看，旅行時間 \(T\) 可以表示為：
\[
T = \int_{\text{ray}} \frac{ds}{v(\mathbf{x})}
\]
其中 \(v(\mathbf{x})\) 是位置 \(\mathbf{x}\) 處的波速，\(ds\) 是路徑長度元素。
b) 逆問題
反過來也是如此：我們擁有觀測資料（走時、走時殘差、面波頻散或波形），並希望估計最佳擬合速度模型。反問題通常是非唯一且不適定的：多個模型可能同樣擬合良好，且資料本身包含不確定性/雜訊。因此，需要進行正規化處理，例如平滑、阻尼或地質約束。
在實際應用中，反演通常以線性化的方式進行：對初始模型進行輕微擾動，然後計算走時差，如下所示：
\[
\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) 波形層析成像/全波形反演（FWI）
理論上，使用完整的波形（而不僅僅是到達時間）可以獲得高分辨率，因為它充分利用了幅度和相位信息，但這需要大量的計算資源和一個良好的初始模型來避免陷入局部最小值陷阱。
4. 處理和反轉的一般階段
1. 數據採集
收集來自永久或臨時地震儀的記錄。在勘探尺度上，震源可能是爆炸或振動；在區域乃至全球尺度上，主要震源是地震。
2. 相識別與選擇
決定P波、S波或面波的到達時間。拾取精度決定模型精度。
3. 初步校正與建模
時間校正（時鐘漂移）、高程校正和參考模型選擇。初始模型可以是1D模型或簡單的3D模型。
4. 光線追蹤/波動模擬
計算合成光線或波形的路徑，以建立靈敏度矩陣。
5. 反演和正規化
求解該系統以獲得速度模型。正規化的選擇取決於目標：是突出尖銳的異常值還是平滑的趨勢。
6. 解析度評估和可靠性測試
例如，使用棋盤格測試、尖峰測試或協方差/點擴散函數分析來查看資料真正解析了哪些部分。
5. 速度模型的解釋
層析成像的結果通常是一張相對速度異常圖：高速帶通常被解釋為溫度較低、密度較高或剛度較大的岩石（例如俯衝板塊）。低速帶可能指示高溫、蝕變岩石、充滿流體的斷裂帶或部分熔融—通常與活火山或熱軟流圈有關。
然而，由於速度受多種因素影響，因此必須謹慎解釋。理想情況下，地震層析成像應與其他資料結合，例如重力、大地電磁（MT）、大地測量、地表地質和岩石學資料。
6. 地震層析成像應用
a) 俯衝帶與板塊動力學研究
全球和區域層析成像技術可以繪製地函俯衝板塊的分佈圖，包括其幾何形狀、深度和分段情況。這些資訊對於理解大地震的震源、板塊耦合機制、區域構造演化至關重要。
b) 火山系統與災害減緩
在火山地區，層析成像技術可以識別與岩漿房、上升流體通道或蝕變岩相關的低速區域。透過定期監測，速度變化（延時層析成像）可以提供與火山爆發早期預警相關的流體/岩漿運動指標。
c) 能源和資源勘探
在淺層地殼尺度上，層析成像技術應用於油氣、地熱和礦產勘探。速度變化有助於繪製岩性、斷層構造、裂縫帶和儲層圖。在地熱勘探中，層析成像與大地電磁法（MT）的結合通常非常有效：層析成像提供彈性訊息，而大地電磁法提供流體電導率資訊。
d) 活動斷層和地震災害的特徵
局部層析成像可以揭示斷層周圍的薄弱帶、損傷帶和非均質性。這有助於了解斷層分段、潛在的閉鎖以及影響波放大（場地效應）的速度變化。
e) 區域地殼與岩石圈結構
利用面波層析成像技術，研究人員可以繪製岩石圈厚度、莫霍面邊界和軟流圈異常圖。這些結果構成了地球動力學模型的基礎，包括盆地形成、造山運動和大陸演化。
7. 限制和挑戰
地震層析成像高度依賴震源和台站的分佈。地震活動少或台站稀疏的區域照明不足，導致解析度低。此外，射線理論假設在強非均勻性區域或某些特定頻率下可能不夠精確。噪音、拾取誤差和震源位置不確定性也會影響反演結果。因此，分辨率評估和多方法融合對於避免過度解釋至關重要。
關閉
地震層析成像技術是研究地球內部結構（從局部到全球尺度）的關鍵工具。其理論基礎建立在岩石彈性性質與波傳播的關係以及求解需要正則化和分辨率檢驗的反演問題之上。地震層析成像技術的應用範圍十分廣泛，涵蓋了俯衝板塊和火山岩漿系統的測繪、能源勘探以及地震災害評估等許多領域。隨著地震台網、計算技術和反演方法（包括全波形反演）的進步，地震層析成像技術能夠產生更清晰、資訊更豐富的影像，從而更好地服務科學研究和災害防治工作。
如果您願意，我可以將這篇文章改編成更學術的版本（附引文和參考書目），或重點介紹其中一個應用（例如火山斷層掃描或地熱能斷層掃描）。
