石琦，劉麗華，倪志康，劉小軍，方廣有

(1.中國(guó)科學(xué)院空天信息創(chuàng)新研究院，北京，100000；2.電磁輻射與探測(cè)技術(shù)院重點(diǎn)實(shí)驗(yàn)室，北京，100000；3.中國(guó)科學(xué)院大學(xué)電子電氣與通信工程學(xué)院，北京，100000)

瞬變電磁法是目前廣泛應(yīng)用的重要電磁探測(cè)方法之一，根據(jù)負(fù)載形式不同，可分為接地長(zhǎng)導(dǎo)線電性源瞬變電磁系統(tǒng)以及多匝回線磁性源瞬變電磁系統(tǒng)。傳統(tǒng)激勵(lì)信號(hào)為雙極性方波，其預(yù)處理方法本質(zhì)上是用觀測(cè)時(shí)間(時(shí)間效率)換取信噪比，使得系統(tǒng)探測(cè)性能受到較大限制。偽隨機(jī)編碼是近年來(lái)應(yīng)用到瞬變電磁領(lǐng)域的一種新的信號(hào)形式，具有良好的相關(guān)特性，可以有效提高系統(tǒng)信噪比[1]。由于電性源負(fù)載為接地長(zhǎng)導(dǎo)線，其負(fù)載特性偏阻性，偽隨機(jī)編碼信號(hào)相對(duì)容易產(chǎn)生，因此，目前只有電性源的偽隨機(jī)編碼電磁探測(cè)系統(tǒng)，典型代表為MTEM系統(tǒng)(multi-transient electromagnetic，MTEM)，成功應(yīng)用于海洋電磁探測(cè)領(lǐng)域，并取得了很好的探測(cè)效果。QUINCY等[2]利用偽噪聲(后來(lái)的偽隨機(jī)編碼)代替低占空比脈沖；并指出該系統(tǒng)使用較小的天線，并且具有較少的發(fā)射機(jī)功率。DUNCAN等[3]利用長(zhǎng)導(dǎo)線雙極性發(fā)射器發(fā)射偽隨機(jī)編碼序列，使用靈敏的單組磁通門(mén)磁力計(jì)或空心線圈來(lái)監(jiān)測(cè)記錄位置處磁場(chǎng)的垂直分量的時(shí)間變化，經(jīng)過(guò)交叉相關(guān)、解卷積和維納濾波得出地球的脈沖響應(yīng)。WRIGHT等[4]在多通道瞬變電磁法中采用m序列作為電性源發(fā)射系統(tǒng)的激勵(lì)源。ZIOLKOWSKI等[5-6]將偽隨機(jī)編碼信號(hào)運(yùn)用在陸地及海洋MTEM中，采用電性源發(fā)射系統(tǒng)進(jìn)行淺?？睖y(cè)，指出偽隨機(jī)編碼信號(hào)比方波具有更強(qiáng)的分辨能力。長(zhǎng)偏移距瞬變電磁法(LOTEM)是一種收發(fā)距大于等于勘探深度的定源瞬變電磁法，主要用于進(jìn)行深部電性分布的研究。HELWIG[7]在LOTEM系統(tǒng)中將傳統(tǒng)方波替換為m序列，獲得了明顯的信噪比增益。國(guó)內(nèi)，趙璧如等[8-9]分析了電性源的偽隨機(jī)編碼在地電阻率測(cè)量的運(yùn)用。湯井田等[10]提出了一種逆重復(fù)m序列的偽隨機(jī)編碼，討論了影響測(cè)量精度與勘探分辨率的信號(hào)參數(shù)設(shè)計(jì)原則。王顯祥等[11]將偽隨機(jī)編碼運(yùn)用到地面多通道電磁探測(cè)(MTEM)電性源系統(tǒng)中，比較了軸向電場(chǎng)和垂直磁場(chǎng)的探測(cè)效果，指出軸向電場(chǎng)探測(cè)效果更加。薛國(guó)強(qiáng)等[12]總結(jié)了MTEM國(guó)內(nèi)外研究進(jìn)展，在新的MTEM系統(tǒng)中采用偽隨機(jī)編碼發(fā)射源，具有更好的抗噪性，介紹了采用MTEM系統(tǒng)在陸地和海洋進(jìn)行勘探的案例。袁哲等[13]對(duì)不同編碼參數(shù)m序列抗噪性能進(jìn)行了定量評(píng)估，給出了m序列編碼參數(shù)的優(yōu)選方案。目前偽隨機(jī)編碼只應(yīng)用于電性源瞬變電磁系統(tǒng)，還未應(yīng)用于磁性源系統(tǒng)。相比電性源，磁性源的負(fù)載為多匝線圈，是一種感性負(fù)載，其感抗較大，產(chǎn)生高質(zhì)量的偽隨機(jī)編碼電流信號(hào)比較困難，需要對(duì)電流脈沖上升沿和下降沿分別進(jìn)行特殊的提升處理，即對(duì)電流波形進(jìn)行高精度控制的技術(shù)難度較大，因此，目前還沒(méi)有成熟的磁性源偽隨機(jī)編碼電磁探測(cè)系統(tǒng)。磁性源偽隨機(jī)編碼電磁響應(yīng)信號(hào)與電性源偽隨機(jī)編碼電磁響應(yīng)信號(hào)的計(jì)算方法不同，分屬不同的理論模型，本文作者對(duì)磁性源瞬變電磁偽隨機(jī)編碼發(fā)射技術(shù)及電磁響應(yīng)相關(guān)理論進(jìn)行深入研究。依托國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目，將偽隨機(jī)編碼應(yīng)用于磁性源瞬變電磁系統(tǒng)中，設(shè)計(jì)完成滿(mǎn)足磁性源偽隨機(jī)編碼應(yīng)用要求的新型瞬變發(fā)射機(jī)電路并進(jìn)行電路動(dòng)態(tài)過(guò)程分析，研究磁性源偽隨機(jī)編碼瞬變電磁響應(yīng)理論分析方法，由發(fā)射激勵(lì)波形可以得到均勻半空間大地模型下的正演波形。運(yùn)用偽隨機(jī)編碼的預(yù)處理方法得到大地脈沖響應(yīng)的估計(jì)值。最后，進(jìn)行偽隨機(jī)編碼電流的發(fā)射實(shí)驗(yàn)，實(shí)測(cè)結(jié)果與預(yù)處理結(jié)果均表明所發(fā)射的偽隨機(jī)編碼電流質(zhì)量較高。

1 偽隨機(jī)編碼特性分析

1.1 編碼類(lèi)型選擇

偽隨機(jī)碼是一種具有類(lèi)似白噪聲性質(zhì)的碼。白噪聲是一種隨機(jī)過(guò)程，它的瞬時(shí)值服從正態(tài)分布，功率譜在很寬頻帶內(nèi)都是均勻的。它具有類(lèi)似白噪聲的性質(zhì)，相關(guān)函數(shù)具有尖銳的特性，功率譜占據(jù)很寬的頻帶，因此，易于從其他信號(hào)或干擾中分離出來(lái)，具有優(yōu)良的抗干擾特性[14-15]。將其作為瞬變電磁探測(cè)發(fā)射波形，通過(guò)預(yù)處理可以得到更精確的大地脈沖響應(yīng)。

常用的偽隨機(jī)編碼有m序列、Golay互補(bǔ)碼、2n序列、Gold碼[16]。從時(shí)域波形和頻譜特征2個(gè)方面對(duì)4種偽隨機(jī)編碼波形進(jìn)行分析介紹，偽隨機(jī)編碼特征如表1所示。偽隨機(jī)編碼自相關(guān)性越好，則最終算法得出的沖激響應(yīng)越接近于理論值，綜合考慮時(shí)間利用率和帶寬內(nèi)能量分布，最終選擇m序列作為偽隨機(jī)編碼的激勵(lì)波形。

1.2 m序列與雙極性方波特性比較

m序列選擇9階，碼元頻率為10 kHz，采樣頻率80 kHz，發(fā)射周期為5，電流I峰值為1A，時(shí)域波形如圖1(a)所示。其非周期自相關(guān)幅度A如圖1(c)所示，m序列功率譜|S(f)|如圖1(e)所示。雙極性方波波形參數(shù)如下：電流峰值為1A，周期為0.06 s，占空比為50%，重復(fù)周期為9。雙極性方波時(shí)域波形如圖1(b)所示，其非周期自相關(guān)如圖1(d)所示，功率譜如圖1(f)所示。

m序列具有很強(qiáng)的自相關(guān)性，能夠極大地抑制相關(guān)噪聲[17]。從圖1(e)和圖1(f)可以看出：在3 dB帶寬內(nèi)，m序列的功率譜比較平坦且頻譜成分豐富。雙極性方波的功率譜不平坦且?guī)挶萴序列窄。m序列與雙極性方波的特性不同，這些區(qū)別決定了二者數(shù)據(jù)預(yù)處理方法不同。

表1 偽隨機(jī)編碼特性比較Table1 Summary of characteristics of pseudo-random sequences

圖1 波形及其相關(guān)特性Fig.1 Waveforms and related characteristics

2 偽隨機(jī)編碼瞬變電磁發(fā)射系統(tǒng)設(shè)計(jì)

2.1 發(fā)射電路設(shè)計(jì)

發(fā)射機(jī)整體結(jié)構(gòu)框圖如圖2所示，主要包括直流電源、有源鉗位電壓源、功率逆變?nèi)珮蚰骐娐?、MATLAB軟件模塊、數(shù)字邏輯控制器FPGA、發(fā)射圈、驅(qū)動(dòng)電路。其中，所述的直流電源為AC-DC電源或者蓄電池組，為發(fā)射機(jī)系統(tǒng)提供能量。有源鉗位電壓源由功率開(kāi)關(guān)控制的直流恒壓源組成，作用是在發(fā)射電流上升沿和下降沿期間，為電感線圈提供穩(wěn)定的大電壓，使發(fā)射電流線性快速上升或下降。功率逆變?nèi)珮螂娐吠ㄟ^(guò)控制功率開(kāi)關(guān)的導(dǎo)通與截止實(shí)現(xiàn)偽隨機(jī)編碼電流波形的產(chǎn)生。MATLAB軟件為偽隨機(jī)編碼以及各時(shí)序邏輯的產(chǎn)生單元。數(shù)字邏輯控制器FPGA為發(fā)射機(jī)的控制核心，調(diào)用存儲(chǔ)在ROM中的時(shí)序邏輯信號(hào)，從而控制功率開(kāi)關(guān)的導(dǎo)通與截止，實(shí)現(xiàn)偽隨機(jī)編碼電流的發(fā)射。發(fā)射回線為不接地方形回線或者圓形回線。驅(qū)動(dòng)電路由MOSFET驅(qū)動(dòng)芯片搭建而成，其主要作用是將主控電路輸出的控制信號(hào)轉(zhuǎn)換為驅(qū)動(dòng)能力更強(qiáng)的MOSFET控制信號(hào)，并實(shí)現(xiàn)數(shù)字邏輯時(shí)序主控電路與MOSFET之間的電氣隔離。

圖2 瞬變電磁發(fā)射機(jī)整體結(jié)構(gòu)圖Fig.2 Transientelectromagnetic transmitter structure diagram

偽隨機(jī)編碼電流發(fā)射機(jī)主拓?fù)潆娐啡鐖D3所示，技術(shù)核心是偽隨機(jī)編碼電流產(chǎn)生技術(shù)。時(shí)序邏輯如圖4所示，M1～M5為MOSFET功率開(kāi)關(guān)。第一路信號(hào)S1控制M1和M4的導(dǎo)通時(shí)刻；第二路信號(hào)S2控制M2和M3的導(dǎo)通時(shí)刻；第三路信號(hào)S3控制M5的導(dǎo)通時(shí)刻，當(dāng)電流上升沿和下降沿到來(lái)時(shí)，M5開(kāi)始導(dǎo)通，可在電感線圈中產(chǎn)生偽隨機(jī)編碼電流。

圖3 發(fā)射機(jī)主拓?fù)潆娐吩韴DFig.3 Transmittermain topology circuitschematic

圖4 邏輯控制信號(hào)示意圖Fig.4 Schematic diagram of logicalcontrolsignal

2.2 發(fā)射電路分析

在主控電路的時(shí)序控制下，偽隨機(jī)編碼電流產(chǎn)生示意圖如圖5所示。從圖5可知：過(guò)程一電流快速提升，過(guò)程二電流緩慢上升并達(dá)到穩(wěn)態(tài)，過(guò)程三電流泄放，過(guò)程四在泄放即將完成時(shí)，電流反向快速提升。

偽隨機(jī)編碼電流產(chǎn)生過(guò)程一如圖6(a)所示?？刂菩盘?hào)S3觸發(fā)M5導(dǎo)通，鉗位電壓Uc供電，由于有源恒壓鉗位電路的恒定電壓Uc遠(yuǎn)大于電源電壓Us，隔離二極管D1不導(dǎo)通，此時(shí)，鉗位電壓Uc為發(fā)射機(jī)供電，并保持電壓恒定不變。同時(shí)控制信號(hào)S1觸發(fā)M1，M4導(dǎo)通，鉗位電壓Uc通過(guò)M5，M1，RL，L，M4構(gòu)成回路，根據(jù)電感兩端感應(yīng)電壓與電流變化：Uc=L×(dI/dt)，此時(shí)，發(fā)射電流快速線性上升，很快達(dá)到穩(wěn)定值。在此過(guò)程中，M2和M3均處于截止?fàn)顟B(tài)。

圖5 電流形成示意圖Fig.5 Schematic diagram of current formation

偽隨機(jī)編碼電流產(chǎn)生過(guò)程二如圖6(b)所示。當(dāng)控制信號(hào)S3觸發(fā)M5截止，有源恒壓鉗位電路停止放電，此時(shí)M1，M4仍處于導(dǎo)通狀態(tài)，電源Us給電感線圈供電，電流緩慢上升至最大值，進(jìn)入穩(wěn)定電流狀態(tài)。

偽隨機(jī)編碼電流產(chǎn)生過(guò)程三如圖6(c)所示。當(dāng)控制信號(hào)S1觸發(fā)M1，M4截止時(shí)，M5處于截止?fàn)顟B(tài)，在M1和M4關(guān)斷瞬間，發(fā)射線圈作為感性負(fù)載，電流的急劇變化使線圈兩端產(chǎn)生很高的感應(yīng)電壓，使M5的反向并聯(lián)體二極管正向?qū)ā４藭r(shí)，發(fā)射線圈中存儲(chǔ)的能量通過(guò)功率開(kāi)關(guān)M5的反向并聯(lián)體二極管向電容組C1，C2，C3，C4，C5和C6泄放，使電容兩端電壓升高至Uc，發(fā)射電流呈指數(shù)規(guī)律緩慢減小。

偽隨機(jī)編碼電流產(chǎn)生過(guò)程四如圖6(d)所示。在電感線圈向電容組泄放結(jié)束，電感線圈中電流接近于0A，控制信號(hào)S3再次觸發(fā)M5導(dǎo)通，鉗位電壓Uc供電，同時(shí)控制信號(hào)S2觸發(fā)M2和M3導(dǎo)通，Uc通過(guò)M5，M2，RL，L，M3構(gòu)成回路，電流在鉗位電壓下快速提升，方向與之前電流方向相反。此過(guò)程中，M1和M4均處于截止?fàn)顟B(tài)。

3 偽隨機(jī)編碼信號(hào)預(yù)處理

3.1 電磁響應(yīng)正演

瞬變電磁系統(tǒng)線圈中的偽隨機(jī)編碼電流，激發(fā)的電磁場(chǎng)作用于大地，大地產(chǎn)生的二次場(chǎng)與一次場(chǎng)疊加后被接收機(jī)接收。瞬變電磁方法基本原理如圖7所示。

圖6 偽隨機(jī)編碼電流產(chǎn)生過(guò)程示意圖Fig.6 Pseudo-random coded currentgeneration process

發(fā)射電流iT(t)為m序列，將大地視為線性時(shí)不變系統(tǒng)，大地脈沖響應(yīng)為g(t)，輸出響應(yīng)為B(t)。經(jīng)過(guò)收機(jī)系統(tǒng)函數(shù)hr(t)，最終觀測(cè)到感應(yīng)電壓u(t)。根據(jù)圖7可以得到發(fā)射端、接收端及感應(yīng)過(guò)程之間的關(guān)系式。假設(shè)電流記錄儀系統(tǒng)函數(shù)和響應(yīng)端接收機(jī)的記錄儀系統(tǒng)函數(shù)一致，均為hl(t)，發(fā)射端測(cè)得的電流波形i(t)為

通過(guò)發(fā)射電流與大地系統(tǒng)脈沖響應(yīng)的卷積可求得：

圖7 瞬變電磁方法基本原理圖Fig.7 Basic schematic diagram of transient electromagnetic method

接收端測(cè)得的感應(yīng)電壓u(t)滿(mǎn)足：

由式(3)可知：若知道大地脈沖響應(yīng)g(t)，則可得到正演波形。下面以均勻半空間為分析模型，納比吉安[18]給出了圓形回線中心處磁場(chǎng)負(fù)階躍響應(yīng)解析表達(dá)式如式(4)所示，磁場(chǎng)負(fù)階躍響應(yīng)對(duì)時(shí)間求偏導(dǎo)后的解析表達(dá)式如式(5)所示。式(5)即為正演所需的大地脈沖響應(yīng)g(t)。

圖8 均勻大地模型下圓形回線中心z軸瞬變電磁響應(yīng)Fig.8 z-axis transient electromagnetic response of circular return line center under uniform earth model

根據(jù)式(3)，引入高斯噪聲n(t)，正演計(jì)算磁性源圓形回線中心的電磁響應(yīng)波形u(t)，其中噪聲峰值為16.7 nV/m2，均值為0 nV/m2，方差為10-9量級(jí)。噪聲電壓U隨時(shí)間t的變化如圖9(a)所示，正演電壓U隨時(shí)間t的變化如圖9(b)所示。

3.2 響應(yīng)波形預(yù)處理方法

預(yù)處理屬于反演范疇，是指由已知的響應(yīng)波形通過(guò)數(shù)據(jù)處理計(jì)算出大地脈沖響應(yīng)估計(jì)值ge(t)[19-20]。利用m序列的自相關(guān)性，將式(3)兩邊對(duì)電流偏導(dǎo)數(shù)求互相關(guān)，得

其中：Rui為電流偏導(dǎo)數(shù)與響應(yīng)波形的互相關(guān)；Rii為電流偏導(dǎo)數(shù)的自相關(guān)；Rni為電流偏導(dǎo)數(shù)與噪聲的互相關(guān)。由于噪聲n(t)是非相干的，Rni被極大抑制，從而實(shí)現(xiàn)抑制非相干噪聲，得

圖9 噪聲與正演波形Fig.9 Noise and forward waveform

對(duì)式(7)卷積展開(kāi)，并離散化，寫(xiě)成矩陣形式[21]有：

可以通過(guò)求偽逆[22]的方式求出g(t)的離散估計(jì)矩陣ge，可將其視為電壓U。

式中：Fs為采樣頻率。

不含噪聲和含高斯噪聲的預(yù)處理結(jié)果U同理論大地脈沖響應(yīng)g(t)隨時(shí)間t的變化分別如圖10(a)和圖10(b)所示，其中，曲線1為理論大地脈沖響應(yīng)，曲線2為預(yù)處理結(jié)果，噪聲參數(shù)與3.1節(jié)中的相同。從圖10可以看出：引入高斯噪聲后，在晚期預(yù)處理結(jié)果與理論結(jié)果出現(xiàn)偏差，通過(guò)改變?cè)肼暤姆?，發(fā)現(xiàn)噪聲越小，預(yù)處理結(jié)果與理論結(jié)果重合度越高，當(dāng)噪聲為0時(shí)，如圖10(a)曲線2所示。當(dāng)系統(tǒng)在高斯噪聲干擾下，采用該預(yù)處理方法仍可以得到與理論值擬合程度較高的估計(jì)值。

圖10 預(yù)處理結(jié)果與理論值Fig.10 Pretreatment results and theoreticalvalues

4 實(shí)測(cè)結(jié)果

4.1 實(shí)測(cè)發(fā)射電流

為探究感性負(fù)載對(duì)偽隨機(jī)編碼電流的阻礙作用和驗(yàn)證基于已有的功率逆變?nèi)珮螂娐肥欠窨砂l(fā)射偽隨機(jī)編碼電流波形，進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)設(shè)備有：編碼產(chǎn)生與控制器，CASTEM-TX41發(fā)射機(jī)，多匝發(fā)射線圈。線圈參數(shù)為1m×1m，10匝，0.3Ω，0.35mH。發(fā)射的偽隨機(jī)編碼為8階，碼元頻率分別為512，1 024，2 882和4 096 Hz。實(shí)測(cè)數(shù)據(jù)參數(shù)如表2所示，實(shí)測(cè)電流I隨時(shí)間t的變化如圖11(a)～11(d)所示，電流峰值約為4.5A。

表2 實(shí)測(cè)數(shù)據(jù)參數(shù)Table2 Measured data parameters

從圖11可知：設(shè)計(jì)的磁性源偽隨機(jī)編碼發(fā)射機(jī)可以產(chǎn)生質(zhì)量較高的偽隨機(jī)編碼電流脈沖信號(hào)。由于線圈為感性負(fù)載，碼元頻率越高，對(duì)電流上升沿和下降沿的變化速度要求越快；當(dāng)發(fā)射碼元頻率為512 Hz和1 024 Hz時(shí)，偽隨機(jī)編碼波形質(zhì)量高；隨著碼元頻率提升至2 882 Hz和4 096 Hz時(shí)，偽隨機(jī)編碼電流信號(hào)波形出現(xiàn)較小失真。

4.2 實(shí)測(cè)波形正演后的預(yù)處理結(jié)果

將實(shí)測(cè)偽隨機(jī)編碼電流數(shù)據(jù)正演得到相應(yīng)的電磁響應(yīng)波形，然后對(duì)正演波形進(jìn)行預(yù)處理，結(jié)合電流自相關(guān)函數(shù)，進(jìn)一步探究電流波形質(zhì)量對(duì)數(shù)據(jù)預(yù)處理的影響，實(shí)測(cè)自相關(guān)函數(shù)如圖12(a)～(d)所示，其中，A為自相關(guān)幅度，t為時(shí)間。預(yù)處理結(jié)果U隨時(shí)間t的變化如圖13(a)～(d)所示。

從圖12可知：4組發(fā)射電流波形具有較好自相關(guān)性，故實(shí)測(cè)電流波形的失真對(duì)其自相關(guān)性影響較小。從圖13可知：4組發(fā)射電流波形的預(yù)處理結(jié)果均與理論值擬合，證明了磁性源偽隨機(jī)編碼發(fā)射系統(tǒng)盡管隨碼元頻率增高電流波形存在一定程度失真，但電流波形失真對(duì)信號(hào)的相關(guān)性影響很小，仍可以獲得與理論值很接近的大地脈沖響應(yīng)估計(jì)值，預(yù)處理結(jié)果不完全依賴(lài)偽隨機(jī)編碼電流波形。該磁性源瞬變電磁系統(tǒng)可以克服感性負(fù)載的阻礙作用，發(fā)射出波形質(zhì)量高的偽隨機(jī)編碼波形。

圖11 實(shí)測(cè)電流波形Fig.11 Measured currentwaveform

圖12 實(shí)測(cè)數(shù)據(jù)自相關(guān)函數(shù)Fig.12 Measured data autocorrelation

圖13 實(shí)測(cè)數(shù)據(jù)預(yù)處理結(jié)果Fig.13 Measured data preprocessing results

5 結(jié)論

1)實(shí)測(cè)電流的時(shí)域波形、自相關(guān)函數(shù)、預(yù)處理結(jié)果均可證明本文設(shè)計(jì)的磁性源偽隨機(jī)編碼瞬變電磁發(fā)射機(jī)系統(tǒng)可以發(fā)射質(zhì)量很高的偽隨機(jī)編碼電流波形，所采用基于恒壓鉗位方法的發(fā)射機(jī)電路設(shè)計(jì)，可以實(shí)現(xiàn)更高碼元頻率的偽隨機(jī)編碼電流波形發(fā)射。

2)預(yù)處理結(jié)果與理論結(jié)果擬合程度高，表明基于系統(tǒng)辨識(shí)理論預(yù)處理方法是正確的，可進(jìn)一步應(yīng)用于真實(shí)電磁響應(yīng)數(shù)據(jù)的預(yù)處理。

3)實(shí)測(cè)數(shù)據(jù)的預(yù)處理結(jié)果與理論結(jié)果擬合程度高，證明了磁性源偽隨機(jī)編碼發(fā)射系統(tǒng)盡管伴隨碼元頻率增高電流波形存在一定程度失真，但電流波形失真對(duì)信號(hào)的相關(guān)性影響很小，仍可以獲得與理論值很接近的大地脈沖響應(yīng)估計(jì)值，預(yù)處理結(jié)果不完全依賴(lài)偽隨機(jī)編碼電流波形。本文的相關(guān)研究可為磁性源偽隨機(jī)編碼瞬變電磁系統(tǒng)電磁響應(yīng)分析和數(shù)據(jù)解釋提供參考。