石加玉,郭鵬,李勇
(1.中國(guó)地質(zhì)大學(xué)(北京)地球物理與信息技術(shù)學(xué)院,北京 100083;2.中國(guó)地質(zhì)調(diào)查局 地球物理調(diào)查中心,河北 廊坊 065000;3.中國(guó)地質(zhì)科學(xué)院 地球物理地球化學(xué)勘查研究所,河北 廊坊 065000)
頻譜激電法(spectrum induced polarization,SIP,簡(jiǎn)稱譜激電)是直流電法勘探的一種,其基于幾何原理實(shí)現(xiàn)測(cè)深,是恒流情況下時(shí)間域激電法的一種擴(kuò)展,為針對(duì)交變電磁場(chǎng)實(shí)施的頻率域激發(fā)極化方法。譜激電測(cè)量的是交流電場(chǎng)作用下巖(礦)石的電化學(xué)極化效應(yīng),即接收電極間總場(chǎng)電位差相對(duì)于發(fā)射電流的相位差—復(fù)電阻率的相位角,它反映的是極化效應(yīng)的“純異?!盵1-3]。近年來(lái),譜激電在資源、環(huán)境和工程等領(lǐng)域得到前所未有的發(fā)展和應(yīng)用,越來(lái)越被工程技術(shù)人員和學(xué)者所重視。利用譜激電可以得到Cole-Cole模型復(fù)電阻率表達(dá)式中4個(gè)表征被探測(cè)目標(biāo)的物性參數(shù):零頻率時(shí)的電阻率、極化率、時(shí)間常數(shù)和頻率相關(guān)系數(shù),這4個(gè)參數(shù)為區(qū)分礦與非礦以及不同類型的礦種提供了重要信息[4-5]。其中,極化率等價(jià)于譜激電的極化相位,這是一個(gè)重要參數(shù):在低頻工作條件下,相位的大小可以反映地下介質(zhì)激電效應(yīng)的強(qiáng)弱,對(duì)于一定的頻率,激電效應(yīng)越強(qiáng),負(fù)相位的絕對(duì)值越大;反之,激電效應(yīng)越弱,負(fù)相位的絕對(duì)值越小。
“十一五”期間,在中國(guó)地質(zhì)調(diào)查局的支持下研制了頻譜激電測(cè)量?jī)x器原型機(jī),并利用原型機(jī)在福建、海南、云南、西藏、新疆、內(nèi)蒙古等礦區(qū)開(kāi)展了方法技術(shù)應(yīng)用示范,取得了良好的應(yīng)用效果,尤其在弱極化區(qū)、接地困難區(qū)和干擾區(qū),頻譜激電法顯示出比常規(guī)激電法更強(qiáng)的異常發(fā)現(xiàn)能力和抗干擾能力。為了充分發(fā)揮頻譜激電法在資源、環(huán)境和工程等領(lǐng)域的作用,“十二五”期間對(duì)頻譜激電原型機(jī)進(jìn)行了完善和改進(jìn),形成可供推廣的實(shí)用化儀器,并將該技術(shù)向全國(guó)進(jìn)行推廣。本文對(duì)頻譜激電測(cè)量中的關(guān)鍵技術(shù),如高精度同步技術(shù)、高精度穩(wěn)流技術(shù)、相位提取技術(shù)和弱信號(hào)同步相關(guān)檢測(cè)技術(shù)等做了充分探討并給出了實(shí)現(xiàn)方法,最后對(duì)實(shí)用化儀器和原型機(jī)在內(nèi)蒙古某礦區(qū)的對(duì)比試驗(yàn)結(jié)果進(jìn)行了分析,結(jié)果表明:基于本文所述技術(shù)方案設(shè)計(jì)的頻譜激電測(cè)量?jī)x器實(shí)用化程度高,相位測(cè)量重現(xiàn)性好。目前,該儀器已向全國(guó)推廣,實(shí)現(xiàn)了預(yù)期目標(biāo)。
Cole-Cole模型復(fù)電阻率表達(dá)式中的4個(gè)物性參數(shù),只有極化相位必須通過(guò)測(cè)量才能獲得,其他幾個(gè)參數(shù)可通過(guò)電位差的幅值、裝置系數(shù)和發(fā)射電流等計(jì)算得到;因此,頻譜激電測(cè)量?jī)x器的關(guān)鍵技術(shù)就集中體現(xiàn)在相位測(cè)量這一環(huán)節(jié)。相位測(cè)量涉及的關(guān)鍵技術(shù)包括高精度同步技術(shù)、恒流供電技術(shù)和相位檢測(cè)技術(shù)等。
譜激電測(cè)量?jī)x器屬于野外工作設(shè)備,發(fā)射機(jī)和接收機(jī)常常需要分開(kāi)工作,某些觀測(cè)裝置甚至要求發(fā)射與接收之間相隔數(shù)千米。本文研究一種基于同步相關(guān)檢測(cè)的絕對(duì)相位測(cè)量技術(shù),發(fā)射與接收之間不受距離限制,接收電路的結(jié)構(gòu)也相對(duì)簡(jiǎn)單。由于采用了數(shù)字式窄帶選頻與相位提取技術(shù),提高了測(cè)量?jī)x器的抗干擾能力和實(shí)用性。該技術(shù)實(shí)現(xiàn)的前提是:①發(fā)射與接收之間在時(shí)間上嚴(yán)格同步;②發(fā)射電流要維持恒定。
常規(guī)的基于硬件電路的相位測(cè)量方式需要在同一塊電路板上對(duì)發(fā)射電流和接收電壓進(jìn)行比較,譜激電測(cè)量?jī)x器屬于野外工作設(shè)備,發(fā)射機(jī)和接收機(jī)相距較遠(yuǎn),這種相位測(cè)量方法不再是最優(yōu)方案。本研究方案采用高精度同步技術(shù),使發(fā)射機(jī)和接收機(jī)在時(shí)間上嚴(yán)格同步,實(shí)現(xiàn)高精度絕對(duì)相位的測(cè)量。
同步精度設(shè)計(jì)規(guī)則是以滿足譜激電法對(duì)相位觀測(cè)的最高精度為原則,計(jì)算公式為
(1)
式中:f為最高工作頻率,單位Hz;n為譜激電觀測(cè)系統(tǒng)要求的相位分辨率,單位mRad;t為同步精度,單位s。本文研究的相位測(cè)量?jī)x器,其相位觀測(cè)精度要求為1 mRad,最高工作頻率為128 Hz,通過(guò)上式計(jì)算出t=1.24 μs。實(shí)際設(shè)計(jì)電路時(shí),同步精度應(yīng)留有足夠的余量,亦即t<1.24 μs。
常用的同步方式包括有線同步、恒溫晶振同步、電臺(tái)同步和GPS同步等。有線同步中的長(zhǎng)導(dǎo)線會(huì)產(chǎn)生延時(shí),距離越遠(yuǎn)延時(shí)越大,同步距離較遠(yuǎn)時(shí)施工很不方便;采用恒溫晶體的晶體同步方式精度較高,但由于各個(gè)晶體之間始終存在頻差,頻差引起的時(shí)間誤差會(huì)不斷累積,所以經(jīng)常需要對(duì)鐘以消除該誤差,當(dāng)同步鐘的數(shù)量超過(guò)3臺(tái)時(shí),同步和對(duì)鐘將很繁瑣;電臺(tái)同步的精度較低,發(fā)射機(jī)和各個(gè)接收機(jī)之間容易受障礙物遮擋,同步距離有限[6-8]。本方案基于GPS+溫補(bǔ)晶振同步設(shè)計(jì),GPS 秒脈沖只具有隨機(jī)誤差、不具有累積誤差的優(yōu)點(diǎn)可以消除同步時(shí)鐘源的誤差,溫補(bǔ)晶振與恒溫晶振相比功耗顯著降低,可有效降低儀器的體積和重量。該設(shè)計(jì)無(wú)啟動(dòng)時(shí)間,可提高儀器的工作效率[9-10];在保證同步精度的同時(shí),可降低儀器設(shè)備的成本,便于推廣應(yīng)用[11]。同步原理框圖如圖1所示,GPS模塊為U-BLOX公司的LEA-5T,授時(shí)精度為30 ns,溫補(bǔ)晶振的頻率準(zhǔn)確度為1×10-7,同步系統(tǒng)的最大誤差為GPS授時(shí)精度與晶振的1個(gè)時(shí)鐘周期之和。本方案中GPS授時(shí)精度為30 ns,晶體頻率為16.384 000 MHz,對(duì)應(yīng)的時(shí)鐘周期為1/163 840 00 s=61 ns,故最大同步誤差為30+61=91 ns,遠(yuǎn)小于測(cè)量?jī)x器要求的1.24 μs,滿足最高工作頻率為128 Hz時(shí)的相位測(cè)量精度。
圖1 高精度同步原理框圖
采用GPS同步時(shí),同步時(shí)鐘由1 PPS的上升沿觸發(fā),主控單片機(jī)讀取GPS模塊的授時(shí)信息,通過(guò)PA口設(shè)置同步時(shí)序控制電路的工作方式,產(chǎn)生發(fā)射機(jī)的工作頻率。采用外同步時(shí),同步時(shí)鐘由外同步脈沖的上升沿觸發(fā)。發(fā)射機(jī)和接收機(jī)采用相同的同步方案,保證發(fā)射電流和接收電壓之間嚴(yán)格同步,如圖2所示。在圖2a同步時(shí)鐘的驅(qū)動(dòng)下,發(fā)射機(jī)以恒流方式供出圖2b所示的理想方波電流。接收機(jī)的接收電壓波形如圖2c所示,與發(fā)射電流波形相比,上升沿變得圓滑,這是地下地質(zhì)體的極化效應(yīng)所致。頻譜激電的觀測(cè)參數(shù)之一就是該極化相位,因此,高精度的同步技術(shù)是相位測(cè)量的重要保證。
圖2 同步時(shí)鐘、發(fā)射電流和接收電壓時(shí)序圖
基于同步相關(guān)檢測(cè)的相位提取技術(shù)中,將發(fā)射電流方波作為參考基準(zhǔn),要求發(fā)射電流為理想方波(如圖2b)。實(shí)際工作中,供電電極深埋于地下,受極化作用的影響接地電阻(發(fā)射機(jī)的負(fù)載)隨時(shí)間發(fā)生變化,為保持發(fā)射電流恒定,使其不受負(fù)載效應(yīng)和源效應(yīng)(電源電壓波動(dòng))的影響,發(fā)射機(jī)必須具有良好的穩(wěn)流機(jī)制。實(shí)踐表明,穩(wěn)流精度應(yīng)達(dá)到0.1% 或更高[12],觀測(cè)到的絕對(duì)相位才有實(shí)際意義。換言之,如果發(fā)射電流不穩(wěn)定,接收端觀測(cè)到的相位可能是發(fā)射電流不穩(wěn)定導(dǎo)致的偽相位,而不是地下目標(biāo)體的真實(shí)極化相位。
圖3是基于ARM單片機(jī)和IGBT功率開(kāi)關(guān)設(shè)計(jì)的高精度穩(wěn)流原理框圖,由穩(wěn)流ARM、隔離驅(qū)動(dòng)、IGBT功率開(kāi)關(guān)、電流隔離檢測(cè)和低通濾波等組成。逆變橋及連接在AB端的大地或假負(fù)載構(gòu)成穩(wěn)流電路的負(fù)載。為了避免調(diào)制脈沖(PWM)與發(fā)射電流波形(AB輸出)之間由于時(shí)鐘源的差異產(chǎn)生系統(tǒng)相位差,本方案中的主控CPU和穩(wěn)流ARM采用同一個(gè)溫補(bǔ)晶振作為工作時(shí)鐘,以消除系統(tǒng)自身相位誤差。該工作時(shí)鐘由同步時(shí)序控制電路產(chǎn)生。
圖3 高精度穩(wěn)流原理框圖
工作時(shí),穩(wěn)流ARM(STM32F103VBT6)根據(jù)主控CPU設(shè)定的電流值啟動(dòng)穩(wěn)流功能。STM32控制內(nèi)部的ADC(12位)進(jìn)行電流采樣,對(duì)電流樣進(jìn)行數(shù)字濾波并與設(shè)定的電流值進(jìn)行比較,若供電輸出電流小于設(shè)定的電流,STM32將增加調(diào)制脈沖的寬度(片內(nèi)16位PWM控制器),使輸出電流增大;反之,若供電輸出電流大于設(shè)定的電流,STM32將會(huì)減小調(diào)制脈沖的寬度,使輸出電流變小。此過(guò)程無(wú)限循環(huán),保持供電電流恒定。
穩(wěn)流精度按照(1/2N)×100%計(jì)算,N為ADC的位數(shù);這里N=12,故穩(wěn)流精度的理論值為0.02%,實(shí)際精度達(dá)到了0.1%。由于電流采樣、電流濾波、脈寬調(diào)制均為數(shù)字形式,因此本穩(wěn)流方案的抗干擾能力強(qiáng)、穩(wěn)流精度高。
常規(guī)相位測(cè)量方法如圖4所示,硬件電路同時(shí)對(duì)圖4a所示發(fā)射電流波形和圖4b所示接收電壓波形做上升沿檢測(cè),得到如圖4d、e所示的脈沖,測(cè)量這兩個(gè)脈沖的時(shí)間差,再將其換算為對(duì)應(yīng)供電周期的相位差,如圖4e中的φ1。圖4b是一種理想的接收電壓波形,實(shí)際工作中,由于接收電極埋設(shè)于大地,電極接收到的是包含天然電磁場(chǎng)和各種人文干擾的混合信號(hào)(如圖4c),這種包含噪聲的信號(hào)經(jīng)放大整形和上升沿檢測(cè)后如圖4f所示;該脈沖與發(fā)射電流脈沖圖4d比較,測(cè)得的相位如圖4f的φ2所示,并非目標(biāo)體的極化相位。由于干擾和噪聲的頻率和幅度不固定,導(dǎo)致相位測(cè)量發(fā)生偏差;在測(cè)量電路中加入工頻陷波器,效果改善也不明顯,而陷波器帶來(lái)的相位移將使處理電路變得更為復(fù)雜,增加調(diào)試工作量和儀器成本。
圖4 基于硬件電路的相位檢測(cè)原理
野外觀測(cè)時(shí),被測(cè)信號(hào)常常淹沒(méi)于噪聲中,為了從干擾環(huán)境中有效提取微弱信號(hào),本文開(kāi)發(fā)了一種基于高精度同步的弱信號(hào)相關(guān)檢測(cè)技術(shù)——同步數(shù)字相關(guān)檢測(cè)技術(shù),該技術(shù)實(shí)現(xiàn)的前提條件是發(fā)射電流恒定、發(fā)射與接收之間在時(shí)間上嚴(yán)格同步,其原理框圖見(jiàn)圖5。
圖5 數(shù)字同步相關(guān)檢測(cè)原理框圖
接收電極MN兩端的信號(hào)經(jīng)過(guò)阻抗變換、信號(hào)放大,直接送入模數(shù)轉(zhuǎn)換器(ADC)進(jìn)行轉(zhuǎn)換。同步時(shí)序控制電路為模數(shù)轉(zhuǎn)換器和參考信號(hào)發(fā)生器提供同步工作時(shí)鐘,該時(shí)鐘與發(fā)射機(jī)端的發(fā)射電流嚴(yán)格同步。ADC轉(zhuǎn)換后的數(shù)字信號(hào)在相關(guān)檢測(cè)器中進(jìn)行同頻正弦變換和余弦變換,由于正弦函數(shù)和余弦函數(shù)是一對(duì)正交基,所以經(jīng)過(guò)變換就可得到被測(cè)信號(hào)的實(shí)分量和虛分量,再利用實(shí)分量和虛分量計(jì)算出被測(cè)信號(hào)的幅值和相位。相關(guān)檢測(cè)器中的乘法器和積分器可以采用硬件實(shí)現(xiàn),也可以采用軟件實(shí)現(xiàn),由于ARM處理器的運(yùn)算速度足夠快,故本方案中采用軟件實(shí)現(xiàn)。接收電壓信號(hào)經(jīng)過(guò)正弦和余弦變換后,若實(shí)部表示為A、虛部表示為B,則該信號(hào)的幅值、相位及電阻率的計(jì)算公式如下:
(2)
信號(hào)相位:φ=1000×arctan(B/A),
(3)
(4)
式中:R為信號(hào)幅值,單位V;φ為接收電壓滯后于發(fā)射電流的絕對(duì)相位,單位mRad;K為裝置系數(shù);I為發(fā)射電流,單位A;ρ為視電阻率,單位Ω·m。
本方案與常規(guī)相位測(cè)量電路相比,接收機(jī)省去了陷波器、信號(hào)整形、時(shí)間測(cè)量、乘法器和積分器等,接收機(jī)電路結(jié)構(gòu)更簡(jiǎn)單、調(diào)試工作量小,采用同步相關(guān)檢測(cè)和窄帶數(shù)字濾波技術(shù)可有效濾除工作頻率之外的干擾信號(hào),具有更強(qiáng)的抗干擾能力。場(chǎng)地試驗(yàn)結(jié)果表明:該技術(shù)可在強(qiáng)干擾背景下有效提取周期性弱信號(hào),效果優(yōu)于基于硬件電路的信號(hào)檢測(cè)方法;在弱極化區(qū)工作時(shí)極化相位的幅度更大,反映異常的效果更明顯[13]。
基于上述方案對(duì)原型儀器進(jìn)行了改進(jìn)和完善,形成實(shí)用化的頻譜激電測(cè)量?jī)x器,并在內(nèi)蒙古西烏珠穆沁旗某礦區(qū)開(kāi)展了實(shí)用化儀器與原型儀器的對(duì)比試驗(yàn)。
試驗(yàn)采用陣列方式的偶極—偶極觀測(cè)裝置,如圖6所示:AB為供電電極,MN為接收電極,接收機(jī)為雙通道,2個(gè)通道對(duì)應(yīng)的接收電極分別為M1-N1和M2-N2。發(fā)射與接收間采用高精度GPS同步方式,一臺(tái)發(fā)射機(jī)供電,多臺(tái)雙通道接收機(jī)按陣列方式同時(shí)測(cè)量,一次供電過(guò)程就可獲得多個(gè)深度的觀測(cè)信息,觀測(cè)深度根據(jù)工作需要選擇相應(yīng)的隔離系數(shù)1,2,…,n;發(fā)射機(jī)和接收機(jī)沿測(cè)線跑極一遍,就可獲得該測(cè)線的相位和電阻率等地電斷面資料。
圖6 野外工作示意
儀器主要指標(biāo):工作頻率128、64、…、2、1、1/2、…、1/128 Hz;最高供電電壓600 V;最大發(fā)射電流1.5A;相位測(cè)量精度優(yōu)于1 mRad;GPS同步精度±30 ns;功率高壓電源采用24 V鋰電池升壓產(chǎn)生。
試驗(yàn)參數(shù):實(shí)用化儀器和原型儀器各投入1臺(tái)發(fā)射機(jī),實(shí)用化接收機(jī)2臺(tái)(編號(hào)1、2),原型接收機(jī)3臺(tái)(編號(hào)16、17、18)。工作方式:軸向偶極—偶極裝置;AB=MN=80 m;點(diǎn)距40 m;工作頻率4、1、0.25 Hz。實(shí)用化儀器和原型儀器的隔離系數(shù)均為n=1、2、3、4。
實(shí)用化儀器和原型儀器的一致性對(duì)比在測(cè)量工作開(kāi)始之前進(jìn)行。在相位異常較大的地段選11個(gè)測(cè)點(diǎn),每臺(tái)儀器均在相同條件下進(jìn)行往返多次觀測(cè),隔離系數(shù)選用1和2,每臺(tái)儀器的第1道與第1道對(duì)比、第2道與第2道對(duì)比。視電阻率采用均方相對(duì)誤差衡量,視相位采用總均方誤差衡量,一致性對(duì)比結(jié)果列于表1。評(píng)判依據(jù)為《陣列相位激電法工作手冊(cè)》(試用稿,2008年9月由中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所編制),該手冊(cè)規(guī)定儀器的相位誤差不應(yīng)大于設(shè)計(jì)總精度的2/3。從表中數(shù)據(jù)可看出,5臺(tái)儀器的一致性檢測(cè)結(jié)果均滿足手冊(cè)要求。
表1 儀器一致性檢測(cè)結(jié)果
實(shí)用化儀器和原型儀器均完成3條剖面的測(cè)量工作,有效剖面長(zhǎng)度4 440 m,總測(cè)點(diǎn)數(shù)114個(gè)。實(shí)用化儀器和原型儀器檢查點(diǎn)數(shù)分別為12個(gè)和10個(gè),檢查點(diǎn)數(shù)占總測(cè)點(diǎn)數(shù)比例分別為10.53%和 8.77%。經(jīng)計(jì)算,實(shí)用化儀器視電阻率均方相對(duì)誤差為2.17%,視相位總均方誤差為0.30 mRad;原型儀器視電阻率均方相對(duì)誤差為3.38%,視相位總均方誤差為0.60 mRad,滿足《陣列相位激電法工作手冊(cè)》檢查觀測(cè)點(diǎn)數(shù)占總觀測(cè)點(diǎn)數(shù)的3%~5%、視電阻率檢查觀測(cè)的均方相對(duì)誤差應(yīng)小于5%、視相位檢查觀測(cè)的總均方誤差應(yīng)小于1 mRad之要求。
為了客觀反映實(shí)用化儀器和原型儀器對(duì)相位異常的發(fā)現(xiàn)能力,兩種儀器開(kāi)展了同測(cè)線、同測(cè)點(diǎn)、同深度的對(duì)比觀測(cè),并繪制成剖面圖,相位和電阻率均為原始數(shù)據(jù),未經(jīng)任何處理或編輯。由于篇幅所限,僅給出124號(hào)線的對(duì)比情況,圖7為2種儀器在124號(hào)線的視電阻率和視相位觀測(cè)結(jié)果對(duì)比,該圖展示了3個(gè)頻率的觀測(cè)結(jié)果,每個(gè)頻率均包括3個(gè)隔離系數(shù)的觀測(cè)數(shù)據(jù)。黑色曲線表示實(shí)用化儀器的測(cè)量結(jié)果,紅色曲線表示原型儀器的測(cè)量結(jié)果,兩種儀器觀測(cè)的曲線形態(tài)相同,數(shù)據(jù)對(duì)應(yīng)關(guān)系較好。
圖7 124線實(shí)用化儀器與原型儀器視電阻率、視相位測(cè)量結(jié)果對(duì)比
對(duì)比試驗(yàn)表明,采用數(shù)字式同步相關(guān)檢測(cè)技術(shù)設(shè)計(jì)的相位測(cè)量?jī)x器,具有輕便、快速、觀測(cè)精度高和抗干擾能力強(qiáng)等特點(diǎn),能夠在接地條件差、人文干擾強(qiáng)等常規(guī)激電難以開(kāi)展工作的地區(qū)開(kāi)展工作。在發(fā)射電流小、接收信號(hào)弱、強(qiáng)干擾等復(fù)雜地質(zhì)環(huán)境下,可實(shí)現(xiàn)對(duì)激電異常的有效探測(cè)[14]。
深入研究了頻譜激電測(cè)量?jī)x器的關(guān)鍵技術(shù),在此基礎(chǔ)上對(duì)原型樣機(jī)做了改進(jìn)和完善,形成實(shí)用化的適合推廣的頻譜激電測(cè)量?jī)x器。全新設(shè)計(jì)的數(shù)字式弱信號(hào)同步相關(guān)檢測(cè)技術(shù)方案,省去了接收機(jī)中的陷波器、信號(hào)整形、時(shí)間測(cè)量、乘法器和積分器等硬件電路,簡(jiǎn)化了電路結(jié)構(gòu),提高了相位檢測(cè)的穩(wěn)定性,降低了儀器的成本,減小了儀器調(diào)試的工作量,更適合推廣應(yīng)用。礦區(qū)對(duì)比試驗(yàn)表明,基于新方案設(shè)計(jì)的實(shí)用化儀器與原型儀器相比,具有同等發(fā)現(xiàn)異常的能力和更強(qiáng)的抗干擾能力。希望本文所述關(guān)鍵技術(shù)及實(shí)現(xiàn)方案對(duì)同行及同類儀器研制有借鑒意義。