王言章，楊蒙蒙，符 磊，林 君

吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院/地球信息探測(cè)儀器教育部重點(diǎn)實(shí)驗(yàn)室，長(zhǎng)春 130026

0 引言

在瞬變電磁法中，大地中存在局部不均勻體時(shí)，野外實(shí)驗(yàn)采集到的數(shù)據(jù)中出現(xiàn)負(fù)響應(yīng)值，不能正確反映地下導(dǎo)體的情況，在后期數(shù)據(jù)處理時(shí)需要引起注意。雖然產(chǎn)生瞬變電磁負(fù)響應(yīng)的原因有若干種，但人們通常認(rèn)為是由激發(fā)極化效應(yīng)（IP）所引起的。

Lee［1］把均勻大地的導(dǎo)電率用 Cole－Cole模型描述，計(jì)算了大回線裝置下的瞬變電磁響應(yīng)，結(jié)果是激電效應(yīng)引起的瞬變曲線最初很快衰減，反向后瞬變衰減變慢；Smith等［2］發(fā)現(xiàn)，單個(gè)球體、半空間或?qū)訝畲蟮匾谥丿B回線TEM響應(yīng)中產(chǎn)生可觀測(cè)到的符號(hào)變化，其極化率一般都必須具有比電磁頻率范圍內(nèi)正常測(cè)得的大得多的值；Hohmann等［3］認(rèn)為，重疊回線或中心回線TEM測(cè)量中測(cè)得的正、負(fù)符號(hào)變化的畸變衰減曲線，是由近地表物質(zhì)如黏土等低極化材料激發(fā)極化引起的。

除了大多數(shù)研究者認(rèn)為T(mén)EM中的負(fù)響應(yīng)與IP效應(yīng)有關(guān)，人們對(duì)其他一些可能的原因也同時(shí)進(jìn)行了研究。在位移電流影響方面，Lee［4］從Morrison等人的有關(guān)公式出發(fā)，經(jīng)過(guò)傅里葉變換，得到重疊回線下無(wú)磁性均勻半空間瞬變電磁響應(yīng)的感應(yīng)電動(dòng)勢(shì)，結(jié)果發(fā)現(xiàn)，位移電流只對(duì)反映高頻響應(yīng)的早期道有影響，而對(duì)反映低頻響應(yīng)的晚期道的影響可以忽略。在集流效應(yīng)方面，McNeil［5］等在印度Gani村用SIROTEM－Ⅱ儀器，在一傾斜導(dǎo)電炭質(zhì)頁(yè)巖上方，垂直巖層走向，在100m×100m的重疊回線的剖面圖中，觀測(cè)到晚延時(shí)刻顯示出了一個(gè)負(fù)響應(yīng)。

由于不均勻體結(jié)構(gòu)的復(fù)雜性，很難準(zhǔn)確地通過(guò)建模來(lái)計(jì)算出礦體的理論衰減曲線，而多匝密繞的閉合線圈（簡(jiǎn)稱異常線圈）電磁信號(hào)衰減規(guī)律與不均勻體的電磁信號(hào)衰減規(guī)律相同，可采用異常線圈來(lái)模擬地下有限導(dǎo)體［6］。筆者采用大定源回線裝置，用不同的異常線圈模型對(duì)回線下方存在的淺層不均勻體進(jìn)行模擬；擬從TEM的理論知識(shí)出發(fā)，通過(guò)matlab仿真，研究發(fā)射線圈、接收線圈采用中心回線方式，當(dāng)發(fā)射回線中的梯形波電流關(guān)斷之后，在電阻率不同的均勻大地中存在異常線圈模型，及位置不同時(shí)產(chǎn)生的響應(yīng)曲線對(duì)接收線圈衰減曲線的影響情況。

1 接收線圈響應(yīng)理論

根據(jù)瞬變電磁法基本原理，當(dāng)發(fā)射回線中的電流突然下降為0A時(shí)，會(huì)在其周圍產(chǎn)生急劇變化的磁場(chǎng)。對(duì)于發(fā)射回線下方存在的異常體，變化的磁場(chǎng)會(huì)在其內(nèi)部激發(fā)產(chǎn)生感應(yīng)電流，形成大致按指數(shù)規(guī)律隨時(shí)間衰減的瞬變磁場(chǎng)。圖1為局部不均勻體位于發(fā)射回線下方的示意圖。

圖1 局部不均勻體位于發(fā)射回線下方示意圖Fig.1 Testing diagram of the local heterogeneity located below the transmit loop

實(shí)際發(fā)射電流波形不可能為理想階躍，大發(fā)射回線采用的是梯形波發(fā)射電流，波形如圖2所示。

發(fā)射回線中電流關(guān)斷之后，在異常線圈中產(chǎn)生的感應(yīng)電動(dòng)勢(shì)為

圖2 梯形波發(fā)射電流Fig.2 Emission current of trapezoidal wave

其中：MTL為T(mén)、A之間的互感系數(shù)；IT（t）為發(fā)射電流。

異常線圈感應(yīng)電流IAL（t）在晚期之后已處于穩(wěn)定狀態(tài)，并按照指數(shù)規(guī)律衰減，可以用等效電感LAL和電阻RAL近似［7］，所以

式（1）、（2）聯(lián)立可求出IAL（t），則接收到的異常線圈感應(yīng)電動(dòng)勢(shì)為

其中，

式中：MRL為R、A之間的互感系數(shù)；I為發(fā)射電流幅值表示0時(shí)刻選在關(guān)斷結(jié)束時(shí)刻；t1為電流上升時(shí)間；t2為電流下降時(shí)間；t3為發(fā)射電流寬度；τ為異常線圈時(shí)間常數(shù)，由異常線圈的幾何形狀、大小和電性參數(shù)所決定。MRL和MTL與發(fā)射、異常及接收線圈的相對(duì)位置、異常線圈的幾何形狀等因素有關(guān)［8－9］。

當(dāng)發(fā)射電流突然斷開(kāi)時(shí)，發(fā)射線圈周圍的磁感線穿過(guò)異常線圈激發(fā)產(chǎn)生二次場(chǎng)，隨著A的幾何形態(tài)以及T、A、R三者之間相對(duì)位置的不同，產(chǎn)生的二次場(chǎng)方向不同，使得A與R之間的互感有正有負(fù)。根據(jù)TEM理論，這會(huì)引起電磁場(chǎng)矢量的疊加方向發(fā)生變化，結(jié)合式（3）、（4），可以得出結(jié)論：接收到的異常線圈感應(yīng)電動(dòng)勢(shì)有正有負(fù)。因此，異常線圈的存在，使得接收線圈中的綜合響應(yīng)可能出現(xiàn)負(fù)響應(yīng)。

對(duì)于梯形波回線源發(fā)射，發(fā)射線圈、接收線圈采用中心回線方式。當(dāng)發(fā)射電流關(guān)斷后，在準(zhǔn)靜態(tài)場(chǎng)中，均勻大地的電磁響應(yīng)為

其中：N為接收線圈匝數(shù)；S為接收線圈面積；B－S（t）為負(fù)階躍電流產(chǎn)生的二次場(chǎng)磁感應(yīng)強(qiáng)度。為了便于matlab仿真計(jì)算，將矩形大發(fā)射回線等效成面積相等的圓形發(fā)射回線［10－12］。根據(jù)納比吉安［13］的做法，單位電流切斷后，回線中心的磁場(chǎng)強(qiáng)度表達(dá)式為

其中：

式中：μ0為真空磁導(dǎo)率，值為4π×10－7H/m；a為圓形發(fā)射回線半徑；ρ為大地電阻率。發(fā)射電流關(guān)斷后，理想接收的一次場(chǎng)響應(yīng)近似為0［14］，接收線圈中的理想電磁響應(yīng)為

在研究接收線圈中的綜合響應(yīng)曲線時(shí)，采用式（8）進(jìn)行仿真計(jì)算。

2 異常線圈模型確定

異常線圈的時(shí)間常數(shù)τ值，與異常線圈本身的材料和電氣特性有關(guān)，即與線圈的自感和電阻有關(guān)。其表達(dá)式為

式中：n為異常線圈匝數(shù)；ρc為銅導(dǎo)線電阻率，值為1.75×10－8Ω·m；rc為導(dǎo)線半徑；lc為導(dǎo)線長(zhǎng)度；sc為導(dǎo)線橫截面積。

從式（9）可以看出：τ值隨n的增加而線性增大；隨rAL的增加而增加；隨rc的增加而增加。

筆者采用矩形發(fā)射回線，發(fā)射邊長(zhǎng)為150m×150m，要研究的是在發(fā)射回線下方或者是框邊下方，反向磁通對(duì)接收線圈綜合響應(yīng)的影響情況。所以，選中的異常線圈的半徑比較大，產(chǎn)生的響應(yīng)值比較大，異常線圈的直徑不能超過(guò)發(fā)射回線邊長(zhǎng)的一半，即rAL不能超過(guò)150/2/2＝37.5m。

由于有經(jīng)濟(jì)價(jià)值的礦體的時(shí)間常數(shù)一般大于2 ms，筆者對(duì)于時(shí)間常數(shù)分別為1、3、5ms的異常線圈進(jìn)行研究，根據(jù)式（9），其對(duì)應(yīng)的模型參數(shù)如表1所示。

表1 不同時(shí)間常數(shù)對(duì)應(yīng)的線圈模型Table 1 Coil model corresponding to the different time constant

根據(jù)表1，用異常線圈對(duì)發(fā)射回線下方存在的局部不均勻體進(jìn)行模擬，仿真計(jì)算不均勻體的電磁響應(yīng)衰減曲線。結(jié)合式（9），選用導(dǎo)線橫截面積為8～10mm2。

為了能夠明顯觀測(cè)到淺層不均勻體對(duì)衰減曲線的影響情況，在異常線圈中心點(diǎn)P坐標(biāo)（x，y，h）的選擇方面，應(yīng)考慮到由于A、R之間相對(duì)位置的變化，引起異常線圈的負(fù)響應(yīng)值幅度的變化。結(jié)合上述異常線圈半徑選擇方面的情況，對(duì)于異常線圈位于發(fā)射回線下方時(shí)，設(shè)x＝40m，y＝0；異常線圈位于發(fā)射回線邊界下方時(shí)，x變化范圍為45～105m，y＝0。

3 均勻大地響應(yīng)

發(fā)射回線采用圖2所示的梯形波，其中發(fā)射電流幅度為10A，上升時(shí)間為1ms，平頂時(shí)間為19 ms，關(guān)斷時(shí)間為0.176ms。接收線圈的半徑為0.3 m，等效面積為0.28m2。

發(fā)射線圈、接收線圈采用中心回線方式。根據(jù)式（5）—（7），不同電阻率的均勻大地產(chǎn)生的響應(yīng)曲線如圖3所示。

由圖3可知，發(fā)射梯形波電流關(guān)斷之后，均勻大地的電阻率越高，大地電磁響應(yīng)幅值越小。無(wú)異常線圈存在時(shí)的均勻大地，產(chǎn)生的響應(yīng)值全為正響應(yīng)。

圖3 不同電阻率的均勻大地響應(yīng)Fig.3 Response of the homogeneous earth when the resistivity is different

4 接收線圈綜合響應(yīng)

不同ρ的均勻大地產(chǎn)生響應(yīng)不同，以及τ、r、h和P點(diǎn)位置中x不同的異常線圈，在接收線圈中產(chǎn)生的響應(yīng)值也不同，根據(jù)式（8）可知，接收到的綜合響應(yīng)衰減情況也不相同。

4.1 異常線圈位于不同電阻率的均勻大地中時(shí)，接收線圈的綜合響應(yīng)

在電阻率ρ不同的均勻大地，當(dāng)rAL＝30m，時(shí)間常數(shù)τ＝3ms的異常線圈模型埋深h＝5m時(shí)，接收線圈的綜合響應(yīng)曲線如圖4所示。

由圖4可知，在電阻率不同的均勻大地中，異常線圈的存在使接收線圈的綜合響應(yīng)中出現(xiàn)了負(fù)值情況，電阻率越大，負(fù)響應(yīng)值出現(xiàn)的越早，早期的衰減曲線過(guò)0情況越明顯。在接下來(lái)的討論中，為了能夠明顯看出因異常線圈存在而導(dǎo)致的衰減曲線與真實(shí)情況的偏離程度，選擇電阻率ρ＝3 000Ω·m的均勻大地進(jìn)行研究。

4.2 時(shí)間常數(shù)τ不同時(shí)，接收線圈的綜合響應(yīng)

當(dāng)rAL半徑＝30m，τ分別為1、3、5ms的異常線圈h＝5m時(shí)，接收線圈的綜合響應(yīng)曲線如圖5所示。

由圖5可知，異常線圈的存在使接收線圈的綜合響應(yīng)值減?。涸陔娏麝P(guān)斷之后2ms范圍內(nèi)，時(shí)間常數(shù)越小的異常線圈，對(duì)綜合響應(yīng)的影響越明顯，導(dǎo)致異常線圈衰減曲線中的過(guò)0幅度越明顯；在2ms之后，時(shí)間常數(shù)越大的異常線圈，對(duì)綜合響應(yīng)的負(fù)值影響越明顯；時(shí)間常數(shù)越小的異常線圈衰減為0的時(shí)間越快。

圖4 ρ不同時(shí)均勻大地接收線圈的綜合響應(yīng)Fig.4 Integrated response of the receiver coil whenρof the homogeneous earth is different

圖5 異常線圈τ時(shí)間常數(shù)不同時(shí)接收線圈的綜合響應(yīng)Fig.5 Integrated response of the receiver coil whenτof the anomaly loop is different

4.3 異常線圈半徑rAL不同時(shí)接收線圈的綜合響應(yīng)

異常線圈rAL半徑分別為35、30、25、15m，τ＝3ms，h＝5m時(shí)，接收線圈的綜合響應(yīng)曲線如圖6所示。

由圖6可知：對(duì)于同一時(shí)間常數(shù)的異常線圈而言，在同一位置，半徑越大，產(chǎn)生的負(fù)值響應(yīng)越明顯，導(dǎo)致接收線圈的響應(yīng)衰減曲線中過(guò)0情況越明顯；隨著半徑的減小，異常線圈產(chǎn)生的響應(yīng)值變小，對(duì)接收線圈綜合響應(yīng)的情況可忽略不計(jì)。無(wú)論異常線圈半徑大小如何變化，綜合響應(yīng)的曲線大致在同一時(shí)間點(diǎn)衰減為0。

圖6 rAL不同時(shí)接收線圈的綜合響應(yīng)Fig.6 Integrated response of the receiver coil when rALof the anomaly loop is different

4.4 異常線圈深度h不同時(shí)接收線圈的綜合響應(yīng)

rAL＝30m，τ＝3ms的異常線圈，h分別為5、10、20、30、40m時(shí)，接收線圈的綜合響應(yīng)曲線如圖7所示。

圖7 異常線圈h不同時(shí)接收線圈的綜合響應(yīng)Fig.7 Integrated response of the receiver coil when h of the anomaly loop is different

由圖7可知：對(duì)于同一時(shí)間常數(shù)的異常線圈，當(dāng)深度不同時(shí)，異常線圈產(chǎn)生的感應(yīng)電動(dòng)勢(shì)強(qiáng)度和方向不同。埋深位置越淺，負(fù)值越小，減小了接收線圈中的綜合響應(yīng)，使衰減曲線出現(xiàn)負(fù)值；隨著埋深深度的增加，感應(yīng)電動(dòng)勢(shì)變?yōu)檎较?，增加了接收線圈中的綜合響應(yīng)，使衰減曲線全為正值；當(dāng)增加到一定深度時(shí)，由于異常線圈、接收線圈之間的距離較大，二者之間的互感越來(lái)越小，此時(shí)，異常線圈的存在對(duì)接收線圈的影響可忽略不計(jì)。從圖7中可以看出，在0～10m的淺層范圍內(nèi)，異常線圈的存在，可導(dǎo)致接收線圈的綜合響應(yīng)中出現(xiàn)過(guò)0情況。

4.5 異常線圈位置x不同時(shí)接收線圈的綜合響應(yīng)

對(duì)于rAL＝30m，τ＝3ms的異常線圈模型，當(dāng)h＝5m，發(fā)射回線下方的位置x分別為0、20、40、70、90m時(shí)，接收線圈的綜合響應(yīng)曲線如圖8所示。

圖8 異常線圈位置x不同時(shí)接收線圈的綜合響應(yīng)Fig.8 Integrated response of the receiver coil when the position（x）of the anomaly loop is different

異常線圈在發(fā)射回線下方，沿與x軸平行的測(cè)線移動(dòng)。由圖8可知：當(dāng)移動(dòng)到發(fā)射回線正下方（x＝0m）時(shí)，異常線圈產(chǎn)生的響應(yīng)值為正，使接收線圈中的綜合響應(yīng)增加，衰減為0的速度變慢；當(dāng)異常線圈逐漸遠(yuǎn)離發(fā)射回線正下方時(shí)，異常線圈的響應(yīng)由正值變?yōu)樨?fù)值，使接收線圈中的綜合響應(yīng)曲線出現(xiàn)過(guò)0情況；當(dāng)異常線圈移動(dòng)到發(fā)射回線框邊下方，并且越來(lái)越遠(yuǎn)離發(fā)射回線時(shí)，異常線圈產(chǎn)生的響應(yīng)正值比較小，可忽略不計(jì)，此時(shí)接收線圈中的綜合響應(yīng)曲線在雙對(duì)數(shù)坐標(biāo)中，大致按線性規(guī)律衰減。

5 結(jié)論與建議

1）在電阻率高的均勻大地中，由于均勻大地產(chǎn)生的響應(yīng)值比較小，不均勻的存在對(duì)接收線圈衰減曲線的影響比較明顯。

2）對(duì)于時(shí)間常數(shù)不同的異常線圈模型來(lái)說(shuō)：早期（如2ms之前），時(shí)間常數(shù)越小，衰減曲線過(guò)0情況越明顯；晚期（如2ms之后），時(shí)間常數(shù)越大，衰減曲線過(guò)0情況越明顯。對(duì)時(shí)間常數(shù)相同的異常線圈模型來(lái)說(shuō)，當(dāng)位于發(fā)射回線下方的位置坐標(biāo)x＝40 m時(shí)，半徑越大、深度越淺，產(chǎn)生的負(fù)值響應(yīng)幅度越大，使接收線圈的綜合響應(yīng)曲線的過(guò)0程度越明顯。

3）對(duì)于在發(fā)射回線下方位置變化的異常線圈模型來(lái)說(shuō)，產(chǎn)生的響應(yīng)值大小和方向發(fā)生變化。當(dāng)位于發(fā)射回線中心點(diǎn)下方位置時(shí)，響應(yīng)值為正，增加了接收線圈中的綜合響應(yīng)；當(dāng)沿測(cè)線向發(fā)射回線框邊移動(dòng)時(shí)，響應(yīng)值變?yōu)樨?fù)值；當(dāng)靠近發(fā)射回線框邊時(shí)，產(chǎn)生的響應(yīng)值為正，與均勻大地的響應(yīng)值相比，幅值較小，可忽略不計(jì)。

所以，均勻大地中存在的淺層不均勻體產(chǎn)生的響應(yīng)值，會(huì)導(dǎo)致接收線圈中的衰減曲線出現(xiàn)負(fù)值，由于這種負(fù)值不能正確的反映地下目標(biāo)體的情況，在后期進(jìn)行數(shù)據(jù)分析時(shí)，應(yīng)對(duì)其進(jìn)行處理。

（References）：

［1］Lee T.The Cole－Cole Model in Time Domain Induced Polarization［J］.Geophysics，1981，46（6）：932－933.

［2］Smith R S，West G F.Inductive Interaction Between Polarizable Conductors：An Explanation of a Negative Coincident－Loop Transient Electromagnetic Response［J］.Geophysics，1988，53：677－690.

［3］Hohmann G，Newman G.Transient Electromagnetic Response of Surficial， Polarizable Patches［J］.Geophysics，1990，55：1098－1100.

［4］Lee T.The Effect of Displacement Currents on Time Domain Electromagnetic Fields［J］.Bull Aust Soc Explor Geophys，1981，12（3）：34－36.

［5］McNeill J D，Edwards R N，Levy M.Approximate Calculations of the Transient Electromagmetic Response from Buried Conducts in a Conductive Half－Space［J］.Geophysies，1984，49（7）：918－924.

［6］嵇艷鞠，李肅義，于生寶，等.基于異常線圈的時(shí)間域AEM系統(tǒng)測(cè)試和標(biāo)定方法研究［J］.地球物理學(xué)報(bào)，2011，54（10）：2690－2697.Ji Yanju，Li Suyi，Yu Shengbao，et al.A Study on Time－Domain AEM Testing and Calibration Method Based on Anomaly Loop［J］.Chinese Journal of Geophysics，2011，54（10）：2690－2697.

［7］牛之璉.時(shí)間域電磁法原理［M］.長(zhǎng)沙：中南大學(xué)出版社，2007.Niu Zhilian.The Theory of Time－Domain Electro－magnetic Methods［M］.Changsha：Central South University of Technology Press，2007.

［8］陳曙東.直升機(jī)瞬變電磁系統(tǒng)響應(yīng)計(jì)算與地回線標(biāo)定［D］.長(zhǎng)春：吉林大學(xué)，2012.Chen Shudong.Response Calculation and Groud－Loop Calibration of Helicopter Transient Electromagnetic System［D］.Changchun：Jilin University，2012.

［9］嵇艷鞠，欒卉，李肅義，等.全波形時(shí)間域航空電磁探測(cè)分辨率［J］.吉林大學(xué)學(xué)報(bào)：地球科學(xué)版，2011，41（3）：885－891.Ji Yanju，Luan Hui，Li Suyi，et al.Resolution of Full－Waveform Airborne TEM［J］.Journal of Jilin University：Earth Science Edition，2011，41（3）：885－891.

［10］許洋鋮.全波形時(shí)間域電磁接收系統(tǒng)及分辨力研究［D］.長(zhǎng)春：吉林大學(xué)，2011.Xu Yangcheng.Research of Resolving Capacity and Development of Full－Waveform Electromagnetic Receiver in Time Domain［D］.Changchun：Jilin University，2011.

［11］劉桂芬.回線源層狀大地航空瞬變電磁場(chǎng)的理論計(jì)算［D］.長(zhǎng)春：吉林大學(xué)，2004.Liu Guifen.Theoretical Calculation of Airborne Transient Electromagnetic Field for Loop Source on Layered Earth［D］.Changchun：Jilin University，2004.

［12］嵇艷鞠，林君，程德福，等.ATEM－Ⅱ瞬變電磁儀數(shù)據(jù)處理軟件的研制與應(yīng)用［J］.吉林大學(xué)學(xué)報(bào)：地球科學(xué)版，2003，33（2）：242－245.Ji Yanju，Lin Jun，Cheng Defu，et al.Development and Application of Data Processing Software of ATEM－ⅡTransient Electromagnetic Instrument［J］.Journal of Jilin University：Earth Science Edition，2003，33（2）：242－245.

［13］米薩克N納比吉安.勘查地球物理電磁法［M］.趙經(jīng)祥，譯.北京：地質(zhì)出版社，1992.Nabighian M N.Electromagnetic Methods in Applied Geophysics［M］.Translated by Zhao Jingxiang.Beijing：Geological Publishing House，1992.

［14］Smith R S，Balch S J.Robust Estimation of the Band－Limited Inductive－Limit Response from Impulse－Response TEM Measurements Taken During the Transmitter Switch－Off and the Transmitter Off－Time：Theory and an Example from Voisey’s Bay，Labrador，Canada［J］.Geophysics，2000，65：476－481.