關(guān)于大地電磁阻抗張量分解技術(shù)之Bahr分解穩(wěn)定性與實用性的簡要分析

（安徽省勘查技術(shù)院，安徽 合肥 230031）

阻抗張量分解一直以來都是大地電磁（MT）資料處理中的關(guān)鍵技術(shù)，在國內(nèi)外得到廣泛的應(yīng)用。張量分解方法眾多，Bahr分解作為常規(guī)方法之一，國內(nèi)有學者使用該方法應(yīng)用于MT實測資料的處理，其目的多是希望能夠消除MT觀測資料的局部畸變效應(yīng)并獲得區(qū)域構(gòu)造的主軸方位角及二維偏離度等相關(guān)參數(shù)。近年來有研究認為Bahr分解法不夠穩(wěn)定，但是并沒有做詳細的說明，這是本文研究的出發(fā)點，作者將依據(jù)已經(jīng)發(fā)表的論文和前人的經(jīng)驗，對Bahr分解技術(shù)做進一步的研究。

1 Bahr分解方法

1.1 Bahr分解基本原理

Bahr分解[2]方法提出于20世紀90年代初，是Bahr在Swift分解[1]的基礎(chǔ)上提出的一種阻抗張量分解方法。

現(xiàn)假設(shè)測量坐標系（x，y，z）中的大地電磁觀測阻抗Z可表示為：

當測量坐標系（x，y，z）旋轉(zhuǎn)至區(qū)域二維構(gòu)造的主軸方位時，在旋轉(zhuǎn)后的坐標系（x'，y'，z'）下，阻抗張量Z可表示為：

巴爾分解法認為在（x'，y'，z'）坐標系中，測量的阻抗張量矩陣的每列，兩個元素的相位差為零[3]。利用阻抗旋轉(zhuǎn)原理和在構(gòu)造主軸方位角處列元素相位相同這一特點，可以將主軸方位角的求解表示成如下形式：

巴爾分解對于區(qū)域二維構(gòu)造中超過二維的部分會被當作局部畸變異常來處理，從而使區(qū)域構(gòu)造和局部構(gòu)造都發(fā)生一定的變化，Bahr（bahr，1991）對原分解公式進行了修正，引入一個相位偏離角δ來補償，從而獲得修正的巴爾分解法。它將引入的偏差角同時加在阻抗張量的兩列上，并將觀測的阻抗張量修改成如下形式：

由（4）、（5）式可知，θ和δ均有兩個值，它們分別對應(yīng)最小相位偏差和最大相位偏差。在實際應(yīng)用過程中，應(yīng)該選擇與最小的相位偏差角δ所對應(yīng)的θ作為主軸方向角，并利用阻抗張量旋轉(zhuǎn)關(guān)系，計算出分解后的視電阻率和相位。

1.2 理論數(shù)據(jù)驗證

利用Python語言編制巴爾分解算法程序，通過對已發(fā)表的論文（Pracser，1999）中的試驗數(shù)據(jù)進行了試算，計算結(jié)果基本一致。

試驗數(shù)據(jù)：

試驗數(shù)據(jù)結(jié)果：

本文計算結(jié)果：

2 實測數(shù)據(jù)計算

利用編制的Bahr分解實用程序?qū)戏逝璧貎?nèi)一條MT剖面的實測數(shù)據(jù)進行了計算，計算的主要參數(shù)包括視電阻率和相位、主軸方位角以及二維偏離度，從以上三個維度對該分解技術(shù)的穩(wěn)定性和有效性進行了簡要的分析。

2.1 視電阻率和相位

視電阻率和相位是初步判斷大地電磁資料優(yōu)劣的一種較為直觀的方式，本次實測MT點位超過100個，采用Bahr分解技術(shù)進行分解后的資料整體上不是很好，0.1Hz以后的低頻段分散性加強（如圖1），在資料比較差的情況下，同時存在一部分MT測點的部分頻點資料無法計算的情況。

2.2 主軸方位角

主軸方位角是確定區(qū)域構(gòu)造走向的關(guān)鍵參數(shù)，該參數(shù)對于資料后期的處理和解釋具有重要的意義。

作者編制了畫圖程序用于繪制MT單測點的方位角玫瑰圖，圖2中MT實測資料的方位角玫瑰圖對應(yīng)134號、196號、250號測點，結(jié)果表明主軸方位角大體上位于北偏東30°±10°左右。

2.3 二維偏離度

二維偏離度是大地電磁法資料處理中衡量測區(qū)電性構(gòu)造特征的主要參數(shù)，通過二維偏離度可以得知該測點局部異常和區(qū)域異常特征及資料受局部畸變影響的程度。

MT剖面位于盆地內(nèi)，圖3大體上能夠反映該剖面的結(jié)構(gòu)維性特性：周期小于10s的數(shù)據(jù)二維偏離度基本上在0.2以下，表明該處電性結(jié)構(gòu)二維性較好，測點受到的局部畸變影響較小；低頻部分（周期大于100s）的二維偏離度明顯增大，反映深部地質(zhì)構(gòu)造三維性較強。

圖1 136號測點資料Bahr分解前（a）、（b）后的對比

圖2 134，196，250號測點方位角玫瑰圖

圖3 Bahr二維偏離度

3 結(jié)論

利用Bahr分解技術(shù)對一條實測MT剖面數(shù)據(jù)進行了初步的計算分析，結(jié)果表明：采用Bahr分解后獲得的視電阻率和相位資料整體上并不理想，并且在原始資料質(zhì)量較差的情況下使用該技術(shù)并不能獲得更好的結(jié)果，因而在數(shù)據(jù)處理過程中不建議直接使用該方法計算的視電阻率和相位參數(shù)。Bahr分解計算的主軸方位角以及二維偏離度參數(shù)可以較好的確定區(qū)域地質(zhì)體的電性走向和構(gòu)造維數(shù)，實用性較強，并能夠為后期的數(shù)據(jù)處理和地質(zhì)解釋提供理論支持。