李敏杰

（中國石化中原油田分公司物探研究院，河南 濮陽 457001）

數(shù)據(jù)規(guī)則化技術(shù)在三維地震老資料疊前偏移中的應用

李敏杰

（中國石化中原油田分公司物探研究院，河南 濮陽 457001）

數(shù)據(jù)規(guī)則化技術(shù)是地震資料處理中的重要技術(shù)，它對改善地震數(shù)據(jù)的面元屬性，提高地震資料處理質(zhì)量具有重要意義。進行疊前時間偏移處理時，更多關注的是偏移速度場精度及偏移參數(shù)對偏移成像效果的影響，忽視了不規(guī)則空間采樣的疊前數(shù)據(jù)對疊前偏移成像效果的影響。文中在分析疊前數(shù)據(jù)規(guī)則化常用技術(shù)方法優(yōu)缺點的基礎上，研究了反假頻傅里葉變換的疊前數(shù)據(jù)規(guī)則化方法，實現(xiàn)了陸地地震勘探資料不規(guī)則空間采樣的疊前數(shù)據(jù)規(guī)則化，滿足了疊前時間偏移技術(shù)對輸入數(shù)據(jù)的要求，在濟源后鄧地區(qū)稀疏三維地震老資料重新處理中取得了較好的應用效果。

三維地震；數(shù)據(jù)規(guī)則化；稀疏采樣；疊前時間偏移

0 引言

隨著油氣勘探開發(fā)程度的深入，地質(zhì)、解釋人員對地震資料的要求不斷提高,常規(guī)疊后時間偏移資料已難以滿足油氣精細勘探的需要。疊前時間偏移［1-3］方法從理論上繞開了水平層狀介質(zhì)和剖面為自激自收記錄的假設條件，與實際地質(zhì)情況相符，偏移結(jié)果更加真實可靠，與疊前深度偏移［4］相比，速度模型的建立和修正相對簡單,偏移結(jié)果對速度模型的依賴性較小，計算效率高，現(xiàn)已廣泛應用于復雜構(gòu)造偏移成像處理中［5］。但在該技術(shù)應用過程中發(fā)現(xiàn)，有些情況下即使進行了精細預處理，建立了較精確的偏移速度場，選取了最佳的偏移處理參數(shù)（如偏移孔徑、偏移傾角和反假頻參數(shù)），也難以取得滿意的偏移成像效果，這在濟源后鄧地區(qū)三維地震老資料疊前偏移處理中得到了印證。

1 濟源后鄧地區(qū)地震資料概況

濟源凹陷后鄧地區(qū)三維地震資料是由中原油田于1989年冬—1990年春采集完成的。野外地震資料采集采用6線4炮線束狀觀測系統(tǒng)，儀器總道數(shù)240道，檢波線距200 m，檢波間距50 m；炮線距100 m，炮間距100 m；覆蓋次數(shù)2×10，面元25 m×50 m。滿覆蓋面積84.7 km2，總炮數(shù)6 688炮，總道數(shù)550 400道。原始資料的特點是炮間距、檢波間距、炮線距及檢波線距較大，覆蓋次數(shù)少（20次），空間采樣稀疏。由于工區(qū)障礙物的存在，炮點、檢波點位置偏離原來采集設計網(wǎng)格點的現(xiàn)象較多，導致偏移距分布、覆蓋次數(shù)分布不均勻，振幅能量空間變化較大。利用這樣的不規(guī)則數(shù)據(jù)進行疊前時間偏移處理，很容易產(chǎn)生偏移噪聲（偏移畫?。?，致使偏移效果不好，振幅能量橫向變化大，難以滿足精細勘探的地質(zhì)要求。通過對原始地震資料的空間屬性分析，研究并應用疊前數(shù)據(jù)規(guī)則化技術(shù)，解決了該區(qū)上述問題，滿足了疊前偏移對輸入數(shù)據(jù)的要求，取得了較好的應用效果。

2 疊前數(shù)據(jù)規(guī)則化主要方法和技術(shù)

數(shù)據(jù)規(guī)則化技術(shù)最初是為了恢復空缺地震道記錄、滿足多道處理算法（多次波衰減、偏移等）的要求而提出來的。隨著疊前時間偏移的廣泛應用，數(shù)據(jù)規(guī)則化技術(shù)也日趨成熟。常用的疊前數(shù)據(jù)規(guī)則化技術(shù)主要有基于拋物線拉冬變換的數(shù)據(jù)規(guī)則化、DMO傾角時差校正的數(shù)據(jù)規(guī)則化、偏移距分布密度的數(shù)據(jù)規(guī)則化和傅里葉變換的數(shù)據(jù)規(guī)則化等。

2.1 基于拋物線拉冬變換的數(shù)據(jù)規(guī)則化方法

基于拋物線拉冬變換的數(shù)據(jù)規(guī)則化方法的特點，都是基于CMP道集的雙曲線走時，采用雙曲線或拋物線拉冬正反變換、多次迭代來實現(xiàn)疊前數(shù)據(jù)規(guī)則化。該方法計算量比較大，而且遠偏移距道的數(shù)據(jù)重構(gòu)效果不理想［6-10］。此外，這種方法是在假設水平層狀介質(zhì)的情況下進行的，當?shù)叵掠袃A斜地層時，共面元內(nèi)的各道不再滿足拋物線時距關系，經(jīng)過數(shù)據(jù)規(guī)則化的疊前數(shù)據(jù)精度受到懷疑。

2.2 基于DMO傾角時差校正的數(shù)據(jù)規(guī)則化方法

該方法是根據(jù)DMO運算可逆的這一特點來實現(xiàn)的。DMO傾角時差校正本來是用來消除地層傾角對疊加成像影響的，但是，它的實現(xiàn)是以波場的傳播理論為基礎，其運算是可逆的，即DMO的效應可以通過反DMO（即DMO-1）來予以消除，利用這一特點便可以實現(xiàn)地震數(shù)據(jù)的規(guī)則化。由于DMO既可以在時間-空間域?qū)崿F(xiàn)，也可以在頻率-波數(shù)域?qū)崿F(xiàn)，因此，地震數(shù)據(jù)規(guī)則化也就可以在這2個域里實現(xiàn)。前者因采用積分法，在積分過程中由于數(shù)據(jù)拉伸容易產(chǎn)生假頻，后者雖然可以避免假頻，但要求規(guī)則的觀測系統(tǒng)。

1991年，D.Hale［11］通過在DMO積分運算中引入抗假頻濾波器，有效避免了積分過程中產(chǎn)生的假頻現(xiàn)象。1993年，J.L.Black等［12］提出的保幅DMO算子修正了1984年D.Hale提出的F-K域DMO算子，實現(xiàn)了保幅反DMO，為疊前數(shù)據(jù)規(guī)則化的振幅保持奠定了良好基礎。2002年，辛可鋒等［13］在借鑒1996年A.Canning等［14］提出的利用DMO實現(xiàn)三維地震數(shù)據(jù)的規(guī)則化和1998年B.Biondi等［15］提出的三維地震數(shù)據(jù)方位角校正的基礎上，從方位角校正的角度入手，利用DMO和反DMO相結(jié)合的方法，即采用抗假頻的積分法DMO，并在F-K域?qū)崿F(xiàn)反DMO，對疊前數(shù)據(jù)進行規(guī)則化，目前該方法具有抗假頻和振幅保持，以及運算效率高等特點。

2.3 基于偏移距分布密度的數(shù)據(jù)規(guī)則化方法

該方法根據(jù)疊前數(shù)據(jù)每個偏移距的空間分布密度，確定適當?shù)恼穹壤蜃樱瑢λ械墓财凭鄶?shù)據(jù)進行振幅加權(quán)處理，即對分布密度大的偏移距數(shù)據(jù)的振幅乘以較小的比例因子，分布密度小的偏移距數(shù)據(jù)的振幅乘以較大的比例因子，從而實現(xiàn)地震數(shù)據(jù)振幅一致性處理［16］。該方法的特點就是消除由于偏移距分布不均而造成的振幅能量差異，疊前數(shù)據(jù)量在規(guī)則化前后沒有發(fā)生任何改變，適合空間高密度采樣的地震數(shù)據(jù)規(guī)則化。

2.4 基于傅里葉變換的數(shù)據(jù)規(guī)則化方法

法國CGG公司GEOVATION 6200地震處理系統(tǒng)中基于二維傅里葉重建的數(shù)據(jù)規(guī)則化方法，是沿著三維數(shù)據(jù)的InLine和CrossLine兩個方向，利用二維正、反傅里葉變換來實現(xiàn)三維疊前數(shù)據(jù)規(guī)則化。步驟為：首先，利用二維傅里葉分解法對輸入數(shù)據(jù)進行變換（輸入數(shù)據(jù)的位置不發(fā)生改變）；然后，進行傅里葉反變換把數(shù)據(jù)投影到一個規(guī)則的網(wǎng)格上。該方法采用標準傅里葉變換和去假頻傅里葉變換，可以按照不同數(shù)據(jù)集（炮集、檢波點集、偏移距域、CMP道集等）對不同類型的數(shù)據(jù)（陸地地震數(shù)據(jù)、海洋地震數(shù)據(jù)、海底電纜OBC數(shù)據(jù)、多分量數(shù)據(jù)等）進行規(guī)則化處理。本文就是采用這種方法對疊前不規(guī)則偏移距數(shù)據(jù)進行規(guī)則化。

3 應用效果

濟源后鄧三維工區(qū)內(nèi)村莊、道路、工業(yè)區(qū)分布密集，這些地面障礙物的存在使得炮點和檢波點位置不能按照采集設計中的規(guī)則網(wǎng)格布設（見圖1），表現(xiàn)為炮點（圖中藍色）、檢波點（圖中草綠色）空間位置分布不均勻，導致了最小、最大偏移距（見圖2），覆蓋次數(shù)（見圖3）和地下CMP點空間位置（見圖4）分布不均勻。采用反假頻二維傅里葉正反變換，實現(xiàn)疊前數(shù)據(jù)規(guī)則化，使野外跨越障礙物而導致的近偏移距資料缺失的現(xiàn)象得到一定改善，面元內(nèi)各反射點處于面元中心位置，覆蓋次數(shù)均勻，振幅能量得到統(tǒng)一。經(jīng)過疊前時間偏移處理后，偏移畫弧現(xiàn)象得以消除，地震反射波振幅屬性趨于一致。

圖5為數(shù)據(jù)規(guī)則化后的地下CMP點空間位置分布，可以看出數(shù)據(jù)規(guī)則化后地下CMP面元分布均勻。圖6為數(shù)據(jù)規(guī)則化前、后偏移距組（500～600 m）內(nèi)的CMP面元空間分布，可以看出規(guī)則化后數(shù)據(jù)空洞現(xiàn)象得以較好消除。

圖7a、圖7b分別為數(shù)據(jù)規(guī)則化前、后的CMP道集，規(guī)則化后的CMP道集數(shù)據(jù)反射波同相軸雙曲線特征更加清楚，信噪比得到提高。圖8為數(shù)據(jù)規(guī)則化前、后的頻譜分析，數(shù)據(jù)規(guī)則化前、后頻譜形態(tài)基本沒有改變，說明該方法具有較高的保真性。

圖7 CMP道集

圖8 數(shù)據(jù)規(guī)則化前、后頻譜

圖9為數(shù)據(jù)規(guī)則化前、后的疊前時間偏移純波剖面，通過疊前數(shù)據(jù)規(guī)則化處理后，較好地消除了由于偏移距分布不均造成的振幅強弱變化明顯的假象。與數(shù)據(jù)規(guī)則化前的疊前時間偏移剖面相比，數(shù)據(jù)規(guī)則化后的偏移剖面斷面波得到較好收斂，斷層更加清楚，角度不整合等地質(zhì)現(xiàn)象更加明顯；反射波同相軸連續(xù)性增強，振幅能量基本趨于一致，剖面上振幅條帶狀變化的現(xiàn)象得以消除。成果資料滿足了地震解釋、綜合地質(zhì)研究和儲層預測對地震資料的要求。

圖9 疊前時間偏移純波剖面

4 結(jié)束語

在分析總結(jié)常用的疊前數(shù)據(jù)規(guī)則化方法優(yōu)缺點的基礎上，研究并應用了反假頻傅里葉變換疊前數(shù)據(jù)規(guī)則化方法，較好地解決了不規(guī)則陸地地震勘探資料的數(shù)據(jù)規(guī)則化難題，滿足了疊前偏移方法對輸入數(shù)據(jù)空間均勻采樣的要求，取得了較好的應用效果。此方法不僅可以解決陸地地震勘探數(shù)據(jù)不規(guī)則化問題，而且對海洋地震資料處理也有借鑒作用。

［1］ 曹延軍，游洪文，林玉英，等.疊前時間偏移技術(shù)及應用［J］.斷塊油氣田，2004，11（5）：21-22.

［2］ 楊春梅，顏軍，陳天鑄，等.克希霍夫疊前時間偏移技術(shù)及應用［J］.斷塊油氣田，2007，14（6）：21-23.

［3］ 朱武.西非深海地區(qū)Kirchhoff偏移技術(shù)的研究與應用［J］.斷塊油氣田，2008，15（2）：37-39.

［4］ 林明玉，匡斌.疊前深度偏移中的旅行時并行計算技術(shù)［J］.斷塊油氣田，2003，10（4）：12-14.

［5］ 郭朝斌，侯春麗，曾紅，等.高精度三維深度域成像技術(shù)在文衛(wèi)濮結(jié)合部的應用［J］.斷塊油氣田，2002，9（4）：23-25.

［6］ Nurul Kabir M M，Verschuur D J.應用“拋物線Radon變換”恢復空缺道記錄：上［J］.趙希剛，譯.國外地質(zhì)勘探技術(shù)，1996（4）：5-10.

［7］ Nurul Kabir M M，Verschuur D J.應用“拋物線Radon變換”恢復空缺道記錄：下［J］.趙希剛，譯.國外地質(zhì)勘探技術(shù)，1996（5）：6-11.

［8］ 王維紅，高紅偉，劉洪道.均衡拋物線Radon變換法地震道重建［J］.石油地球物理勘探，2005，40（5）：518-522，560.

［9］ 王維紅，劉洪.拋物Radon變換法近偏移距波場外推［J］.地球物理學進展，2005，2（2）：289-293.

［10］王維紅，裴江云，張劍鋒.加權(quán)拋物Radon變換疊前地震數(shù)據(jù)重建［J］.地球物理學報，2007，50（3）：851-859.

［11］Hale D.A nonaliases integral method for dipmoveout［J］.Geophysics，1991，56（6）：795-805.

［12］Black J L，Schleicher K L，Zhang L.True-amplitude imaging and dip moveout［J］.Geophysics，1993，58（1）：47-66.

［13］辛可鋒，王華忠，王成禮，等.疊前地震數(shù)據(jù)的規(guī)則化［J］.石油地球物理勘探，2002，37（4）：311-317.

［14］Canning A，Gerald G H F.Regularizing 3-D data sets with DMO［J］. Geophysics，1996，61（4）：1103-1114.

［15］Biondi B，F(xiàn)omel S，Chemingui N.Azimuth moveout for 3-D perstack imaging［J］.Geophysics，1998，63（5）：574-588.

［16］王棣，馬秀紅，崔興福，等.偏移距規(guī)則化技術(shù)在疊前時間偏移中的應用［J］.勘探地球物理進展，2009，32（1）：44-47.

(編輯 朱麗)

Application of data regularization in prestack migration of old 3D seismic data

Li Minjie
(Geophysical Research Institute,Zhongyuan Oilfield Company,SINOPEC,Puyang 457001,China)

Data regularization is a key technology in seismic data processing which has great importance to improve the attribute of bin size and to enhance the processing quality of seismic data.Technician pay more attention to the effect of migration velocity field accuracy and migration parameters on migration imaging result when prestack time migration is carried on,but the effect of prestack data on prestack migration imaging result in irregular space sampling is neglected.Based on the analysis of the advantages and disadvantages of common prestack data regularization technology,the prestack data regularization method of anti-aliasing Fourier transform is studied and the prestack data regularization of irregular space sampling in land seismic exploration is realized,which meets the requirement of prestack time migration on input data and obtains a good application result in the reprocessing of old 3D seismic data in Houdeng Area of Jiyuan.

3D seismics;data regularization;sparse sampling;prestack time migration

TE132.1；P631.4+43

A

10.6056/dkyqt201204006

2012-02-15；改回日期：2012-05-19。

李敏杰，男，1964年生，高級工程師，1988年畢業(yè)于石油大學（華東）物探專業(yè)，長期從事地震資料處理方法研究及應用工作。E-mail：zycclmj@163.com。

李敏杰.數(shù)據(jù)規(guī)則化技術(shù)在三維地震老資料疊前偏移中的應用［J］.斷塊油氣田，2012，19（4）：430-433.

Li Minjie.Application of data regularization in prestack migration of old 3D seismic data［J］.Fault-Block Oil&Gas Field，2012，19（4）：430-433.