MT資料處理與解釋技術
介紹
大地電磁法（MT）是一種被動式地球物理勘探技術，它利用地球電場和磁場的自然變化來繪製地下電阻率分佈圖。電阻率與岩石類型、流體含量、溫度以及斷層或蝕變帶等地質構造密切相關。因此，大地電磁法廣泛應用於地熱勘探、礦化作用、沉積盆地（油氣）勘探以及構造和地殼研究。
然而，原始大地電磁資料很少能直接使用。它需要一系列處理步驟來提高訊號雜訊比，然後基於地質一致的電阻率模型進行解釋。本文討論了大地電磁資料的處理和解釋技術，從資料收集和品質控製到阻抗處理，再到反演和地質解釋。
MT資料庫：電場、磁場和阻抗張量
在大地電磁測深（MT）測量中，感測器記錄電場成分（Ex、Ey）和磁場成分（Hx、Hy，有時也會記錄Hz）隨時間的變化。根據頻域中電場和磁場的關係，可以得到阻抗張量：
\[
\begin{bmatrix} E_x \\ E_y \end{bmatrix}
=
\begin{bmatrix} Z_{xx} & Z_{xy} \\ Z_{yx} & Z_{yy} \end{bmatrix}
\begin{bmatrix} H_x \\ H_y \end{bmatrix}
\]
這些張量元素隨後被導出為關鍵的解釋參數，例如視電阻率 (ρa) 和相位 (φ)。通常，高品質的數據在不同頻率下，ρa 和 φ 的變化趨勢應平滑且一致，且誤差範圍合理。
第一階段：品質控制與預處理
良好的大氣壓力微波輻射處理始於品質控制，而品質控制則從現場數據開始。此階段包括：
1. 時間序列檢查：找出尖峰、漂移、感測器飽和或週期性幹擾。
2. 文化噪音：常見的噪音源有 50/60 Hz 的電力線、火車、工業活動、電圍籬和電信。
3. 電極和接地接點的狀況：接觸電阻過高會降低 Ex/Ey 的品質。
4. 感測器方向和位置：感測器方位角誤差會導致解釋失真，尤其是在定向結構研究中。
在這個階段，通常會切斷壞資料段，進行偏移校正，如果使用遠端參考站，也會進行時間同步。
第二階段：頻域變換和頻譜估計
由於MT分析的是頻域中E和H之間的關係，因此需要使用傅立葉變換等頻譜技術對時間序列進行變換。然後，為了確保統計穩定性，將資料分割成視窗（片段）。一些常用的技術包括：
– 透過加窗和加錐（例如漢寧窗）來減少光譜洩漏。
– 將視窗平均，以獲得更穩健的頻譜估計。
– 利用交叉功率譜建立 E 與 H 的關係。
目標是在從高頻（淺層）到低頻（深層）的寬頻率範圍內獲得準確的阻抗估計值。
第三階段：穩健處理與遠程參考
MT面臨的挑戰之一是特定頻道內噪音往往具有相關性。因此，人們開發了穩健性處理技術來降低異常值和不良資料段的影響。魯棒性方法採用迭代加權來抑制不一致資料的影響。
另一種非常重要的技術是遠端參考（RR）。其原理是在遠離本地噪音源的測站測量磁場。透過將主測站的電場強度（E）和磁場強度（H）與遠程參考站的磁場強度（H）進行關聯，可以降低本地雜訊對磁場強度（H）的影響。遠端參考技術在雜訊較大的區域，例如居民區或基礎設施附近，尤其有效。
此階段的輸出通常為：
– Zxy 和 Zyx 的 ρa 和 φ 曲線
– 誤差線（標準差）
– 初級訊號之間的相干性
第四階段：降噪與失真處理（靜態偏移）
除了雜訊之外，大地電磁法也經常受到靜態偏移的影響，靜態偏移是指由於淺層非均質性（例如薄黏土層、礫石層或電極條件）導致的視電阻率曲線的垂直偏移。靜態偏移不會顯著改變相位，但會使視電阻率ρa發生倍增偏移。
靜態輪班處理可以透過以下方式完成：
1. 基於TDEM/CSAMT資料的校正作為淺層電阻率控制。
2. MT-TDEM聯合反演，使淺層模型更具約束性。
3. 穩健的移位反演方法，例如允許每個台站設定移位參數。
此外，還進行了頻率異常值檢測：曲線上某些偏離較大的點通常會被丟棄或賦予較大的誤差。
第五階段：維度與打擊分析
在反演之前，確定地下結構是1D、2D或3D至關重要。這會影響反演方法的選擇和解釋。常見的分析方法包括：
– 偏度參數（例如 Bahr 偏度）用於評估 3D 水準。
– 相位張量用於觀察結構的主要方向，而不受靜態位移的影響。
– 走向分析，確定二維構造的主要走向方向。
如果資料呈現明顯的二維特徵，通常會進行沿走向的阻抗張量旋轉，以明確主成分（TE/TM）。如果三維特徵明顯，則優選進行三維反演。
步驟 6：MT 資料反演（1D、2D、3D）
大地電磁（MT）的定量解釋通常透過反演實現，即找到最符合資料響應的電阻率模型。 MT反演是非線性的且不適定的，需要正規化以防止模型變得「失控」。反演類型：
1. 一維反演：適用於層狀區域（例如水平沉積物）。速度快，但限制較大。
2. 二維反演：適用於斷層、地塹或定向地熱系統等狹長結構。
3. 三維反演：對於複雜的地質情況來說最現實，但需要密集的數據、大量的計算和嚴格的品質控制。
逆目標函數通常考慮資料擬合誤差和模型粗糙度：
\[
\Phi = \Phi_d + \lambda \Phi_m
\]
其中 λ 為正規化參數。 λ 的選擇至關重要：太小會導致模型過於粗糙（過度擬合），而太大則會導致模型過於平滑，從而遺漏重要的地質特徵。
步驟 7：反演結果和靈敏度評估
反演之後，必須對模型進行評估，而不僅僅是目測。評估內容包括：
– 不匹配度（RMS）：是否與目標相符（例如，RMS ~ 1–2，取決於誤差的定義）。
– 各站點資料曲線與模型響應的比較。
– 解析度測試：例如棋盤格測試或靈敏度圖。
– 先前模型的影響：嘗試幾個初始模型，以檢驗結果的穩定性。
如果一個小的參數發生變化時，模型發生劇烈的變化，這意味著在解釋模型結果時應該更加謹慎，並且可能需要更多的數據。
地質解釋：電阻率與地下系統的關係
電阻率並非直接的“岩石類型”，而是一種綜合的物理響應。然而，一些通用模式通常可作為指導原則：
– 低電阻率（導電性）：黏土、熱液蝕變（黏土蓋）、鹽流體帶、石墨或硫化物礦物。
– 高電阻率（電阻性）：塊狀火成岩、乾岩、矽化帶或結晶基底。
例如，在地熱勘探中，經典模型通常顯示：
1. 頂部覆蓋導電黏土帽
2. 下方儲層的阻力較大
3. 構造控制上升流區（斷層）
4. 熱源有時會根據岩性和溫度表現為電阻異常或複雜異常。
好的解釋總是與其他數據連結起來的：地表地質、熱表現、地球化學、重力、地震或井資料。
關閉
大地電磁（MT）資料處理與解釋技術環環相扣：從時間序列品質控制（QC）、頻譜估計、穩健/遠端參考處理、靜態偏移校正、維度分析，到二維/三維反演和解析度評估。 MT的成功不僅取決於反演軟體，還取決於測量品質、雜訊控制以及與地質背景的整合。
憑藉其嚴謹的工作流程和基於多數據的解釋，MT 是一種強大的工具，可用於繪製地下電阻率結構圖，並輔助複雜的勘探和地球科學研究中的決策。
