📊 觀測數據分析：震測剖面圖 (Record Section)
本專案透過自動化 Python 腳本 seismicwave_sort.py 讀取 CWASN SAC 格式的垂直分量 (Z) 地震資料，進行帶通濾波 (0.5-5.0 Hz) 與振幅正規化後，繪製出以下之專業地震波形觀測剖面圖：
💡 圖片詳細說明與地震學解讀：
1. 軸向定義：本圖的 X 軸為震央距離 (Distance, km)，Y 軸為走時 (Travel Time, seconds)。圖中的每一條震盪波形線代表不同觀測測站所接收到的垂直分量（Z分量）實測震波訊號。
2. 波形排列與走時趨勢：所有觀測站依照與震央的距離由近到遠（由左至右）依序排列。隨著震央距離的增加，地震波抵達測站所需的時間也隨之增加，在空間中清晰地勾勒出地震學中經典的走時曲線 (Travel-Time Curve) 軌跡。
3. P波與S波辨識：在剖面圖中，最左側各通道中最先抵達（走時較短）、振幅較小且頻率較高的初達波相為 P 波 (Primary Wave)；而隨後抵達、振幅明顯放大、週期較長且能量顯著的波相則為 S 波 (Secondary Wave)。藉由觀測兩者的抵達時間差（S-P 時差），可進一步反演該區域地殼與上部地函的速度構造。
1. 軸向定義：本圖的 X 軸為震央距離 (Distance, km)，Y 軸為走時 (Travel Time, seconds)。圖中的每一條震盪波形線代表不同觀測測站所接收到的垂直分量（Z分量）實測震波訊號。
2. 波形排列與走時趨勢：所有觀測站依照與震央的距離由近到遠（由左至右）依序排列。隨著震央距離的增加，地震波抵達測站所需的時間也隨之增加，在空間中清晰地勾勒出地震學中經典的走時曲線 (Travel-Time Curve) 軌跡。
3. P波與S波辨識：在剖面圖中，最左側各通道中最先抵達（走時較短）、振幅較小且頻率較高的初達波相為 P 波 (Primary Wave)；而隨後抵達、振幅明顯放大、週期較長且能量顯著的波相則為 S 波 (Secondary Wave)。藉由觀測兩者的抵達時間差（S-P 時差），可進一步反演該區域地殼與上部地函的速度構造。
📚 第三章 課程講義與核心主題
以下為本課程第三章的核心單元規劃，重點在於利用地震波探測地殼結構與地球內部的物理性質：
| 章節 | 主題名稱 | 核心內容概述 |
|---|---|---|
| CH 3.1 | 前言 Introduction | 介紹如何利用地震波走時、波形作為探針來研究地球內部結構的基本原理，以及走時逆問題的初探。 |
| CH 3.2 | 折射地震學 Refraction Seismology | 探討地震波在不同岩層界面產生的折射現象與首波 (Head Waves)。包含平坦地層與傾斜地層的走時計算方法，以及如何定出康拉德面與莫霍面。 |
| CH 3.3 | 反射地震學 Reflection Seismology | 詳細介紹地震波在垂直/傾斜界面的反射機制、走時曲線（雙曲線特徵）及其在震測勘探與地殼深部構造研究中的應用。 |
| CH 3.4 | 球體地球中的地震波傳播 | 考量地球球體幾何效應下的射線參數 (Ray Parameter, p)、射線路徑與走時方程。 |
| CH 3.5 | 體波走時與地球內部構造 | 分析體波在全球尺度下的走時曲線與核心震相命名（如 PcP, PKP, SKS），探討一維地球速度模型（PREM）。 |
| CH 3.6 | 上部地函結構 Upper Mantle | 解構低速帶（LVZ）以及由礦物晶格相變引發的 410 km 與 660 km 地函不連續面。 |
| CH 3.7 | 下部地函與地核結構 | 探討核幔邊界（CMB）的 D" 層、超低速帶（ULVZ）、液態外核影帶以及固態內核的各向異性。 |
| CH 3.8 | 地函與地核的組成 Composition | 結合高溫高壓礦物物理學與地震波速，反演地球內部的化學組成、密度分布與三維層析成像。 |
🌍 CH 3.1 前言 Introduction
利用地表測站接收到的震波訊號，解構並反演（Inversion）震波前進時穿透的深部介質構造。
📊 圖 3-1-1 說明【地震學探針與基本射線幾何】
展示了實體震波在非均勻地球介質中傳播時的射線軌跡。隨著深度增加、波速改變，波動前進方向發生連續彎曲，建立起地表觀測站走時與地下深處速度構造的非線性光學映射關係。
展示了實體震波在非均勻地球介質中傳播時的射線軌跡。隨著深度增加、波速改變，波動前進方向發生連續彎曲，建立起地表觀測站走時與地下深處速度構造的非線性光學映射關係。
📐 CH 3.2 折射地震學 Refraction Seismology
當入射角達到臨界角時，地震波在高速界面滑行並激發出「首波 (Head Wave)」，這是精密估算莫霍面深度與地殼厚度的經典物理手段。
📊 圖 3-2-1 說明【臨界折射與首波產生幾何】
依據惠更斯原理推導臨界折射。當入射角 $i$ 滿足 $\sin i_c = v_1 / v_2$ 時，震波能量沿速度不連續面邊界滑行，並向上一層不斷激發出具備線性平面波前的「首波」。
依據惠更斯原理推導臨界折射。當入射角 $i$ 滿足 $\sin i_c = v_1 / v_2$ 時，震波能量沿速度不連續面邊界滑行，並向上一層不斷激發出具備線性平面波前的「首波」。
📊 圖 3-2-2 說明【直接波與折射波線性走時曲線 (T-X 圖)】
在距離-走時關係圖中，直接波為通過原點、斜率為 $1/v_1$ 的直線；折射波則是不過原點（具截距時間 $t_i$）、斜率為 $1/v_2$ 的直線。兩線交點即為「交叉距離」，可用以計算上層地殼厚度。
在距離-走時關係圖中，直接波為通過原點、斜率為 $1/v_1$ 的直線；折射波則是不過原點（具截距時間 $t_i$）、斜率為 $1/v_2$ 的直線。兩線交點即為「交叉距離」，可用以計算上層地殼厚度。
📊 圖 3-2-3 說明【多層平坦介面的階梯式走時方程】
當地下構造推廣至多層速度模型（如沉積層、上地殼、下地殼）時，各界面激發的首波會形成多條具備不同截距時間與斜率的走時直線，使地震學家得以由上至下逐層反演波速。
當地下構造推廣至多層速度模型（如沉積層、上地殼、下地殼）時，各界面激發的首波會形成多條具備不同截距時間與斜率的走時直線，使地震學家得以由上至下逐層反演波速。
📊 圖 3-2-4 說明【傾斜地層界面與雙向放炮觀測】
若地下界面存在傾角，朝迎風（Updip）與朝順風（Downdip）兩個反方向觀測到的折射波斜率將不對稱。圖中解構了如何利用雙向視速度（Apparent Velocity）解出地層真實傾角與厚度。
若地下界面存在傾角，朝迎風（Updip）與朝順風（Downdip）兩個反方向觀測到的折射波斜率將不對稱。圖中解構了如何利用雙向視速度（Apparent Velocity）解出地層真實傾角與厚度。
📊 圖 3-2-6 說明【地函頂部 Pn 折射波實測震幅特徵】
展示實際震測觀測紀錄。隨著震央距拉遠，直接波（Pg）能量急速衰減，而沿莫霍面高速滑行的 Pn 折射波則超越直接波，在遠距離觀測站率先抵達，成為清晰的初達波分支。
展示實際震測觀測紀錄。隨著震央距拉遠，直接波（Pg）能量急速衰減，而沿莫霍面高速滑行的 Pn 折射波則超越直接波，在遠距離觀測站率先抵達，成為清晰的初達波分支。
📊 圖 3-2-9 說明【大陸與海洋地殼折射構造剖面對比】
利用大尺度折射地震勘探所得出的莫霍面深部構造圖，清晰展現了海洋地殼（較薄，~5-7 km）與大陸地殼（較厚，~30-50 km）在速度層序與厚度上的巨大板塊構造不連續性。
利用大尺度折射地震勘探所得出的莫霍面深部構造圖，清晰展現了海洋地殼（較薄，~5-7 km）與大陸地殼（較厚，~30-50 km）在速度層序與厚度上的巨大板塊構造不連續性。
🔍 CH 3.3 反射地震學 Reflection Seismology
利用近垂直入射的震波能量，捕捉地下地層界面因波阻抗（Acoustic Impedance）突變所產生的反射同相軸。
📊 圖 3-3-1 說明【單層平坦界面反射射線路徑】
展示震源發射後經地下介面反射至地表接收點的幾何路徑。利用幾何直角三角形，其走時方程在距離-走時圖（T-X圖）上呈現特徵的雙曲線軌跡。
展示震源發射後經地下介面反射至地表接收點的幾何路徑。利用幾何直角三角形，其走時方程在距離-走時圖（T-X圖）上呈現特徵的雙曲線軌跡。
📊 圖 3-3-2 說明【反射雙曲線與正常時差 (NMO) 物理概念】
反射波走時隨偏移距增加而產生的時間延遲，稱為正常時差（Normal Moveout, NMO）。本圖展示如何透過 NMO 修正（將雙曲線波形拉平），將非垂直入射走時精準轉換為垂直入射走時。
反射波走時隨偏移距增加而產生的時間延遲，稱為正常時差（Normal Moveout, NMO）。本圖展示如何透過 NMO 修正（將雙曲線波形拉平），將非垂直入射走時精準轉換為垂直入射走時。
📊 圖 3-3-6 說明【共同中間點 (CMP) 資料採集與多道疊加】
展示工業震測核心技術。利用多對震源與接收器重複抽樣地下同一個「共同中間點（CMP）」，經 NMO 修正後予以重合疊加（Stacking），可有效抑止隨機雜訊，倍增地下界面的訊噪比。
展示工業震測核心技術。利用多對震源與接收器重複抽樣地下同一個「共同中間點（CMP）」，經 NMO 修正後予以重合疊加（Stacking），可有效抑止隨機雜訊，倍增地下界面的訊噪比。
📊 圖 3-3-7 說明【實務工業反射地震時間剖面成像】
這是一張經過完整信號處理與偏移成像（Migration）後的反射時間剖面。圖中高連續性的明暗反射同相軸，細緻地描繪了地下沉積層序的物理波阻抗界面與微細的斷層滑移錯動。
這是一張經過完整信號處理與偏移成像（Migration）後的反射時間剖面。圖中高連續性的明暗反射同相軸，細緻地描繪了地下沉積層序的物理波阻抗界面與微細的斷層滑移錯動。
🌐 CH 3.4 球體地球中的地震波傳播
考量真實地球的球面幾何幾何效應時，地震射線參數與走時方程必須引入半徑因子進行修正。
📊 圖 3-4-1 說明【球體坐標系下的斯乃爾定律守恆】
在圓弧球層模型中，地震射線守恆參數修正為球體斯乃爾定律：$p = (r \sin i) / v = \text{常數}$。隨著射線向內深入、半徑 $r$ 減小，入射角 $i$ 會持續增大，直至在最低轉折點折返回地表。
在圓弧球層模型中，地震射線守恆參數修正為球體斯乃爾定律：$p = (r \sin i) / v = \text{常數}$。隨著射線向內深入、半徑 $r$ 減小，入射角 $i$ 會持續增大，直至在最低轉折點折返回地表。
📊 圖 3-4-3 說明【速度連續遞增介面下的彎曲射線軌跡】
當地心速度隨深度連續漸變（非斷崖突變）時，射線軌跡不會形成尖銳的折角，而是呈現圓滑、凹口朝上的弧形曲線。圖中繪製了不同出射角射線的下探深度與傳播弧長關係。
當地心速度隨深度連續漸變（非斷崖突變）時，射線軌跡不會形成尖銳的折角，而是呈現圓滑、凹口朝上的弧形曲線。圖中繪製了不同出射角射線的下探深度與傳播弧長關係。
📊 圖 3-4-4 說明【韋謝特-黑格羅茲 (Wiechert-Herglotz) 走時逆問題積分】
展示經典的一維速度反演幾何物理。只要在地表觀測到完整的角度走時變化率 $dT/d\Delta$（即射線參數 $p$ 隨震央距的分布），便可利用此積分方程直接推算地下任意深度處的絕對波速。
展示經典的一維速度反演幾何物理。只要在地表觀測到完整的角度走時變化率 $dT/d\Delta$（即射線參數 $p$ 隨震央距的分布），便可利用此積分方程直接推算地下任意深度處的絕對波速。
🌌 CH 3.5 體波走時與地球內部構造
實體波（P波與S波）在穿透全球尺度的核幔邊界（CMB）與內核邊界（ICB）時，會產生高度複雜的核心透鏡折射與反射震相。
📊 圖 3-5-1 說明【全球主要核心震相傳播射線全景圖】
這是一幅經典的全球射線幾何光學圖譜。詳細描繪了穿過地函的 P/S 波、穿過液態外核的 K 波（如 PKP）、穿過固態內核的 I 波（如 PKIKP）以及在核幔邊界反射的 PcP 射線路徑。
這是一幅經典的全球射線幾何光學圖譜。詳細描繪了穿過地函的 P/S 波、穿過液態外核的 K 波（如 PKP）、穿過固態內核的 I 波（如 PKIKP）以及在核幔邊界反射的 PcP 射線路徑。
📊 圖 3-5-2 說明【PREM 標準一維地球速度與密度分佈模型】
展示了初步參考地球模型（PREM）隨深度（或半徑）的變化。圖中可見在外核（2891 km以下）橫波速度（Vs）驟降至 0，這是外核為流體狀態最直接的物理鐵證。
展示了初步參考地球模型（PREM）隨深度（或半徑）的變化。圖中可見在外核（2891 km以下）橫波速度（Vs）驟降至 0，這是外核為流體狀態最直接的物理鐵證。
📊 圖 3-5-3 說明【全球大尺度走時主曲線分支圖 (Jeffreys-Bullen Map)】
橫軸為震央距角度（$0^\circ - 180^\circ$），縱軸為總走時（分鐘）。此圖彙整了數十年來全球測站的核心震相走時交織曲線，重疊的多個分支對應了深部大界面引起的震波聚焦現象。
橫軸為震央距角度（$0^\circ - 180^\circ$），縱軸為總走時（分鐘）。此圖彙整了數十年來全球測站的核心震相走時交織曲線，重疊的多個分支對應了深部大界面引起的震波聚焦現象。
📊 圖 3-5-4 說明【液態低速外核造成的 P 波地表觀測影帶 (Shadow Zone)】
由於地函基底波速顯著高於外核頂部波速，當 P 波入射核幔邊界時會產生強烈的向內折射（凹透鏡聚焦效應），導致地表在震央距 $103^\circ - 143^\circ$ 之間形成直達 P 波無法到達的「影帶」。
由於地函基底波速顯著高於外核頂部波速，當 P 波入射核幔邊界時會產生強烈的向內折射（凹透鏡聚焦效應），導致地表在震央距 $103^\circ - 143^\circ$ 之間形成直達 P 波無法到達的「影帶」。
📊 圖 3-5-7 說明【核幔邊界反射波 PcP 走時與 CMB 半徑估算】
本圖聚焦於在核幔邊界（CMB）垂直反射回地表的 PcP 震相。藉由高精密觀測 PcP 與直達 P 波的時間差（PcP - P），科學家得以將地球外核的半徑精確鎖定在約 3480 公里。
本圖聚焦於在核幔邊界（CMB）垂直反射回地表的 PcP 震相。藉由高精密觀測 PcP 與直達 P 波的時間差（PcP - P），科學家得以將地球外核的半徑精確鎖定在約 3480 公里。
🌋 CH 3.6 上部地函結構 Upper Mantle Structure
包含地函動力學最活躍的軟流圈低速帶，以及由橄欖石高溫高壓晶格相變所主導的過渡帶界面。
📊 圖 3-6-1 說明【軟流圈低速帶 (LVZ) 速度局部凹陷特徵】
此速度剖面圖放大展示了深度在 100-200 km 之間，縱波與橫波波速同時出現明顯凹陷。該區域對應「低速帶（Low-Velocity Zone）」，通常解釋為高溫高壓下岩石發生微量部分熔融（~1%）的軟流圈。
此速度剖面圖放大展示了深度在 100-200 km 之間，縱波與橫波波速同時出現明顯凹陷。該區域對應「低速帶（Low-Velocity Zone）」，通常解釋為高溫高壓下岩石發生微量部分熔融（~1%）的軟流圈。
📊 圖 3-6-2 說明【地函過渡帶 410km 與 660km 速度階梯跳躍】
展示地函過渡帶內速度隨深度呈梯形階躍跳升的特徵。在 410 km 處，橄欖石相變為瓦茲利石；在 660 km 處，林伍德石分解為橋石與方鎂石，此兩大突變面會在走時上引發顯著的三叉分支（Triplication）。
展示地函過渡帶內速度隨深度呈梯形階躍跳升的特徵。在 410 km 處，橄欖石相變為瓦茲利石；在 660 km 處，林伍德石分解為橋石與方鎂石，此兩大突變面會在走時上引發顯著的三叉分支（Triplication）。
🔩 CH 3.7 下部地函與地核結構
下部地函基底存在化學與熱不均勻的 D" 層，核幔邊界（CMB）為固液過渡面，而內核邊界（ICB）則重新過渡為固態結晶核。
📊 圖 3-7-1 說明【下部地函基底 D" 層之超低速帶 (ULVZ) 橫向不均勻性】
描繪緊鄰核幔邊界上方約 200-300 km 厚的特殊 D" 區域。近代層析成像在此處發現了局部波速下降達 10-30% 的「超低速帶（ULVZ）」，被視為主導全球地函熱柱湧升的發源根基。
描繪緊鄰核幔邊界上方約 200-300 km 厚的特殊 D" 區域。近代層析成像在此處發現了局部波速下降達 10-30% 的「超低速帶（ULVZ）」，被視為主導全球地函熱柱湧升的發源根基。
📊 圖 3-7-2 說明【核幔邊界 (CMB) 固液介面兩側物理參數斷崖】
以數據直觀呈現 2891 km 深度（CMB）兩側物理性質的劇烈跳躍：密度由 5.5 暴增至 10 $g/cm^3$；縱波波速自 13.7 驟降至 8.1 $km/s$，橫波完全歸零，為地球內部最顯著的物理化學邊界。
以數據直觀呈現 2891 km 深度（CMB）兩側物理性質的劇烈跳躍：密度由 5.5 暴增至 10 $g/cm^3$；縱波波速自 13.7 驟降至 8.1 $km/s$，橫波完全歸零，為地球內部最顯著的物理化學邊界。
📊 圖 3-7-3 說明【外核液態介面中 PKP 震相分支幾何計算】
詳解進入液態外核後再度折射而出的 PKP 射線路徑。由於外核流體的強烈聚焦效應，原本在影帶中消失的波動會在震央距大於 $143^\circ$ 的遠距離觀測站再次強烈現身，並形成 AB、BC、DF 多條走時分支。
詳解進入液態外核後再度折射而出的 PKP 射線路徑。由於外核流體的強烈聚焦效應，原本在影帶中消失的波動會在震央距大於 $143^\circ$ 的遠距離觀測站再次強烈現身，並形成 AB、BC、DF 多條走時分支。
📊 圖 3-7-4 說明【內核邊界 (ICB) 波動轉換與固態內核剪切波】
展示 5150 公里深內核邊界（ICB）的波相轉換。當 P 波自液態外核擊中固態內核時，除反射出 PKiKP 外，亦能局部轉換為內核中的剪切橫波（J波，震相為 PKJKP），直接證實地核中心為固態。
展示 5150 公里深內核邊界（ICB）的波相轉換。當 P 波自液態外核擊中固態內核時，除反射出 PKiKP 外，亦能局部轉換為內核中的剪切橫波（J波，震相為 PKJKP），直接證實地核中心為固態。
📊 圖 3-7-7 說明【固態內核地震波速各向異性 (Anisotropy) 觀測證據】
呈現穿過內核的 PKIKP 波速方向性差異：沿地球自轉軸（南北向）穿過內核的震波波速，比沿赤道面（東西向）穿過快了約 1-2%，揭示內核鐵晶格在高壓下具備方向性定向排列。
呈現穿過內核的 PKIKP 波速方向性差異：沿地球自轉軸（南北向）穿過內核的震波波速，比沿赤道面（東西向）穿過快了約 1-2%，揭示內核鐵晶格在高壓下具備方向性定向排列。
💎 CH 3.8 地函與地核的組成 Composition of Mantle and Core
結合高溫高壓實驗室岩石力學量測，將觀測到的地震波速、密度曲線，轉譯為真實地球內部的化學組成與礦物相成分。
📊 圖 3-8-2 說明【高溫高壓礦物物理相變與實測速度比對】
將金剛石砧壓機模擬之高溫高壓實驗數據疊加於標準 PREM 曲線。比對結果高度支持地函主要由富含鐵鎂的矽酸鹽岩石（如 Pyrolite 橄欖岩模型）組成，相變深度完美對應速度不連續面。
將金剛石砧壓機模擬之高溫高壓實驗數據疊加於標準 PREM 曲線。比對結果高度支持地函主要由富含鐵鎂的矽酸鹽岩石（如 Pyrolite 橄欖岩模型）組成，相變深度完美對應速度不連續面。
📊 圖 3-8-3 說明【純鐵與鐵合金在核心環境下的熔融溫度關係】
展示鐵鎳合金在數百萬大氣壓下的熔融相圖。由於地核實測密度略低於純鐵，說明核心富含硫、矽、氧等輕元素；外核因溫度高於熔點呈流體，內核則因巨大的靜水壓力壓制，提升熔點而維持固態結晶。
展示鐵鎳合金在數百萬大氣壓下的熔融相圖。由於地核實測密度略低於純鐵，說明核心富含硫、矽、氧等輕元素；外核因溫度高於熔點呈流體，內核則因巨大的靜水壓力壓制，提升熔點而維持固態結晶。
📊 圖 3-8-5 說明【地球全尺度品質密度與重力加速度 (g) 深度模型】
展示密度隨深度呈階梯狀上升，而重力加速度（g）因物質分布反而在地函中大致守恆於 10 $m/s^2$，直至進入地核內部才因核心品質高度向內集中，呈線性遞減至地心處的 0。
展示密度隨深度呈階梯狀上升，而重力加速度（g）因物質分布反而在地函中大致守恆於 10 $m/s^2$，直至進入地核內部才因核心品質高度向內集中，呈線性遞減至地心處的 0。
📊 圖 3-8-6 說明【地震波非彈性衰減品質因子 (Q值) 全球分布】
呈現描述介面吸收衰減波動能量內能消耗的指標——品質因子 $Q$ 值。軟流圈低速帶因部分熔融與晶格蠕變，表現出極低的 $Q$ 值（高衰減特徵）；而下部地函與內核則具備高 $Q$ 值。
呈現描述介面吸收衰減波動能量內能消耗的指標——品質因子 $Q$ 值。軟流圈低速帶因部分熔融與晶格蠕變，表現出極低的 $Q$ 值（高衰減特徵）；而下部地函與內核則具備高 $Q$ 值。
📊 圖 3-8-7 說明【橫波與縱波波速比值 (Vp/Vs) 與岩石岩性判別】
探討 $Vp/Vs$ 比值（與帕松比密切相關）作為地下岩石分類的關鍵指標。長英質花崗岩與基性玄武岩/輝長岩在比值分布上具備顯著不連續性，是無開挖狀態下鑑別深部地殼化學成分的利器。
探討 $Vp/Vs$ 比值（與帕松比密切相關）作為地下岩石分類的關鍵指標。長英質花崗岩與基性玄武岩/輝長岩在比值分布上具備顯著不連續性，是無開挖狀態下鑑別深部地殼化學成分的利器。
📊 圖 3-8-9 說明【三維全球地震層析成像 (Tomography) 與地函熱對流】
打破一維層狀模型的物理侷限，展示現代三維地震層析成像全景。圖中藍色偏快代表低溫（古老板塊俯衝殘骸），紅色偏慢代表高溫（地函熱柱或湧升流），具體勾勒出地球內部動態熱對流景象。
打破一維層狀模型的物理侷限，展示現代三維地震層析成像全景。圖中藍色偏快代表低溫（古老板塊俯衝殘骸），紅色偏慢代表高溫（地函熱柱或湧升流），具體勾勒出地球內部動態熱對流景象。