智慶全，武軍杰，王興春，楊 毅，張 杰，鄧曉紅

(1.中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所，河北 廊坊 065000；2.自然資源部地球物理電磁法探測(cè)技術(shù)重點(diǎn)實(shí)驗(yàn)室，河北 廊坊 065000；3.長(zhǎng)安大學(xué) 地質(zhì)工程與測(cè)繪學(xué)院，陜西 西安 710054)

地-井瞬變電磁法(Transient Electromagnetic Method，TEM)是一種將發(fā)射裝置布設(shè)在井孔上方或附近平面上，在鉆孔中逐點(diǎn)測(cè)量的電磁探測(cè)方法。相比于傳統(tǒng)的地面瞬變電磁法，其接收系統(tǒng)置于鉆孔中，大地的濾波作用能夠有效壓制高頻干擾，提高觀(guān)測(cè)數(shù)據(jù)的信噪比；接收裝置更加靠近地質(zhì)體，能夠近距離獲取異常響應(yīng)，對(duì)薄層、井旁地質(zhì)體等探測(cè)目標(biāo)具有更強(qiáng)的分辨能力。多年來(lái)在固體礦產(chǎn)資源勘查[1-3]、煤田采空區(qū)探測(cè)[4-5]巷道超前預(yù)報(bào)[6-7]等方面得到了廣泛的應(yīng)用并取得了良好的勘查效果。

作為一種地下電磁探測(cè)方法，地-井瞬變電磁的響應(yīng)特征和解釋理論相比地面瞬變電磁法更為復(fù)雜。在地-井瞬變電磁響應(yīng)的正演模擬方面，J.R.Wait 等[8-10]較早地開(kāi)展了理論分析工作，推導(dǎo)了均勻半空間情況下地-井瞬變電磁響應(yīng)的積分表達(dá)式和特殊情況下的解析表達(dá)式，奠定了響應(yīng)特征分析的理論基礎(chǔ)。之后，學(xué)者們采用數(shù)字濾波方法實(shí)現(xiàn)了地下電磁場(chǎng)的一維正演[11-13]，通過(guò)有限元[14-15]、有限差分[16-18]、有限體積[19]等方法實(shí)現(xiàn)了復(fù)雜模型的全空間電磁場(chǎng)三維模擬。目前，地-井瞬變電磁的響應(yīng)模擬理論已較為成熟，關(guān)于其響應(yīng)特征和規(guī)律的總結(jié)也漸趨完善。

在地-井瞬變電磁的數(shù)據(jù)解釋方面，M.N.Nabighian[20]首先分析了瞬變電磁場(chǎng)在地下的擴(kuò)散特征，指出地表電磁場(chǎng)可以通過(guò)地下的一個(gè)電流環(huán)進(jìn)行等效，電流環(huán)在大地-空氣界面的作用下向深部和外部擴(kuò)散，即著名的“煙圈”理論，從物理上闡明了瞬變電磁場(chǎng)在地下的擴(kuò)散規(guī)律，同時(shí)也提供了一種較為快捷的數(shù)據(jù)解釋方法。但該理論的應(yīng)用范圍局限于地表以上的電磁場(chǎng)，對(duì)于地-井瞬變電磁數(shù)據(jù)并不適用。A.Ziolkowski[21]、羅維斌[22]、湯井田[23]等通過(guò)對(duì)海洋瞬變電磁法的模擬分析，發(fā)現(xiàn)不同偏移距的電場(chǎng)響應(yīng)在時(shí)間上存在一個(gè)“峰值”，進(jìn)而提出了瞬變沖激時(shí)刻[21-23]的概念，并指出沖激時(shí)刻與介質(zhì)電阻率相關(guān)，通過(guò)沖激時(shí)刻可計(jì)算出介質(zhì)電阻率值。部分學(xué)者總結(jié)提出了一套實(shí)用的地-井TEM 數(shù)據(jù)采集技術(shù)，以及從定性分析到半定量的解釋技術(shù)，采用矢量交會(huì)技術(shù)實(shí)現(xiàn)了地-井TEM 井旁異常的快速定位[24]，采用電流環(huán)等效實(shí)現(xiàn)了純異常近似三維反演[25]，并初步研究了地-井TEM 數(shù)據(jù)的視電阻率定義方法[26]。

綜合以上研究可以發(fā)現(xiàn)，當(dāng)前已有的地-井TEM解釋方法中，基本上以定性分析和近似反演等半定量解釋方法為主，以視電阻率定義、一維反演為代表的快速定量解釋技術(shù)研究較少且尚不成熟，而解釋效果最佳的三維反演技術(shù)時(shí)間成本過(guò)高，難以實(shí)現(xiàn)快速定量解釋。最終導(dǎo)致地-井瞬變電磁測(cè)量成果往往僅能以響應(yīng)曲線(xiàn)、交會(huì)矢量、電流環(huán)位置等形式呈現(xiàn)，不能直觀(guān)地反映地下介質(zhì)的電阻率分布。鑒于此，筆者基于前人關(guān)于瞬變沖激時(shí)刻的研究，提出一種地-井TEM 數(shù)據(jù)的快速定量解釋方法，利用地下瞬變電磁場(chǎng)的沖激時(shí)刻定義大地視電阻率，通過(guò)研究瞬變電磁場(chǎng)擴(kuò)散速度與深度的關(guān)系，導(dǎo)出地層真電阻率的求取方法。該方法理論基礎(chǔ)可靠，計(jì)算速度快、對(duì)大地真實(shí)電阻率恢復(fù)程度高，具有較高的實(shí)用價(jià)值。

1 地-井瞬變電磁響應(yīng)的沖激時(shí)刻

關(guān)于瞬變電磁煙圈效應(yīng)和沖激時(shí)刻的已有研究結(jié)果表明，當(dāng)發(fā)射源和接收點(diǎn)的位置存在偏移量時(shí)，電場(chǎng)分量在早期較小，隨時(shí)間推移逐漸增大并達(dá)到峰值，然后開(kāi)始衰減，晚期呈現(xiàn)單調(diào)的衰減特征，電場(chǎng)達(dá)到峰值的時(shí)刻稱(chēng)為瞬變沖激時(shí)刻。沖激時(shí)刻的早晚取決于瞬變電磁場(chǎng)在空間上的擴(kuò)散速度，大地電阻率越高，電磁場(chǎng)擴(kuò)散速度越快，沖激時(shí)刻越早。對(duì)于地-井瞬變電磁法而言，一般觀(guān)測(cè)地下的磁場(chǎng)分量(dB/dt)，因此，首先需要對(duì)磁場(chǎng)響應(yīng)的沖激時(shí)刻特征進(jìn)行分析。

J.R.Wait 等[9]首先推導(dǎo)了均勻半空間的地-井瞬變電磁響應(yīng)表達(dá)式。對(duì)于電導(dǎo)率為σ的均勻大地，半徑為a的發(fā)射回線(xiàn)位于地表，其中供有大小為I0的電流，接收點(diǎn)P位于地下，如圖1 建立柱坐標(biāo)系，取向上為Z軸正方向。在t=0 時(shí)刻瞬間撤掉發(fā)射回線(xiàn)中的電流后，P點(diǎn)磁場(chǎng)響應(yīng)表達(dá)式為：

圖1 均勻半空間模型和地-井TEM 測(cè)量Fig.1 Diagram of the half-space model and surface-borehole TEM

式中：Br和Bz分別為磁場(chǎng)的徑向和垂向分量；μ為大地磁導(dǎo)率；h為接收點(diǎn)深度；r為接收點(diǎn)與坐標(biāo)原點(diǎn)的水平偏移距離；A為發(fā)射回線(xiàn)的相對(duì)半徑，A=a/h；m(t)和q(t)分別為包含有第一類(lèi)貝塞爾函數(shù)J0和J1的無(wú)窮積分；D=r/h；T=t/(σμh2)；Z=z0/h；z0為發(fā)射回線(xiàn)在地表以上的高度；f(x)=2J1(x)/x。式中的無(wú)窮積分可通過(guò)快速漢克爾變換法[27]求取。

為了直觀(guān)地分析地-井瞬變電磁響應(yīng)的特征和規(guī)律，分別計(jì)算了均勻半空間模型在不同深度接收點(diǎn)和不同電阻率、相同接收點(diǎn)2 種情況下的dBz/dt響應(yīng)，結(jié)果如圖2 所示。在圖2a 中，均勻半空間電阻率ρ=100 Ω·m，發(fā)射回線(xiàn)位于地表z=0 m 位置，半徑為100 m，發(fā)射電流1 A，接收點(diǎn)位于地下(-10 m,10 m,-h)位置，深度h的范圍設(shè)置為0～600 m，間隔100 m。圖2b 中，接收點(diǎn)深度固定為300 m，模型電阻率分別設(shè)置為10、30、100、300、1 000 Ω·m。由圖2 可知，瞬變電磁場(chǎng)在地下傳播時(shí)，dBz/dt響應(yīng)在時(shí)間上存在一個(gè)峰值，其對(duì)應(yīng)時(shí)刻即為dBz/dt響應(yīng)的沖激時(shí)刻，且沖激時(shí)刻的早晚取決于大地電阻率和接收點(diǎn)深度：深度越大、大地電阻率越低，峰值時(shí)刻出現(xiàn)越晚。

圖2 均勻半空間的地-井TEM 響應(yīng)Fig.2 Surface-borehole TEM responses of uniform half-space

根據(jù)沖激時(shí)刻的定義，分別采用電阻率為10、100、1 000 Ω·m 的均勻半空間模型，繪制在不同井位、不同深度上dBz/dt響應(yīng)的沖激時(shí)刻tp，如圖3 所示，并進(jìn)行多項(xiàng)式曲線(xiàn)擬合。其中設(shè)計(jì)井位Hole 1 的水平坐標(biāo)為(-10 m,10 m)，Hole 2 的水平坐標(biāo)為(80 m，50 m)。由圖3 可知，沖激時(shí)刻隨大地電阻率的增加而變?cè)?，二者之間存在反比例關(guān)系；沖激時(shí)刻隨深度的增大而變晚，其增加速率高于線(xiàn)性函數(shù)。結(jié)合不同電阻率大地的沖激時(shí)刻曲線(xiàn)擬合結(jié)果，進(jìn)行分析可以發(fā)現(xiàn)，沖激時(shí)刻隨深度變化的函數(shù)呈現(xiàn)二次型特征，近似有函數(shù)關(guān)系：

圖3 沖激時(shí)刻隨深度變化的函數(shù)圖Fig.3 Function diagram of impulse time varying with depth

式中：tp為沖激時(shí)刻，ms；ρ為大地電阻率，Ω·m；K為一個(gè)常量系數(shù)，ms·Ω·m-1，對(duì)于Hole1有K≈ 1.23×10-4ms·Ω·m-1，對(duì)于Hole2 有K≈ 1.30×10-4ms·Ω·m-1，二者之間的差距約為5.38%。

2 視電阻率定義和大地電阻率求取方法

關(guān)于準(zhǔn)靜態(tài)條件下電磁場(chǎng)擴(kuò)散規(guī)律的研究表明，電磁場(chǎng)在均勻大地中的極大值位置是僅僅關(guān)于大地電導(dǎo)率和時(shí)間的函數(shù)，與發(fā)射回線(xiàn)幾何形狀、發(fā)射電流等無(wú)關(guān)[20]。進(jìn)一步研究表明，地下瞬變電磁響應(yīng)的dBz/dt分量的極大值深度[28]可以近似表示為：

式中：b為一個(gè)常數(shù)，在1 階至3 階級(jí)數(shù)近似條件下，分別有b(1)=8，b(2)=12，b(3)≈ 9.112。簡(jiǎn)單整理可以發(fā)現(xiàn)，該式的計(jì)算結(jié)果與式(3)十分接近，表明二次函數(shù)完全可以作為沖激時(shí)刻隨深度變化函數(shù)的一個(gè)較好的近似。

在已知接收點(diǎn)深度h和沖激時(shí)刻tp的情況下，根據(jù)式(4)即可求得大地電阻率值，即有基于瞬變沖激時(shí)刻的視電阻率定義：

根據(jù)此式即可直接計(jì)算不同深度接收點(diǎn)的視電阻率?？紤]到級(jí)數(shù)近似所固有的誤差，本文在實(shí)際計(jì)算時(shí)直接采用式(2)，基于反函數(shù)理論[29]，通過(guò)迭代來(lái)獲取更為精確的視電阻率值。

以上視電阻率定義方法直接基于瞬變電磁場(chǎng)在大地中的沖激時(shí)刻，主要反映接收點(diǎn)所在深度以淺大地介質(zhì)的總體電性特征，具有較強(qiáng)的“體積效應(yīng)”。為更加精細(xì)地恢復(fù)大地電阻率的分布，本文基于電磁場(chǎng)擴(kuò)散速度，給出一種改進(jìn)的大地電阻率計(jì)算方法。假設(shè)大地由多個(gè)電阻率均勻的薄層組成，根據(jù)式(4)有均勻介質(zhì)中瞬變電磁場(chǎng)的擴(kuò)散速度：

對(duì)于第i個(gè)薄層，其頂、底界面深度分別為hi和hi+1，沖激時(shí)刻分別為ti和ti+1，對(duì)應(yīng)的視電阻率分別為和，則有：

又根據(jù)速度的物理定義有：

聯(lián)立式(6)和式(8)，將式(7)代入并整理可得第i薄層的電阻率：

式(9)即為各層大地電阻率的計(jì)算公式，該式不包含系數(shù)項(xiàng)b，與漸近級(jí)數(shù)的階數(shù)無(wú)關(guān)。這種電阻率計(jì)算方法實(shí)際上基于瞬變電磁場(chǎng)在大地中的擴(kuò)散速度，相比基于沖激時(shí)刻的視電阻率而言，通過(guò)式(9)計(jì)算得到的電阻率受體積效應(yīng)影響更小，更接近真實(shí)大地電阻率。由于無(wú)需進(jìn)行復(fù)雜模型的響應(yīng)擬合計(jì)算，該算法具有計(jì)算簡(jiǎn)單、時(shí)間成本低、無(wú)多解性的特點(diǎn)，適用于地-井TEM 數(shù)據(jù)的快速定量解釋。

3 模型算例

為檢驗(yàn)上述基于沖激時(shí)刻的地-井TEM 數(shù)據(jù)快速解釋方法的正確性和有效性，選擇了一維模型系列進(jìn)行數(shù)值模擬計(jì)算。一維模型主要選擇D、G、H、K型4 種典型地電斷面，各層的電阻率ρ和厚度H見(jiàn)表1。觀(guān)測(cè)參數(shù)參照實(shí)際工作中常用的工作裝置進(jìn)行設(shè)置，發(fā)射回線(xiàn)位于地表z=0 m 位置，回線(xiàn)邊長(zhǎng)為100 m，發(fā)射電流1 A，井孔位于(-10 m,10 m)位置，鉆孔為直孔，接收點(diǎn)深度范圍為0～600 m，間隔10 m。

表1 一維層狀模型參數(shù)Table 1 Parameters of 1D models

圖4a 為均勻半空間情況下基于沖激時(shí)刻定義的視電阻率ρs和基于擴(kuò)散速度恢復(fù)的大地電阻率值ρd。由圖4a 可知，在不同深度的接收點(diǎn)上，基于沖激時(shí)刻定義的視電阻率和基于擴(kuò)散速度恢復(fù)的大地電阻率基本一致，且與設(shè)計(jì)模型吻合較好，證明了所提出算法的正確性。圖4b 為計(jì)算電阻率數(shù)據(jù)與設(shè)計(jì)模型的相對(duì)誤差，ρs和ρd與設(shè)計(jì)半空間電阻率的最大相對(duì)誤差分別為0.28%和2.06%。這是由于式(9)中包含了對(duì)大地電導(dǎo)率的差分計(jì)算，使得ρd在數(shù)值穩(wěn)定性上弱于基于沖激時(shí)刻的視電阻率ρs，導(dǎo)致其相對(duì)誤差較大。

圖4 均勻半空間計(jì)算電阻率隨深度的變化Fig.4 Variation of calculated resistivity with depth in uniform half-space

圖5 為兩層D 型模型和G 型模型的ρs和ρd計(jì)算結(jié)果，并與設(shè)計(jì)模型電阻率(ρtrue)進(jìn)行了對(duì)比。由圖5 可知，ρs和ρd均能較好地反映出兩層模型的電阻率變化趨勢(shì)，但ρs曲線(xiàn)較為寬緩，ρd曲線(xiàn)變化較為迅速且存在振蕩現(xiàn)象。隨深度增大，ρd總能更早地體現(xiàn)出大地的電阻率變化，且與真實(shí)大地電阻率更為接近。這表明，雖然ρd的數(shù)值穩(wěn)定性較弱，但其對(duì)大地電阻率變化的靈敏度高于ρs，對(duì)大地的電性界面位置反映也更為精細(xì)。

圖5 兩層模型計(jì)算電阻率隨深度的變化Fig.5 Variation of calculated resistivity with depth in 2-layer models

注意的一點(diǎn)是，雖然ρd對(duì)大地電阻率變化更為靈敏，但其電阻率曲線(xiàn)存在與傳統(tǒng)視電阻率定義類(lèi)似的“overshoot”和“undershoot”現(xiàn)象，即在低阻地層中會(huì)伴隨產(chǎn)生一個(gè)略弱的“假”高阻層，反之在高阻地層中也同樣會(huì)伴隨產(chǎn)生一個(gè)略弱的“假”低阻層。在進(jìn)行資料解釋工作時(shí)，需對(duì)這種現(xiàn)象加以考慮，避免造成與真實(shí)電阻率變化特征相反的錯(cuò)誤解釋。

圖6 為三層H 型模型和K 型模型的ρs和ρd及與設(shè)計(jì)模型的對(duì)比結(jié)果。由圖6 可知，ρs和ρd對(duì)低阻薄層反映較好，ρs完整體現(xiàn)出了大地電阻率高-低-高的變化趨勢(shì)，ρd基本上準(zhǔn)確地恢復(fù)了薄層和背景大地的電阻率值。對(duì)比薄層頂、底界面的恢復(fù)電阻率結(jié)果可以發(fā)現(xiàn)，ρd對(duì)頂界面反映較為清晰，而底界面較為模糊，低阻范圍延續(xù)較長(zhǎng)。這是由于隨深度增加，瞬變電磁場(chǎng)對(duì)地層的分辨能力迅速降低。ρs和ρd對(duì)高阻薄層反映效果均欠佳，僅能體現(xiàn)出高阻薄層的存在和大致位置，未能準(zhǔn)確恢復(fù)薄層電阻率值。在一定程度上而言，這是由瞬變電磁法對(duì)于高阻層的不靈敏特性所造成的。高阻層的存在本身對(duì)于瞬變電磁場(chǎng)分布和傳播特征的影響，相對(duì)低阻層而言較弱。

圖6 三層模型計(jì)算電阻率隨深度的變化Fig.6 Variation of calculated resistivity with depth in 3-layer models

4 體積效應(yīng)探討

瞬變電磁法作為一種基于擴(kuò)散場(chǎng)理論的感應(yīng)類(lèi)電磁勘探方法，本身固有體積效應(yīng)。式(9)基于電磁場(chǎng)的垂向擴(kuò)散速度給出了大地電阻率的改進(jìn)計(jì)算方法，在一定程度上削弱了垂向體積效應(yīng)的影響，所獲得的大地電阻率可以近似看作井孔中相鄰接收點(diǎn)間介質(zhì)的綜合效應(yīng)，但其對(duì)應(yīng)不同深度的徑向體積效應(yīng)仍需進(jìn)一步探討。

在均勻半空間條件下，瞬變電磁感生渦流的極大值同時(shí)向下和向外傳播[20,28]，其包絡(luò)面為一個(gè)與地面夾角約為25°的錐形斜面。據(jù)此可導(dǎo)出包絡(luò)面半徑(R)與深度的關(guān)系：

式中：a為發(fā)射回線(xiàn)半徑。

根據(jù)渦流電場(chǎng)所產(chǎn)生的磁場(chǎng)分布規(guī)律，在該半徑范圍內(nèi)，瞬變電磁磁場(chǎng)響應(yīng)較為均勻且具有較大的強(qiáng)度，并在超出該半徑范圍后快速衰減，此時(shí)可將井孔中接收點(diǎn)接收到的瞬變電磁信號(hào)看作主要是渦流半徑范圍內(nèi)介質(zhì)的綜合效應(yīng)。由式(10)可以發(fā)現(xiàn)，通過(guò)前文算法所得出大地電阻率的體積效應(yīng)取決于發(fā)射回線(xiàn)的半徑和接收點(diǎn)深度：發(fā)射回線(xiàn)半徑越大、接收點(diǎn)深度越深，體積效應(yīng)越強(qiáng)。

5 實(shí)測(cè)數(shù)據(jù)試算

為進(jìn)一步檢驗(yàn)本文算法的有效性和實(shí)用性，選取了某礦區(qū)的地-井TEM 實(shí)測(cè)數(shù)據(jù)進(jìn)行了試算。所測(cè)量鉆孔的終孔深度為2 000 m，受鉆孔套管和井況限制，地-井TEM 實(shí)測(cè)深度范圍為950～1 700 m，該井段同步開(kāi)展了0.5 m 普通電阻率測(cè)井工作。地-井TEM 工作采用的觀(guān)測(cè)參數(shù)為：方形發(fā)射回線(xiàn)，邊長(zhǎng)300 m，發(fā)射電流20 A，發(fā)射基頻12.5 Hz，觀(guān)測(cè)時(shí)窗0.052～15.96 ms，接收點(diǎn)垂直間隔5 m。

圖7a 為不同深度上的dBz/dt響應(yīng)曲線(xiàn)。由圖7a可知，隨深度增大，dBz/dt響應(yīng)曲線(xiàn)的峰值時(shí)刻逐漸后移。圖7b 為根據(jù)dBz/dt響應(yīng)曲線(xiàn)拾取的沖激時(shí)刻曲線(xiàn)，可以發(fā)現(xiàn)沖激時(shí)刻總體隨深度增加而變晚。但受拾取誤差、噪聲和測(cè)量誤差等因素的影響，沖激時(shí)刻曲線(xiàn)存在振蕩。沖激時(shí)刻曲線(xiàn)的變化率反映了電磁場(chǎng)的擴(kuò)散速度特征，整體呈現(xiàn)“慢-快-慢”趨勢(shì)，表明地下介質(zhì)電阻率呈“低-高-低”趨勢(shì)變化。采用本文算法分別計(jì)算了ρs和ρd，并與電阻率測(cè)井結(jié)果(ρl)進(jìn)行了對(duì)比，如圖7c 所示。由圖7c 可知：ρs較好地反映出了大地電阻率“低-高-低”的變化趨勢(shì)，與電阻率測(cè)井結(jié)果的總體趨勢(shì)一致；ρd更加靈敏地反映出了大地電阻率值的變化，且與電阻率測(cè)井結(jié)果基本吻合，特別是對(duì)于1 500 m 左右深度的低阻界面位置反映較為精準(zhǔn)。該實(shí)測(cè)數(shù)據(jù)算例表明，基于瞬變沖激時(shí)刻的視電阻率定義算法能夠較好地反映出大地電阻率的整體變化趨勢(shì)，基于瞬變電磁場(chǎng)擴(kuò)散速度的改進(jìn)算法能較為準(zhǔn)確地恢復(fù)大地電阻率值，對(duì)電性界面反映較好，具有較強(qiáng)的實(shí)用性。

圖7 實(shí)測(cè)數(shù)據(jù)計(jì)算結(jié)果Fig.7 Calculated results of the field data

6 結(jié)論

a.借鑒關(guān)于TEM“煙圈”效應(yīng)和沖激時(shí)刻的已有研究成果，探討了地-井瞬變電磁響應(yīng)的沖激時(shí)刻特征，給出了dBz/dt沖激時(shí)刻與大地電導(dǎo)率和深度之間的函數(shù)關(guān)系。在此基礎(chǔ)上，提出了一種基于瞬變沖激時(shí)刻的地-井TEM 視電阻率定義方法，并改進(jìn)大地電阻率的精確恢復(fù)算法。模型計(jì)算和實(shí)測(cè)數(shù)據(jù)試算結(jié)果表明，基于沖激時(shí)刻的視電阻率定義算法能夠較好地反映設(shè)計(jì)模型的電性變化趨勢(shì)，基于擴(kuò)散速度的大地電阻率計(jì)算方法能夠更為精確地恢復(fù)出大地電阻率值和電性界面，可用于地-井TEM 數(shù)據(jù)的初步定量解釋。

b.基于瞬變沖激時(shí)刻的地-井TEM 定量解釋方法相比傳統(tǒng)的定性和半定量解釋方法，具有直接獲取大地電阻率參數(shù)的優(yōu)勢(shì)。相比反演解釋方法，不需要進(jìn)行復(fù)雜模型的正演迭代，具有較高的計(jì)算效率。

c.模型計(jì)算結(jié)果表明，基于沖激時(shí)刻的視電阻率定義方法仍然存在“overshoot”和“undershoot”現(xiàn)象，在實(shí)際資料處理中，應(yīng)加以考慮，避免造成錯(cuò)誤解釋。

d.鑒于瞬變電磁電場(chǎng)E和磁場(chǎng)B的優(yōu)良特性，利用E場(chǎng)和B場(chǎng)沖激時(shí)刻的視電阻率定義方法或能避免“overshoot”和“undershoot”現(xiàn)象，具備進(jìn)一步研究的價(jià)值。

e.基于瞬變沖激時(shí)刻的地-井TEM 定量解釋方法屬于一種快速、初步的定量方法，僅利用了地-井TEM 數(shù)據(jù)的沖激時(shí)刻特征，其余大量的有效信息并未參與計(jì)算。而地-井TEM 方法深入地下展開(kāi)測(cè)量，既包含有隨深度增加的“幾何測(cè)深”信息，又包含有隨時(shí)間增加的“時(shí)間測(cè)深”信息，獲取的地電信息十分豐富。先利用沖激時(shí)刻、視電阻率定義等方法提供初始地電模型，再利用高維反演或近似高維反演獲取鉆孔周邊和深部的高分辨電性分布信息，是地-井瞬變電磁法數(shù)據(jù)處理解釋的一條可取途徑。