程輝， 崔峻卿， 付國(guó)紅， 傅崧原

1 湖南科技大學(xué)先進(jìn)礦山裝備教育部工程研究中心， 湘潭 411201 2 桂林理工大學(xué)廣西隱伏金屬礦產(chǎn)勘查重點(diǎn)實(shí)驗(yàn)室， 桂林 541006

0 引言

巖礦石電性參數(shù)研究是電法(電磁法)領(lǐng)域的一項(xiàng)重要基礎(chǔ)性工作，其研究程度直接影響到方法理論研究與勘查實(shí)踐工作．眾多學(xué)者進(jìn)行了卓有成效的研究工作．Wait(1959a，b)、Madden 和Cantwell(1967)、Ward和Fraser(1976)、Zonge等(1972)、Zonge和Wynn(1975)、Pelton(1977)、Pelton等(1978)、Dias(2000)、Gurin等(2015)、Zhdanov等(2018)、Gurin等(2018)進(jìn)行了大量露頭、標(biāo)本以及人造標(biāo)本的測(cè)定工作，先后提出了不同的導(dǎo)電機(jī)理與導(dǎo)電模型以及相應(yīng)的參數(shù)定義，如阻容模型、德拜模型、瓦爾堡模型、Cole-Cole模型、Dias模型、Gemtip模型等；何繼善(2006)、張賽珍等(1984)、王慶乙和徐立忠(2012)、底青云等(2008，2009)、薛國(guó)強(qiáng)等(2020)、陳儒軍等(2004)、黃理善等(2014)、柯式鎮(zhèn)等(2006)、李建軍等(2005)、趙云山等(2015)、何蘭芳(2014)、程輝等(2010，2014)、付國(guó)紅等(2020)進(jìn)行了不同的實(shí)驗(yàn)與理論研究工作．

以上研究多是基于變頻法進(jìn)行的頻率特性測(cè)量，然而，使用變頻法觀測(cè)有著明顯不足(何繼善，2010)：(1)頻率改變前后的外部條件與干擾因素各不相同，抗干擾能力相對(duì)較弱；(2)各個(gè)頻點(diǎn)的觀測(cè)時(shí)間不同，觀測(cè)效率較低且增加了觀測(cè)數(shù)據(jù)一致性降低的可能性；(3)為了保證必要的穩(wěn)流條件，會(huì)增加發(fā)送端的復(fù)雜性．

為了達(dá)到實(shí)驗(yàn)?zāi)康?，研究人員使用了各類(lèi)儀器設(shè)備進(jìn)行觀測(cè)，包括阻抗分析儀(Solartron 1250、1260，Agilent4294A，HP4263B等)、LCR測(cè)試儀(數(shù)字電橋測(cè)試儀)(ZL5，Hioki 3522-50等)、巖心測(cè)試設(shè)備(LDT-10發(fā)射機(jī)與接收機(jī)組合形式、GDD公司生產(chǎn)的樣本巖芯IP測(cè)量?jī)xSCIP等)、激電儀、電化學(xué)工作站等．雖然這些觀測(cè)裝置類(lèi)型不一，但大多是基于變頻法思想逐個(gè)頻點(diǎn)觀測(cè)．

何繼善(2010)以2n序列偽隨機(jī)波形作為信號(hào)源應(yīng)用于電法勘探，提出了雙頻激電法與廣域電磁法．近年來(lái)基于相關(guān)辨識(shí)的檢測(cè)技術(shù)應(yīng)用到了電磁法勘探中，趙璧如等(2006)、萬(wàn)健如等(2005)、羅維斌等(2012)、武欣等(2015)、羅延鐘等(2015，2020)、王顯祥等(2016)、王若等(2016)將m序列(或逆重復(fù)m序列)應(yīng)用于電法勘探中，并進(jìn)行相應(yīng)的理論、實(shí)驗(yàn)與野外對(duì)比觀測(cè)研究．若將相關(guān)辨識(shí)技術(shù)應(yīng)用于巖礦石標(biāo)本電性參數(shù)實(shí)驗(yàn)，觀測(cè)數(shù)據(jù)的抗干擾能力將顯著提高，有利于提高測(cè)量精度與測(cè)量效率．

在針對(duì)巖礦石標(biāo)本的室內(nèi)電性參數(shù)觀測(cè)工作中，巖礦石標(biāo)本受其巖性、結(jié)構(gòu)與構(gòu)造、導(dǎo)電礦物類(lèi)型與含量及其導(dǎo)電礦物空間分布與接觸關(guān)系、導(dǎo)電礦物與圍巖的接觸關(guān)系等多方面因素共同影響，即使同一類(lèi)巖性的標(biāo)本，在電性參數(shù)表現(xiàn)上也可能有數(shù)量級(jí)的區(qū)別．針對(duì)觀測(cè)對(duì)象的特殊性與復(fù)雜性，一般地，選用頻率域信號(hào)為激勵(lì)源進(jìn)行實(shí)驗(yàn)時(shí)，有以下要求：(1)實(shí)驗(yàn)頻帶要足夠?qū)挘?2)為保證標(biāo)本的電性參數(shù)為線性響應(yīng)，流過(guò)標(biāo)本的電流激勵(lì)信號(hào)幅值最小僅需維持在10-6A量級(jí)；(3)信號(hào)源要有較強(qiáng)的負(fù)載驅(qū)動(dòng)能力；(4)為了保證觀測(cè)數(shù)據(jù)的可靠，要求發(fā)送端的電流輸出模塊必須保證具有較高的輸出穩(wěn)定度．

本文針對(duì)以上問(wèn)題，提出基于逆重復(fù)m序列偽隨機(jī)電流信號(hào)的巖礦石標(biāo)本電性參數(shù)觀測(cè)方案.研發(fā)了μA級(jí)編碼雙極性電流源，產(chǎn)生實(shí)驗(yàn)所需的高精度小電流，測(cè)量頻帶寬度設(shè)置為10-3～104Hz，通過(guò)合理改變m序列的編碼參數(shù)，覆蓋整個(gè)測(cè)量頻帶，能有效抑制周期性干擾，滿足巖礦石標(biāo)本電性參數(shù)觀測(cè)過(guò)程中的特殊要求．文中介紹了逆重復(fù)m序列的產(chǎn)生方法與頻譜特征，分析了不同階數(shù)序列的測(cè)量效率，針對(duì)工頻干擾提出編碼參數(shù)選擇要求；圍繞巖礦石標(biāo)本電性特征，從觀測(cè)裝置角度論述了發(fā)送端驅(qū)動(dòng)能力及接收端阻抗匹配、相關(guān)辨識(shí)法校正測(cè)試系統(tǒng)頻率響應(yīng)以及高頻微弱電流信號(hào)的產(chǎn)生與拾取等關(guān)鍵技術(shù)問(wèn)題，并對(duì)所研發(fā)的電流源進(jìn)行了實(shí)驗(yàn)驗(yàn)證．

1 逆重復(fù)m序列的產(chǎn)生及其頻譜特征

1.1 逆重復(fù)m序列偽隨機(jī)信號(hào)的產(chǎn)生

逆重復(fù)m序列是將 2 個(gè)周期的m序列隔位取反后得到的一種偽隨機(jī)序列，并且與m序列互不相關(guān)．因此由m序列可方便得到逆重復(fù)m序列．

m序列偽隨機(jī)信號(hào)在電法勘探中采用二進(jìn)制、雙極性的信號(hào)形式，n階m序列的周期長(zhǎng)度為N=2n-1．m序列可由n級(jí)線性反饋移位寄存器生成，也可使用線性遞推公式產(chǎn)生：

ai=p1ai-1⊕p2ai-2⊕…⊕pnai-n,i≥0，

(1)

式中，⊕表示異或運(yùn)算，pi為n階m序列的本原多項(xiàng)式系數(shù)．表1列出部分階數(shù)m序列的遞推公式．

表1 遞推公式表Table 1 Recursive formula

令初始條件a-1,a-2,…,a-n不全為 0 (一般取初始條件全為 1)，則m序列{ai} 可通過(guò)式(1)獲得．通過(guò)映射關(guān)系：

(2)

使用周期為Δt的時(shí)鐘控制ui輸出，則可獲得雙極性m序列偽隨機(jī)信號(hào)，周期為T(mén)=NΔt=N/fc(fc=1/Δt，為碼元頻率)．同樣地對(duì)逆重復(fù)m序列進(jìn)行處理，可得到周期為2T的逆重復(fù)m序列偽隨機(jī)信號(hào)L(t)，圖1為fc=3.1 Hz的5階L(t)波形圖．

圖1 5階L(t)偽隨機(jī)信號(hào)波形(fc=3.1 Hz)Fig.1 Waveforms of 5th-order L(t) pseudo-random singal (fc=3.1 Hz)

1.2 逆重復(fù)m序列的頻譜特征

L(t)屬于循環(huán)平穩(wěn)隨機(jī)過(guò)程，其自相關(guān)函數(shù)為(鐘延炯，1979)：

RLL(τ)=

(3)

其中：

式(3)中第一項(xiàng)為三角脈沖函數(shù)，第二項(xiàng)為周期三角波函數(shù)，其自相關(guān)函數(shù)為三角脈沖函數(shù)與周期三角波之和，歸一化自相關(guān)函數(shù)圖如圖2a所示．由式(3)中各項(xiàng)的系數(shù)可知，階數(shù)n越大，周期長(zhǎng)度N越大，周期三角波的幅值越小，自相關(guān)函數(shù)越接近脈沖函數(shù)，辨識(shí)誤差越?。虼诉x取適當(dāng)大的階數(shù)n，可以提高辨識(shí)精度．

根據(jù)維納—辛欽定理，將L(t)的自相關(guān)函數(shù)RLL(τ)進(jìn)行傅里葉變換即可得到L(t)的功率譜：

(4)

將式(3)代入式(4)，求得：

(5)

考察函數(shù) sinc(f/fc)，令sinc(f/fc)=1/2，求得f≈0.44fc，即 -3 dB帶寬截止頻率．

(6)

圖2 5階L(t)自相關(guān)函數(shù)圖(a)及理論功率譜圖(b)Fig.2 The 5th-order L(t)-signal autocorrelation function (a) and theoretical power spectrum (b)

kmax=0.44(2n-1)+0.5，

(7)

式中，kmax向左取整.

(8)

L(t)的周期公式為：

(9)

因此兩種方法所用時(shí)間之比為：

(10)

根據(jù)式(10)可知，測(cè)量效率與階數(shù)n有關(guān)，n越大，測(cè)量效率提升越明顯，表2列出了部分階數(shù)序列對(duì)應(yīng)的測(cè)量頻點(diǎn)個(gè)數(shù)以及測(cè)量效率提高的倍數(shù)．

表2 不同階數(shù)L序列測(cè)量效率表Table 2 Efficiency of L-sequence measurement with different orders

2 抗工頻干擾參數(shù)選取條件

巖礦石頻譜測(cè)量工作中觀測(cè)數(shù)據(jù)受自身固有特性、外部因素、溶液條件、儀器和實(shí)驗(yàn)條件的共同作用，從干擾類(lèi)型與強(qiáng)度分析，工頻干擾最為明顯且難以屏蔽．由于L(t)具有良好的自相關(guān)性，通過(guò)合理選擇編碼參數(shù)，可以使工頻干擾信號(hào)正交于L(t)，從而通過(guò)相關(guān)辨識(shí)的方法壓制工頻干擾．

工頻干擾為典型的周期性干擾信號(hào)，在接收端表現(xiàn)為L(zhǎng)(t)經(jīng)過(guò)巖礦石標(biāo)本后的響應(yīng)信號(hào)與干擾信號(hào)疊加．設(shè)周期干擾信號(hào)n(t)=Acos(ω0t+φ)，周期為T(mén)0=2π/ω0，其自相關(guān)函數(shù)為：

(11)

功率譜為：

Sn n(ω)=F[Rn n(τ)]

(12)

(13)

式中，f0=1/T0，其自功率譜也是離散譜，只在±f0處取值．設(shè)兩信號(hào)的互功率譜為SLn(f)，則有不等式：

(14)

RLn(τ)=F-1[SLn(f)]=0.

(15)

因此，當(dāng)L(t)信號(hào)中不包含工頻頻率分量時(shí)，其與工頻干擾的互相關(guān)函數(shù)為0，可以有效壓制工頻干擾．考慮到實(shí)際測(cè)量過(guò)程中存在工頻干擾頻率偏移等問(wèn)題，在選取碼元頻率fc時(shí)，可根據(jù)實(shí)際工作條件，避開(kāi)工頻干擾的頻率．

圖3為抑制工頻干擾效果仿真圖，L(t)參數(shù)為：fc=2000 Hz，n=7，在50 Hz附近的頻點(diǎn)為39.37 Hz和55.12 Hz，波形如圖3a．經(jīng)過(guò)圖10所示的系統(tǒng)后，輸出響應(yīng)為圖3b．引入信噪比為-20 dB的50 Hz干擾信號(hào)后，輸出響應(yīng)如圖3c．未被干擾信號(hào)和被干擾信號(hào)的辨識(shí)結(jié)果如圖3d，兩曲線基本吻合，50 Hz附近無(wú)明顯波動(dòng)，按照此要求選擇頻組，具有良好的工頻干擾抑制效果．

圖3 工頻干擾壓制仿真(a) 系統(tǒng)輸入信號(hào); (b) 系統(tǒng)輸出響應(yīng); (c) 被50 Hz干擾的輸出響應(yīng); (d) 測(cè)得系統(tǒng)頻率響應(yīng)對(duì)比.Fig.3 Suppression simulation of power frequency interference(a) System input signal; (b) System output response; (c) Output response disturbed by 50 Hz; (d) Comparison of measured system frequency responses.

3 觀測(cè)系統(tǒng)精細(xì)測(cè)量的解決方法

在實(shí)現(xiàn)精準(zhǔn)的測(cè)量過(guò)程中，面臨電阻率分布跨度大、頻率響應(yīng)存在觀測(cè)誤差大、頻率分布范圍寬等問(wèn)題，在開(kāi)展被測(cè)對(duì)象特征的研究基礎(chǔ)上，確定觀測(cè)系統(tǒng)精細(xì)測(cè)量的解決方法．

針對(duì)巖礦石標(biāo)本電阻率分布跨度大(n×10～106Ωm)，采用高驅(qū)動(dòng)能力設(shè)計(jì)及阻抗匹配技術(shù)；針對(duì)頻率響應(yīng)存在觀測(cè)誤差，采用系統(tǒng)頻率響應(yīng)校正技術(shù)；針對(duì)頻率分布范圍寬(10-3～104Hz)以及非線性效應(yīng)問(wèn)題，采用高頻微弱信號(hào)產(chǎn)生與拾取的相關(guān)技術(shù)．圖4為觀測(cè)系統(tǒng)示意圖，由發(fā)送機(jī)產(chǎn)生逆重復(fù)m序列電流信號(hào)，作為待測(cè)標(biāo)本的激勵(lì)源，使用數(shù)據(jù)采集卡采集標(biāo)本響應(yīng)．

圖4 觀測(cè)系統(tǒng)示意圖Fig.4 Schematic diagram of the observation system

圖5 激發(fā)極化效應(yīng)等效電路(a)及恒流源供電示意圖(b)Fig.5 Equivalent circuit of IP effect (a) and power supply of constant current source (b)

3.1 發(fā)送端驅(qū)動(dòng)能力及接收端阻抗匹配問(wèn)題

巖礦石標(biāo)本電阻率可達(dá)n×106Ωm，以長(zhǎng)l為5 cm，半徑r為2 cm的標(biāo)準(zhǔn)圓柱體標(biāo)本為例，其直流阻值約為n×39.8 MΩ.W.H.Pelton經(jīng)過(guò)對(duì)巖礦石的大量實(shí)驗(yàn)和研究，提出巖礦石激發(fā)極化效應(yīng)可以等效為圖5a的等效電路．

此模型復(fù)阻抗可由式(16)計(jì)算得到：

(16)

恒流信號(hào)源通常采用壓控電流源的方式實(shí)現(xiàn)，其驅(qū)動(dòng)能力主要受限于最高輸出電壓，原理如圖5b所示．恒流源輸出恒定電流，當(dāng)負(fù)載RL發(fā)生變化時(shí)，輸出電壓UO亦隨負(fù)載變化．實(shí)際器件輸出電壓受供電電源限制，當(dāng)負(fù)載足夠大時(shí)，輸出電壓達(dá)到最大值，此時(shí)恒流源將輸出恒定電壓，反饋回路失效，類(lèi)似于恒壓源，電流將隨負(fù)載變大而減小，輸出信號(hào)將嚴(yán)重失真．另一方面，在頻率較高時(shí)，過(guò)大的負(fù)載將與線路間的雜散分布電容CP形成RC網(wǎng)絡(luò)，使τ值變大，導(dǎo)致輸出波形上升沿平緩，增大測(cè)量誤差．

圖6 測(cè)量回路阻抗等效電路Fig.6 Measurement circuit impedance equivalent circuit

對(duì)于接收端而言，測(cè)量回路阻抗等效電路如圖6所示，實(shí)際測(cè)得阻抗為負(fù)載RL與測(cè)量系統(tǒng)輸入阻抗Rin的并聯(lián)阻抗，測(cè)量時(shí)必須加以考慮．

因此在系統(tǒng)設(shè)計(jì)時(shí)，發(fā)送端要保證較高的輸出電壓范圍，選用高帶寬、高壓擺率的運(yùn)放器件，并盡量降低回路分布電容，以擴(kuò)大測(cè)量負(fù)載范圍，提高頻率響應(yīng)范圍．接收端應(yīng)采用差分輸入模式，并設(shè)置幾組不同的輸入阻抗，根據(jù)負(fù)載阻抗大小合理切換，在保證測(cè)量精度的前提下，保持較高的信噪比．

3.2 系統(tǒng)頻率響應(yīng)校正

在實(shí)際工作中，盡管測(cè)試系統(tǒng)采用高性能器件，但系統(tǒng)頻率響應(yīng)誤差仍然存在且不可避免．下面根據(jù)相關(guān)辨識(shí)測(cè)量原理，給出系統(tǒng)頻率響應(yīng)誤差校正方法．相關(guān)辨識(shí)測(cè)量原理如圖7所示，L(t) 經(jīng)過(guò)測(cè)試系統(tǒng)H′(ω)后，引入干擾信號(hào)e(t)，最終測(cè)得輸出信號(hào)y(t)．然后將y(t)與L(t)通過(guò)前文第一小節(jié)的方法，分別計(jì)算出互功率譜和自功率譜并相除，即可求得系統(tǒng)響應(yīng)H′(f)．

圖7 相關(guān)辨識(shí)原理圖Fig.7 Schematic diagram of correlation identification

測(cè)試系統(tǒng)H′(ω)簡(jiǎn)化框圖如圖8所示．信號(hào)L(t)經(jīng)過(guò)發(fā)送通道的V-I轉(zhuǎn)換后變?yōu)殡娏餍图?lì)信號(hào)，激發(fā)標(biāo)本產(chǎn)生響應(yīng)，后經(jīng)接收端測(cè)量通道進(jìn)行采集，記錄為時(shí)間序列y[n]．

圖8 測(cè)量系統(tǒng)信號(hào)傳遞框圖Fig.8 Signal transmission block diagram of the measuring system

對(duì)測(cè)得數(shù)據(jù)處理后，得到的系統(tǒng)響應(yīng)為發(fā)送通道響應(yīng)h0(t)、標(biāo)本響應(yīng)h(t)和測(cè)量通道響應(yīng)h1(t)共同作用的結(jié)果，最終測(cè)得系統(tǒng)頻率響應(yīng)H′(ω)可用式(17)表示：

(17)

其中，H0(ω)、H(ω)、H1(ω)分別對(duì)應(yīng)圖7中各部分的頻率響應(yīng)，X(ω)、Y(ω)為系統(tǒng)輸入輸出采樣信號(hào)x[n]、y[n]的傅里葉變換．由式(18)可得標(biāo)本真實(shí)頻率響應(yīng)的修正公式：

(18)

根據(jù)式(18)可知，采用標(biāo)準(zhǔn)電阻作為負(fù)載，測(cè)得系統(tǒng)頻率響應(yīng)即可修正誤差．

3.3 避免非線性效應(yīng)的高頻弱信號(hào)產(chǎn)生及拾取

對(duì)于巖礦石的非線性效應(yīng)，何繼善等(1995)提出不同礦物出現(xiàn)非線性特征所需的電流密度不同，最小僅1 μA/cm2，因此發(fā)送信號(hào)的電流幅值需實(shí)現(xiàn)μA級(jí)穩(wěn)定輸出．本文提出的觀測(cè)方法，所用信號(hào)采用雙極性形式，且碼元頻率fc是測(cè)量帶寬上限頻率的2.23倍，逆變頻率至少應(yīng)達(dá)到22.3 kHz．微弱電流的高速逆變，不僅要保證電流的直流輸出精度，還要保證高頻時(shí)電流不會(huì)因寄生感抗而衰減，實(shí)現(xiàn)較為困難．經(jīng)過(guò)研究與多次試驗(yàn)，采用圖9所示的運(yùn)放級(jí)聯(lián)技術(shù)，實(shí)現(xiàn)了設(shè)計(jì)要求．輸入級(jí)采用帶寬盡量高的精密運(yùn)放(對(duì)應(yīng)圖9中的A1、A3)，輸出級(jí)采用高速運(yùn)放(對(duì)應(yīng)圖9中的A2)，保證精度與帶寬．輸入信號(hào)Uin為電壓信號(hào)L(t)，經(jīng)過(guò)運(yùn)放電路實(shí)現(xiàn)電壓-電流轉(zhuǎn)換，將電壓信號(hào)轉(zhuǎn)換為電流信號(hào)輸出．輸出電流信號(hào)幅值由輸入信號(hào)幅值決定，通過(guò)調(diào)節(jié)輸入端電壓信號(hào)幅值獲得不同輸出電流幅值．

圖9 恒流源級(jí)聯(lián)原理示意圖Fig.9 Schematic diagram of the flow source cascade

接收端采集的信號(hào)是電壓信號(hào)，在測(cè)量低阻標(biāo)本時(shí)，電壓通常只有n×100 μV，因此測(cè)量端至少需要保證0.1 μV測(cè)量分辨率．根據(jù)采樣定理，采樣頻率至少為44.6 kHz才能保證頻率不失真．在室內(nèi)環(huán)境中，工頻干擾耦合在測(cè)量回路中的幅值可達(dá)幾十毫伏，為了提高測(cè)量信噪比，在分辨率足夠的前提下，可在測(cè)量回路中串聯(lián)精密電阻提高輸出響應(yīng)幅值，再根據(jù)基爾霍夫電壓定理(KVL)校正為標(biāo)本響應(yīng)．

4 實(shí)驗(yàn)測(cè)試

為驗(yàn)證方法的可行性與設(shè)計(jì)儀器的性能，實(shí)驗(yàn)分別對(duì)所研發(fā)電流源樣機(jī)的直流恒流性能和交流測(cè)量精度進(jìn)行測(cè)試，根據(jù)本文提出的測(cè)量方法，研制了信號(hào)發(fā)送樣機(jī)，電路板實(shí)物圖如圖10所示，采用NI公司的NI-PXIe4497數(shù)據(jù)采集卡采集數(shù)據(jù)，最大采樣率為204.8 kS/s，具有24位A/D轉(zhuǎn)換精度．

圖10 樣機(jī)電路板實(shí)物圖Fig.10 Real object diagram of prototype circuit board

直流恒流性能選取輸出電流幅值為2 μA，依次使用阻值為100 Ω、1 kΩ、10 kΩ、100 kΩ、1 MΩ的電阻作為負(fù)載，采用電流-電壓(I-V)轉(zhuǎn)換增益為104的I-V轉(zhuǎn)換電路，測(cè)量I-V轉(zhuǎn)換電路的輸出電壓，測(cè)量結(jié)果如圖11所示．在不同負(fù)載下，電流輸出穩(wěn)定度較高，整體波動(dòng)范圍小于0.5%．從圖11中可以看到，在高阻負(fù)載情況下輸出電流偏低，是因?yàn)殡娮璐嬖谥圃旃睿瑢?dǎo)致恒流源V-I轉(zhuǎn)換電路的反饋回路電阻阻值不嚴(yán)格匹配，可經(jīng)過(guò)調(diào)試進(jìn)一步降低誤差，獲取更高的輸出穩(wěn)定度．

圖11 不同負(fù)載的穩(wěn)定性Fig.11 Stability under different loads

交流精度測(cè)試采用阻容模型模擬巖礦石標(biāo)本，圖12給出了測(cè)量所用阻容模型及其幅頻特性曲線．

圖12 阻容模型及其幅頻特性Fig.12 Resistive and capacitive model and amplitude-frequency characteristics

圖13 轉(zhuǎn)折頻率附近頻組響應(yīng)波形及功率譜(a) (e) fc=62 Hz； (b) (f) fc=625 Hz； (c) (g) fc=6.25 kHz； (d) (h) fc=62 kHz.Fig.13 Response waveforms of the frequency group near the turning frequency and power spectrum

圖14 頻率響應(yīng)對(duì)比圖Fig.14 Comparison of frequency response

發(fā)送信號(hào)詳細(xì)參數(shù)為I=50 μA，n=5，fc=0.062 Hz、0.62 Hz、6.2 Hz、62 Hz、625 Hz、6.25 kHz、62 kHz．測(cè)量頻帶范圍為0.001～11 kHz，按十倍頻設(shè)置頻組，兩相鄰頻組間有重復(fù)頻段，共包含90個(gè)頻點(diǎn)如表4所示，實(shí)現(xiàn)單次單周期覆蓋測(cè)量頻帶用時(shí)約18.5 min，采用變頻法單次單周期測(cè)量以上頻點(diǎn)需要耗時(shí)約42.6 min，與變頻法相比測(cè)量效率提高2.3倍．圖13展示了轉(zhuǎn)折頻率附近頻組的波形及其功率譜，圖13a—d依次為fc=62 Hz、625 Hz、6.25 kHz、62 kHz的時(shí)域響應(yīng)波形圖，圖13e—h為對(duì)應(yīng)的功率譜圖．

圖14為測(cè)得頻率響應(yīng)與理論值對(duì)比圖，從圖中可知測(cè)得幅頻曲線在工頻附近無(wú)明顯變化．頻率較大時(shí)存在系統(tǒng)誤差，按照式(18)的修正方法，計(jì)算出系統(tǒng)頻率響應(yīng)并標(biāo)定后，另取fc=25 kHz的頻組重新測(cè)量，測(cè)得響應(yīng)呈現(xiàn)于圖14中，測(cè)量誤差如表3所示，標(biāo)定后的測(cè)量誤差小于0.5%．

表3 測(cè)量誤差表Table 3 Measurement errors

5 結(jié)論與展望

本文基于逆重復(fù)m序列的編碼原理，論述了逆重復(fù)m信號(hào)的功率譜特點(diǎn)，推導(dǎo)了抗工頻干擾的參數(shù)選擇方法．從觀測(cè)裝置角度詳細(xì)分析了驅(qū)動(dòng)能力及阻抗匹配、測(cè)試系統(tǒng)的頻率響應(yīng)、以及高頻微弱信號(hào)的產(chǎn)生與拾取等關(guān)鍵問(wèn)題．以5階序列為例的實(shí)驗(yàn)結(jié)果表明：

(1)該方法可在10-3～104kHz頻帶內(nèi)測(cè)量90個(gè)有效頻點(diǎn)，實(shí)現(xiàn)寬頻帶、多頻點(diǎn)的觀測(cè)目的；同時(shí)，可根據(jù)需要改變逆重復(fù)m信號(hào)參數(shù)加密頻點(diǎn)，獲得更高的頻率分辨率．

(2)抗工頻干擾能力強(qiáng)，通過(guò)合理選擇信號(hào)參數(shù)，使發(fā)送信號(hào)不包含工頻分量，可以有效壓制工頻干擾．

(3)測(cè)量精度高，通過(guò)系統(tǒng)響應(yīng)修正公式對(duì)系統(tǒng)進(jìn)行標(biāo)定后，測(cè)量誤差小于0.5%．

表4 5階信號(hào)在測(cè)量頻帶內(nèi)包含的頻點(diǎn)Table 4 Frequency points of 5th-order signal contained in the measurement frequency band

(4)測(cè)量效率高，測(cè)量效率由逆重復(fù)m序列階數(shù)決定，5階序列的測(cè)量效率相比變頻法測(cè)量提高2.3倍，在實(shí)際觀測(cè)工作中可根據(jù)需要適當(dāng)提高階數(shù)以獲得更高測(cè)量效率．

(5)恒流源輸出采用分檔連續(xù)可調(diào)方案，在μA 級(jí)精度上可實(shí)現(xiàn)以1 μA為單位的輸出電流調(diào)整.

(6)觀測(cè)數(shù)據(jù)中有效頻點(diǎn)增加，為巖礦石電性參數(shù)計(jì)算與分析提供了基礎(chǔ)數(shù)據(jù)保證.

通過(guò)理論分析與實(shí)驗(yàn)結(jié)果，可以展望下一步工作重點(diǎn)：

(1)接收系統(tǒng)通過(guò)MN拾取大于3 dB有效信號(hào)后利用虛擬儀器技術(shù)上傳至上位機(jī)(PC機(jī))對(duì)時(shí)間波形保存，以便于進(jìn)行后續(xù)計(jì)算與分析.

(2)通過(guò)時(shí)間波形的時(shí)頻計(jì)算，可在一次測(cè)量中，分別求取到時(shí)間域與頻率域電性參數(shù)的各類(lèi)曲線，同時(shí)可利用導(dǎo)電模型，得到模型的各二級(jí)譜參數(shù)，進(jìn)一步增強(qiáng)觀測(cè)數(shù)據(jù)的通用性.

感謝在文章撰寫(xiě)過(guò)程，得到中國(guó)科學(xué)院地質(zhì)與地球物理研究所底青云研究員的指導(dǎo)，在此表示感謝.