崔永福, 苗　 青， 黨青寧， 吳國忱， 裴廣平， 袁　源

(1．中國石油大學(xué)　地球科學(xué)與技術(shù)學(xué)院，青島　266580；

2．中國石油塔里木油田分公司　勘探開發(fā)研究院，庫爾勒　841000)

?

基于OMP算法五維數(shù)據(jù)規(guī)則化技術(shù)

崔永福1,2, 苗 青2， 黨青寧2， 吳國忱1， 裴廣平2， 袁源2

(1．中國石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，青島266580；

2．中國石油塔里木油田分公司勘探開發(fā)研究院，庫爾勒841000)

摘要：數(shù)據(jù)規(guī)則化技術(shù)是地震資料處理中的重要技術(shù)，它對改善地震數(shù)據(jù)的面元屬性、提高地震資料的信噪比和成像質(zhì)量有很大的優(yōu)勢。基于OMP算法五維數(shù)據(jù)規(guī)則化技術(shù)能綜合利用三維地震數(shù)據(jù)的“縱向、橫向、時(shí)間、偏移距、方位角”五個(gè)維度的信息，利用基于正交匹配追蹤算法(OMP)實(shí)現(xiàn)數(shù)據(jù)規(guī)則化。與常規(guī)數(shù)據(jù)規(guī)則化方法相比，該方法規(guī)則化和插值更精確，解決了連片偏移處理中塊與塊之間覆蓋次數(shù)不均造成的能量不均和偏移畫弧問題。該技術(shù)在塔北隆起LG地震工區(qū)得到成功應(yīng)用，新資料清楚地刻畫出了奧陶系碳酸鹽巖的潛山面，并能夠更加準(zhǔn)確地識(shí)別出奧陶系內(nèi)幕孔洞型“串珠”儲(chǔ)層特征。

關(guān)鍵詞：數(shù)據(jù)規(guī)則化； 五維插值； OMP算法； 孔洞型儲(chǔ)層； 碳酸鹽巖

0引言

隨著勘探開發(fā)的不斷深入，對地震成像的精度要求也越來越高。精度更高的偏移算法要求原始資料的觀測方式比較規(guī)則，覆蓋次數(shù)較高，方位角以及偏移距分布比較均勻。然而在連片資料拼接處理時(shí)，在多塊資料的交界區(qū)域，覆蓋次數(shù)往往高達(dá)上百次，而正常覆蓋只有幾十次，有些區(qū)域因地表?xiàng)l件限制、面元不同或采集方位不同導(dǎo)致出現(xiàn)空洞。這些差異都會(huì)導(dǎo)致用于疊前偏移的地震數(shù)據(jù)空間采樣不均勻[1-5]，偏移時(shí)容易產(chǎn)生假頻，帶來地震成像質(zhì)量的下降(比如成像腳印問題、成像噪音問題及偏移畫弧問題[3-4])。更重要的是，不規(guī)則的數(shù)據(jù)使得精度高的偏移方法的優(yōu)勢發(fā)揮不出來，導(dǎo)致成像振幅不保真，影響了成像的精度。因此，進(jìn)行偏前數(shù)據(jù)規(guī)則化研究勢在必行[5-6]。

實(shí)際生產(chǎn)主要應(yīng)用的數(shù)據(jù)規(guī)則化方法有面元均化技術(shù)和加權(quán)補(bǔ)償方法。面元均化技術(shù)是丟掉冗余道或者向相鄰面元借道，在一定距離內(nèi)補(bǔ)齊缺失道，其優(yōu)點(diǎn)是可以消除高、低覆蓋次數(shù)導(dǎo)致的能量差異，可以使得覆蓋次數(shù)達(dá)到基本一致，面元內(nèi)炮檢距分布均勻；缺點(diǎn)是該方法修改了原始地震道的真實(shí)坐標(biāo)，會(huì)導(dǎo)致疊前偏移出現(xiàn)空間假頻，對疊前偏移效果帶來很大的副作用。加權(quán)補(bǔ)償方法主要采用基于覆蓋次數(shù)的補(bǔ)償方法，根據(jù)面元中覆蓋次數(shù)，求取比例因子，然后對同一CMP道集應(yīng)用相同的比例因子。優(yōu)點(diǎn)是不會(huì)改變連片各工區(qū)原始覆蓋次數(shù)和地震道坐標(biāo)，能使疊前偏移的純波數(shù)據(jù)體在能量上保持一致性。此方法較面元均化方法更易實(shí)現(xiàn)，缺點(diǎn)是該方法并不能彌補(bǔ)所有缺失的偏移距[7]。

這里針對目前數(shù)據(jù)規(guī)則化方法的特點(diǎn)，研究了基于OMP算法五維數(shù)據(jù)規(guī)則化技術(shù)，將空間上不規(guī)則的樣點(diǎn)插值成規(guī)則均勻樣點(diǎn)，在一定范圍內(nèi)補(bǔ)全缺失的地震道，從而消除地震資料空間采樣不規(guī)則對成像效果造成的影響。同時(shí)內(nèi)插后的地震數(shù)據(jù)具有完整的方位角、偏移距信息。該數(shù)據(jù)規(guī)則化處理前，要對地震資料進(jìn)行預(yù)處理，如一致性拼接處理、近地表靜校正、疊前去噪、地表一致性處理等，其中一致性拼接處理為數(shù)據(jù)規(guī)則化前比較關(guān)鍵的步驟。之后采用五維數(shù)據(jù)規(guī)則化技術(shù)，以得到高質(zhì)量的道集，從而能夠進(jìn)行精確的速度拾取、準(zhǔn)確的偏移成像。

1規(guī)則化前預(yù)處理

LG地區(qū)共有A、B、C、D四塊地震工區(qū)，各次采集年代不同、觀測系統(tǒng)也不同，但互有重疊，工區(qū)的炮點(diǎn)空間分布不規(guī)則(圖1)。A、B、C三個(gè)工區(qū)的有效偏移距較小、覆蓋次數(shù)低；D工區(qū)有效偏移距較大，覆蓋次數(shù)高；A、B、D三個(gè)工區(qū)采集方位角基本一致；D工區(qū)縱橫比最大；C工區(qū)采集方位與其他工區(qū)是垂直關(guān)系(圖2)。不同區(qū)塊本身覆蓋次數(shù)不同，炮點(diǎn)分布不規(guī)則和不同區(qū)塊重疊，導(dǎo)致覆蓋次數(shù)空間變化大，不同面元內(nèi)的道集偏移距分布也不規(guī)則。在偏移之前，必須做規(guī)則化處理。

圖1　LG地區(qū)三維地震炮、檢點(diǎn)分布圖Fig．1　Source and receiver points distribution of 3D seismic acquisition in LG district

圖2　LG地區(qū)三維地震分布圖Fig．2　Offset distribution of 3D seismic data in LG district(a)偏移距分布圖；(b)方位角分布圖

不同年代采集的地震資料，采集參數(shù)不同，導(dǎo)致相位、時(shí)差、振幅能量及頻率上都會(huì)有差異[5]。不一致的相位和時(shí)差會(huì)導(dǎo)致地表一致性靜校正計(jì)算不準(zhǔn)確以及動(dòng)校正和偏移成像時(shí)不能同相疊加精確成像，不一致的振幅和子波會(huì)帶來后續(xù)保幅和精準(zhǔn)偏移成像困難，因此在進(jìn)行數(shù)據(jù)規(guī)則化之前，須先解決不同工區(qū)地震子波的空間不一致性。

根據(jù)資料特點(diǎn)，設(shè)計(jì)針對性處理流程(圖3)。通過相位一致性處理技術(shù)解決了相位和時(shí)差不一致問題(圖4(a))；地表一致性反褶積解決了子波的一致性問題，橫向分辨率趨于一致(圖5)。采用地表一致性振幅補(bǔ)償解決了相同地震工區(qū)內(nèi)振幅能量差異，但疊加剖面上仍表現(xiàn)出振幅能量差異，其振幅能量差異來源于不同工區(qū)因覆蓋次數(shù)不同導(dǎo)致。此外在野外采集施工時(shí)，由于地表變化、地表障礙等因素影響，造成炮、檢點(diǎn)分布不規(guī)則，而在克希霍夫疊前偏移成像中不僅需要滿足規(guī)則化網(wǎng)格點(diǎn)假設(shè)，而且需要滿足覆蓋次數(shù)保持一致，否則偏移成像中會(huì)出現(xiàn)較嚴(yán)重的偏移畫弧現(xiàn)象，很難做到保幅成像。因此，解決能量不均和偏移畫弧[7]等問題，需要進(jìn)行數(shù)據(jù)規(guī)則化處理。

圖3　地震數(shù)據(jù)預(yù)處理及五維數(shù)據(jù)規(guī)則化處理流程Fig．3　The flow of seismic routine processing and 5D regularization

2五維數(shù)據(jù)規(guī)則化基本原理

一般來說，數(shù)據(jù)規(guī)則化和插值的實(shí)現(xiàn)是，先將不規(guī)則網(wǎng)格點(diǎn)的數(shù)據(jù)從時(shí)間域進(jìn)行傅里葉變換到空間

頻率域，從而進(jìn)行規(guī)則網(wǎng)格點(diǎn)的數(shù)據(jù)規(guī)則化和插值，再經(jīng)過傅里葉反變換到時(shí)間域，進(jìn)行數(shù)據(jù)規(guī)則化和插值時(shí)數(shù)據(jù)和維度選擇[8]，即如何選擇和排列數(shù)據(jù)，在什么樣的維度上進(jìn)行數(shù)據(jù)規(guī)則化和插值，同樣是數(shù)據(jù)規(guī)則化和插值的關(guān)鍵。

正交匹配追蹤(OMP)算法是基于傅立葉展開式表示輸入數(shù)據(jù)以及用展開式的因子，在任何期望的空間位置來進(jìn)行數(shù)據(jù)重建的。該方法在空間頻率域中沿頻率切片計(jì)算，并且能在2、3、4、5個(gè)空間維度中應(yīng)用。偏移距的擴(kuò)展(如坐標(biāo))是由不規(guī)則離散傅立葉變換計(jì)算出。在4維和5維插值中，共中心點(diǎn)或炮點(diǎn)坐標(biāo)的擴(kuò)展，使用分級的規(guī)則離散傅立葉變換計(jì)算出。為了清晰地了解其原理，這部分插值的擴(kuò)展是用二維地震數(shù)據(jù)來說明。

假定f(xj)，(j=0…Nx-1)表示對于給定頻率的空間頻率域地震數(shù)據(jù)，其中Nx是地震道數(shù)，0≤xj≤xmax是每一道的空間坐標(biāo)，則輸入數(shù)據(jù)的擴(kuò)展形式如下：

(1)

反變換如下：

(2)

公式(1)與公式(2)的右邊能通過不規(guī)則快速傅立葉變換得到。

圖4　相位一致性疊加對比圖Fig．4　Seismic stack section phase consistency match(a)匹配前；(b)匹配后

圖5　地表一致性疊加對比Fig．5　Seismic stack section surface consistent deconvolution(a)振幅補(bǔ)償前；(b)振幅補(bǔ)償后

空間波數(shù)kl與指數(shù)l有關(guān)，kl=2πl(wèi)/xmax，(l=-Nk/2…Nk/2)，其中Nk表示空間頻率的個(gè)數(shù)。假定最大波數(shù)為固定值，增加xmax的數(shù)值會(huì)導(dǎo)致波數(shù)域更加精細(xì)的采樣率以及冗余擴(kuò)展空間。同理，反過來會(huì)提高插值計(jì)算效率。通常xmax是數(shù)據(jù)中最大偏移距的4倍。

正交匹配追蹤算法在公式(1)擴(kuò)展時(shí)每次增加一個(gè)分量，經(jīng)過m步后，誤差表示為式(3)。

(3)

其中Pl是在第l步選取的系數(shù)的指數(shù)。算法的步驟描述如下：

初始化：設(shè)定m=0,R0(xj)=f(xj) 。

迭代過程：

2)通過一個(gè)最小二乘最小化公式重新計(jì)算所有的系數(shù)。

3)通過公式(3)重新計(jì)算一個(gè)剩余值Rm+1(xj)。

4)當(dāng)‖Rm+1‖變得足夠小時(shí)或者當(dāng)達(dá)到所允許系數(shù)的最大個(gè)數(shù)時(shí)終止。

這里步驟2)是通過共軛梯度方法計(jì)算，由于共軛梯度的解需要大量正向和反向的不規(guī)則傅立葉變換，這就導(dǎo)致增加計(jì)算量，但是減少了共軛梯度迭代次數(shù)(通常只需要5到10步)，因此，執(zhí)行步驟2)是通過減少迭代次數(shù)減少了計(jì)算量。

選取數(shù)據(jù)規(guī)則化的維度也非常關(guān)鍵。常規(guī)的數(shù)據(jù)選擇有共偏移距體(相同偏移距組的數(shù)據(jù))、分方位角扇區(qū)共偏移距體(不同方位角相同偏移距的數(shù)據(jù))、共偏移距矢量體(具有相近X和Y兩個(gè)方向偏移距的數(shù)據(jù))[6]，如圖6所示。

用上述三種數(shù)據(jù)進(jìn)行數(shù)據(jù)規(guī)則化，通常的做法是分別在不同部分的數(shù)據(jù)中分別進(jìn)行縱向、橫向和時(shí)間三個(gè)維度的數(shù)據(jù)規(guī)則化和插值，這忽略了數(shù)據(jù)在不同偏移距和不同方位角之間的聯(lián)系，數(shù)據(jù)規(guī)則化和插值之后的數(shù)據(jù)不能很好地保真、保幅以及保持偏移距和方位角信息。將所有數(shù)據(jù)放在一起，在不同面元中同時(shí)考慮數(shù)據(jù)本身所包含的各種信息(縱向、橫向、時(shí)間、偏移距、方位角)[8-10]，在不同區(qū)塊、不同偏移距、不同方位角數(shù)據(jù)適當(dāng)重疊以增加更多信息，進(jìn)而進(jìn)行數(shù)據(jù)規(guī)則化和插值必然更為精確，即五維數(shù)據(jù)規(guī)則化。

五維插值數(shù)據(jù)規(guī)則化是在頻率域中運(yùn)算，一旦建立了傅里葉系數(shù)，即可運(yùn)用相應(yīng)的算法重建規(guī)則數(shù)據(jù)輸出。獲得不規(guī)則采樣數(shù)據(jù)的傅里葉系數(shù)并不是很直接，需要解決的一個(gè)問題是從振幅強(qiáng)的傅里葉系數(shù)向振幅弱的傅里葉系數(shù)漸變的過程中會(huì)有漏失。另外，為了獲得有意義的系數(shù)，輸入的主測線和聯(lián)絡(luò)測線必須有典型的偏移距，這是需要采用“超道集”的原因?！俺兰毙枰銐虼螅园璧钠凭嗄軌蛑亟?，但是也不能太大而包含太多的遠(yuǎn)數(shù)據(jù)點(diǎn)。

基于OMP算法五維數(shù)據(jù)規(guī)則化技術(shù)與常規(guī)技術(shù)的區(qū)別是：①是綜合考慮了“縱向、橫向、時(shí)間、偏移距、方位角”五個(gè)維度；②是定義輸出觀測系統(tǒng)，重新構(gòu)建道集；③是實(shí)現(xiàn)了不同面元內(nèi)炮檢距和方位角規(guī)則分布。而常規(guī)技術(shù)僅僅考慮“縱向、橫向、時(shí)間”三個(gè)維度，能夠彌補(bǔ)部分缺失的炮檢距，但不能實(shí)現(xiàn)炮檢距均勻分布，容易導(dǎo)致疊前偏移出現(xiàn)空間假頻。

3實(shí)例效果分析

圖6　數(shù)據(jù)規(guī)則化和插值時(shí)不同數(shù)據(jù)的選擇和排列示意圖Fig．6　The choice and interpolation arrange of the different seismic data in data regularization and interpolation (a)共偏移距體；(b)分方位角扇區(qū)共偏移距體；(c)共偏移距

圖7　五維數(shù)據(jù)CMP道集對比Fig．7　The contrast of CMP gathers before and after 5D regularization(a)規(guī)則化前；(b)規(guī)則化后

為了驗(yàn)證該技術(shù)在數(shù)據(jù)規(guī)則化方面的有效性，對塔里木盆地塔北隆起LG地區(qū)地震資料進(jìn)行了方法驗(yàn)證和效果測試。測試取得了較好的實(shí)際應(yīng)用效果。五維數(shù)據(jù)規(guī)則化前、后的CMP道集，五維插值數(shù)據(jù)規(guī)則化技術(shù)之后，道集偏移距分布更規(guī)則、信噪比更高(圖7)。采取五維插值數(shù)據(jù)規(guī)則化技術(shù)，不僅實(shí)現(xiàn)了整個(gè)工區(qū)具有比較一致的覆蓋次數(shù)，解決了區(qū)塊間能量不一致問題，避免了疊前偏移的畫弧現(xiàn)象，而且目的層潛山面較原偏移成像改善明顯，能量更均衡、連續(xù)性變好(圖8)，能夠準(zhǔn)確識(shí)別內(nèi)幕溶洞“串珠”[11]。

圖8　五維插值數(shù)據(jù)規(guī)則化技術(shù)對比Fig．8　The contrast of the seismic migration section(a)偏移剖面；(b)老剖面

4結(jié)論與認(rèn)識(shí)

陸上地震資料，由于受野外采集環(huán)境的限制、人為施工的誤差，不同年度采集地震資料的拼接都可能導(dǎo)致用于疊前偏移的地震數(shù)據(jù)空間采樣不均勻。作者正是基于這方面的考慮，提出了基于疊前五維數(shù)據(jù)規(guī)則化的一種新思路，充分利用地震數(shù)據(jù)的五個(gè)維度，來構(gòu)建規(guī)則道集，實(shí)際資料的應(yīng)用效果證明了該技術(shù)的實(shí)用性。

基于正交匹配追蹤算法的五維數(shù)據(jù)規(guī)則化技術(shù)，能夠有效消除偏移距、覆蓋次數(shù)、方位角等分布不規(guī)則對能量和成像效果的影響，對改善面元屬性、炮檢距分組、提高數(shù)據(jù)信噪比及偏移成像質(zhì)量優(yōu)勢明顯，能夠解決不同區(qū)塊和拼接區(qū)之間偏前覆蓋次數(shù)不均，偏移剖面上振幅能量不均和成像失真等問題；同時(shí)該方法為陸上采集炮、檢點(diǎn)靈活布置提供了條件。

該技術(shù)同樣適用于3維或4維數(shù)據(jù)規(guī)則化插值，即可以去掉偏移距和方位角信息。但是更多空間維度的插值算法有很多的優(yōu)點(diǎn)，特別是同時(shí)在5個(gè)地震數(shù)據(jù)維度進(jìn)行插值，能夠很好地通過修正振幅和相位變量來預(yù)測缺失的地震數(shù)據(jù)。

參考文獻(xiàn)：

[1]張紅梅，劉洪. 基于稀疏離散τ-p變換的疊后地震道外推 [J]. 物探化探計(jì)算技術(shù)，2007，29(1)：1-6.

ZHANG H M，LIU H. Sparseness discrete τ-p transform in poststack seismic trace extrapolaion[J].Computing techniques for geophysical and geochemical exploration，2007，29(1)：1-6.(In Chinese )

[2]徐輝，張軍華. 面元細(xì)分處理資料的分辨率定量評價(jià)[J]. 物探化探計(jì)算技術(shù)，2010，32(1)：59-63.

XUE H，ZHANG J H.The quantitative evaluation of resolution for seimic data processing by bin-dividing[J].Computing techniques for geophysical and geochemical exploration，2010，32(1)：59-63.(In Chinese)

[3]石穎，張振，王建民，等. 地震數(shù)據(jù)反假頻規(guī)則化方法研究[J]. 地球物理學(xué)進(jìn)展，2013，28(1)：0250-0256.

SHI Y，ZHANG Z，WAGN J M，et al. Investigation on de-alias regularization approach for seismic data[J]. Progress in Geophys，2013，28(1):0250-0256.(In Chinese )

[4]劉禮農(nóng)，周立宏，朱振宇，等. 復(fù)雜觀測系統(tǒng)下的三維波動(dòng)方程疊前深度偏移[J]. 地球物理學(xué)進(jìn)展，2002，17(4)：658-662.

LIU L N，ZHOU L H，ZHU Z Y，et al.3D wave equation prestack depth migration with complex Layouts [J]. Progress in Geophys，2002，17(4)：658-662.(In Chinese )

[5]高建軍，陳小宏，王芳芳，等.不規(guī)則地震道數(shù)據(jù)規(guī)則化重建方法研究[J]. 地球物理學(xué)進(jìn)展，2011，26(3)：983-991.

GAO J J，CHEN X H，WANG F F，et al. Study on regularized reconstruction of uneven seismic traces. [J]. Progress in Geophys，2011，26(3)：983-991.(In Chinese )

[6]王棣，馬秀紅，崔興福，等. 偏移距規(guī)則化技術(shù)在疊前時(shí)間偏移中的應(yīng)用[J]. 勘探地球物理進(jìn)展，2009，32(1)：44-47.

WANG D，MA X H，CUI X F，et al. Application of prestack offset regularization prestack time migration[J]. Progress in Exploration Geophysics ，2009，32 ( 1)： 44- 47. (In Chinese)

[7]趙忠華，初海紅，陳志德. 頻率-空間域疊前數(shù)據(jù)規(guī)則化技術(shù)研究及其應(yīng)用[J]. 油氣藏評價(jià)與開發(fā)，2013，3(2)：15-20.

ZHAO Z H，CHU H H，CHEN Z D. Research on frequency-spatial domain prestack data regularization and its application[J].Reservoir Evaluation and Development，2013，3(2)：15-20.(In Chinese)

[8]KAWIN NIMSAILA， PING WANG.Fast progressive sparse Tau-P transform for regularization of spatially aliased seismic data[C]. The 76thEAGE conference and exhibition，2014：230-234.

[9]G.POOLE， P.HERMANN.Multi-dimensional data regularization for modern acquisition geometries[C]. The 69thEAGE conference and exhibition incorporating SPE EUROPEC， 2007：70-75.

[10]SHENG XU， YU ZHANG， GILLES LAMBARE.Anti-leakage Fourier transform for seismic data regularization practice in higher dimensions[J]. Geophysics， 2010，70(4)：87-95.

[11]高彩霞.波動(dòng)方程疊前成像數(shù)據(jù)規(guī)則化技術(shù)研究與應(yīng)用[J]. 石油天然氣學(xué)報(bào)，2010，32(6)：271-273.

GAO C X. Technical Study and Application of Wave Equation Prestack Imaging Data Rule[J]. Journal of Oil and Gas Technology，2010，32(6)：271-273.(In Chinese)

5D regularization processing based on OMP arithmetic in seismic data CUI Yong-fu1,2， MIAO Qing2， DANG Qing-ning2

WU Guo-chen1， PEI Guang-ping2， YUAN Yuan2

(1. School of Geosciences，China University of Petroleum(Huadong), Qingdao286580， China；2. Research Institute of Exploration and Development，Tarim Oilfield Company， PetroChina， Korla841000， China)

Abstract:Regularization of seismic data is one of the important technologies in seismic data processing. It can improve bin properties, seismic S/N ratio and seismic migration imaging. It is different from the conventional data regularization methods, 5D regularization based on OPM arithmetic comprehensively utilizes 3D seismic data in five dimensions including longitudinal coordinate, transverse coordinate, time, offset, and azimuth angle, leading to the regularization and interpolation more precise. This method is successfully applied in LunGu 7 district of Tabei lift in Tarim basin. The results show that this technology solves the problem of uneven amplitude and migrating arcing caused by fold uneven between contiguous blocks in PSTM processing. The new seismic data can describe clearly the surface of burial hill of Ordovician dolomite and the "string" characteristics in interior burial hill.

Key words:regularization； 5D interpolation； orthogonal matching pursuit； cavity reservoirs； carbonate

中圖分類號：P 631.4

文獻(xiàn)標(biāo)志碼：A

DOI：10.3969/j.issn.1001-1749.2016.01.10

文章編號：1001-1749(2016)01-0067-07

作者簡介：崔永福(1978-)，男，博士，高級工程師，主要從事石油物探工作， E-mail：cuiyongfu-tlm@petrochina.com.cn。

基金項(xiàng)目：國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2013CB228604)； 國家油氣重大專項(xiàng)(2011ZX05046, 2011ZX05019-003, 2011ZX05009-003-004)

收稿日期：2014-19-09改回日期：2015-05-15