王飛，孫亞杰，裴金梅，宋建國，李文建

(東方公司 青海物探處,甘肅 敦煌 736202)

0 引言

當前，地震勘探已進入面向地質(zhì)目標體的綜合地震勘探階段[1]，地質(zhì)目標趨于更復雜、地質(zhì)需求越來越高，勘探開發(fā)對地震勘探數(shù)據(jù)的信噪比、分辨率、保真度和成像精度等提出了更高的要求，最為迫切需要就是提高地震勘探精度[2-4]。在構(gòu)造勘探時期奉為經(jīng)典的大組合接收隨著高密度采集技術的推廣應用逐漸不再適用，另外隨著野外施工物理點數(shù)增多，大組合接收施工難度大，勘探成本過大[5-8]。

近年來青海探區(qū)進行了大量的單點檢波器采集試驗生產(chǎn)工作，截至2019年青海探區(qū)共完成了18個單點地震采集項目，單點數(shù)字檢波器小道距、高覆蓋、不組合的野外采集方式逐漸得到應用。

在地震數(shù)據(jù)處理中使用地震道組合的方法, 把單點檢波器輸出的信號組合疊加到一起，和采集時使用檢波器面積組合相比，能夠避免組合的低通濾波特性和大組合基距帶來的地表靜校正量致使信號發(fā)生畸變。采集時使用單點接收，每個單獨的檢波器接收到的信號單獨記錄，這種方法具有明顯優(yōu)勢[9-10]。

1 理論分析

1.1 道距與假頻

地震波的傳播，既有空間的變化，也有時間的變化[9],用y(t,x)表示地震信號，其中t為時間變量，x為空間變量。則二維傅里葉變換為：

(1)

式中：f表示頻率，代表時間信號y(t)周期T的倒數(shù)，即f=1/T；k表示波數(shù)，代表信號y(x)波長λ的倒數(shù)，即k=1/λ。Y(f,k)即稱為地震信號y(t,x)的頻率—波數(shù)譜，簡稱頻波譜或FK譜。通過FK譜轉(zhuǎn)換，一些在時間—空間域無法分離的地震信號可以在頻率—波數(shù)域區(qū)分開，由此為去噪提供了基礎。

在實際地震采集生產(chǎn)過程中，有效反射由于視速度高，很容易滿足不出現(xiàn)假頻現(xiàn)象。但是低速的干擾波，目前的地震勘探方法很難控制假頻出現(xiàn)[10-14]。如圖1所示，從柴達木盆地尖頂山地區(qū)某條測線單炮FK譜可以看出，即使是10 m小道距，低速干擾波也存在假頻折疊，隨著道距擴大，這種假頻折疊逐漸嚴重，道距達到30 m時，這種折疊與有效反射混淆在一起，大道距產(chǎn)生的假頻折疊現(xiàn)象對處理中的噪聲壓制非常不利。

a、b、c、d、e—道間距分別為10,20,30,40,50 m

1.2 道距與疊加的關系

道組合其原理與水平疊加是相同的，參照多次疊加的特性做如下假設，需要組合的每一組數(shù)據(jù)組內(nèi)有效波在波形和能量上是相同的,只考慮反射時間的差別，公式推導如下[15-17]。

設動校正后剩余時差為零的振動函數(shù)為f(t), 頻譜為g(jω), 則動校正后剩余時差為Δti的振動函數(shù)為f(t-Δti)。就是說動校正后各道波形分別為:f(t-Δt1),f(t-Δt2),…,f(t-Δtn), 則水平疊加輸出為：

(2)

其道組合后的頻譜為：

G(jω)=g(jω)e-jωΔt1+g(jω)e-jωΔt2+

如果令

則有

G(jω)=K(jω)g(jω)。

(3)

上式表明道組合相當于一個線性濾波器,對波形的改造作用可以用K(jω)反映出來, 即K(jω)的模K(ω)反映了多次疊加的振幅頻率特性：

對反射波來說, 最理想的情況是它的剩余時Δti差為0, 則K(ω)為n, 表明組合后反射波增強了n倍; 對于其他干擾波,Δti不為0,K(ω)一定小于n,則組合對于干擾波起到相對削弱的作用。

為了便于分析不同組合道數(shù)的垂直疊加效果,把K(ω)除以組合次數(shù)n, 即得:

進行組合時,Δti可用正常時差表示：

(5)

式(5)中，有效波反射時間t0取0.8 s，動校正速度vd取2 000 m/s，道間距離xi和時差Δti關系如圖2所示;有效波波至時間t0取0.8 s，道間距離xi取值為30 m，則道間時差Δti與動校正速度vd關系如圖3所示。

圖2 道間時差與道間距的關系

圖3 不同動校速度與時差關系

觀察圖2、圖3,可以得出以下兩點結(jié)論:①隨著組合道數(shù)的增加,正常時差對組合的影響也增大；②動校正速度減小,正常時差對組合的影響增大。

2 處理流程

通過前面分析得知，小道距在壓制視速度較低的噪聲有明顯優(yōu)勢。在做道組合時，每一組的道內(nèi)時差也對組合效果有很大的影響[18]。

高密度單點接收地震采集數(shù)據(jù)處理的基本思路是：首先要分離線性噪聲和消除野外靜校正對數(shù)據(jù)的影響；其次，做好動校正；最后，進行道組合處理。這樣就能夠充分利用小道距壓制低速噪聲，有效壓制線性噪聲，消除組合道道間時差對信噪比的影響，達到最優(yōu)組合，其具體流程如圖4。

圖4 處理流程

3 處理效果

柴達木盆地尖頂山地區(qū)開展了單點檢波器和常規(guī)檢波器的對比試驗，采用了單點SG5、單點SN5、常規(guī)SN8-10 Hz這3種型號的檢波器，其中SG5、SN5檢波器道距10 m，SN8-10 Hz檢波器“口”字組合，道距30 m，3種檢波器并列布設。采用可控震源激發(fā)，3種檢波器同時接收。本文選擇SN8-10 Hz和SG5兩種檢波器接收的數(shù)據(jù)處理分析。

從單炮數(shù)據(jù)上看，常規(guī)SN8-10 Hz 10個檢波器面積組合壓制噪聲效果好，比SG5單個檢波器接收的數(shù)據(jù)信噪比高(圖5)。

a—常規(guī)檢波器數(shù)據(jù);b—單點檢波器數(shù)據(jù)

該條測線設計面元為15 m，故在處理過程中需要對10 m道距數(shù)據(jù)進行道組合，使組合后道距為30 m。由于小道距有利于干擾波的分離與壓制，在處理過程中，做了先組合再去噪和先去噪后組合試驗。通過試驗結(jié)果看，先去噪再組合去噪效果明顯優(yōu)于先組合再去噪，特別是對于低速干擾壓制尤為明顯(圖6)。

a—先組合后去噪;b—先去噪后組合

通過前面分析，低降速帶厚度或橫向速度變化都會造成微小的靜態(tài)時移, 該時移會影響到處理中道組合的效果，因此組合前要做靜校正處理。

通過圖7可以看出，靜校正前后道集同相軸連續(xù)性有一定影響，靜校正后同相軸更連續(xù)，一致性也明顯變好，更有利于道組合時同相疊加。

a—靜校正前;b—靜校正后

道集動校正，會最大限度減小道間時差(圖8),能夠提高組合的同向疊加效果。通過試驗對比，先做動校正再組合，其同相軸更連續(xù)，信噪比較高, 組合效果比先做組合再做動校正好(圖9)。

圖8 動校正前(b)和動校正后(b)道集

a—先做動校正再組合;b—先組合再做動校正后

由疊加剖面對比分析，先做靜校正和動校正再進行道組合，剖面連續(xù)性明顯提高(圖10)。

a—動靜校正后組合;b—組合后做動靜校正

常規(guī)SN8-10 Hz 10個檢波器面積組合接收的數(shù)據(jù)和單點SG5檢波器接收的數(shù)據(jù)從單炮品質(zhì)上看，信噪比存在很大差別，使用常規(guī)檢波器面積組合，會有一定的壓制噪聲作用。常規(guī)10個檢波器面積組合接收的數(shù)據(jù)和單點檢波器接收的數(shù)據(jù)對比，信噪比明顯偏高。但是通過優(yōu)化處理流程，在適合的步驟中增加道組合處理，從最終疊加剖面看，兩者剖面品質(zhì)基本相當，在某些細節(jié)上，單點檢波器數(shù)據(jù)剖面比常規(guī)檢波器數(shù)據(jù)對細節(jié)刻畫更清晰，如圖11。

a—常規(guī)流程常規(guī)檢波器數(shù)據(jù);b—本文流程單點檢波器數(shù)據(jù)

4 結(jié)論

高分辨率勘探需要深度域和空間域都有高精度采樣的數(shù)據(jù), 小道距、不組合、高覆蓋的野外采集方式順應了高分辨率勘探數(shù)據(jù)采集的特點。高密度單點檢波器地震數(shù)據(jù)道組合在青海柴達木盆地大部分地區(qū)實施效果良好。由于道距減小, 在數(shù)據(jù)處理過程中通常要進行道組合處理。

1)地震數(shù)據(jù)處理中的道組合最好是在做好噪聲壓制、靜校正、動校正之后進行。

2)單點接收野外采集方式靈活, 數(shù)據(jù)處理時可根據(jù)干擾波分布特征靈活地選擇組合方式和組合參數(shù)，這也是道組合下一步需要嘗試的工作。