E-Ex廣域電磁法靜態校正方法探討

陸柏樹

(湖北省地質局 水文地質工程地質大隊,湖北 武漢 430051)

受地形和地表電性不均勻體的影響，電磁測深數據在靜態效應的影響下，導致視電阻率曲線產生嚴重的靜態畸變，對后期數據反演與解釋的準確性有較大影響。

為消除靜態效應的影響，可控源音頻大地電磁法(CSAMT)常用的靜態校正方法有相位空間濾波法、空間濾波法、小波多尺度分析法等。

大地電磁測深法(MT)一般采用空間濾波法進行靜態校正。

和CSAMT、MT一樣，E-Ex廣域電磁法(WFEM)測試數據也存在靜態效應。E-Ex廣域電磁法測試的是單分量數據，沒有測量相位曲線，因此常用的靜態校正方法有空間濾波法、小波多尺度分析法等。

以上幾種靜態校正方法都是針對視電阻率進行靜態校正，其成立的前提條件是視電阻率與電場成正比。但實際上這個前提條件不一定都成立。因此，目前所用的靜態校正方法是存在問題的。

本文從電磁測深法中的靜態效應產生的原因出發，分析目前常用靜態校正方法所導致的問題，試圖給出新的靜態校正方法，并以E-Ex廣域電磁法為例，進行理論計算和模擬對比。

1 電磁測深法中的靜態效應[1]

根據電磁場傳播規律，在不均勻體的界面上，所有穿過邊界的場和位都是連續的，只有電感應強度的法向分量不連續。如果已知不均勻體的形狀、導電電性等，可以據此初步計算出不均勻體表面的面電荷密度的大小。這個面電荷密度雖然小，但它對電場的作用卻是不可忽略的。

當趨膚深度比不均勻體的尺寸大得多時，便能察覺到這種表面電荷的影響。這就意味著，在地表或地表附近小的二維或三維不均勻體，可能對整個電場測量都有影響。當然，較深的物體也能引起靜態位移，但大地電磁場在觀測過程中，靜態位移現象主要是由淺部電阻率的橫向變化引起的。

在地形平坦的二維介質條件下，電磁場建立后，地表淺層的不均勻體表面上的電荷分布可能使電場數據向上或向下移動一個數值，表現為在lgE-lgf曲線上發生向上或向下的平移[2-3]，這個平移的數值與頻率無關。由于地下介質和空氣的磁導率相差很小，淺地表的不均勻體幾乎不會對磁場的分量產生靜態影響，對磁場觀測數據的影響不大。

根據靜態位移后的電場值計算出的視電阻率也發生了移動，但相位曲線所受影響不大。

目前頻率域電磁法基本上都采用針對視電阻率的靜態校正方法。這種方法成立的前提是根據電場值計算出的視電阻率，在lgρ-lgf曲線上表現為整體的平移。

2 常用靜態校正方法及存在的問題

2.1 常用靜態校正方法

2.1.1 相位空間濾波法

相位空間濾波法是根據無靜態效應的阻抗相位數據積分來求得視電阻率靜態校正數據的方法。其理論依據是觀測的電場雖然存在靜態效應，但其相位不存在靜態效應，以此出發推導出一套靜態校正方法。

2.1.2 空間濾波法

利用空間濾波法做靜校正的基本思想是假定地下電性異常體或地質構造引起的視電阻率沿測線是平緩漸變的；而地表局部電性不均勻體或局部地形不平則會引起視電阻率的急劇變化。若設計一種低通濾波器沿測線做空間濾波，則可壓制高頻的靜態效應。

空間濾波法包括平衡移動平均空間濾波法、定長滑動平均空間濾波法、自適應空間濾波法、中值空間濾波法等等。這些方法都是對計算出來的視電阻率曲線組，選取合適的頻點、頻段視電阻率曲線組進行低通濾波，從而獲取靜態校正系數，對視電阻率曲線進行平移。

2.1.3 小波多尺度分析法

小波多尺度分析法是采用新的數學手段，對視電阻率曲線組進行低通濾波，其本質也是空間濾波方法。

以上較為常用的三種方法，都適用于CSAMT的靜態校正處理。由于E-Ex廣域電磁法測試的是單分量測試數據，沒有測量相位曲線，只適用于空間濾波法、小波多尺度分析法。

2.2 靜態校正方法存在的問題

目前靜態校正方法基本上都是針對視電阻率曲線進行靜態校正。而視電阻率曲線是根據存在的靜態效應的電場曲線計算而來的。

根據前面所述，電場的靜態效應在lgE-lgf曲線上表現為產生向上或向下的平移，即電場的靜態效應基本可以看成所有觀測頻率上的電場都放大了一個系數。

那么就存在一個問題，如果根據存在靜態效應的電場計算視電阻率，在lgρ-lgf曲線上是否也表現為整體的平移呢？

2.2.1 CSAMT

CSAMT的視電阻率計算公式如下：

ρs=|Ex/Hy|2/ω/μ

(1)

從式中可以看出,CSAMT的視電阻率與電場呈平方關系。進一步分析可知，針對視電阻率的靜態校正，從理論上講是成立的，因為ρ=f(a*Ex)=a*a*f(Ex)。若針對電場的校正系數是a，則針對電阻率的校正系數是a*a，也表現為在lgρ-lgf曲線上的整體平移。

那么，CSAMT的觀測數據，采用空間濾波法和小波多尺度分析法，均可以較好地消除靜態效應。

2.2.2 E-Ex廣域電磁法

廣域電磁法是由何繼善院士提出并發展起來的一種勘探方法[4]。廣域電磁法英文全稱是：Wide Field Electromagnetic Methods，簡寫為WFEM。

均勻半空間表面水平電偶源的電磁場Ex分量表達式為[4-6]:

(2)

E-Ex廣域電磁法假定觀測到的電場Ex是均勻半空間表面水平電偶源的電場Ex，然后根據式(2)可求視電阻率。

式(2)等號右邊中的k也含有大地電阻率參數，因此，不能直接從式(2)中求解出視電阻率，需要采取迭代的方法進行求解。

也就是說，E-Ex廣域電磁法視電阻率與電場呈復雜的、非線性關系，那么針對視電阻率所進行的靜態校正，從理論上講是不成立的。

綜上所述，在MT、CSAMT的數據預處理中，針對視電阻率的靜態校正法，從理論上講是成立的；但對廣域電磁法(WFEM)采用針對視電阻率的靜態校正法是不成立的。

3 新的電磁法靜態校正方法

根據前面的描述，目前廣域電磁法針對視電阻率所采用的靜態校正方法是不成立的。因產生靜態效應的根源在于電場受到了影響，所以，新的靜態校正方法應該是針對電場進行靜態校正，然后利用靜態校正后的電場計算視電阻率。

至于具體采用何種數學方法進行靜態校正，則不在本文討論范圍。下面就以E-Ex廣域電磁法為例，通過理論計算對比，看看兩種靜態校正方法所帶來的差別。

4 E-Ex廣域電磁法靜態校正理論計算

E-Ex是指采用水平電偶源測量Ex分量方法。以下理論計算的E-Ex廣域電磁測深場源參數均為收發距8 000 m，偶極源赤道上觀測。

4.1 模型1

本模型參數如下:層數5層，h1=20 m、h2=50 m、h3=300、h4=600、h5=∞，ρ1=15、ρ2=35、ρ3=7、ρ4=35、ρ5=100 Ω·m。

該模型的廣域電磁法正演視電阻率ρs-視深度D曲線和模型ρ-H曲線見圖1所示。視深度D取用視電阻率和頻率計算的視波長1/10，其具體含義是該深度至地表的綜合電性反映,就是該頻率所測得的視電阻率。

圖1 模型1的WFEM正演ρs-D曲線和模型Fig.1 The model 1 and WFEM forward ρs-D curve

將正演數據中的Ex分別放大30、10、5、2、1、1/2、1/5、1/10、1/30倍，計算出相應的視電阻率曲線，將這9條視電阻率曲線再相應地除以對應的放大倍數，并繪制ρs-D曲線在圖2中進行對比。

圖2 模型1針對電阻率靜態校正的情況Fig.2 In view of the static correction for resistivity of the model 1

圖2中，從上至下依次為Ex放大30、10、5、2、1、1/2、1/5、1/10、1/30倍后，將電阻率縮小相應倍率的ρs-D曲線。很明顯，這種數據處理和靜態校正方法，已經使曲線形態發生明顯的改變。

4.2 模型2

本模型參數如下:層數5層，h1=58 m、h2=25 m、h3=85、h4=546、h5=∞，ρ1=218、ρ2=138、ρ3=67.5、ρ4=369、ρ5=1 000 Ω·m。

該模型的廣域電磁法正演視電阻率ρs-視深度D曲線和模型ρ-H曲線見圖3所示。

將正演數據中的Ex分別放大20、10、5、2、1、1/2、1/5、1/10、1/20倍，計算出相應的視電阻率曲線，將這9條視電阻率曲線再相應地除以對應的放大倍數，并繪制ρs-D曲線在圖4中進行對比。

圖4中，從上至下依次為Ex放大20、10、5、2、1、1/2、1/5、1/10、1/20倍后，將電阻率縮小相應倍率的ρs-D曲線。同樣的，這種數據處理和靜態校正方法，已經使曲線形態發生明顯的改變。

圖3 模型2的WFEM正演ρs-D曲線和模型Fig.3 The model 2 and WFEM forward ρs-D curve

圖4 模型2針對電阻率靜態校正的情況Fig.4 In view of the static correction for resistivity of the model 2

為進一步說明其曲線的變化特征，給出模型2靜態校正的ρs-f曲線組,見圖5所示。

圖5 模型2針對電阻率靜態校正的情況Fig.5 In view of the static correction for resistivity of the model 2

4.3 對比情況初步分析

顯而易見，以上模型的計算結果,采用針對Ex進行靜態校正的方法,對電場放大縮小同樣的比例,將完全還原成計算模型的結果。

以上對Ex進行放大再計算視電阻率的模擬情況，類似于通過直流電法或瞬變電磁法獲得淺部地層的電阻率，并將淺部曲線校正一致。

對于E-Ex廣域電磁法曲線，可以看出，如果用帶有靜態效應的電場，計算出的視電阻率曲線形態將發生改變。理論計算對比時,可以較清晰地發現這一種現象，但實際上，在實測數據中，一旦采用針對視電阻率的空間濾波法，選定的頻點已經將非線性結果帶入，則這條曲線將可能與實際情況相差較遠。

根據理論計算和實際測試數據對比分析，測線中靜態校正系數>1.5或<0.7的實測曲線，都將受到較大影響，一旦靜態校正選擇的頻點位于畸變段，則難以得到接近真實情況的反演結果。

對E-Ex廣域電磁法曲線來說，如果在勘探深度內對應的頻段出現了非線性現象，則采用針對電阻率進行靜態校正的方法，電阻率曲線將發生畸變，導致反演的厚度和電阻率都不正確，甚至無法較好地進行反演擬合。若選擇校正的頻點位于畸變頻段時，還可能使淺部的曲線形態發生錯誤的平移，進而影響解譯結果。

另外，如果野外觀測采用的MN較大，則以上現象將有所緩解。分析原因可能是MN較大后，地形影響變小，局部不均勻體也相對影響變小。

一個比較有趣的現象是，如果模擬計算的頻率進一步降低，針對電阻率進行靜態校正的曲線尾枝也是重合的，在ρs-D曲線中，影響較大的不重合區域大致在D=800～80 000 m之間。因D與電阻率和頻率都有關，所以，電阻率高的地層，出現這一現象的頻率會相對較高，反之亦然。這一現象還需要進一步進行研究。

5 生產實例

靜態校正不當的例子，在實際工作中可以找到很多，尤其在高阻地層中，這種現象很普遍。

在筆者前期編制的廣域電磁法反演軟件中，也采用的是針對視電阻率的靜態校正方法。實際工作發現，在高阻地層中觀測的曲線往往差別極大，是因為MN極距較小后，靜態效應顯著增大所致。而且高阻地段在相對較高的頻段就出現廣域電阻率計算的非線性現象，所以，針對視電阻率進行靜態校正的方法的不合理性就顯得相對突出。

筆者在認真分析了原因后，把軟件修改為針對電場進行靜態校正的方法，解決了實際生產問題。

湖北省羅田縣白廟河地熱田分布于片麻巖地區，完整的片麻巖電阻率極高，可>10 000 Ω·m，表層風化層極不均勻，導致很多測點的靜態效應較大。采用廣域電磁法進行勘探，收發距7 800 m，觀測頻率8 192～0.75 Hz，供電電流0.5～10 A，MN=20 m。

圖6是一個測點的校正情況，電場的靜態校正系數為0.27，兩種靜態校正方法所得的ρs-D曲線差別較大。如果采用針對電阻率的校正方法，曲線難以反演擬合上，但采用針對電場的校正方法，則可以較好地進行反演擬合。

圖6 實測數據采用不同靜態校正方法對比情況(校正系數=0.27)Fig.6 Comparison of different static correction methods for measuring points

圖7 實測數據采用不同靜態校正方法對比情況(校正系數=2.26)Fig.7 Comparison of different static correction methods for measuring points

圖8 針對電場靜態校正后的反演實例Fig.8 An inversion example which adopt the static correction method for electric field

圖7是一個測點的校正情況，電場的靜態校正系數為2.26，同樣，兩種靜態校正方法所得的ρs-D曲線差別較大。如果采用針對電阻率的校正方法，曲線難以反演擬合上，但采用針對電場的校正方法，則可以較好地進行反演擬合。

采用針對電場的靜態校正方法，反演所得的湖北省羅田縣白廟河地熱田1線的地電斷面圖，見圖8。解譯結果與實際情況吻合較好。在230 m樁號附近成功打出溫度57 ℃地熱，水量>2 000 m3/d，水位高出地表5 m。

6 結語

(1) 目前廣域電磁法靜態校正方法一般是針對視電阻率曲線進行的，但電場強度和根據其計算得來的視電阻率之間呈復雜的非線性關系，理論分析和模擬計算結果表明，針對視電阻率進行靜態校正的方法存在一定程度的不合理性。

(2) 本文針對這一情況，提出了新的針對電場進行靜態校正的方法，即先對電場進行靜態校正，然后采用靜態校正后的電場計算視電阻率。并在筆者編制的廣域電磁法處理軟件中，采用了這一種靜態校正方法。

(3) 生產實踐表明，本文提出的靜態校正方法較好地解決了實際生產問題。

(4) 為減小地形和淺地表不均勻體的影響，在實際觀測過程中，MN不宜取得過小。

[1] 何繼善.可控源音頻大地電磁法[M].長沙:中南大學出版社,1999.

[2] 李金銘.地電場與電法勘探[M].北京:地質出版社,2005.

[3] 程志平.電法勘探教程[M].北京:冶金工業出版社,2007.

[4] 樸化榮.電磁測深法原理[M].北京:地質出版社,1990.

[5] 何繼善.廣域電磁法和偽隨機信號電法[M].北京:高等教育出版社,2010.

[6] 何繼善.頻率域電法的新進展[J].地球物理學進展,2007,22(4):1250-1254.