相位激電和時(shí)域激電對(duì)激電效應(yīng)響應(yīng)關(guān)系研究

郭 鵬 ， 肖 都 ，石福升 ， 林品榮

(中國(guó)地質(zhì)科學(xué)院 地球物理地球化學(xué)勘查研究所, 廊坊 065000)

0 引言

近年來(lái)，相位激電法以其快速、輕便、對(duì)弱異常發(fā)現(xiàn)能力強(qiáng)等特點(diǎn)廣受關(guān)注，并在礦產(chǎn)普查與詳查中得到廣泛應(yīng)用[1-4]。相位激電法測(cè)量交流電場(chǎng)作用下的巖(礦)石的電化學(xué)極化效應(yīng),即測(cè)量電極間總場(chǎng)電位差相對(duì)供電電流所產(chǎn)生的相位差,即復(fù)電阻率的相位角[5]，它反映的是極化效應(yīng)的“純異常”。在低頻工作條件下,相位移的大小可以反映地下介質(zhì)激電效應(yīng)的強(qiáng)弱。對(duì)于一定的頻率,激電效應(yīng)越強(qiáng),負(fù)相位移的絕對(duì)值越大；反之,激電效應(yīng)越弱,負(fù)相位移的絕對(duì)值越小。激電相位與時(shí)間域激電測(cè)量的極化率同樣是反映巖(礦)石激電效應(yīng)強(qiáng)弱的參數(shù)，二者之間是否存在某種聯(lián)系，學(xué)者們?cè)谶@方面做了一些工作，分別得到了不同的結(jié)果[6-7]。作者結(jié)合柯?tīng)?柯?tīng)柲Ｐ团c實(shí)測(cè)數(shù)據(jù)對(duì)激電相位和極化率之間的聯(lián)系做了分析。

1 數(shù)值模擬激電相位與極化率的變化規(guī)律

在激電的頻譜特性研究中，國(guó)內(nèi)外學(xué)者通常引用柯?tīng)?柯?tīng)柲Ｐ停浔磉_(dá)式如式(1)。

(1)

其中：ρ0為零頻時(shí)的電阻率；m為極化率；c為頻率相關(guān)系數(shù)；τ為時(shí)間常數(shù)。

柯?tīng)?柯?tīng)柲Ｐ兔枋隽思る娦?yīng)引起的復(fù)電阻率隨頻率的變化。根據(jù)式(1)即可計(jì)算出用柯?tīng)?柯?tīng)柲Ｐ捅磉_(dá)式表示的復(fù)電阻率的虛、實(shí)分量和振幅、相位頻譜。

從柯?tīng)?柯?tīng)柲Ｐ偷玫较辔槐磉_(dá)式

φ(ω)=

(2)

式(2)中同時(shí)包含了相位和極化率參數(shù)，但其關(guān)系不是簡(jiǎn)單的對(duì)應(yīng)關(guān)系，相位是一個(gè)包含極化率、時(shí)間常數(shù)、頻率和頻率相關(guān)系數(shù)的綜合參數(shù)，很難根據(jù)式(2)簡(jiǎn)單定義相位和極化率的相互關(guān)系。為此考慮使用數(shù)值模擬的方法，通過(guò)繪制極化率與相位關(guān)系圖，來(lái)分析極化率和相位的關(guān)系。數(shù)值模擬中分別考慮了不同頻率相關(guān)系數(shù)、時(shí)間常數(shù)和不同頻率下的相位隨極化率的變化情況，并在不同坐標(biāo)系下進(jìn)行繪圖如圖1所示。

圖1 相位隨極化率變化圖Fig．1 Phase change with polarization

從圖1中可以看到，在算術(shù)坐標(biāo)下，極化率從0 %～100 %范圍內(nèi)(圖1中(a)、(e)、(i))，在不同的c值、τ值和f值時(shí)，相位都隨極化率的增大而增大，呈現(xiàn)非線性變化；而在極化率較小時(shí)(≤10%)，相位和極化率在算術(shù)坐標(biāo)下呈現(xiàn)出近似線性變化的關(guān)系(圖1中(c)、(g)、(k))，只是由于c值、τ值和f值的差異，該曲線的斜率存在不同。由此我們給出在極化率較小時(shí)，相位與極化率近似滿足關(guān)系式

φ≈k·m+b

(3)

其中：k為曲線的斜率，其大小與頻率、時(shí)間常數(shù)和頻率相關(guān)系數(shù)有關(guān)；b為線段在y軸上的截距。

在雙對(duì)數(shù)坐標(biāo)下，相位隨極化率的增大而單調(diào)增加，極化率從0 %～100 %范圍內(nèi)呈現(xiàn)非線性變化(圖1中(b)、(f)、(j))；在極化率較小時(shí)(≤10%)相位與極化率呈現(xiàn)出近似線性變化的規(guī)律(圖1中(d)、(h)、(l))。曲線斜率在不同頻率、時(shí)間常數(shù)和頻率相關(guān)系數(shù)下基本保持不變。

從上述相位隨極化率的變化規(guī)律給出在雙對(duì)數(shù)坐標(biāo)下，極化率較小時(shí)(≤10 %)相位隨極化率變化的近似關(guān)系式

logφ≈k1·logm+b1

(4)

其中：k1為線段的斜率，為一常數(shù)；b1為線段在y軸上的截距，其大小與頻率、時(shí)間常數(shù)和頻率相關(guān)系數(shù)有關(guān)。

由圖1結(jié)合上面關(guān)系式(3)與關(guān)系式(4)分析可得，在極化率≤10%的情況下(這種情況包括了絕大多數(shù)時(shí)間域激電實(shí)測(cè)視極化率數(shù)值)，在算術(shù)坐標(biāo)和雙對(duì)數(shù)坐標(biāo)下，相位與極化率都呈現(xiàn)出近似線性關(guān)系。不同點(diǎn)是，在算術(shù)坐標(biāo)下，頻率、時(shí)間常數(shù)和頻率相關(guān)系數(shù)的變化，影響曲線斜率k的變化；在雙對(duì)數(shù)坐標(biāo)下頻率、時(shí)間常數(shù)和頻率相關(guān)系數(shù)的變化，則是影響到線段在y軸上的截距b1的變化。

2 實(shí)測(cè)激電相位隨視極化率變化規(guī)律

在得到算術(shù)坐標(biāo)和雙對(duì)數(shù)坐標(biāo)系下相位與極化率的關(guān)系式(3)與關(guān)系式(4)后，究竟哪一個(gè)關(guān)系式更貼近實(shí)際情況，為此在內(nèi)蒙某礦區(qū)進(jìn)行了野外同裝置、同測(cè)點(diǎn)的相位激電和時(shí)間域激電數(shù)據(jù)采集工作，將野外采集數(shù)據(jù)進(jìn)行了算術(shù)坐標(biāo)系和雙對(duì)數(shù)坐標(biāo)系下的對(duì)比。

工作中分別采集了中梯裝置和偶極-偶極裝置下的時(shí)間域視極化率和頻率域激電相位數(shù)據(jù)。時(shí)間域激電選用的供電周期為16 S，供電電流為7 A，數(shù)據(jù)采集斷電延時(shí)的時(shí)間為100 ms，數(shù)據(jù)的采樣寬度為160 ms。相位激電觀測(cè)頻率為4 S,供電的電流為5 A。中梯裝置的供電極距為1 500 m，接收極距為40 m；偶極-偶極裝置AB=MN=80 m，點(diǎn)距為40 m，隔離系數(shù)n= 3～8。

圖2是中梯裝置下，將采集的65個(gè)測(cè)點(diǎn)的相位激電和時(shí)間域激電數(shù)據(jù)繪制成圖。從圖2中可見(jiàn),算術(shù)坐標(biāo)系下數(shù)據(jù)點(diǎn)的分布較雙對(duì)數(shù)坐標(biāo)系下的離散度要小。算術(shù)坐標(biāo)系下的數(shù)據(jù)更接近于線性關(guān)系。

圖3是偶極-偶極裝置下采集的513個(gè)測(cè)點(diǎn)的實(shí)測(cè)相位和視極化率數(shù)據(jù)繪制成圖。從圖3中可見(jiàn),算術(shù)坐標(biāo)系下，極化率數(shù)據(jù)點(diǎn)分布，較雙對(duì)數(shù)坐標(biāo)系下的離散度要小的多。從數(shù)據(jù)點(diǎn)總體分布來(lái)看，算術(shù)坐標(biāo)系下的數(shù)據(jù)點(diǎn)分布更接近于線性關(guān)系。

從野外實(shí)測(cè)數(shù)據(jù)分析發(fā)現(xiàn)，激電相位和視極化率在算術(shù)坐標(biāo)系下更接近于線性關(guān)系，因此式(3)更接近于實(shí)測(cè)相位和極化率的分布規(guī)律。

根據(jù)關(guān)系式(3)，使用直線來(lái)擬合這些數(shù)據(jù)點(diǎn)，確定k值和b值如圖4所示。計(jì)算得中梯裝置下k值為6.86，b值為 -3.18,激電相位與視極化率相互轉(zhuǎn)換關(guān)系式為：

φ=6.86·η-3.18；η=(φ+3.18)/6.86

(5)

偶極-偶極裝置下k值為6.84，b值為-0.75，激電相位與視極化率相互轉(zhuǎn)換關(guān)系式為：

φ=6.84·η-0.75;η=(φ+0.75)/6.84

(6)

3 驗(yàn)證關(guān)系式

在各個(gè)參數(shù)確定后，進(jìn)行了視極化率和激電相位的相互轉(zhuǎn)換對(duì)比，驗(yàn)證轉(zhuǎn)換公式的適用性。圖5和圖6是該礦區(qū)，中梯裝置觀測(cè)的激電相位和視極化率相互轉(zhuǎn)換對(duì)比圖。從圖中可見(jiàn)，由視極化率轉(zhuǎn)換的激電相位與實(shí)測(cè)激電相位數(shù)值大小非常相近，曲線形態(tài)完全一致；由激電相位轉(zhuǎn)換的視極化率，與實(shí)測(cè)視極化率的數(shù)值和曲線形態(tài)，同樣具有良好的一致性。

圖2 中梯裝置不同坐標(biāo)系下視相位與視極化率關(guān)系圖

圖3 偶極-偶極裝置不同坐標(biāo)系下相位與視極化率關(guān)系圖Fig．3 Relation of apparent phase and apparent polarization under dipole-dipole device in different coordinate system

圖7和圖8是該礦區(qū)122線偶極-偶極裝置下，繪制的視極化率和激電相位相互轉(zhuǎn)換斷面對(duì)比圖。如圖中所示，視極化率和激電相位相互轉(zhuǎn)換繪制的斷面圖，與實(shí)測(cè)視極化率和激電相位斷面圖異常形態(tài)、異常位置相同，實(shí)測(cè)數(shù)值與轉(zhuǎn)換數(shù)值大小接近，具有很好的一致性。

4結(jié)論

通過(guò)數(shù)值模擬和實(shí)測(cè)數(shù)據(jù)分析表明，激電相位和時(shí)間域激電極化率對(duì)激電效應(yīng)的反映是相同的，在一定條件下，激電相位和視極化率可進(jìn)行相互轉(zhuǎn)換。

圖4 直線擬合數(shù)據(jù)點(diǎn)圖Fig．4 Fitting a straight line

圖5 內(nèi)蒙某礦區(qū)122線激電相位與視極化率相互轉(zhuǎn)換對(duì)比圖

圖6 內(nèi)蒙某礦區(qū)126線激電相位與視極化率相互轉(zhuǎn)換對(duì)比圖

圖7 內(nèi)蒙古某礦區(qū)122線激電相位轉(zhuǎn)換視極化率測(cè)深斷面對(duì)比圖

圖8 內(nèi)蒙古某礦區(qū)122線視極化率轉(zhuǎn)換激電相位測(cè)深斷面對(duì)比圖

參考文獻(xiàn)：

[1] 林品榮．電磁法綜合探測(cè)系統(tǒng)研究[J]．地質(zhì)學(xué)報(bào)，2006，80(10)：1539-1548.

[2] 郭鵬，肖都，石福升．陣列相位激電法在弱極化異常區(qū)的試驗(yàn)效果[J] 物探與化探，2012，36(5)：772-774.

[3] 劉瑞德，陸桂福，張強(qiáng)，等．相位激發(fā)極化在寧波油氣勘探中的應(yīng)用[J]．地質(zhì)與勘探，2004，40(增刊)：14-16.

[4] 李建華，林品榮，徐寶利，等．陣列相位激電法與時(shí)域激電法的對(duì)比分析[J]．物探化探計(jì)算技術(shù)，2007，29(增刊)：98-100.

[5] 劉崧．譜激電法[M]．北京：中國(guó)地質(zhì)大學(xué)出版社，1998．

[6] 王書(shū)民，雷達(dá)．相位激電法(偶極一偶極)單頻電磁耦合校正方法[J]物探與化探，2002，26(1)：57-59.

[7] 徐寶利，林品榮．巖礦石標(biāo)本激電相位和極化率測(cè)定及其關(guān)系探討[J]．物探化探計(jì)算技術(shù)，2006，28(增刊):81-83.