一種電磁大地沖激響應(yīng)的辨識(shí)方法
【專利摘要】本發(fā)明公開了一種電磁大地沖激響應(yīng)的辨識(shí)方法，按照Wiener-Hopf方程設(shè)計(jì)電磁大地沖激響應(yīng)辨識(shí)系統(tǒng)，利用偽隨機(jī)序列對人工源電磁方法發(fā)射波形進(jìn)行編碼并發(fā)射，在觀測中對編碼發(fā)射波形與響應(yīng)觀測信號(hào)進(jìn)行同步采集，之后采用基于互相關(guān)辨識(shí)原理的方法，由收發(fā)信號(hào)中消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響，實(shí)現(xiàn)對大地沖激響應(yīng)的高精度辨識(shí)。與傳統(tǒng)階躍源激勵(lì)方式相比，本發(fā)明的方法由收發(fā)互相關(guān)中去除發(fā)射波形自相關(guān)旁瓣影響，從而顯著提高了電磁大地沖激響應(yīng)的辨識(shí)精度。
【專利說明】-種電磁大地沖激響應(yīng)的辨識(shí)方法

【技術(shù)領(lǐng)域】
[0001] 本發(fā)明設(shè)及地球物理勘探領(lǐng)域，更具體地，設(shè)及一種電磁大地沖激響應(yīng)的辨識(shí)方 法，該方法適用于由地球物理電磁勘探系統(tǒng)的觀測數(shù)據(jù)中高精度地辨識(shí)出大地沖激響應(yīng)， 尤其適用于基于編碼發(fā)射電流的地球物理電磁勘探系統(tǒng)。

【背景技術(shù)】
[0002] 隨著科學(xué)技術(shù)地發(fā)展，研究人員面對的對象越發(fā)復(fù)雜，比如化學(xué)化工過程、生物醫(yī) 學(xué)系統(tǒng)、社會(huì)經(jīng)濟(jì)系統(tǒng)、環(huán)境系統(tǒng)等，該些復(fù)雜的對象通常很難由理論分析的方法得到其數(shù) 學(xué)模型。對該類對象的研究就提出了一個(gè)問題；如何對該類對象的數(shù)學(xué)模型與參數(shù)進(jìn)行辨 識(shí)？
[0003] 正是基于上述現(xiàn)實(shí)需求而產(chǎn)生了系統(tǒng)辨識(shí)理論與方法。系統(tǒng)辨識(shí)是現(xiàn)代控制論的 一個(gè)分支，它與狀態(tài)估計(jì)、控制理論構(gòu)成了現(xiàn)代控制論的=大支柱。從控制理論的角度出 發(fā)，系統(tǒng)辨識(shí)的通俗定義是根據(jù)被控對象或被辨識(shí)系統(tǒng)的輸入、輸出觀測信息來估計(jì)它的 數(shù)學(xué)模型。凡是需要通過對系統(tǒng)的輸入輸出進(jìn)行觀測，利用觀測數(shù)據(jù)對系統(tǒng)的數(shù)學(xué)模型與 重要參數(shù)進(jìn)行研究的場合都屬于系統(tǒng)辨識(shí)過程。
[0004] 電磁法作為地球物理勘探方法的一種，依賴外界源（包括人工源與天然源）對大 地媒質(zhì)進(jìn)行激勵(lì)，通過觀測大地媒質(zhì)受到激勵(lì)后產(chǎn)生的響應(yīng)（傳統(tǒng)上稱作二次場）來建立 大地媒質(zhì)電性參數(shù)隨空間分布變化的模型。該個(gè)過程從現(xiàn)代控制論的角度看實(shí)際上就是一 個(gè)系統(tǒng)辨識(shí)的過程，其包括4個(gè)方面；激勵(lì)源、信號(hào)觀測、辨識(shí)模型的建立與系統(tǒng)辨識(shí)方法。
[0005] 對于常見的頻率域方法，其激勵(lì)源包括天然源和人工源（電流波形常見為占空比 為100%的雙極性方波）；信號(hào)觀測包括觀測數(shù)據(jù)的采集與觀測數(shù)據(jù)的處理，觀測數(shù)據(jù)的 采集一般包括正交分布的電場觀測和磁場觀測，觀測數(shù)據(jù)的處理一般包括濾波等預(yù)處理過 程；辨識(shí)模型的建立依據(jù)遠(yuǎn)區(qū)假設(shè)，使發(fā)射源在測點(diǎn)滿足平面波條件，之后按照波阻抗模型 建立映射關(guān)系；系統(tǒng)辨識(shí)方法一般指基于最小二乘法的反演過程。
[0006] 對于常見的時(shí)間域方法，其激勵(lì)源一般為人工源（電流波形常見為占空比非 100%的雙極性方波、雙極性半正弦W及占空比為100%的雙極性梯形波、雙極性S角波）； 信號(hào)觀測包括觀測數(shù)據(jù)的采集與觀測數(shù)據(jù)的處理，觀測數(shù)據(jù)的采集一般包括單分量或多分 量磁場觀測和/或單分量或多分量電場觀測，觀測數(shù)據(jù)的處理一般包括平滑等預(yù)處理過 程；辨識(shí)模型的建立一般分為遠(yuǎn)區(qū)法和近區(qū)法，通過特定的收發(fā)及觀測條件建立相對簡化 的映射關(guān)系；系統(tǒng)辨識(shí)方法一般包括基于最小二乘法的反演過程W及基于簡化算法的成像 過程。
[0007] 無論是頻率域方法還是時(shí)間域方法，就整個(gè)辨識(shí)過程而言實(shí)際都可W被區(qū)分為兩 個(gè)環(huán)節(jié)：第一個(gè)環(huán)節(jié)實(shí)現(xiàn)對電磁大地沖激響應(yīng)的辨識(shí)，第二個(gè)環(huán)節(jié)實(shí)現(xiàn)由電磁大地沖激響 應(yīng)對地下媒質(zhì)電性參數(shù)分布模型的辨識(shí)。第一個(gè)辨識(shí)環(huán)節(jié)完全屬于"黑箱"辨識(shí)，因?yàn)樗?能通過測試系統(tǒng)的響應(yīng)數(shù)據(jù)進(jìn)行辨識(shí)；第二個(gè)辨識(shí)環(huán)節(jié)屬于"灰箱"辨識(shí)，因?yàn)殡姶糯蟮叵?統(tǒng)響應(yīng)與地下媒質(zhì)電性參數(shù)分布之間的映射基本規(guī)律是掌握的（基于麥克斯韋方程組）， 但一些具體的機(jī)理尚未清楚，因此可稱之為"灰箱"辨識(shí)。因?yàn)榈诙€(gè)辨識(shí)環(huán)節(jié)的輸入信號(hào) 是第一個(gè)辨識(shí)環(huán)節(jié)的結(jié)果，因此第二個(gè)環(huán)節(jié)的辨識(shí)精度將很大程度上取決于第一個(gè)環(huán)節(jié)的 辨識(shí)精度。故此，如何提高第一個(gè)環(huán)節(jié)的辨識(shí)精度就對于整個(gè)辨識(shí)的精度具有重大意義。 [000引電磁法的辨識(shí)過程一般不具有實(shí)時(shí)性，按照系統(tǒng)辨識(shí)術(shù)語，屬于離線辨識(shí)。離線辨 識(shí)對于輸入信號(hào)具有一定要求，其最低要求是在整個(gè)觀測周期內(nèi)，系統(tǒng)的所有模態(tài)必須被 輸入信號(hào)持續(xù)激勵(lì)。該就是要求輸入信號(hào)的帶寬要能夠覆蓋辨識(shí)系統(tǒng)的帶寬，輸入信號(hào)帶 寬相對于辨識(shí)系統(tǒng)帶寬足夠大，則對辨識(shí)系統(tǒng)的辨識(shí)也將更加精確。因此，為了能夠更高精 度地實(shí)現(xiàn)系統(tǒng)辨識(shí)，就必須對輸入信號(hào)進(jìn)行設(shè)計(jì)，其包括輸入信號(hào)的類型選擇、幅值與帶寬 等參數(shù)的選擇。從該個(gè)角度觀察常見的人工源電磁法激勵(lì)波形，其在帶寬、幅值等方面均有 欠缺之處。正是因?yàn)槿绱耍絹碓蕉嗟难芯空呖紤]使用W具有最長循環(huán)周期的偽隨機(jī)二進(jìn) 制序列為代表的編碼波形作為激勵(lì)源波形。
[0009] 偽隨機(jī)二進(jìn)制序列（Pseudo Random Binary Sequence,簡稱PRB巧是對多種二進(jìn) 制偽隨機(jī)序列的總稱，所謂二進(jìn)制意味著序列中的每一個(gè)隨機(jī)量只有1和0兩個(gè)邏輯狀態(tài)。 在所有PRBS中，有一類稱作具有最長循環(huán)周期的線性移位寄存器序列（Maximal Length Sequence,簡稱m序列），因具具有近似白噪聲的性質(zhì)，工程上也易于實(shí)現(xiàn)，被廣泛地用作系 統(tǒng)辨識(shí)的輸入信號(hào)。
[0010] 在人工源電磁勘探中使用編碼波形作為發(fā)射電流波形并非新鮮事物。早在二十世 紀(jì)^；：十年代，Quinw、^vewort、Lindsay 1^及Duncan就已經(jīng)展開了基于編碼電流波形的 人工源電磁勘探系統(tǒng)及方法的研究。在國內(nèi)，何繼善于1982年提出了《^^序列偽隨機(jī)信號(hào)電 法，并于21世紀(jì)初對系統(tǒng)進(jìn)行了進(jìn)一步完善。1985年至1986年，中國地質(zhì)大學(xué)電法科研組 在羅延鐘教授的主持下，開展偽隨機(jī)信號(hào)寬帶激電儀的研制工作。上世紀(jì)九十年代，Strack 在其著作中也同樣將偽隨機(jī)編碼作為L0TEM方法發(fā)射波形進(jìn)行了研究，該項(xiàng)研究啟發(fā)蘇格 蘭愛了堡大學(xué)的科研工作者于2004年左右推出了 MTEM系統(tǒng)。在過去十年時(shí)間里，國內(nèi)基 于編碼發(fā)射電流的勘探系統(tǒng)研究也方興未艾，如趙碧如團(tuán)隊(duì)于2007年至2009年研制了適 用于礦產(chǎn)勘查的偽隨機(jī)信號(hào)電阻率法和激電法儀器KGR。2013年啟動(dòng)的國家重大科研裝備 研制項(xiàng)目"深部資源探測核屯、裝備研發(fā)"，其子項(xiàng)目"多通道大功率電法勘探儀"展開基于m 序列編碼電流波形的電法勘探裝置研發(fā)。
[001U 利用m序列進(jìn)行系統(tǒng)辨識(shí)的基本原理基于Wiener-Hopf方程。對于如圖2所示的 一個(gè)線性時(shí)不變系統(tǒng)：
[0012] y (t) = z (t)+n(t) (13)
[0013] z (t) = g(t)*u(t) (14)
[0014] 其中u(t)為輸入信號(hào)，g(t)為系統(tǒng)的沖激響應(yīng)，z(t)為u(t)經(jīng)過系統(tǒng)的輸出信 號(hào)，n(t)為噪聲，y(t)是包含噪聲的輸出信號(hào)。基于圖2的系統(tǒng)描述，首先僅考慮u(t)與 Z (t)，Wiener-Hopf 方程可寫為；
[0015] CR(z，U) = g*AR(u) (15)
[0016] 其中，CR(1，m)表示信號(hào)1與m的互相關(guān)，ARa)表示信號(hào)1的自相關(guān)。式（巧） 表示輸入信號(hào)的自相關(guān)與系統(tǒng)沖激響應(yīng)的卷積為輸入信號(hào)與輸出信號(hào)的互相關(guān)?？紤]噪聲 n，則有：
[0017] CR(y,u) = CR(z, u)+CR(n, u) = g*AR(u)+CR(n, u) (16)
[0018] 式（16)即為基本相關(guān)辨識(shí)方法的數(shù)學(xué)描述。w式（16)為基礎(chǔ)選擇輸入信號(hào)，其 要求包括；首先，輸入信號(hào)應(yīng)具有隨機(jī)性，使式中的CR(n，u)可忽略；其次，輸入信號(hào)自相關(guān) 應(yīng)具有類似5 (t)函數(shù)的形態(tài)，則可使g*AR(u)近似等于g(t)。
[0019] m序列能夠較好地滿足上述要求，因此成為一種常見的系統(tǒng)辨識(shí)輸入信號(hào)。然而， 相對于5 (t)函數(shù)僅在一點(diǎn)上具有非零值的特性，即便m序列選擇較高的編碼階數(shù)、較短的 碼元寬度，其自相關(guān)序列尖峰兩側(cè)的非零旁瓣依然存在。該些旁瓣的影響復(fù)雜且無法忽略， 因此對于辨識(shí)精度有較高要求的EM方法，則不能簡單引用式（16)進(jìn)行大地沖激響應(yīng)辨識(shí)。
[0020] 綜上所述，使用m序列作為激勵(lì)源波形對電磁大地脈沖響應(yīng)進(jìn)行的辨識(shí)，具有一 系列優(yōu)勢，比如高抗干擾性、更寬的帶寬等。然而，對于使用基于m序列作為激勵(lì)源波形的 系統(tǒng)，由于如何克服發(fā)射信號(hào)自相關(guān)旁瓣影響的方法尚不成熟，導(dǎo)致辨識(shí)精度不足。本發(fā)明 正是針對此問題，提出一種去除自相關(guān)旁瓣影響的方法，從而顯著提高了對大地沖激響應(yīng) 的辨識(shí)精度。


【發(fā)明內(nèi)容】

[0021] 有鑒于此，本發(fā)明公開了一種電磁大地沖激響應(yīng)的精確辨識(shí)方法，W從收發(fā)結(jié)果 中更加精確地辨識(shí)出大地電磁脈沖響應(yīng)，
[0022] 為了實(shí)現(xiàn)上述目的，作為本發(fā)明的一個(gè)方面，本發(fā)明提供了一種電磁大地沖激響 應(yīng)的辨識(shí)方法，其特征在于，按照Wiener-Hopf方程設(shè)計(jì)電磁大地沖激響應(yīng)辨識(shí)系統(tǒng)，利用 偽隨機(jī)序列對人工源電磁方法發(fā)射波形進(jìn)行編碼并發(fā)射，在觀測中對編碼發(fā)射波形與響應(yīng) 觀測信號(hào)進(jìn)行同步采集，之后采用基于互相關(guān)辨識(shí)原理的方法，由收發(fā)信號(hào)中消除發(fā)射信 號(hào)自相關(guān)旁瓣的復(fù)雜影響，實(shí)現(xiàn)對大地沖激響應(yīng)的高精度辨識(shí)。
[0023] 作為本發(fā)明的另一個(gè)方面，本發(fā)明還提供了一種電磁大地沖激響應(yīng)的辨識(shí)方法， 包括W下步驟：
[0024] 對人工源電磁方法發(fā)射波形進(jìn)行編碼；
[0025] 按照所述編碼生成發(fā)射驅(qū)動(dòng)信號(hào)；
[0026] W所述發(fā)射驅(qū)動(dòng)信號(hào)驅(qū)動(dòng)發(fā)射機(jī)進(jìn)行發(fā)射；
[0027] 使用具有相同系統(tǒng)響應(yīng)的記錄裝置或各裝置間系統(tǒng)響應(yīng)函數(shù)關(guān)系已知的記錄裝 置同時(shí)對所述發(fā)射機(jī)實(shí)際發(fā)射電流波形與觀測到的響應(yīng)信號(hào)進(jìn)行記錄存儲(chǔ)；
[002引根據(jù)記錄的所述發(fā)射機(jī)實(shí)際發(fā)射電流波形計(jì)算所述發(fā)射機(jī)實(shí)際發(fā)射電流波形的 自相關(guān)函數(shù)；
[0029] 根據(jù)記錄的所述觀測到的響應(yīng)信號(hào)計(jì)算所述觀測到的響應(yīng)信號(hào)與所述發(fā)射機(jī)實(shí) 際發(fā)射電流波形的互相關(guān)函數(shù)；
[0030] 基于互相關(guān)辨識(shí)原理，通過數(shù)學(xué)方法消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響，實(shí)現(xiàn) 對大地沖激響應(yīng)的高精度辨識(shí)。
[0031] 其中，所述對人工源電磁方法發(fā)射波形進(jìn)行編碼的步驟是按照Wiener-Hopf方程 進(jìn)行編碼的。
[0032] 其中，所述對人工源電磁方法發(fā)射波形進(jìn)行編碼的步驟是采用偽隨機(jī)二進(jìn)制序 列，即m序列進(jìn)行編碼的。
[0033] 其中，所述計(jì)算發(fā)射機(jī)實(shí)際發(fā)射電流波形的自相關(guān)函數(shù)的步驟包括：
[0034] 計(jì)算AR(T,(t))，其中T,(t)為實(shí)際發(fā)射電流波形，ARa)表示信號(hào)1的自相關(guān)。
[0035] 其中，所述計(jì)算觀測到的響應(yīng)信號(hào)與發(fā)射機(jī)實(shí)際發(fā)射電流波形的互相關(guān)函數(shù)的步 驟包括：
[0036] 對馬(t)與Tx W做互相關(guān)計(jì)算：
[0037] CR (Rx (t)，Tx (t)) = CR (g (t) *1, (t)曲,(t)，Tx (t)) +CR (V，Tx (t)) (3)
[003引其中，實(shí)際發(fā)射電流波形為T,(t)、觀測到的響應(yīng)信號(hào)為R,(t) W及電磁大地沖激 響應(yīng)為g(t)，I，(t)為發(fā)射機(jī)實(shí)際輸出的發(fā)射電流波形，htf(t)為用于記錄I，(t)的接收機(jī) 的系統(tǒng)響應(yīng)，hf(t)為用于響應(yīng)信號(hào)觀測的接收機(jī)的系統(tǒng)響應(yīng)，V為噪聲，CR(l，m)表示信號(hào) 1與m的互相關(guān)。
[0039] 其中，所述具有相同系統(tǒng)響應(yīng)的記錄裝置或各裝置系統(tǒng)響應(yīng)函數(shù)關(guān)系已知的記錄 裝置滿足W下關(guān)系：
[0040] htr (t) = hr (t)妨
[0041] 或Mt)與hr(t)滿足w下關(guān)系：
[00創(chuàng) htr (t) = f (t) *hr (t)化）
[00創(chuàng)其中，f(t)為已知的、t(t)和hfU)的關(guān)系函數(shù)。
[0044] 其中，所述基于互相關(guān)辨識(shí)原理的方法，消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響的 步驟包括：
[0045]

【權(quán)利要求】
1. 一種電磁大地沖激響應(yīng)的辨識(shí)方法，其特征在于，按照Wiener-Hopf方程設(shè)計(jì)電磁 大地沖激響應(yīng)辨識(shí)系統(tǒng)，利用偽隨機(jī)序列對人工源電磁方法發(fā)射波形進(jìn)行編碼并發(fā)射，在 觀測中對編碼發(fā)射波形與響應(yīng)觀測信號(hào)進(jìn)行同步采集，之后采用基于互相關(guān)辨識(shí)原理的方 法，由收發(fā)信號(hào)中消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響，實(shí)現(xiàn)對大地沖激響應(yīng)的高精度辨 識(shí)。
2. -種電磁大地沖激響應(yīng)的辨識(shí)方法，包括以下步驟： 對人工源電磁方法發(fā)射波形進(jìn)行編碼； 按照所述編碼生成發(fā)射驅(qū)動(dòng)信號(hào)； 以所述發(fā)射驅(qū)動(dòng)信號(hào)驅(qū)動(dòng)發(fā)射機(jī)進(jìn)行發(fā)射； 使用具有相同系統(tǒng)響應(yīng)的記錄裝置或各裝置間系統(tǒng)響應(yīng)函數(shù)關(guān)系已知的記錄裝置同 時(shí)對所述發(fā)射機(jī)實(shí)際發(fā)射電流波形與觀測到的響應(yīng)信號(hào)進(jìn)行記錄存儲(chǔ)； 根據(jù)記錄的所述發(fā)射機(jī)實(shí)際發(fā)射電流波形計(jì)算所述發(fā)射機(jī)實(shí)際發(fā)射電流波形的自相 關(guān)函數(shù)； 根據(jù)記錄的所述觀測到的響應(yīng)信號(hào)計(jì)算所述觀測到的響應(yīng)信號(hào)與所述發(fā)射機(jī)實(shí)際發(fā) 射電流波形的互相關(guān)函數(shù)； 基于互相關(guān)辨識(shí)原理，通過數(shù)學(xué)方法消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響，實(shí)現(xiàn)對大 地沖激響應(yīng)的高精度辨識(shí)。
3. 根據(jù)權(quán)利要求2所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中所述對人工源電磁方法 發(fā)射波形進(jìn)行編碼的步驟是按照Wiener-Hopf?方程進(jìn)行編碼的。
4. 根據(jù)權(quán)利要求3所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中所述對人工源電磁方法 發(fā)射波形進(jìn)行編碼的步驟是采用偽隨機(jī)二進(jìn)制序列，即m序列進(jìn)行編碼的。
5. 根據(jù)權(quán)利要求2所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中所述計(jì)算發(fā)射機(jī)實(shí)際發(fā) 射電流波形的自相關(guān)函數(shù)的步驟包括： 計(jì)算AR(Tx(t))，其中Tx(t)為實(shí)際發(fā)射電流波形，AR(I)表示信號(hào)1的自相關(guān)。
6. 根據(jù)權(quán)利要求2所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中所述計(jì)算觀測到的響應(yīng) 信號(hào)與發(fā)射機(jī)實(shí)際發(fā)射電流波形的互相關(guān)函數(shù)的步驟包括： 對艮(0與1；(〇做互相關(guān)計(jì)算： CR(Rx (t)，Tx (t)) =CR(g(t) *IW (t) *hr (t)，Tx (t))+CR(v，Tx (t)) (3) 其中，實(shí)際發(fā)射電流波形為Tx(t)、觀測到的響應(yīng)信號(hào)為Rx(t)以及電磁大地沖激響應(yīng) 為g(t)，Iw(t)為發(fā)射機(jī)實(shí)際輸出的發(fā)射電流波形，ht,(t)為用于記錄Iw(t)的接收機(jī)的系 統(tǒng)響應(yīng)，Mt)為用于響應(yīng)信號(hào)觀測的接收機(jī)的系統(tǒng)響應(yīng)，V為噪聲，CR(l，m)表示信號(hào)1與 m的互相關(guān)。
7. 根據(jù)權(quán)利要求2所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中所述具有相同系統(tǒng)響應(yīng) 的記錄裝置或各裝置系統(tǒng)響應(yīng)函數(shù)關(guān)系已知的記錄裝置滿足以下關(guān)系： htr(t) =hr(t) (5) 或Mt)與^⑴滿足以下關(guān)系：hr(t) =f(t)*hr(t) (6) 其中，f(t)為已知的ht,(t)和h,(t)的關(guān)系函數(shù)。
8. 根據(jù)權(quán)利要求7所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中所述基于互相關(guān)辨識(shí)原 理的方法，通過數(shù)學(xué)方法消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響的步驟包括：
且當(dāng)htt(t)與匕⑴之間的關(guān)系如權(quán)利要求7所述時(shí)， 則式（4)可被矩陣化為： A=BG+cv (11) 其中：
Ng為辨識(shí)大地沖激響應(yīng)的采樣點(diǎn)數(shù)；nJPn2分別為a(n)和b(n)序列中最大值采樣點(diǎn) 的序列號(hào)； 通過最小二乘法計(jì)算，將矢量G從矢量A中分離出來： G= (BtBK1BtA(12)， 由此得到的G實(shí)現(xiàn)了將發(fā)射信號(hào)自相關(guān)旁瓣的影響從收發(fā)互相關(guān)中去除，實(shí)現(xiàn)了對大 地沖激響應(yīng)的高精度辨識(shí)。
9. 根據(jù)權(quán)利要求8所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中在所述基于互相關(guān)辨識(shí) 原理的方法，消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響的步驟中，對所述式（12)進(jìn)行多次迭代 計(jì)算，以進(jìn)一步提高辨識(shí)精度。
10. 根據(jù)權(quán)利要求8所述的電磁大地沖激響應(yīng)的辨識(shí)方法，其中在所述基于互相關(guān)辨 識(shí)原理，消除發(fā)射信號(hào)自相關(guān)旁瓣的復(fù)雜影響的步驟之前還包括對數(shù)據(jù)進(jìn)行去偏置處理， 即使用低通濾波將信號(hào)帶寬限制在相應(yīng)的范圍內(nèi)的步驟；以及 在完成辨識(shí)之后，還包括根據(jù)接收機(jī)采樣頻率對G的幅度進(jìn)行修正，并對G曲線進(jìn)行平 滑處理的步驟。
【文檔編號(hào)】G01V3/08GK104502980SQ201410743525
【公開日】2015年4月8日 申請日期:2014年12月8日 優(yōu)先權(quán)日:2014年12月8日
【發(fā)明者】武欣, 方廣有, 薛國強(qiáng), 底青云, 張一鳴 申請人:中國科學(xué)院電子學(xué)研究所, 中國科學(xué)院地質(zhì)與地球物理研究所, 北京工業(yè)大學(xué)