賈海青,姜 弢,徐學(xué)純,林 君

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

基于全方向波束定向的地震記錄信噪比改善方法

賈海青1,2,姜 弢1,2,徐學(xué)純3,林 君1,2

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

常規(guī)波束定向方法利用多炮地震記錄、采用恒定延時(shí)差實(shí)現(xiàn)目標(biāo)方向上地震波場的加強(qiáng),所用的地震記錄是由相同震源在不同位置激發(fā)得到的,由于不同地震記錄的震源與大地的耦合程度不同,即震源的一致性存在差異,因而降低了波束定向的質(zhì)量,同時(shí),該方法未能實(shí)現(xiàn)全方向波束加強(qiáng),波束定向結(jié)果的整體信噪比仍較低。為此,提出了一種全方向波束定向方法(Omnidirectional Beam-forming,OB),該方法基于波場疊加原理,通過不同方向有效波場的疊加,實(shí)現(xiàn)地震記錄全方向有效信號(hào)的增強(qiáng)。其具體實(shí)現(xiàn)包括5個(gè)步驟:單炮記錄的復(fù)制、角度域劃分、延時(shí)參數(shù)計(jì)算及波束定向、有效波場的提取和有效波場的合并。利用數(shù)值模擬的含隨機(jī)噪聲單炮記錄,驗(yàn)證了OB方法的正確性和有效性。將OB方法應(yīng)用于某地區(qū)實(shí)際地震資料處理,使地震資料中的隨機(jī)噪聲得到明顯壓制,信噪比提高了6.75dB。

全方向波束定向;單炮記錄;信噪比;隨機(jī)噪聲壓制

波束定向技術(shù)最早源自于雷達(dá)天線領(lǐng)域,其核心是波的干涉加強(qiáng)。作為一種陣列處理技術(shù),波束定向方法合成的定向波束實(shí)現(xiàn)了目標(biāo)區(qū)域的信號(hào)成倍增強(qiáng),大大提高了對(duì)目標(biāo)的探測、識(shí)別和追蹤能力。由于地震波與電磁波具有波的共性,因此近幾年來通過借鑒雷達(dá)天線領(lǐng)域的波束形成理論,波束定向技術(shù)在地震勘探領(lǐng)域得到了快速的發(fā)展。

20世紀(jì)70年代,美國石油公司在Oklahoma州Tulsa地區(qū)采用單頻控制信號(hào)進(jìn)行了地震波束形成(beam-forming)實(shí)驗(yàn),驗(yàn)證了定向地震波的存在[1]。2005年,日本東京大學(xué)地震研究所利用4個(gè)精密可控常時(shí)震源系統(tǒng)(ACROSS)組成的相控陣列實(shí)現(xiàn)了對(duì)深部地層構(gòu)造的研究[2-3]。在國內(nèi),姜弢等[4-8]從2002年開始研制第一套電磁式可控震源相控系統(tǒng),并開始對(duì)相控震源地震波定向方法進(jìn)行研究。2009年,姜弢等[9]利用研制的電磁式可控震源開展了相控震源地震勘探實(shí)驗(yàn),證明了相控震源改善地震記錄信噪比的有效性。2006年,王忠仁等[10]利用有限差分法對(duì)相控陣震源單頻和變頻激發(fā)方式的地震波場進(jìn)行了數(shù)值模擬,研究分析了這兩種方式下陣列的方向特性。2012年,葛麗華等[11]利用有限差分方法數(shù)值模擬了不同礦區(qū)模型下相控震源產(chǎn)生的定向地震波場對(duì)礦區(qū)的加強(qiáng)效果。

地震波束的形成是指多個(gè)地震波場通過延時(shí)、疊加,在特定方向上實(shí)現(xiàn)波場的加強(qiáng)。由地震勘探互易性原理可知,地震波束的形成既可以在源端,也可以在接收端。源端地震波束的合成是在地震勘探的采集階段實(shí)現(xiàn)的,即通過控制震源陣列中每個(gè)震源的初始相位,在目標(biāo)區(qū)域激發(fā)定向地震波,因此該方法被稱為相控震源方法。相控震源勘探方法存在需要的震源數(shù)多、震源一致性難以保障以及采集成本高等問題,因而發(fā)展較為緩慢。2009年,劉暢[12]提出了單震源定向照明方法,該方法不僅有效避免了相控震源不一致性的問題,而且能夠方便地實(shí)現(xiàn)不同方向上的波束定向。2012年,姜弢等[13]提出了基于接收陣列的時(shí)域地震波束形成方法,并將其成功應(yīng)用于野外地震數(shù)據(jù)?；诮邮斩说卣鸩ㄊ暮铣煞椒?TSBBRA)是在地震勘探的數(shù)據(jù)預(yù)處理階段實(shí)現(xiàn)的,處理的數(shù)據(jù)是同一個(gè)可控震源按照傳統(tǒng)采集方式采集的多炮數(shù)據(jù),基本解決了相控震源勘探方法存在的問題。

但是,利用TSBBRA方法進(jìn)行波束定向存在以下3個(gè)問題:第一,在實(shí)際勘探中,由于震源在不同位置與大地的耦合程度不一致等原因,即使用同一個(gè)震源進(jìn)行地震數(shù)據(jù)采集,震源的一致性仍比較低,導(dǎo)致對(duì)采集的數(shù)據(jù)進(jìn)行波束定向會(huì)存在主波束方向上信號(hào)被加強(qiáng)的程度降低問題;第二,TSBBRA方法通常將多炮數(shù)據(jù)合成一炮數(shù)據(jù),因而存在合成后的單炮數(shù)據(jù)總炮數(shù)減少的問題;第三,TSBBRA方法僅能加強(qiáng)單個(gè)主波束方向上的有效信號(hào),主波束方向外的區(qū)域信號(hào)會(huì)發(fā)生畸變，信噪比很低。針對(duì)以上問題,本文提出了一種全方向波束定向方法(OB),該方法利用單個(gè)地震記錄進(jìn)行全方向波束定向,因而更大程度地保證了震源的一致性,實(shí)現(xiàn)了地震記錄全方向的信號(hào)加強(qiáng),從而能有效改善原始資料整體的質(zhì)量。模擬數(shù)據(jù)和實(shí)際地震資料處理結(jié)果均表明,OB方法能夠有效提高單個(gè)地震記錄的信噪比。

1 理論基礎(chǔ)

1.1 TSBBRA波束定向原理

波束定向?qū)嶋H上是對(duì)具有恒定相位差的波場進(jìn)行疊加。TSBBRA方法主要利用采集的多炮地震記錄在目標(biāo)點(diǎn)進(jìn)行波束定向和波束合成,其基本原理如圖1a 和圖1b所示。同一個(gè)震源分別在a,b和c三個(gè)點(diǎn)處激發(fā)得到地震記錄Sa,Sb和Sc,震源間距為d。對(duì)地震記錄Sa,Sb和Sc分別引入延時(shí)-τ,0和τ,得到延時(shí)后的地震記錄Sa(t-τ),Sb(t)和Sc(t+τ),將延時(shí)后的地震記錄疊加,則相當(dāng)于對(duì)震源在點(diǎn)a,b和c處同時(shí)激發(fā)的地震波場的疊加。由于延時(shí)τ補(bǔ)充了震源在不同位置激發(fā)時(shí)信號(hào)到達(dá)e點(diǎn)的波程差,因而主波束方向目標(biāo)點(diǎn)的有效信號(hào)得到加強(qiáng)。同時(shí),由于TSBBRA是在接收端進(jìn)行處理,因而能夠方便快速地實(shí)現(xiàn)對(duì)目標(biāo)體的波束定向。

1.2 OB波束定向原理

與TSBBRA波束定向方法類似,OB波束定向方法同樣是通過波場疊加來實(shí)現(xiàn)有效信號(hào)的加強(qiáng)(圖1)。由圖1a可知,TSBBRA合成定向波束的過程相當(dāng)于同一個(gè)震源在a,b和c三個(gè)點(diǎn)處激發(fā)波場的延時(shí)疊加。設(shè)圖1a中e點(diǎn)處接收的信號(hào)分別為Sae,Sbe和Sce,對(duì)其依次做-τ,0和τ延時(shí),得到Sae(t-τ),Sbe(t)和Sce(t+τ),對(duì)延時(shí)后的信號(hào)疊加,得：

(1)

圖1 TSBBRA和OB波束定向原理及波束合成方法

由地震波的傳播規(guī)律可知,在震源點(diǎn)處激發(fā)的地震波經(jīng)過地下介質(zhì)到達(dá)檢波點(diǎn)的路徑與檢波點(diǎn)處激發(fā)的地震波到達(dá)震源點(diǎn)處的旅行路徑完全一樣,因此,可以將圖1a和圖1b中的檢波器和震源進(jìn)行等效對(duì)換,結(jié)果如圖1c和圖1d所示。圖1c表明,檢波器和震源等效對(duì)換后相當(dāng)于震源在e點(diǎn)激發(fā),在a,b和c三個(gè)點(diǎn)接收信號(hào),得到一個(gè)單炮記錄。設(shè)該單炮記錄為Re,a,b和c處接收的信號(hào)分別為Rea,Reb和Rec,同樣依次做-τ,0和τ延時(shí),得到Rea(t-τ),Reb(t)和Rec(t+τ),對(duì)延時(shí)后的信號(hào)疊加得：

(2)

由波場的一致性分析可知,在不考慮噪聲的情況下,有：

(3)

對(duì)比(1)式、(2)式和(3)式可知,對(duì)單炮記錄Re的波束定向過程相當(dāng)于將多個(gè)相同單炮記錄Re以道間距d為震源間距的波場延時(shí)、疊加的過程。因此,一個(gè)單炮記錄可以利用TSBBRA波束定向方法合成定向地震波束。由上述分析可知,在不考慮環(huán)境噪聲的情況下,對(duì)單個(gè)記錄進(jìn)行n元波束定向的總波場為[13]：

(4)

式中:E1為單個(gè)震源產(chǎn)生的場強(qiáng);φ為相鄰震源的相位差;k為自由空間波數(shù);d為陣元間距;α為地下任一觀測點(diǎn)方向的軸向夾角;αmax為主波束方向。

同樣可得到均勻介質(zhì)條件下的地震波場方向因子[13]:

(5)

根據(jù)前人研究成果[14]可得：

(6)

式中:v為均勻介質(zhì)速度;τ為波束定向過程中的延時(shí)參數(shù)。

由(6)式可知,延時(shí)參數(shù)τ與主波束方向αmax是一一對(duì)應(yīng)的,單一延時(shí)參數(shù)僅能加強(qiáng)一個(gè)主波束方向上的地震信號(hào)。因此,要得到地下半空間全方向加強(qiáng)的波束定向結(jié)果,就需要對(duì)地下半空間所有方向同時(shí)實(shí)現(xiàn)波束的加強(qiáng)。當(dāng)?shù)叵陆橘|(zhì)為非均勻介質(zhì)時(shí),(6)式中的速度則為等效速度,對(duì)于某一延時(shí)參數(shù),主波束方向會(huì)隨著速度的變化而變化。

由于不同的延時(shí)參數(shù)τ可以加強(qiáng)不同方向的波場,因此,q個(gè)不同的延時(shí)參數(shù)(圖1d)可以得到q個(gè)不同的加強(qiáng)波場。根據(jù)波場的疊加原理,將地下半空間全方向加強(qiáng)后的有效波場進(jìn)行疊加,就可以得到全方向加強(qiáng)的波束定向結(jié)果。

因此,地下半空間所有方向有效波場疊加后的總波場為：

(7)

式中:E(τi)表示延時(shí)為τi合成的有效地震波場。

2 OB波束定向方法

由OB波束定向原理(圖1c)可知,對(duì)一個(gè)單炮記錄進(jìn)行波束定向的前提條件是滿足定向中多元數(shù)據(jù)的需要,以便對(duì)單炮記錄進(jìn)行常規(guī)波束定向,這種多元數(shù)據(jù)可通過復(fù)制單炮記錄的方法得到。由(7)式可知,要實(shí)現(xiàn)真正的全方向波束定向,關(guān)鍵在于對(duì)定向結(jié)果的有效波場提取和合并。因此,OB波束定向方法主要包括5個(gè)步驟:單炮記錄復(fù)制、角度域劃分、延時(shí)參數(shù)計(jì)算及波束形成、有效波場提取和有效波場合并。

2.1 單炮記錄復(fù)制

設(shè)一個(gè)單炮記錄為R,對(duì)R進(jìn)行n元波束定向,需要將R復(fù)制為n個(gè)地震記錄R1,R2,…,Rn?？蓪⑵渥鳛橛蒼個(gè)震源的線性等間隔陣列激發(fā)所得到的n個(gè)地震記錄,陣元間距等于道間距d,其對(duì)應(yīng)的地震波場分別為E1,E2,…,En,這樣就給R的定向提供了基礎(chǔ)條件。

2.2 角度域劃分

理論上,要實(shí)現(xiàn)對(duì)地下全方向的波束定向,需要對(duì)每一個(gè)方向進(jìn)行波束定向,但在實(shí)際操作中這是無法實(shí)現(xiàn)的。由于波束定向后的主波束方向是有一定波束寬度的,在波束寬度范圍內(nèi),地震波場的加強(qiáng)程度大體相當(dāng),因此,為了等效實(shí)現(xiàn)對(duì)地下全方向的波束定向,可以對(duì)地下半空間的方向進(jìn)行角度域劃分,將劃分后的角度區(qū)間的中心角度作為主波束方向。角度域劃分要以所有主波束方向的寬度能夠覆蓋角度區(qū)間為原則。

(8)

由(8)式可知,在地下介質(zhì)、定向元數(shù)以及陣元間距一定的情況下,主波束寬度主要與主波束方向有關(guān)。為了直觀地認(rèn)識(shí)主波束寬度與主波束方向的關(guān)系,令n=25,v=1000m/s,f=80Hz,d=1m,由(5)式可得到對(duì)應(yīng)的定向地震波束方向,如圖2所示。圖2a至圖2d的主波束方向αmax分別為90°,60°,45°和30°,其對(duì)應(yīng)的延時(shí)τ分別為0，0.87，0.71，0.50ms,主波束寬度分別為16.99°,19.75°,24.75°和44.10°。由圖2可以看出,對(duì)不同方向進(jìn)行定向的主波束寬度并非一成不變,而是隨著主波束方向αmax(αmax∈[0,90°])的增大,波束寬度減小,當(dāng)αmax=90°時(shí),波束寬度最小。

地下半空間角度域的劃分有多種方式,本文采用均勻等間隔角度劃分方法,即將α∈[0,180°]從α=90°開始向兩個(gè)方向劃分為2m個(gè)區(qū)域,角度間隔為Δα,劃分后的2m+1個(gè)角度分界分別為(90°,90°±Δα,90°±2Δα,…,90°±m(xù)Δα),令mΔα=90°,m為正整數(shù)。由于主波束寬度隨主波束方向的變化而變化,為了實(shí)現(xiàn)多個(gè)主波束完全覆蓋的目的,本文根據(jù)最小主波束寬度確定Δα。設(shè)αmax=90°,主波束寬度為2α0.5(90°),則0<Δα≤α0.5(90°)。

因此,對(duì)于均勻等間隔角度區(qū)間劃分,有:

(9)

2.3 延時(shí)參數(shù)計(jì)算及波束形成

2.4 有效波場提取

源點(diǎn)的波場在接收端以地震記錄的形式被表現(xiàn)出來,由于地下介質(zhì)的未知,難以根據(jù)定向地震波波束的方向來確定地震記錄中被加強(qiáng)的有效波場。為此,本文利用局部相關(guān)系數(shù)譜來提取波束定向結(jié)果中的有效波場,具體方法如下:

3) 滑動(dòng)時(shí)窗T,計(jì)算得到所有道的局部最大相關(guān)系數(shù)譜C(τp)。

2.5 有效波場合并

圖2 不同主波束方向的地震波束方向

(10)

由此可得到最終的疊加數(shù)據(jù):

(11)

3 模擬數(shù)據(jù)驗(yàn)證

為了驗(yàn)證OB方法的正確性,本文采用3層水平層狀介質(zhì)模型(第1層速度為800m/s,厚度為40m;第2層速度為1000m/s,厚度為75m;第3層速度為1400m/s,厚度為175m)進(jìn)行數(shù)值模擬,信號(hào)主頻為45Hz;采用5m道間距,道數(shù)201，中間激發(fā)、兩邊接收方式得到包含3個(gè)反射波信號(hào)的單炮記錄(圖3a),其采樣間隔為0.1ms,記錄時(shí)長0.7s。對(duì)該單炮記錄加入-23.23dB隨機(jī)噪聲，結(jié)果如圖3b所示。

圖3 模擬單炮地震記錄

圖4 不同延時(shí)9元波束定向后的單炮記錄

利用局部相關(guān)系數(shù)譜法,對(duì)不同延時(shí)9元波束定向后的結(jié)果進(jìn)行畸變切除,提取有效波場,圖5為對(duì)圖4提取有效波場后的地震記錄。由圖5可見,主波束方向和主波束方向外區(qū)域的隨機(jī)噪聲得到有效壓制,實(shí)現(xiàn)了有效波場的準(zhǔn)確提取。

圖6為18個(gè)方向的有效波場疊加后的地震記錄和原始地震記錄的對(duì)比圖,其中圖6a為原始含噪單炮記錄,圖6b為全方向加強(qiáng)后的地震記錄,圖6c 為圖6a與圖6b的差值剖面。對(duì)比圖6a和圖6b可知,經(jīng)過OB方法進(jìn)行全方向波束定向后,單炮記錄的全局質(zhì)量得到明顯改善,信噪比明顯提高；由圖6c可見,差值剖面中有效信號(hào)殘存較少,隨機(jī)噪聲在很大程度上被消除,表明OB方法能夠有效地去除單炮記錄中的隨機(jī)噪聲。為了定量分析OB方法處理前后對(duì)單炮記錄質(zhì)量的改善,計(jì)算了圖6a和圖6b中3個(gè)反射層的信噪比曲線,圖7a至圖7c分別為第1層、第2層和第3層反射波的信噪比,其中虛線為OB方法處理前的信噪比,實(shí)線為OB方法處理后的信噪比。由圖7可知,經(jīng)過OB方法處理后,不僅單炮記錄中所有道反射波的信噪比都得到不同程度的提高,而且對(duì)于同一反射波,所有道的信噪比提高程度基本一致。經(jīng)計(jì)算,3層反射波信號(hào)平均提高了7.14dB,4.19dB和5.51dB,單炮記錄整體質(zhì)量提高了13.64dB。

圖5 不同延時(shí)9元波束定向結(jié)果利用局部互相關(guān)系數(shù)譜提取的有效數(shù)據(jù)

圖6 OB方法處理前(a)、后(b)的單炮記錄及兩者之差(c)

圖7 OB方法處理前(虛線)、后(實(shí)線)反射波的信噪比曲線

4 實(shí)際數(shù)據(jù)處理

圖8a是某地區(qū)一個(gè)被隨機(jī)噪聲嚴(yán)重干擾的地震記錄,道間距為10m,在0.8s左右可以看到由多個(gè)相位組成的模糊的反射同相軸。為了分析OB方法改善地震記錄信噪比的能力,對(duì)該地震記錄進(jìn)行了9元波束定向(即將地震數(shù)據(jù)復(fù)制9次)處理。由公式(8)可知,波束寬度隨地震波長的增大而增大,為了全方向波束定向結(jié)果能夠完全覆蓋地下半空間,選取最小地震波長為10m(有效信號(hào)的頻率范圍為25～100Hz,地下介質(zhì)最小速度為1000m/s),主波束方向?yàn)?0°的波束寬度,進(jìn)行地下半空間的角度域劃分。經(jīng)計(jì)算,α0.5(90°)=11.37°/2≈5.685°,因此,取Δα=90°/17≈5.294°,則把地下半空間劃分為34個(gè)角度區(qū)間,依據(jù)公式(6)求取對(duì)應(yīng)的延時(shí)參數(shù),進(jìn)行9元波束定向。將波束定向后的結(jié)果用局部相關(guān)系數(shù)譜提取有效波場,通過加權(quán)疊加得到最終的全方向定向結(jié)果,如圖8b所示。

對(duì)比圖8a和圖8b可以發(fā)現(xiàn),經(jīng)過OB方法處理后的單炮記錄中隨機(jī)噪聲得到了很好的壓制,0.8s左右的有效反射信號(hào)得到了很大程度的加強(qiáng)和改善,信噪比得到明顯的提高。

圖9為圖8中L1矩形框內(nèi)放大顯示(210～249道,0.6～0.9s),圖10為圖8中L2矩形框內(nèi)放大顯示(50～100道,0.75～0.85s)。從圖9可以看出,處理前單炮記錄中相對(duì)較強(qiáng)的同相軸經(jīng)過OB方法處理后,有效信號(hào)的連續(xù)性和幅度都得到明顯的提高;從圖10可以看出,處理前單炮記錄中局部微弱、幾乎難以分辨的同相軸經(jīng)過OB方法處理后可以清晰地分辨出有效信號(hào)。為了定量分析L1和L2矩形框內(nèi)的信噪比變化,本文采用時(shí)域奇異值(SVD)分解的方法進(jìn)行信噪比的估算[15]。經(jīng)計(jì)算,L1和L2內(nèi)信噪比在應(yīng)用OB方法后分別提高了6.75dB和3.03dB。實(shí)際數(shù)據(jù)處理效果證明了OB方法的有效性。

圖8 OB方法處理前(a)、后(b)的野外單炮記錄

圖9 L1局部放大

圖10 L2局部放大

5 結(jié)論與討論

本文針對(duì)傳統(tǒng)波束定向方法存在的問題，提出了一種全方向波束定向(OB)方法。該方法通過復(fù)制單炮記錄的方法得到多元數(shù)據(jù)，通過對(duì)地下半空間角度域劃分等效實(shí)現(xiàn)地下全方向的波束定向，利用局部相關(guān)系數(shù)譜提取波束定向結(jié)果中的有效波場并進(jìn)行疊加。模擬數(shù)據(jù)和實(shí)際單炮記錄處理結(jié)果表明，該方法可以有效壓制單炮記錄全方向的隨機(jī)噪聲，明顯改善地震記錄的質(zhì)量。

波束定向元數(shù)越多，主波束方向波場強(qiáng)度越強(qiáng)，波束寬度越小，地下全方向角度域劃分的角度區(qū)間越多，計(jì)算量相應(yīng)越大，因此，在主波束方向波場強(qiáng)度滿足要求時(shí)，波束定向元數(shù)應(yīng)盡量減少，以便提高OB方法的計(jì)算效率。

[1] Arnold M E.Beam forming with vibrator arrays[J].Geophysics,1977,42(7):1321-1338

[2] 廖成旺,莊燦濤,梁鴻森.精密可控常時(shí)震源系統(tǒng)(ACROSS)的初步實(shí)驗(yàn)[J].中國地震,2003,19(1):89-96 Liao C W,Zhuang C T,Liang H S.Initial experimental result of accurately controlled routinely operated signal system (ACROSS)[J].Earthquake Research in China,2003,19(1):89-96

[3] Yoko H,Hiroyuki K,Mineo K.Autoregressive modeling of transfer functions in frequency domain to determine complex travel times[J].Earth Planets Space,2001,53 (1):3-11

[4] 姜弢,林君,陳祖斌,等.相控震源地震波定向技術(shù)[J].吉林大學(xué)學(xué)報(bào)(信息科學(xué)版),2004,22(3):181-184Jiang T,Lin J,Chen Z B,et al.Seismic beam-forming in phased array of vibroseis[J].Journal of Jilin University (Information Science Edition),2004,22(3):181-184

[5] 姜弢,林君,陳祖斌,等.可控震源相控系統(tǒng)的信號(hào)檢測與時(shí)延估計(jì)[J].儀器儀表學(xué)報(bào),2005,26(4):336-339 Jiang T,Lin J,Chen Z B,et al.The correlation detection method of signals in phased array of vibroseis system [J].Chinese Journal of Scientific Instrument,2005,26(4):336-339

[6] 姜弢,林君,李桐林,等.相控震源對(duì)地震信號(hào)信噪比的改善研究[J].地球物理學(xué)報(bào),2006,49(6):1819-1825 Jiang T,Lin J,Li T L,et al.Boosting signal to noise ratio of seismic signals using the phased-array vibrator system[J].Chinese Journal of Geophysics,2006,49(6):1819-1825

[7] 姜弢,林君,陳祖斌,等.相控震源對(duì)水平層狀地下介質(zhì)的高信噪比檢測[J].儀器儀表學(xué)報(bào),2006,27(11):1369-1372 Jiang T,Lin J,Chen Z B,et al.High signal-to-noise ratio detection of horizontal layer underground medium with phased-array vibroseis[J].Chinese Journal of Scientific Instrument,2006,27(11):1369-1372

[8] 姜弢,林君,楊冬,等.相控震源定向地震波信號(hào)分析[J].地球物理學(xué)報(bào),2008,51(5):1551-1556 Jiang T,Lin J,Yang D,et al.Analysis of directional seismic signal based on phased-array vibrator system[J].Chinese Journal of Geophysics,2008,51(5):1551-1556

[9] 姜弢,林君,陳祖斌.延時(shí)參數(shù)對(duì)地震信號(hào)信噪比的改善[J].地球物理學(xué)報(bào),2009,52(6):1651-1656 Jiang T,Lin J,Chen Z B.SNR improvement of seismic signals with different time delays[J].Chinese Journal of Geophysics,2009,52(6):1651-1656

[10] 王忠仁,林君,馮聲涯.地震勘探中相控陣震源的方向特性研究[J].地球物理學(xué)報(bào),2006,49(4):1191-1197 Wang Z R,Lin J,Feng S Y.Directivity of phase array vibrators in seismic exploration [J].Chinese Journal of Geophysics,2006,49(4):1191-1197

[11] 葛利華,姜弢,林君,等.相控震源在礦產(chǎn)勘探中應(yīng)用的數(shù)值模擬研究[J].地球物理學(xué)報(bào),2012,55(12):4266-4276 Ge L H,Jiang T,Lin J,et al.Numerical simulation of application of the phased-array vibrator system in mineral exploration [J].Chinese Journal of Geophysics,2012,55(12):4266-4276

[12] 劉暢.定向照明地震數(shù)據(jù)合成的仿真算法設(shè)計(jì)與實(shí)現(xiàn)[D].長春:吉林大學(xué),2009 Liu C.Design and realization of the simulation algorithm on the synthesis of directional illumination data[D].Changchun:Jilin University,2009

[13] 姜弢,林君,賈海青,等.基于接收陣列的時(shí)域地震波束形成方法[J].地球物理學(xué)報(bào),2012,55(12):4277-4287 Jiang T,Lin J,Jia H Q,et al.Time-domain seismic beam-forming based on receiver arrays [J].Chinese Journal of Geophysics,2012,55(12):4277-4287

[14] 姜弢.基于相控震源的地震波定向方法研究[D].長春:吉林大學(xué),2006 Jiang T.Study on seismic beam-forming method based on phased-array vibrator system[D].Changchun:University of Jilin,2006

[15] 劉洋,李承楚.地震資料信噪比估計(jì)的幾種方法[J].石油地球物理勘探,1997,32(2):257-262 Liu Y,Li C C.Some methods for estimating the signal noise ratio of seismic data[J].Oil Geophysical Prospecting,1997,32(2):257-262

(編輯:戴春秋)

The SNR improvement method based on omnidirectional beam-forming for vibrator seismic record

Jia Haiqing1,2,Jiang Tao1,2,Xu Xuechun3,Lin Jun1,2

(1.KeyLaboratoryofGeo-explorationInstrumentationofMinistryofEducation,JilinUniversity,Changchun130026,China; 2.CollegeofInstrumentScienceandElectricalEngineering,JilinUniversity,Changchun130026,China; 3.CollegeofEarthSciences,JilinUniversity,Changchun130026,China)

Conventional beam-forming method is applied to strengthen the seismic wavefield in target directions by adopting multi seismic records and constant delay difference.Since the multi vibrator seismic records were achieved by the same vibrator vibrating at different positions,the coupling between the vibrator and the ground is different,which causes the difference on the consistency of vibrators and eventually decreases the quality of beam-forming.Meanwhile,the method cannot realize beam strengthened in all directions and the overall SNR of the result is low.Therefore,we propose an omnidirectional beam-forming (OB) method,which is based on wavefield superposition principle by the superposition of effective wavefield in different directions to strengthen the effective signals in all directions.The workflow includes single shot record copying,angle-domain division,time-lapse parameters calculation and omnidirectional beam-forming,effective wavefield extraction and merge.The correctness and effectiveness of the OB method are proved by the processing on the single shot with random noise.Moreover,the OB method is applied on the actual seismic data,and the random noise is obviously suppressed and the SNR is increased by 6.75dB.

omnidirectional beam-forming,seismic record,signal-to-noise ratio (SNR),random noise suppression

2015-01-15;改回日期:2015-04-09。

賈海青(1987—),男,博士在讀,研究方向?yàn)榭煽卣鹪醇暗卣鹦盘?hào)處理。

姜弢(1969—),男,教授,博士生導(dǎo)師,主要研究方向?yàn)榭煽卣鹪醇暗卣鹦盘?hào)處理。

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

P631

A

1000-1441(2015)05-0521-10

10.3969/j.issn.1000-1441.2015.05.005