王 益，王 鵬

(1.中煤科工集團(tuán)西安研究院有限公司，西安 710077; 2.長安大學(xué)環(huán)境科學(xué)與工程學(xué)院，西安 710054)

?

地面瞬變電磁數(shù)據(jù)的早期畸變與校正方法

王 益1，2，王 鵬1

(1.中煤科工集團(tuán)西安研究院有限公司，西安 710077; 2.長安大學(xué)環(huán)境科學(xué)與工程學(xué)院，西安 710054)

地面瞬變電磁勘探實(shí)踐中發(fā)現(xiàn)，不同邊長外框采集的數(shù)據(jù)在早期存在規(guī)律性畸變,推斷關(guān)斷時間和視電阻率轉(zhuǎn)換算法均有可能引起畸變,為此采用試驗方法進(jìn)行驗證。通過計算證實(shí)：實(shí)際關(guān)斷時間并不足以引發(fā)電磁數(shù)據(jù)畸變；而在常規(guī)情況下的早晚期時間節(jié)點(diǎn)均存在數(shù)據(jù)畸變，且畸變主要發(fā)生在節(jié)點(diǎn)前，畸變程度與回線大小正相關(guān)；采用晚期公式計算是導(dǎo)致畸變的主要原因。使用二分搜索法對試驗數(shù)據(jù)進(jìn)行視電阻率轉(zhuǎn)換，可以有效的克服電磁數(shù)據(jù)的早期畸變，獲得與實(shí)際地層更吻合的視電阻率曲線。對不同外框試驗數(shù)據(jù)的計算表明校正效果隨外框邊長的增大而減弱。

瞬變電磁早期畸變；關(guān)斷時間；早晚期節(jié)點(diǎn)；全區(qū)視電阻率；二分搜索法

0 引言

時間域瞬變電磁法由于其探測深度大、高阻穿透能力強(qiáng)、受地形影響相對較小和施工效率高的特點(diǎn)，已廣泛應(yīng)用于金屬礦勘探[1]、油氣田勘探[2]和煤田水文地質(zhì)勘探[3～6]中。目前常用的裝置介于大定源和中心回線之間，即在地面敷設(shè)較大的正方形發(fā)射回線，在回線中心1/3邊長范圍內(nèi)進(jìn)行數(shù)據(jù)采集，近似認(rèn)為該區(qū)域為均勻場，電阻率轉(zhuǎn)換按中心回線裝置進(jìn)行計算。

地面瞬變電磁數(shù)據(jù)采集中一個重要的參數(shù)是發(fā)射回線的大小，一般根據(jù)最大探測深度予以確定，但鮮見關(guān)于不同大小回線信號的研究。瞬變電磁觀測信號分為早期和晚期，對早期信號響應(yīng)特征的研究相對較少，楊海燕等[7]對多匝回線的自感影響早期信號的特征進(jìn)行了研究，F(xiàn).Kamenetsky等[8]從計算線圈自感的影響上分析早期信號的響應(yīng)特性，李創(chuàng)社等[9]從儀器線路設(shè)計上研究早期信號的檢測方法。常規(guī)地面瞬變電磁信號的視電阻率轉(zhuǎn)換一般按晚期公式進(jìn)行計算，這導(dǎo)致得到的早期視電阻率曲線畸變，與實(shí)際值偏差較大。由于瞬變場的表達(dá)式是電阻率和時間的隱函數(shù)，不能通過解析式從觀測值得到巖層的視電阻率。在全區(qū)視電阻率計算方面，已有眾多學(xué)者做了大量研究。白登海等[10]給出了一種時間域瞬變電磁法視電阻率的數(shù)值計算方法, 該方法根據(jù)中心方式磁場垂直分量時間變化率的核函數(shù)的表現(xiàn)特征, 把整個瞬變過程分為早期階段、早期到晚期的轉(zhuǎn)折點(diǎn)和晚期階段。分別得到早期視電阻率和晚期視電阻率, 然后通過轉(zhuǎn)折點(diǎn)構(gòu)成一條完整的全程視電阻率曲線。楊生等[11]給出了中心回線裝置發(fā)射電流為斜階躍波形條件下全區(qū)視電阻率迭代反演計算方法。李建平等[12]把回線分解為水平電偶極子,然后給出利用電偶極子求取全區(qū)視電阻率的方法。王華軍等[13]依據(jù)均勻半空間瞬變響應(yīng)曲線隨地下電導(dǎo)率、發(fā)射回線邊長與觀測時間具有平移伸縮特性, 提出了一種直接計算全區(qū)視電阻率的方法。付志紅等[14]研究了斜階躍場源瞬變電磁法的全程視電阻率數(shù)值計算。陳清禮基于[15]均勻半空間中心回線觀測方式瞬變電磁測深法的感應(yīng)電動勢隨電阻率的增大而單調(diào)下降的特性，設(shè)計出計算全區(qū)視電阻率的一種快速且精確的算法，張成范等[16]通過對理論模型和實(shí)際數(shù)據(jù)的計算證明了該方法的有效性和正確性。

地面瞬變電磁實(shí)際工作中發(fā)現(xiàn)，當(dāng)采用不同邊長的外框進(jìn)行相同測點(diǎn)的測量時，相應(yīng)的衰減值曲線存在有規(guī)律性的畸變。對邊長越大的外框，轉(zhuǎn)換的早期視電阻率值越高，但晚期均趨于一致。針對這種現(xiàn)象，在陜北某地進(jìn)行針對性試驗，以驗證這種現(xiàn)象。從關(guān)斷時間、早晚期節(jié)點(diǎn)與全區(qū)算法方面，對畸變現(xiàn)象進(jìn)行了研究與校正。

1 早期畸變現(xiàn)象

大量實(shí)測結(jié)果表明，回線邊長對視電阻率的影響較大。對同一測點(diǎn)進(jìn)行瞬變電磁測深，其它參數(shù)不變，只改變回線大小，采集的衰減值曲線在早期相交，而視電阻率曲線在早期分叉。針對此現(xiàn)象，在陜北某煤田進(jìn)行試驗，圖1為120、240、360m回線在相同測點(diǎn)觀測的感應(yīng)曲線和換算的視電阻率曲線，采取雙對數(shù)坐標(biāo)。

圖1 不同大小回線觀測結(jié)果Figure 1 Observed results from different sized loops

從圖1可以看出，改變發(fā)射回線邊長后，采集的感應(yīng)曲線和換算的視電阻率曲線均存在明顯不同。對感應(yīng)曲線，表現(xiàn)為早期時間內(nèi)大回線的感應(yīng)值相對低，小回線的感應(yīng)值相對高，即感應(yīng)值與回線大小負(fù)相關(guān)；晚期時間差距明顯縮小。對視電阻率曲線，晚期時間各曲線基本重合，而早期時間分叉加大，表現(xiàn)為回線越長則電阻率值越高。由瞬變電磁方法的特點(diǎn)可知，時間-視電阻率曲線反映地層深度上的導(dǎo)電性變化，而相同測點(diǎn)下地層導(dǎo)電性是唯一的，因此上面不同大框得到的結(jié)果并不統(tǒng)一這個特點(diǎn)與實(shí)際不符。

2 關(guān)斷時間的影響

關(guān)斷時間是影響早期信號的一個重要因素。由于完全斷掉大框電流需要一定的時間，因此很容易地聯(lián)想到是否因為外框越大越難關(guān)斷，使其二次場受到影響而產(chǎn)生畸變。試驗設(shè)備為V8工作站，關(guān)斷時間的估算公式：

L=0.005×S1.127×N2

TΓ=I×L

式中L——電感，mH；S——發(fā)射框邊長，m；N——發(fā)射框匝數(shù)；T?！P(guān)斷時間，μs；l——發(fā)射電流，A。

試驗中120、240、360m大框?qū)?yīng)的電流分別是16A、13.7A、12.9A，按照上面公式計算關(guān)斷時間分別為17.6μs、33.0μs、49.0μs。試驗數(shù)據(jù)第一道時間為79.0μs，大于計算的關(guān)斷時間，在忽略接收系統(tǒng)過渡過程對響應(yīng)影響的條件下，可認(rèn)為原始數(shù)據(jù)并沒有受到關(guān)斷時間的影響。

3 早晚期節(jié)點(diǎn)

薛國強(qiáng)[17]在對瞬變電磁視電阻率數(shù)據(jù)向平面波場測深視電阻率數(shù)據(jù)轉(zhuǎn)換時，發(fā)現(xiàn)由于瞬變電磁使用晚期計算公式及裝置問題, 使測深曲線早期數(shù)據(jù)發(fā)生畸變。在建立一個兩層介質(zhì)模型，使用不同大小的框進(jìn)行正演計算，按晚期公式轉(zhuǎn)換的視電阻率曲線所表現(xiàn)出來的規(guī)律，與試驗數(shù)據(jù)相同。該項研究對這種現(xiàn)象的解釋是：不同裝置的早晚期時間界限不同，小框的早期時間相對較短，更早地進(jìn)入晚期時間；大框的早期時間較長，進(jìn)入晚期時間相對較晚。地表電阻率的相對高低也同樣影響進(jìn)入晚期時間的界限。對均勻半空間的正演結(jié)果表明，當(dāng)接收時間段滿足下式時，按晚期公式計算的視電阻率值與真值的偏差小于10%。

式中 t——晚期開始時間，ms； S——大框邊長，m； ρ——地表電阻率，Ω·m。

針對目前常用的120、240、360、480、600、720、840、960m大框邊長，分別計算出晚期開始時間t1、t2、t3、t4、t5、t6、t7、t8(表1)。

試驗區(qū)地表電阻率約為100Ω·m，對試驗數(shù)據(jù)進(jìn)行晚期節(jié)點(diǎn)的計算，可得到120、240、360m大框在試驗區(qū)進(jìn)入晚期的時間節(jié)點(diǎn)，分別為0.1856、0.7425、1.6705ms。為驗證晚期電阻率與真值的誤差，分別將240m和360m晚期視電阻率曲線與120m曲線取比值，見圖2。圖中已將晚期節(jié)點(diǎn)標(biāo)注在對應(yīng)的曲線上，120m大框最早進(jìn)入晚期，其次是240m外框和360m外框。由于120m外框最早進(jìn)入晚期時間，節(jié)點(diǎn)后的曲線可近似視為真值。240m外框與360m外框在節(jié)點(diǎn)處的視電阻率比值均小于10%，說明實(shí)際數(shù)據(jù)表現(xiàn)的特征與理論計算一致。

以上分析表明，用晚期公式計算的瞬變電磁視電阻率曲線，在晚期節(jié)點(diǎn)之前畸變，使之與真值存在明顯偏差。這種偏差在理論和實(shí)測數(shù)據(jù)上均表現(xiàn)與大框邊長正相關(guān)的規(guī)律，且大框越長，畸變時間越長。

表1 晚期時間節(jié)點(diǎn)統(tǒng)計表

圖2 不同大框視電阻率比值曲線與晚期節(jié)點(diǎn)Figure 2 Apparent resistivity ratio curves and late phase nodes from different sized loops

4 全區(qū)視電阻率計算實(shí)例

根據(jù)鉆孔資料，試驗區(qū)淺部為風(fēng)積沙，中部為新近系紅土，深部為延安組含煤地層(主要成分為中、粗砂巖)和三疊系永坪組粉、細(xì)砂巖。將鉆孔測井曲線數(shù)字化，可得試驗區(qū)巖層隨深度的電阻率變化，見圖3(a)。圖3(a)為側(cè)向電阻率測井曲線及其4階擬合結(jié)果，該圖顯示由淺至深巖層電阻率逐漸升高，隨后出現(xiàn)下降趨勢，整體表現(xiàn)為“低-高-低”的地電斷面。

采用陳清禮的二分搜索法[15]對試驗數(shù)據(jù)進(jìn)行

全區(qū)視電阻率計算。圖3(b)(c)(d)分別為120、240、360m外框的晚期曲線和全區(qū)曲線。觀察該圖可以發(fā)現(xiàn)，120m外框數(shù)據(jù)在早期整體呈現(xiàn)最小電阻率；240m外框的早期曲線大幅下降，接近于最低阻；360m外框數(shù)據(jù)在早期也大幅下降，但整體仍為高阻，與實(shí)際值具有一定的偏差?？傮w而言，經(jīng)過全區(qū)計算后，結(jié)果曲線在早期時間均有一定程度的下降，與真實(shí)值更趨于吻合。

需要注意的是，盡管使用二分搜索法得到的全區(qū)視電阻率曲線與實(shí)際更加吻合，但對不同的外框，得到的全區(qū)曲線并不一致，這說明該算法并不能完全解決全區(qū)視電阻率的轉(zhuǎn)換問題。全區(qū)計算對早期電阻率計算不能完全恢復(fù)的原因在于，各全區(qū)計算方法均采用數(shù)值算法，公式推導(dǎo)中為便于計算，忽略一些影響因素(如過渡過程)，使得計算結(jié)果與真實(shí)值存在一定偏差，且不同算法的偏差并不相同。此處采用效果更好的二分搜索法，已對真實(shí)電阻率進(jìn)行最大程度恢復(fù)。

圖3 測井曲線與不同大框電阻率曲線Figure 3 Logging traces and resistivity curves of different sized loops

5 結(jié)論

在發(fā)現(xiàn)不同邊長外框得到的視電阻率曲線存在畸變后，展開針對性試驗，在對試驗數(shù)據(jù)分析、試算的基礎(chǔ)上，對畸變原因、特征及校正方法進(jìn)行總結(jié)，得出以下主要結(jié)論：

1)隨著大框邊長的加大，早期信號變?nèi)酢０赐砥诠接嬎愕碾娮杪是€在早期明顯畸變，表現(xiàn)為外框越大，早期電阻率值越高。

2)對關(guān)斷時間的計算結(jié)果表明，各邊長大框的關(guān)斷時間均小于第一道采樣時間。在忽略過渡過程影響的條件下，這種畸變現(xiàn)象與關(guān)斷時間無關(guān)。

3)產(chǎn)生畸變現(xiàn)象的原因是在不適用的早期段使用晚期公式，早晚期時間節(jié)點(diǎn)與大框邊長和地表電阻率有關(guān)。框越大、地表電阻率越低，早晚期節(jié)點(diǎn)越靠后。

4)使用全區(qū)算法能改善早期畸變現(xiàn)象。計算實(shí)例表明，二分搜索法轉(zhuǎn)換的全區(qū)曲線與實(shí)際結(jié)果更吻合，但修正效果隨框的加大而減弱。

[1]柳建新,劉春明,童艷光.強(qiáng)場源瞬變電磁儀及在某銅礦的試驗性應(yīng)用[J].物探與化探,2004,28(2):130-132.

[2]唐新功,胡文寶,嚴(yán)良俊，等.瞬變電磁法油藏動態(tài)監(jiān)測模擬[J].石油物探,2004,43(2):192-196.

[3]周韜,韓自豪,周建雄，等. 瞬變電磁法在煤礦水害防治中的應(yīng)用[J].中國煤田地質(zhì),2004,16(2):44-46.

[4]葉劍湘.瞬變電磁測深法找煤效果研究[J].中國煤田地質(zhì),2000,12(4):71-74.

[5]王惠亮.應(yīng)用瞬變電磁法探測煤層風(fēng)氧化帶[J].采礦技術(shù),2004,4(2):61-62.

[6]路軍臣,蘇維濤,張濟(jì)懷.瞬變電磁法在探測小窯采空區(qū)中的應(yīng)用[J].河北煤炭,2002,2:39-40.

[7]楊海燕,岳建華,胡文武，等.多匝回線的自感對瞬變電磁早期信號的影響特征[J].物探化探計算技術(shù),2007.3:96-99.

[8]F.Kamenetsky and Chr.Oelsner. Distortions of EM transients in coincident loops at short-time-dalay[J]. Geophysical Prospecting, 2000, 48(6):983-993.

[9]李創(chuàng)社,李實(shí),宋建平，等.瞬變電磁勘探法中早期信號的檢測方法研究[J].煤田地質(zhì)與勘探,1994,27(2):58.

[10]白登海,Meju M A,盧健，等.時間域瞬變電磁法中心方式全程視電阻率的數(shù)值計算[J].地球物理學(xué)報,2003,46(5):697-704.

[11]楊生.TEM中心回線法計算考慮關(guān)斷時間的全區(qū)視電阻率[J].物探與化探,2003,32(6): 647-6511.

[12]李建平,李桐林,趙雪峰，等.層狀介質(zhì)任意形狀回線源瞬變電磁全區(qū)視電阻率的研究[J]球物理學(xué)進(jìn)展, 2007,22(6):1777-1780.

[13]王華軍.時間域瞬變電磁法全區(qū)視電阻率的平移算法[J].地球物理學(xué)報,2008,51(6):1936-1942.

[14]付志紅,孫天,陳清禮，等.斜階躍場源瞬變電磁法的全程視電阻率數(shù)值計算[J].電工技術(shù)學(xué)報,2008,23 (11):15-211.

[15]陳清禮.瞬變電磁法全區(qū)視電阻率的二分搜索算法[J].石油天然氣學(xué)報(江漢石油學(xué)院學(xué)報)，2009.4,31(2):45-49.

[16]張成范,翁愛華,孫世棟，等.計算矩形大定源回線瞬變電磁測深全區(qū)視電阻率[J].吉林大學(xué)學(xué)報,2009.7,39(4):755-758.

[17]薛國強(qiáng),李貅.瞬變電磁測深早期數(shù)據(jù)的修正[J].地球物理學(xué)進(jìn)展,2008.2,23(1):220-224

Surface TEM Data Early Phase Distortion and Calibration Method

Wang Yi1, 2, Wang Peng1

(1. Xi’an Research Institute, China Coal Technology and Engineering Group Corp, Xi’an, Shaanxi 710077;2. School of Environmental Science and Engineering, Chang’an University, Xi’an, Shaanxi 710054)

During the surface TEM prospecting practices, the early phase distortion can be found in data acquired from different sized loops have existed early phase regular distortions. It is inferred that the turn-off time and apparent resistivity conversion algorithm all can cause distortion, therefore, test procedure has been used to verify. Through computation has justified that the practical turn-off time is not enough to trigger electromagnetic data distortion. While under normal circumstances, early and late phase time nodes have data distortions, and mostly happened before the nodes; distortion degree is positively correlated with loop size; using late phase formula computation is the main reason to cause distortion. Through dichotomizing search carried out tested data apparent resistivity conversion, can effectively overcome early phase electromagnetic data distortion, obtained apparent resistivity curves are more accorded with actual strata. The computation of different sized loops tested data has shown that the calibrated effect is attenuating along with the increasing of loop size.

TEM data early phase distortion; turn-off time; early and late phase time nodes; region-dependent apparent resistivity; dichotomizing search

10.3969/j.issn.1674-1803.2017.06.15

1674-1803(2017)06-0074-05

王益(1982—)，男，四川人，工程師，在讀碩士研究生，長安大學(xué)環(huán)境科學(xué)與土木工程學(xué)院地質(zhì)工程專業(yè)，從事煤田電法實(shí)踐與研究工作。

2017-04-01

文獻(xiàn)標(biāo)識碼：A

責(zé)任編輯：孫常長