考慮IP效應的可控源電磁法電磁響應畸變分析

2020-06-10 01:41:12潘威

物探化探計算技術(shù) 2020年2期

潘威

(中鐵一院甘肅鐵道綜合工程勘察院有限公司，蘭州 730000)

0 引言

巖礦石是復雜的多相介質(zhì)，電阻率、極化率等參數(shù)是巖礦石本身的固有物理屬性，以激發(fā)極化效應理論為基礎(chǔ)的頻譜激電法，是一種有效的地球物理勘探方法，對此方法國內(nèi)、外皆有研究[1—9]。可控源電磁法(CSEM)現(xiàn)已廣泛應用于油氣、礦場等資源勘探中，目前，國內(nèi)、外學者對CSEM的研究主要只考慮純電磁感應[10—15]，但是當?shù)叵陆橘|(zhì)中含有極化層時，會對CSEM的電磁響應產(chǎn)生影響，而用常規(guī)的CSEM理論已無法進行有效解釋。因此筆者首先建立不同的三層地電模型(僅考慮中間層極化)，將極化參數(shù)加入到有限長導線頻率域電磁法一維正演中去，對含IP效應、純電磁感應的Ex、Hy、Hz電磁響應進行理論分析，為含IP效應的人工源頻率域測深方法提供一定的理論參考。

1 不同地電模型CSEM一維正演

1.1 CSEM基本理論

設發(fā)射導線長度為L，鋪設于地表，取坐標原點O，接收點A。如圖1所示，圖中r為收發(fā)距。

圖1 均勻半空間下的有限長導線源示意圖Fig.1 Schematic diagram of finite-length traverse source in Geodetic medium

則純電磁感應的有限長導線源的電場x分量、磁場y分量以及磁場z分量表達式如下：

J1(mr)cos2θ/r]dm

(1)

m2J0(mr)cos2θ]dm

(2)

(3)

(4)

本次測試將設置導線長為2 km，發(fā)射點坐標為(0,0),接收點坐標為(5 000，9 000)，電流大小為1 A。

1.2 含激電數(shù)值模擬

描述巖礦石中的激發(fā)極化效應的激電模型有多種，筆者僅對應用最為廣泛的Cole-Cole模型進行討論。

(5)

式中：ω= 2πf為圓頻率；ρ(ω)為頻率為ω時的巖、礦石復電阻率(Ω · m)；ρ0為零頻電阻率(Ω ·m)；m為極化率；τ為時間常數(shù)(s)；C為頻率相關(guān)系數(shù)。經(jīng)過推導[16]可得Cole-Cole模型的振幅表達式為式(6)。

(6)

將式(6)代替式(1)、式(2)、式(3)中的地層電阻率ρ，從而進行含IP效應的可控源頻率域一維正演計算，并利用式(7)進行電磁場響應畸變分析。

(7)

式中：|F|為含IP效應的頻率域電磁場響應的振幅值；|FB|為純電磁感應的電磁場響應的振幅值；δ為最大相對異常值或最大相對畸變值。

1.3 均勻半空間正演模擬

為驗證本文所用數(shù)值算法的正確性，采用均勻半空間情況下的水平電偶極子的電磁場Ex，Hy以及Hz分量解析解表達式：

(8)

式中：PE=Idl，K1、K0、I0、I1表示虛宗量貝塞爾函數(shù)。然后計算出純電磁感應以及考慮IP效應的解析解[16]，如圖2、圖3中虛線所示。均勻半空間地電參數(shù)為ρ0= 150 Ω·m，m= 0.5，C= 0.45，τ= 150 s。

在層狀介質(zhì)情況下，式(1)、式(2)、式(3)一般是沒有解析解的，因此可采用快速漢克爾變換方法對其進行計算，進而得到各電磁場分量的數(shù)值表達式，其具體算法可參考文獻[16]。計算數(shù)值解時每層的地電參數(shù)皆為均勻半空間地電參數(shù)。根據(jù)圖2、圖3可知，均勻半空間情況下含IP效應以及純電磁感應的數(shù)值解與解析解在頻率10-2Hz ～103Hz 的范圍內(nèi)擬合很好，說明該數(shù)值算法的正確性。

圖2 純電磁感應的數(shù)值解與解析解Fig.2 Numerical solution and analytical solution only considering electrical variation

圖3 考慮IP效應的數(shù)值解與解析解Fig.3 Numerical solution and analytical solution with the IP effect

圖4 M-200的Ex響應曲線Fig.4 The thickness of 200m of Ex curve

圖5 M-80的Ex響應曲線Fig.5 The thickness of 80m of Ex curve

1.4 層狀模型正演模擬

為進行含IP效應的三層層狀介質(zhì)模型正演，設置中間層為極化層，層厚分別為200 m、80 m。第一層厚度等于中間層的厚度除以厚度埋深比(極化層的厚度與其埋深的比值)，每層極化參數(shù)如表1所示。

圖4～圖9是不同厚度埋深比情況下考慮IP效應和純電磁感應的一維正演模擬結(jié)果。為了進行畸變分析，計算了考慮IP效應和純電磁感應的最大相對畸變值，見表2、表4；同時也計算了層狀介質(zhì)情況下考慮IP效應的電磁場響應與均勻半空間情況下純電磁感應的電磁響應的最大相對異常，見表3、表5。

圖4、圖6、圖8為極化層厚度為200 m時電場分量Ex，磁場分Hz響應曲線。對比純電磁感應的響應曲線，加入激電效應后的響應曲線形態(tài)有所改變，水平磁場Hy以及垂直磁場Hz幅值的極大值與極小值之間的差值變大，這是由于極化層的電阻率隨頻率的變化發(fā)生改變，而隨著厚度埋深比的減小以及電磁波傳播距離的增加，響應畸變出現(xiàn)在更低的頻率處。當厚度埋深比達1:10時，根據(jù)表2可知，由激電效應引起的響應畸變逐漸減小。如表3所示，根據(jù)最大相對異常值可知，Hz對淺部極化層更為靈敏，其次是電場分量Ex，而水平磁場Hy最弱。對比于其他分量，Ex對深部極化層依然有一定的反映能力。

圖5、圖7、圖9表示目標層厚度為80 m時電磁場各場值的理論響應曲線。同樣隨著埋深的加深，激電效應對電磁場分量的幅值的影響減弱。根據(jù)表4，當厚度埋深比為1:20時，考慮IP效應的響應與純電磁感應的響應之間的差別減小。