基于重構(gòu)激發(fā)信號的電磁式可控震源相關(guān)檢測參考信號方法研究

楊志超， 姜弢*， 徐學(xué)純， 林君， 賈海青， 馬瑞琪， 陳振振

1 吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院， 長春　130061 2 地球信息探測儀器教育部重點(diǎn)實(shí)驗(yàn)室(吉林大學(xué))， 長春　130061 3 吉林大學(xué)地球科學(xué)學(xué)院， 長春　130061

?

基于重構(gòu)激發(fā)信號的電磁式可控震源相關(guān)檢測參考信號方法研究

楊志超1,2， 姜弢1,2*， 徐學(xué)純3， 林君1,2， 賈海青1,2， 馬瑞琪1,2， 陳振振1,2

1 吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院， 長春130061 2 地球信息探測儀器教育部重點(diǎn)實(shí)驗(yàn)室(吉林大學(xué))， 長春130061 3 吉林大學(xué)地球科學(xué)學(xué)院， 長春130061

摘要針對電磁式可控震源地震數(shù)據(jù)的相關(guān)檢測，研究發(fā)現(xiàn)，在地下結(jié)構(gòu)復(fù)雜、基板-大地耦合不佳時(shí)，常規(guī)方法——基于震源控制信號或基板附近信號作為參考信號檢測得到的地震記錄中，存在子波到時(shí)誤差和虛假多次波問題.本文分析了上述問題的理論原因，并提出基于重構(gòu)激發(fā)信號的相關(guān)檢測參考信號方法(Correlation Detection Reference Signal Based on the Reconstructed Excitation Signal，CDRSBRES).首先，利用直達(dá)波與其他地震波到時(shí)不一致的特點(diǎn)，從震源基板附近信號中分離、提取直達(dá)波.然后，利用直達(dá)波重構(gòu)震源激發(fā)信號并作為參考信號對地震數(shù)據(jù)進(jìn)行相關(guān)檢測.最后，應(yīng)用譜白化技術(shù)提高檢測結(jié)果質(zhì)量.數(shù)值模擬研究表明，重構(gòu)激發(fā)信號與理想激發(fā)信號的相關(guān)系數(shù)為0.9869，達(dá)到高度線性相關(guān)，CDRSBRES方法檢測的地震記錄在子波到時(shí)和波形特征上均與模型相符.隨后，在某金屬礦區(qū)開展了可控震源對比實(shí)驗(yàn).與液壓式可控震源MiniVib T15000檢測結(jié)果相比，電磁式可控震源PHVS 500的檢測結(jié)果中：基于震源控制信號的檢測結(jié)果存在子波到時(shí)誤差約0.012 s，對應(yīng)垂向精度誤差約11.16 m；基于基板附近信號的檢測結(jié)果部分區(qū)域出現(xiàn)虛假多次波，信噪比降低；而CDRSBRES方法的檢測結(jié)果子波到時(shí)誤差約0.001 s，對應(yīng)垂向精度誤差約0.93 m，波形特征一致，相同區(qū)域無虛假多次波.綜上，本方法適用于電磁式可控震源地震數(shù)據(jù)的高精度檢測，尤其對于地下結(jié)構(gòu)復(fù)雜區(qū)域的高分辨率地震勘探具有重要意義.

關(guān)鍵詞電磁式可控震源； 相關(guān)檢測； 參考信號； 激發(fā)信號； 重構(gòu)

This paper analyzes the errors above systematically and then proposes a method called correlation detection reference signal based on the reconstructed excitation signal (CDRSBRES). Firstly, the direct wave is separated and extracted from the near base-plate signal according to different arrival time among the direct wave and other seismic waves. Secondly, applying the extracted direct wave to reconstruct the vibrator excitation signal as the reference signal to detect the seismic data. Finally, the spectral whitening is used to improve the quality of detected result, such as the SNR and resolution.

This paper designs a horizontal layer model to confirm the CDRSBRES method. The simulation data demonstrates that the correlation coefficient between the reconstructed excitation signal and the ideal excitation signal is 0.9869, up to high linear correlation. The correlation detected record by CDRSBRES method is consistent with the designed model not only on wavelet arrival time but also on waveform characters. Then a comparison experiment is carried out at one metal deposit in West China. It uses the correlation detected record of the hydraulic seismic vibrator MiniVib T15000 as a standard. Compared to MiniVib T15000, there are different results for the electromagnetic seismic vibrator PHVS 500. The 0.012 s error of wavelet arrival time is detected based on the control-signal method, corresponding error of vertical position is about 11.16 m. The false multiple waves are detected based on the near base-plate signal method and it drops the SNR of detected record. But the error of wavelet arrival time is only 0.001 s, corresponding error on vertical position is down to 0.93 m. The waveform characteristic is consistent, no false multiple waves detected at the same zone using CDRSBRES method.

In summary, the CDRSBRES method can be used for high precision detection of electromagnetic seismic vibrator data, especially for the region where the subsurface media is complex. It is of great significance for high resolution seismic exploration.

1引言

根據(jù)可控震源地震勘探原理(林君，2004；佟訓(xùn)乾等，2012；趙春蕾，2013)，震源激發(fā)信號(Ziolkowski，2010；Krohn et al.，2010)應(yīng)作為理想的參考信號與地震數(shù)據(jù)進(jìn)行相關(guān)檢測(姜弢等，2005)，得到類似脈沖震源的地震子波.但由于震源結(jié)構(gòu)和工作模式等限制，使得震源激發(fā)信號難以直接測量.Castenet和Lavergne(1965)提出“加權(quán)和”方法，即將基板加速度與質(zhì)量乘積和反應(yīng)體加速度與質(zhì)量乘積進(jìn)行加權(quán)，估計(jì)震源激發(fā)力Ground Force.該方法假設(shè)震源基板是剛性的，且忽略了基板的彎曲力(Baeten and Ziolkowski，1990).在此基礎(chǔ)上，Sallas(1984)繼續(xù)驗(yàn)證了方法的有效性，并提出應(yīng)用鎖相等技術(shù)反饋控制 “加權(quán)和”方法估計(jì)結(jié)果，使激發(fā)信號更加逼近震源控制信號.進(jìn)一步地，Wei(2010)和Wei等(2010)完善了基板與大地耦合模型，研究了基板材料的剛性程度和大地耦合程度等因素對“加權(quán)和”方法估計(jì)結(jié)果的影響，主要是諧波干擾影響.目前，“加權(quán)和”方法及其基礎(chǔ)上的研究成為估計(jì)震源激發(fā)信號的主要方法，但其模型和實(shí)驗(yàn)主要基于液壓式可控震源(Boucard and Ollivrin,2010；Sallas，2010).電磁式可控震源(林君，2004；Haines,2006)研究的起步時(shí)間較晚，但擁有更寬的信號頻帶范圍，在高分辨率地震勘探領(lǐng)域具有廣泛的應(yīng)用前景(葛利華等，2012；王俊秋等，2012).“加權(quán)和”方法對于電磁式可控震源是否適用？Veen等(1999)針對OYO輕量電磁式可控震源，將“加權(quán)和”估計(jì)結(jié)果與放置在基板上的力傳感器測量結(jié)果進(jìn)行對比，指出該方法對激發(fā)信號高頻段的估計(jì)存在誤差.同時(shí)，通過對便攜式高頻可控震源系統(tǒng)(Portable High-frequency Vibrator System，PHVS，吉林大學(xué)自主研發(fā)的電磁式可控震源)(陳祖斌等，2003；林 君，2004；梁鐵成等，2005)的測試研究發(fā)現(xiàn)，該方法的另一個(gè)限制在于對可控震源硬件要求較高，需要在反應(yīng)體和基板上增加加速度傳感器.因此，針對具有5～1000 Hz頻率應(yīng)用范圍的電磁式可控震源，需要尋找更適宜的方法估計(jì)激發(fā)信號，作為相關(guān)檢測的參考信號.在當(dāng)前的電磁式可控震源應(yīng)用中，常規(guī)方法之一是參照液壓式可控震源，基于震源控制信號(Nigel，1991)作為參考信號.該方法應(yīng)用簡單，但對震源質(zhì)量控制(QC)及基板-大地耦合條件要求較高.此外，陳祖斌等(2005)通過對電磁式可控震源基板-大地耦合過程的建模分析，提出在震源基板附近放置一道檢波器記錄震動(dòng)信號，基于基板附近信號作為參考信號.該方法的優(yōu)勢在于考慮了基板-大地耦合情況對震源控制信號產(chǎn)生的影響.

但是，通過分析電磁式可控震源PHVS 500與液壓式可控震源MiniVib T15000的對比實(shí)驗(yàn)結(jié)果發(fā)現(xiàn)，在地下結(jié)構(gòu)復(fù)雜、基板-大地耦合不佳的情況下，PHVS 500常規(guī)方法，即基于震源控制信號和基板附近信號作為參考信號相關(guān)檢測得到的地震記錄均存在問題.因此，本文提出基于重構(gòu)激發(fā)信號的相關(guān)檢測參考信號方法(Correlation Detection Reference Signal Based on the Reconstructed Excitation Signal，CDRSBRES).通過數(shù)值模擬和實(shí)測數(shù)據(jù)研究結(jié)果表明，本方法不僅可有效重構(gòu)震源激發(fā)信號，而且較常規(guī)相關(guān)檢測方法更好地解決了復(fù)雜地質(zhì)條件下可控震源地震數(shù)據(jù)的高精度檢測問題.

2問題描述

2010年在中國西部某金屬礦區(qū)開展了電磁式可控震源PHVS 500與液壓式可控震源MiniVib T15000的對比實(shí)驗(yàn).該工區(qū)地下結(jié)構(gòu)復(fù)雜：銅鎳礦床規(guī)模小，層位較淺，連續(xù)性差(王俊秋等，2011).地表為干枯河道，起伏不平，基板-大地耦合效果不佳.兩震源工作參數(shù)一致：控制信號為chirp信號(林君，2004)，頻帶范圍為10～120 Hz，掃描時(shí)間為10 s.但輸出力差距較大，MiniVib T15000輸出力為30000 N，而PHVS 500輸出力為300 N.兩震源共用一套觀測系統(tǒng)：無纜自定位地震儀(董樹文等，2012)采集，偏移距為16 m，192道檢波器全排列接收，道間距為2 m，采樣率為1000 Hz，采集時(shí)間為12 s.

圖1對比展示了本次實(shí)驗(yàn)同一震源點(diǎn)不同震源的地震記錄.由于MiniVib T15000在生產(chǎn)中已廣泛應(yīng)用，現(xiàn)階段其數(shù)據(jù)可視為對比參照標(biāo)準(zhǔn).與MiniVib T15000地震記錄圖1a相比，PHVS 500存在不一致問題：基于震源控制信號作為參考信號檢測得到的地震記錄圖1b中，直觀上與圖1a的波形形態(tài)較相似(標(biāo)記實(shí)線框區(qū)域).但精確對比發(fā)現(xiàn)，圖1b與圖1a的子波到時(shí)存在偏差，計(jì)算得知，對于箭頭指示的同一子波，兩者偏差約0.012 s；基于基板附近信號作為參考信號檢測得到的地震記錄圖1c中，標(biāo)記實(shí)線框區(qū)域與圖1a相比較，子波的波形形態(tài)出現(xiàn)畸變，產(chǎn)生了虛假多次波信號.

圖1　西部某金屬礦區(qū)可控震源地震記錄對比 (a) MiniVib T15000地震記錄； (b) PHVS 500地震記錄——基于震源控制信號作為參考信號；(c) PHVS 500地震記錄——基于基板附近信號作為參考信號.

根據(jù)可控震源地震勘探原理，相關(guān)檢測方法在理論和實(shí)踐上都是有效的，上述問題在地下結(jié)構(gòu)簡單、震源基板-大地耦合良好的情況下也很少出現(xiàn).為此，需要首先分析出現(xiàn)上述問題理論上的原因.

3問題分析

3.1震源激發(fā)信號

電磁式可控震源由信號發(fā)生器產(chǎn)生控制信號s(t)，控制震源激發(fā)信號掃描形式.控制信號經(jīng)過電控箱系統(tǒng)，傳遞到電磁激振器系統(tǒng)，經(jīng)基板-大地耦合系統(tǒng)向地下輸出，最終形成震源激發(fā)信號m(t)，過程如圖2所示.其中，電控箱的電子系統(tǒng)和電磁激振器的機(jī)械系統(tǒng)均可視為是線性系統(tǒng)，信號經(jīng)過此過程只產(chǎn)生線性變化.而基板-大地耦合系統(tǒng)則為非線性系統(tǒng)(Lebedev and Beresnev，2004； Saragiotis et al.,2010)，對信號的幅度及相位均產(chǎn)生非線性變化(陶知非，2007; Boucard and Ollivrin，2010).因此，震源激發(fā)信號m(t)實(shí)際上是控制信號s(t)經(jīng)過上述線性系統(tǒng)和非線性系統(tǒng)后產(chǎn)生一定畸變的結(jié)果，可表示為

(1)

公式(1)中K為幅度畸變，t1為相位畸變等效的時(shí)間延遲.激發(fā)信號經(jīng)地下介質(zhì)作用后由檢波器接收形成地震信號.設(shè)一道檢波器信號為y(t)，其中的1個(gè)反射波為f(t)，可表示為

(2)

式中B為地下介質(zhì)作用產(chǎn)生的幅度畸變，t2為該反射波的到時(shí).針對該反射波f(t)，基于震源激發(fā)信號作為參考信號，則相關(guān)檢測結(jié)果為

(3)

當(dāng)τ=t2時(shí)，rfm(τ)達(dá)到最大值，即對于f(t)的到時(shí)檢測為t2，與準(zhǔn)確值相符.

但在實(shí)際應(yīng)用中，由于震源激振器質(zhì)量較大，且需要基板與大地盡可能緊密耦合，才能實(shí)現(xiàn)震源能量有效傳入地下，因此無法在基板與大地之間放置傳感器直接測量激發(fā)信號.目前針對電磁式可控震源，常規(guī)方法基于震源控制信號或基板附近信號作為相關(guān)檢測的參考信號.

3.2震源控制信號

震源控制信號為信號發(fā)生器產(chǎn)生的標(biāo)準(zhǔn)掃描信號，但通過對圖2的分析，由于基板-大地耦合系統(tǒng)的非線性作用(Wei and Phillips，2012)，尤其當(dāng)耦合效果不佳時(shí)，震源激發(fā)信號與控制信號相比，幅度與相位均產(chǎn)生較大畸變.

針對3.1節(jié)中的反射波信號f(t)，基于震源控制信號s(t)作為參考信號，則相關(guān)檢測結(jié)果為

(4)

式中當(dāng)τ=t1+t2時(shí)，rfs(τ)達(dá)到最大值.因此，基于控制信號s(t)作為參考信號時(shí)，對于反射波f(t)的到時(shí)檢測結(jié)果為t1+t2，與準(zhǔn)確值t2不符，使到時(shí)檢測出現(xiàn)誤差.

3.3基板附近信號

如圖2所示，震源激發(fā)信號經(jīng)地表作用進(jìn)入檢波器形成直達(dá)波.同時(shí)，由地震波傳播原理，震源激發(fā)信號也會(huì)向地下傳播，經(jīng)地下各層介質(zhì)反射、透射、折射等作用進(jìn)入檢波器，形成反射波和干擾波，此處干擾波的含義包括折射波、面波等以及環(huán)境噪聲等(姜弢等，2008).因此，每一道檢波器采集的信號同時(shí)包含直達(dá)波、反射波和干擾波，是震源激發(fā)信號經(jīng)過一系列系統(tǒng)作用后結(jié)果的累加.基板附近信號與其他檢波器信號并無差別，只是位置上與震源更近，直達(dá)波能量更強(qiáng).

(5)

由公式(5)，基板附近信號為

(6)

基于基板附近信號y1(t)作為參考信號，則第i道檢波器記錄進(jìn)行相關(guān)檢測之后的結(jié)果為

(7)

(8)

設(shè)di(t)=ri,0(t)?r1,0(t)，

(9)

公式(7)—(8)中，?代表相關(guān).公式(9)中，di(t)為直達(dá)波，Ri(t)為相關(guān)檢測出的n個(gè)反射波，而mi(t)為實(shí)際不存在但卻被檢測出的虛假多次波信號.實(shí)測數(shù)據(jù)中，由于地下介質(zhì)對地震波的吸收和衰減(葛利華等，2012)，深層反射波能量很弱，不經(jīng)過自動(dòng)增益控制(姜弢等，2012)，可識別的虛假多次波雖不明顯，但仍會(huì)對后續(xù)的處理和解釋造成困難.

4CDRSBRES原理

通過上述實(shí)驗(yàn)與理論分析，在地下結(jié)構(gòu)復(fù)雜、基板-大地耦合不佳時(shí)，基于震源控制信號或基板附近信號作為相關(guān)檢測參考信號的方法存在到時(shí)檢測誤差和虛假多次波問題.因此，考慮重新構(gòu)建震源激發(fā)信號.結(jié)合圖2分析基板附近檢波器，實(shí)際工作中，該檢波器位置與震源基板距離極近，通常小于1 m.震源激發(fā)信號m(t)首先在地表傳播了極短距離后即被其接收，形成為直達(dá)波g(t).若忽略地表對激發(fā)信號在此極短距離內(nèi)的作用，則直達(dá)波g(t)的波形與激發(fā)信號可等效為一致.因此，在無法直接測量震源激發(fā)信號的情況下，考慮利用直達(dá)波與其他地震波到時(shí)不一致的特點(diǎn)，通過合理設(shè)計(jì)匹配濾波器，從基板附近信號中有效分離、提取直達(dá)波來重構(gòu)激發(fā)信號，將更有效地實(shí)現(xiàn)電磁式可控震源地震數(shù)據(jù)的高精度檢測，這就是本文提出的重構(gòu)激發(fā)信號原理，簡要示意如圖3所示.

圖3　基板附近信號中提取直達(dá)波過程

5CDRSBRES方法

5.1重構(gòu)震源激發(fā)信號

(1)由于震源激發(fā)信號為調(diào)頻信號，基板附近信號中的直達(dá)波、各層反射波、干擾波等信號雖然到時(shí)不同，但由于持續(xù)時(shí)間長，在時(shí)域內(nèi)混合在一起，無法直接提取直達(dá)波.因此，首先構(gòu)建控制信號的匹配濾波器.

(10)

公式(10)中，A為信號幅度，f0為信號起始頻率，f1為信號終止頻率，T為掃描時(shí)間.

設(shè)計(jì)控制信號s(t)的匹配濾波器(段云卿，2006)，濾波器長度與s(t)相同，延時(shí)為0，其沖激響應(yīng)h(t)為

(11)

即控制信號s(t)的匹配濾波器沖激響應(yīng)h(t)是原信號的共軛映像.

(2)使基板附近信號y1(t)通過該匹配濾波器，結(jié)果為

(12)

公式(12)中，*代表卷積.從卷積結(jié)果v1(t)中，選取適當(dāng)長度窗口，分離、提取初至信號d1(t)(圖3中，右側(cè)矩形框內(nèi)兩條黑線標(biāo)記部分)，作為直達(dá)波.其中，窗口長度的選擇要求保證盡量完整地分離初至信號，同時(shí)不包含反射波、折射波等強(qiáng)干擾信號.

此時(shí)分離的直達(dá)波仍為地震子波形式，尚需變換為直達(dá)波本來的調(diào)頻形式g(t).

(3) 由圖3可得

(13)

根據(jù)卷積定理，兩信號在時(shí)域的卷積等效于兩信號在頻域的乘積，因此，對公式(13)兩邊進(jìn)行傅里葉變換，得到

(14)

通過傅里葉反變換，得到

(15)

公式(14)—(15)中，F(xiàn)和F-1代表傅里葉變換及其反變換.由于h(t)已知，d1(t)已由基板附近信號中分離、提取得到，所以可求得直達(dá)波g(t)，即對震源激發(fā)信號的估計(jì)，稱為重構(gòu)震源激發(fā)信號.

5.2相關(guān)檢測

基于重構(gòu)震源激發(fā)信號作為參考信號，與地震記錄進(jìn)行相關(guān)檢測.其過程與常規(guī)相關(guān)檢測一致.

5.3譜白化處理

由于震源激發(fā)信號相對于控制信號產(chǎn)生了一定畸變，在相關(guān)檢測中會(huì)導(dǎo)致檢測結(jié)果的子波相關(guān)噪聲較強(qiáng)，即旁瓣干擾明顯.同理，基于重構(gòu)激發(fā)信號作為參考信號依然存在此問題，且與激發(fā)信號的相似程度越大，相關(guān)噪聲越明顯.為壓制相關(guān)噪聲，對相關(guān)檢測后的結(jié)果再進(jìn)行譜白化處理(邊國柱和張立群，1986)，補(bǔ)償信號的高頻部分能量，其過程可由下列步驟實(shí)現(xiàn).

(1) 針對一道信號，在其有效頻帶范圍內(nèi)，進(jìn)行若干窄帶通分頻濾波.其中，濾波器個(gè)數(shù)及頻帶范圍根據(jù)信號頻譜分析確定.

(2) 對每個(gè)窄帶通分頻濾波結(jié)果進(jìn)行時(shí)變增益.首先將信號分成若干時(shí)窗，分別求出各自時(shí)窗內(nèi)振幅均方根值，再用時(shí)窗內(nèi)各值除以該均方根值，最后乘以常數(shù)因子.所分時(shí)窗長度根據(jù)記錄中各層有效信號的到時(shí)確定.

(3) 將所有分頻濾波結(jié)果進(jìn)行時(shí)變增益后，在時(shí)域內(nèi)疊加，則完成一道信號的譜白化處理.

(4) 最后，對每一道信號都進(jìn)行上述分頻濾波、時(shí)變增益、疊加的操作，則完成整個(gè)地震記錄的譜白化處理.

綜上，基于重構(gòu)激發(fā)信號的相關(guān)檢測參考信號方法CDRSBRES包括：重構(gòu)震源激發(fā)信號，與地震記錄進(jìn)行相關(guān)檢測和譜白化處理.

6數(shù)值模擬與分析

數(shù)值模擬中，采用5層水平層狀地下模型，其形態(tài)及其各層深度、速度如圖4所示.激發(fā)端，震源控制信號為線性掃頻的chirp信號，掃描頻帶為10～500 Hz，掃描時(shí)間為4 s，波形如圖5a所示.基板-大地耦合系統(tǒng)可看作帶通濾波器(林君，2004)，因此理想激發(fā)信號可認(rèn)為是控制信號通過非線性帶通濾波器的結(jié)果.其中，采用凱澤窗函數(shù)的FIR濾波器，通頻帶為20～250 Hz，相位滯后等效產(chǎn)生的延時(shí)為0.02 s，輸出結(jié)果如圖5b所示.接收端，基板附近檢波器與震源點(diǎn)位置為1 m，布置192道檢波器全排列接收，偏移距為4 m，檢波器道間距為4 m，記錄采樣率為4000 Hz，采集時(shí)間為6 s.這里僅考慮檢波器接收直達(dá)波、反射波和隨機(jī)噪聲，忽略折射波等影響，采用射線追蹤方法(姜弢等，2006)合成地震記錄.

圖4　5層水平層狀地下模型

圖5　數(shù)值模擬中的信號(a)震源控制信號；(b)理想激發(fā)信號.

6.1理想激發(fā)信號作為參考信號

首先針對由圖4模型模擬生成的地震記錄，基于圖5b所示的理想激發(fā)信號作為參考信號進(jìn)行相關(guān)檢測，得到理論上正確的檢測結(jié)果，如圖6a所示，作為對比其他方法檢測結(jié)果的參照標(biāo)準(zhǔn).

6.2震源控制信號作為參考信號

針對由圖4模型模擬生成的地震記錄，基于震源控制信號作為參考信號進(jìn)行相關(guān)檢測，結(jié)果如圖6b所示.總體上比較，波形形態(tài)與圖5a一致.分析到時(shí)檢測精度，同時(shí)抽取圖6b，6a的第1道記錄，截取0～0.7 s在圖7中對比展示.其中，準(zhǔn)確到時(shí)由豎線標(biāo)出.因此，在地下結(jié)構(gòu)復(fù)雜、基板-大地耦合不佳時(shí)，基于震源控制信號作為參考信號的檢測結(jié)果雖然波形特征與模型一致，但對子波到時(shí)的檢測存在誤差，易對地下目標(biāo)體深度的解釋造成困難.

6.3基板附近信號作為參考信號

針對由圖4模型模擬生成的地震記錄，采用基板附近信號作為參考信號進(jìn)行相關(guān)檢測，結(jié)果如圖6c所示，其中初至之前信號被切除.與圖6a相比，直達(dá)波下方(標(biāo)記實(shí)線框區(qū)域)檢測到重復(fù)出現(xiàn)的相似信號，即虛假多次波信號，嚴(yán)重干擾了反射波信號.產(chǎn)生虛假多次波信號的原因可通過3.3節(jié)中的公式推導(dǎo)得出，即地震記錄與基板附近信號中的各層反射波重復(fù)進(jìn)行了相關(guān)檢測.虛假多次波信號成組出現(xiàn)，出現(xiàn)的個(gè)數(shù)等于基板附近信號中的反射波個(gè)數(shù)，出現(xiàn)的間隔等于各反射波的到時(shí)差.尤其在近偏移距端或地下層位較淺時(shí)，虛假多次波信號出現(xiàn)的間隔較短，造成的干擾更為明顯.而在遠(yuǎn)偏移距端或地下層位較深時(shí)，各反射波到時(shí)差較長，虛假多次波信號造成的干擾降低(標(biāo)記虛線框區(qū)域).

因此，基于基板附近信號作為參考信號，雖然考慮到基板-大地耦合系統(tǒng)對控制信號產(chǎn)生的非線性變化，但在地下結(jié)構(gòu)復(fù)雜、基板-大地耦合不佳時(shí)，將使檢測結(jié)果中出現(xiàn)虛假多次波，降低了地震記錄的信噪比(姜弢等，2009)，易使處理結(jié)果出現(xiàn)虛假層位.

6.4CDRSBRES方法

針對由圖4模型模擬生成的地震記錄，應(yīng)用CDRSBRES方法進(jìn)行檢測，過程及結(jié)果如下.

6.4.1重構(gòu)震源激發(fā)信號

首先，由公式(10)—(15)，從基板附近信號中分離、提取直達(dá)波以重構(gòu)震源激發(fā)信號.其中，針對匹配濾波結(jié)果，在初至信號峰值點(diǎn)處選取0.03 s的窗口分離初至，作為直達(dá)波信號，經(jīng)過頻域與時(shí)域處理，其波形如圖8所示.

圖6　數(shù)值模擬可控震源地震記錄(a) 基于理想激發(fā)信號作為參考信號； (b) 基于震源控制信號作為參考信號； (c) 基于基板附近信號作為參考信號.

圖7　地震信號單道記錄波形(豎線標(biāo)明準(zhǔn)確到時(shí))(a) 圖6b的第1道信號(0～0.7 s)；(b) 圖6a的第1道信號(0～0.7 s).

通過計(jì)算，圖8所示重構(gòu)激發(fā)信號與圖5b所示理想震源激發(fā)信號的相關(guān)系數(shù)(崔子健等，2012)為0.9869，相似程度為高度線性相關(guān)，說明該方法可以實(shí)現(xiàn)對震源激發(fā)信號的有效重構(gòu).

6.4.2相關(guān)檢測

基于重構(gòu)激發(fā)信號作為參考信號，與圖4模擬得到的地震記錄相關(guān)檢測，結(jié)果如圖9a所示.

6.4.3譜白化處理

分析圖9a的波形，子波相關(guān)噪聲較大，對深層微弱反射波造成干擾(標(biāo)記實(shí)線框區(qū)域)，因此繼續(xù)進(jìn)行譜白化處理.由于信號頻帶范圍為10～500 Hz，因此，分頻窄帶通濾波選擇的頻帶寬度為10 Hz.根據(jù)圖6a中各層反射波的到時(shí)，確定將長度為1 s的每個(gè)窄帶通濾波結(jié)果分成10段，每段0.1 s進(jìn)行時(shí)變增益，常數(shù)因子選擇1.譜白化處理后的結(jié)果如圖9b所示，圖10a，10b分別為譜白化處理前后的單道信號的幅頻特性.可看出，譜白化處理有效地補(bǔ)償了信號高頻部分幅度，雖然信噪比有所下降，但相關(guān)噪聲減小，消除了對深層微弱反射波的干擾.

圖9b即為應(yīng)用本文提出的CDRSBRES方法的檢測結(jié)果，其中直達(dá)波和5層反射波清晰可見，與基于基板附近信號作為參考信號的檢測結(jié)果圖6b相比，同樣位置并沒有虛假多次波出現(xiàn)，證明本方法檢測結(jié)果波形特征準(zhǔn)確，與模型相符.

圖9　基于重構(gòu)激發(fā)信號作為參考信號的相關(guān)檢測結(jié)果(a) 譜白化處理前； (b) 譜白化處理后.

圖10　地震信號幅頻特性(a) 圖9a的第10道信號； (b) 圖9b的第10道信號.

圖11　地震信號單道波形(豎線標(biāo)明準(zhǔn)確到時(shí)) (a) 圖9b的第1道信號(0～0.7 s)； (b) 圖6a的第1道信號(0～0.7 s).

分析到時(shí)檢測精度，抽取圖9b與圖6a的第1道信號，如圖11對比展示.準(zhǔn)確到時(shí)由豎線標(biāo)出，證明本方法子波到時(shí)檢測準(zhǔn)確.

圖12　PHVS 500數(shù)據(jù)重構(gòu)的震源激發(fā)信號

7實(shí)測數(shù)據(jù)結(jié)果與分析

數(shù)值模擬應(yīng)用證明了本方法的有效性.現(xiàn)將該方法應(yīng)用于問題描述中的西部某金屬礦區(qū)可控震源對比實(shí)驗(yàn)的PHVS 500地震數(shù)據(jù)，實(shí)驗(yàn)參數(shù)見問題描述.參照數(shù)值模擬中的步驟，針對PHVS 500地震數(shù)據(jù)進(jìn)行震源激發(fā)信號的重構(gòu)，結(jié)果如圖12所示，可看出信號低頻和高頻段均產(chǎn)生了畸變.

隨后對該P(yáng)HVS 500地震數(shù)據(jù)基于重構(gòu)激發(fā)信號作為參考信號進(jìn)行相關(guān)檢測，并進(jìn)行譜白化處理，得到的地震記錄如圖13a所示，并與相同激發(fā)點(diǎn)MiniVib T15000地震記錄圖1a進(jìn)行對比分析，此處為方便對比，如圖13b重新所示.

圖13a中，箭頭指示的有效信號及其下方其他波組特征，以及面波、聲波等信號，直觀上與圖13b表現(xiàn)一致.抽取圖13a，13b的第155道信號，截取0.15～0.4 s在圖14中對比展示，其中虛線為PHVS 500數(shù)據(jù)，實(shí)線為MiniVib T15000數(shù)據(jù).針對標(biāo)記點(diǎn)劃線窗口范圍內(nèi)，由箭頭標(biāo)注的同一子波，兩地震記錄的到時(shí)誤差約為0.001 s，兩波形相關(guān)系數(shù)為0.90，為高度線性相關(guān).

由于PHVS 500輸出力為300 N，而MiniVib T15000輸出力為30000 N，在相同噪聲條件下兩震源地震記錄信噪比相差約100倍，即40 dB，導(dǎo)致子波到時(shí)檢測存在微小偏差.同樣由于輸出力差距，遠(yuǎn)偏移距端(約400 m，圖13a，13b標(biāo)記點(diǎn)劃線框區(qū)域)，PHVS 500雖能夠檢測到有效信號，但能量相對較弱，信噪比低于MiniVib T15000.兩震源在輸出力相差100倍情況下，PHVS 500基于CDRSBRES方法檢測得到的地震記錄在波形特征和子波到時(shí)檢測上與大輸出力液壓震源MiniVib T15000地震記錄一致，由此可見CDRSBRES方法在實(shí)際應(yīng)用中效果明顯.

此外，由于驅(qū)動(dòng)形式不同，雖然兩臺(tái)震源采用相同的控制信號頻帶范圍(10～120 Hz)，但PHVS 500的實(shí)際激發(fā)信號頻帶范圍更寬(王俊秋等，2011)，比較近偏移距端淺層信號(圖13a，13b標(biāo)記虛線框區(qū)域)，圖13a的信噪比和分辨率較圖13b都得到了明顯改善，同相軸更加清晰.針對該P(yáng)HVS 500地震數(shù)據(jù)，基于常規(guī)兩種方法檢測得到的記錄已于圖1b，1c展示，現(xiàn)將CDRSBRES方法檢測結(jié)果圖13a與其進(jìn)行對比.針對子波到時(shí)檢測，圖1a，1b中，箭頭指示的同一子波到時(shí)檢測相差約0.012 s.該子波對應(yīng)圖14中箭頭所示的子波，精確對比得知，圖13a，13b對于該子波到時(shí)檢測相差約0.001 s.根據(jù)數(shù)據(jù)處理結(jié)果得到的均方根速度約為1860 m·s-1，以MiniVib T15000檢測結(jié)果為標(biāo)準(zhǔn)，基于震源控制信號作為參考信號的檢測結(jié)果對應(yīng)垂向精度誤差約11.16 m，而基于CDRSBRES方法的檢測結(jié)果誤差僅為0.93 m，后者檢測精度大幅提高.

圖13　西部某金屬礦區(qū)可控震源實(shí)驗(yàn)地震記錄(a) PHVS 500地震記錄，CDRSBRES方法； (b) MiniVib T15000地震記錄.

針對波形形態(tài)檢測，分別將圖13a與圖1c中標(biāo)記實(shí)線框區(qū)域放大，如圖15a，15b對比展示.圖15a中同樣區(qū)域并沒有類似的虛假多次波信號出現(xiàn).因此，在該測區(qū)地下結(jié)構(gòu)條件下，CDRSBRES方法檢測得到的結(jié)果解決了基于基板附近信號作為參考信號時(shí)存在的虛假多次波問題，提高了地震記錄的信噪比.

圖14　不同震源地震記錄同道信號波形比對(0.15～0.4 s)，虛線為PHVS 500的CDRSBRES檢測結(jié)果，實(shí)線為MiniVib T15000檢測結(jié)果

圖15　PHVS 500基于不同參考信號方法得到的地震記錄局部區(qū)域放大圖(a) 圖13a中標(biāo)記實(shí)線框區(qū)域； (b) 圖1c中標(biāo)記實(shí)線框區(qū)域.

8結(jié)論與討論

通過對本文的總結(jié)，可得出如下結(jié)論.

(1) 基于重構(gòu)激發(fā)信號的相關(guān)檢測參考信號方法適用于地下結(jié)構(gòu)復(fù)雜、基板-大地耦合不佳時(shí)電磁式可控震源地震數(shù)據(jù)的相關(guān)檢測.有利于解決基于基板附近信號作為參考信號引起的虛假多次波問題和基于震源控制信號作為參考信號導(dǎo)致的子波到時(shí)檢測誤差問題.

(2) CDRSBRES方法為微弱可控震源地震數(shù)據(jù)的高分辨率檢測提供了有效方法.電磁式可控震源PHVS 500由于設(shè)計(jì)上對便攜性的考慮(PHVS 500可兩人抬行)，輸出力較小，抗干擾能力弱，但是采用合適的檢測方法，如CDRSBRES方法，仍可獲得與大型液壓式可控震源MiniVib T15000質(zhì)量相當(dāng)?shù)牡卣鹩涗?

(3) 合理利用譜白化處理，可有效去除相關(guān)噪聲，進(jìn)一步提高可控震源地震記錄的分辨率，使電磁式可控震源更好地應(yīng)用于淺層高分辨率地震勘探.

本論文提出的方法以直達(dá)波提取為基礎(chǔ)，但未贅述直達(dá)波提取細(xì)節(jié).這里簡單進(jìn)行討論.直達(dá)波分離過程中，提取直達(dá)波應(yīng)以保證盡量完整分離子波，同時(shí)減少混入其他干擾信號為基本原則，具體地要考慮震源信號帶寬、地震子波能量、地表介質(zhì)特性等方面.

1) 震源控制信號頻帶.頻帶范圍越寬，相關(guān)檢測后的地震子波寬度越窄，此時(shí)子波提取窗口長度不應(yīng)過寬，避免混入其他干擾信號.

2) 提取子波能量要求.當(dāng)直達(dá)波與淺層反射信號、地表各類干擾信號混合得較嚴(yán)重時(shí)，需要縮短直達(dá)波提取窗口，以減少干擾信號，提高信噪比，此時(shí)可設(shè)定當(dāng)窗口內(nèi)波形能量大于全部直達(dá)波能量的80%～90%時(shí)達(dá)到要求，本文仿真和實(shí)際數(shù)據(jù)選取90%.

3) 地表介質(zhì)物性差異.直達(dá)波傳播過程中受地表介質(zhì)物性的影響，會(huì)被吸收和衰減.介質(zhì)成分不同，吸收和衰減影響存在較大差異，如果直達(dá)波衰減嚴(yán)重，為保證信噪比，需減小直達(dá)波窗口，反之增大窗口長度.

綜上，基于重構(gòu)激發(fā)信號的相關(guān)檢測參考信號方法能更好地實(shí)現(xiàn)地下結(jié)構(gòu)復(fù)雜區(qū)域的電磁式可控震源地震數(shù)據(jù)的高精度檢測，對于便攜式可控震源地震勘探具有重要意義.

致謝特別感謝吉林大學(xué)儀器科學(xué)與電氣工程學(xué)院教育部地球信息探測重點(diǎn)實(shí)驗(yàn)室的陳祖斌教授、張林行、孫峰、鄭凡、王俊秋副教授等和葛利華博士、馮博碩士等在西部某金屬礦區(qū)可控震源對比實(shí)驗(yàn)中提供的幫助！特別感謝張雪昂博士的全程幫助！

References

Baeten G, Ziolkowski A. 1990. The Vibroseis Source. Oxford: Elsevier.

Bian G Z, Zhang L Q. 1986. Spectral whitening of seismic data.GeophysicalProspectingforPetroleum(in Chinese), 25(2): 26-33.

Boucard D, Ollivrin G. 2010. Developments in vibrator control.GeophysicalProspecting, 58(1): 33-40.

Castenet A, Lavergne M. 1965. Vibrator controlling system. U. S. Patent No. 3208550.

Chen Z B, Lin J, Yu S B, et al. 2003. The development of portable shallow seismic vibrator.ChineseJournalofScientificInstrument(in Chinese), 24(3): 311-314.Chen Z B, Teng J W, Lin J, et al. 2005. Nonlinear analysis on the coupling process of electromagnetic vibrator and earth.ScienceinChina(SeriesD:EarthSciences) (in Chinese), 48(8): 1175-1182.Cui Z J, Li Z X, Chen Z L, et al. 2012. A study on the new method for determining small earthquake sequence type—Correlation analysis of spectral amplitude.ChineseJ.Geophys. (in Chinese), 55(5): 1718-1724, doi: 10.6038/j.issn.0001-5733.2012.02.028.Dong S W, Li T D, Chen X H, et al. 2012. Progress of deep exploration in mainland China: A review.ChineseJ.Geophys. (in Chinese), 55(12): 3884-3901, doi: 10.6038/j.issn.0001-5733.2012.12.002.

Duan Y Q. 2006. Matched filtering and wavelet shaping.OGP(in Chinese), 41(2): 156-159.

Ge L H, Jiang T, Lin J, et al. 2012. Numerical simulation of application of the Phased-array Vibrator System in mineral exploration.ChineseJ.Geophys. (in Chinese), 55(12): 4266-4276, doi: 10.6038/j.issn.0001-5733.2012.12.038.

Haines S S. 2006. Design and application of an electromagnetic vibrator seismic source.JournalofEnvironmentalandEngineeringGeophysics, 11(1): 9.Jiang T, Lin J, Chen Z B, et al. 2005. The correlation detection method of signals in phased array of vibroseis system.ChineseJournalofScientificInstrument(in Chinese), 26(4): 336-339. Jiang T, Lin J, Li T L, et al. 2006. Boosting signal to noise ratio of seismic signals using the phased-array vibrator system.ChineseJ.Geophys. (in Chinese), 49(6): 1819-1825.

Jiang T, Lin J, Yang D, et al. 2008. Analysis of directional seismic signal based on phased-array vibrator system.ChineseJ.Geophys. (in Chinese), 51(5): 1551-1556.Jiang T, Lin J, Chen Z B. 2009. SNR improvement of seismic signals with different time delays.ChineseJ.Geophys. (in Chinese), 52(6): 1651-1656. Jiang T, Lin J, Jia H Q, et al. 2012. Time-domain seismic beam-forming based on receiver arrays.ChineseJ.Geophys. (in Chinese), 55(12): 4277-4287, doi: 10.6038/j.issn.0001-5733.2012.12.039. Krohn C, Johnson M, Ho R, et al. 2010. Vibroseis productivity: shake and go.GeophysicalProspecting，58(1)：101-122.

Lebedev A V, Beresnev I A. 2004. Nonlinear distortion of signals radiated by vibroseis sources.Geophysics, 69(4): 968-977.

Liang T C, Lin J, Chen Z B. 2005. Explorating tunnel using electromagnetic vibrator.JournalofJilinUniversity(EarthScienceEdition) (in Chinese), 35(S1): 87-89.

Lin J. 2004. Seismic Exploration and Application of Vibrator Driven by Electromagnetic (in Chinese). Beijing: Science Press.

Nigel A. 1991. Vibroseis. New Jersey: Prentice Hall PTR.

Sallas J J. 1984. Seismic vibrator control and the downgoing P-wave.Geophysics, 49(6): 732-740.

Sallas J J. 2010. How do hydraulic vibrators work? A look inside the black box.GeophysicalProspecting, 58(1): 3-18.

Saragiotis C, Scholtz P, Bagaini C. 2010. On the accuracy of the ground force estimated in vibroseis acquisition.GeophysicalProspecting, 58(1): 69-80.

Tao Z F. 2007. The study on vibrator output distortion control.EGP(in Chinese), 17(1): 1-13.

Tong X Q, Lin J, Jiang T, et al. 2012. Summary of development of land vibrator.ProgressinGeophysics(in Chinese), 27(5): 1912-1921, doi: 10.6038/j.issn.1004-2903.2012.05.011.

van der Veen M, Brouwer J, Helbig K. 1999. Weighted sum method for calculating ground force: an evaluation by using a portable vibrator system.GeophysicalProspecting, 47(3): 251-267.

Wang J Q, Lin J, Jiang T, et al. 2011. Experiment and application of controlled vibrator seismic method at Jinchang copper-nickel metal deposit.JournalofJilinUniversity(EarthScienceEdition) (in Chinese), 41(5): 1617-1622.

Wang J Q, Lin J, Chen Z B, et al. 2012. Experiment and application of mini vibrators seismic exploration in oil shale.JournalofJilinUniversity(EarthScienceEdition) (in Chinese), 42(S3): 265-270.Wei Z H. 2010. Modelling and modal analysis of seismic vibrator baseplate.GeophysicalProspecting, 58(1): 19-32.

Wei Z H, Philips T F, Hall M A. 2010. Fundamental discussions on seismic vibrators.Geophysics, 75(6): W13-W25.

Wei Z H, Phillips T F. 2012. Break through the limits of vibroseis data quality.GeophysicalProspecting, 60(2): 373-393.

Zhao C L, Lu C, Hao T Y, et al. 2013. A study of the high-precision modular lightweight small vibrator.ChineseJ.Geophys. (in Chinese), 56(11): 3690-3698, doi: 10.6038/cjg20131110.Ziolkowski A. 2010. Review of vibroseis data acquisition and processing for better amplitudes: adjusting the sweep and deconvolving for the time-derivative of the true groundforce.GeophysicalProspecting, 58(1): 41-54.

附中文參考文獻(xiàn)

邊國柱, 張立群. 1986. 地震數(shù)據(jù)的譜白化處理. 石油物探, 25(2): 26-33.

陳祖斌, 林君, 于生寶等. 2003. 輕便淺層地震可控震源的研制. 儀器儀表學(xué)報(bào), 24(3): 311-314.

陳祖斌, 騰吉文, 林君等. 2005. 電磁式可控震源-大地耦合過程中的非線性特征分析. 中國科學(xué)(D輯: 地球科學(xué)), 35(4): 333-339.

崔子健, 李志雄, 陳章立等. 2012. 判別小震群序列類型的新方法研究——振幅譜相關(guān)分析法. 地球物理學(xué)報(bào), 55(5): 1718-1724, doi: 10.6038/j.issn.0001-5733.2012.02.028.

董樹文, 李廷棟, 陳宣華等. 2012. 我國深部探測技術(shù)與實(shí)驗(yàn)研究進(jìn)展綜述. 地球物理學(xué)報(bào), 55(12): 3884-3901, doi: 10.6038/j.issn.0001-5733.2012.12.002.

段云卿. 2006. 匹配濾波與子波整形. 石油地球物理勘探, 41(2): 156-159.

葛利華, 姜弢, 林君等. 2012. 相控震源在礦產(chǎn)勘探中應(yīng)用的數(shù)值模擬研究. 地球物理學(xué)報(bào), 55(12): 4266-4276, doi: 10.6038/j.issn.0001-5733.2012.12.038.

姜弢, 林君, 陳祖斌等. 2005. 可控震源相控地震的相關(guān)檢測技術(shù). 儀器儀表學(xué)報(bào), 26(4): 336-339.

姜弢, 林君, 李桐林等. 2006. 相控震源對地震信號信噪比的改善研究. 地球物理學(xué)報(bào), 49(6): 1819-1825.

姜弢, 林君, 楊冬等. 2008. 相控震源定向地震波信號分析. 地球物理學(xué)報(bào), 51(5): 1551-1556.

姜弢, 林君, 陳祖斌. 2009. 延時(shí)參數(shù)對地震信號信噪比的改善. 地球物理學(xué)報(bào), 52(6): 1651-1656.

姜弢, 林君, 賈海青等. 2012. 基于接收陣列的時(shí)域地震波束形成方法. 地球物理學(xué)報(bào), 55(12): 4277-4287, doi: 10.6038/j.issn.0001-5733.2012.12.039.

梁鐵成, 林君, 陳祖斌. 2005. 電磁驅(qū)動(dòng)可控震源在隧道探測中的應(yīng)用. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 35(S1): 87-89.

林君. 2004. 電磁驅(qū)動(dòng)可控震源地震勘探原理及應(yīng)用. 北京: 科學(xué)出版社.

陶知非. 2007. 可控震源輸出信號畸變控制問題的研究. 物探裝備, 17(1): 1-13.

佟訓(xùn)乾, 林君, 姜弢等. 2012. 陸地可控震源發(fā)展綜述. 地球物理學(xué)進(jìn)展, 27(5): 1912-1921, doi: 10.6038/j.issn.1004-2903.2012.05.011.

王俊秋, 林君, 姜弢等. 2011. 可控震源地震方法在金昌銅鎳礦區(qū)的應(yīng)用實(shí)驗(yàn). 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 41(5): 1617-1622.

王俊秋, 林君, 陳祖斌等. 2012. 小型可控震源在油頁巖地震勘探中的應(yīng)用試驗(yàn).吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 42(S3): 265-270.

趙春蕾, 盧川, 郝天珧等. 2013. 高精度組合式輕便小型可控震源的研究. 地球物理學(xué)報(bào), 56(11): 3690-3698, doi: 10.6038/cjg20131110.

(本文編輯胡素芳)

基金項(xiàng)目深部探測技術(shù)與實(shí)驗(yàn)研究專項(xiàng)SinoProbe-09-06 (201311197；201011083) ，SinoProbe-09-04 (201011081)和高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(20110061110053)聯(lián)合資助.

作者簡介楊志超，男，1986年生，博士，主要研究方向?yàn)榭煽卣鹪吹卣鹦盘柼幚?E-mail: zcyang12@mails.jlu.edu.cn *通訊作者姜弢，女，1969年生，博士，教授，博士生導(dǎo)師，主要研究方向?yàn)榈卣饍x器及信號處理.E-mail:jiang_t@jlu.edu.cn

doi:10.6038/cjg20160127 中圖分類號P631

收稿日期2014-12-09，2015-10-09收修定稿

Research on correlation detection reference signal based on the reconstructed excitation signal for electromagnetic seismic vibrator

YANG Zhi-Chao1,2， JIANG Tao1,2*， XU Xue-Chun3， LIN Jun1,2，JIA Hai-Qing1,2， MA Rui-Qi1,2， CHEN Zhen-Zhen1,2

1CollegeofInstrumentScienceandElectricalEngineering,JilinUniversity,Changchun130061,China2KeyLaboratoryofGeo-explorationInstruments,MinistryofEducationofChina(JilinUniversity),Changchun130061,China3CollegeofEarthScience,JilinUniversity,Changchun130061,China

AbstractThe electromagnetic seismic vibrator has wide application prospect in high resolution seismic exploration. However, in searching of correlation detection reference signal for electromagnetic seismic vibrator, two conventional methods both have problems compared to hydraulic seismic vibrator in the condition that the subsurface media is complex or the base-plate is poorly coupled with ground. One method that based on the control-signal will cause error of wavelet arrival time. The other one that based on the near base-plate signal will cause error of waveform characters, called false multiple waves in this paper.

KeywordsElectromagnetic seismic vibrator； Correlation detection； Reference signal； Excitation signal； Reconstruction

楊志超， 姜弢， 徐學(xué)純等. 2016. 基于重構(gòu)激發(fā)信號的電磁式可控震源相關(guān)檢測參考信號方法研究.地球物理學(xué)報(bào)，59(1):318-329,doi:10.6038/cjg20160127.

Yang Z C, Jiang T, Xu X C, et al. 2016. Research on correlation detection reference signal based on the reconstructed excitation signal for electromagnetic seismic vibrator.ChineseJ.Geophys. (in Chinese),59(1):318-329,doi:10.6038/cjg20160127.