起伏地表采集數(shù)據(jù)的三維直接疊前時(shí)間偏移方法

張 浩，張劍鋒

1中國科學(xué)院地質(zhì)與地球物理研究所中國科學(xué)院油氣資源研究重點(diǎn)實(shí)驗(yàn)室，北京 100029 2中國科學(xué)院研究生院，北京 100049

起伏地表采集數(shù)據(jù)的三維直接疊前時(shí)間偏移方法

張 浩1，2，張劍鋒1

1中國科學(xué)院地質(zhì)與地球物理研究所中國科學(xué)院油氣資源研究重點(diǎn)實(shí)驗(yàn)室，北京 100029 2中國科學(xué)院研究生院，北京 100049

提出一種可對起伏地表采集的三維地震資料直接進(jìn)行偏移成像的疊前時(shí)間偏移方法和流程．它用兩個(gè)等效速度描述近地表和上覆層對地震波傳播的影響，可對炮、檢點(diǎn)不在同一水平面的三維地震資料直接進(jìn)行疊前時(shí)間偏移處理．該方法不對近地表地震波傳播做垂直出、入射假定，因此可適應(yīng)高速層出露等不存在明顯低、降速帶情況．描述近地表和上覆層的兩個(gè)等效速度參數(shù)可依據(jù)偏移道集的同相軸是否平直來確定，避免了確定近地表速度的困難；而對已知近地表速度的情況，則可進(jìn)一步修正近地表速度，獲得更好的成像效果．用三維起伏地表的理論數(shù)據(jù)和中國東部某工區(qū)實(shí)際數(shù)據(jù)驗(yàn)證了所發(fā)展方法和處理流程的有效性和實(shí)用性．

起伏地表，三維成像，近地表速度，疊前時(shí)間偏移

1 引 言

隨著我國陸上油氣勘探向西部和南方地區(qū)的深入，復(fù)雜地表已成為制約地震資料處理的關(guān)鍵因素．由于地表起伏較大、高速的老地層出露，現(xiàn)行的以靜校正處理為基礎(chǔ)的陸上地震資料處理方法和流程遇到困難．浮動基準(zhǔn)面技術(shù)［1－2］可解決地表起伏較大問題，但高速層出露導(dǎo)致垂直入射和出射的靜校正應(yīng)用條件不成立，繼續(xù)使用靜校正方法將帶來較大的誤差［3］．基于波場延拓的基準(zhǔn)面重建方法［4］可同時(shí)解決地表起伏和高速層出露問題，但這一方法是以獲得準(zhǔn)確的近地表速度模型為基礎(chǔ)的；由于起伏地表采集的地震資料信噪比較低，獲得準(zhǔn)確的近地表速度模型是一個(gè)相當(dāng)困難的任務(wù)［5］．此外，基于波場延拓的基準(zhǔn)面重建方法導(dǎo)致重建數(shù)據(jù)的地震道數(shù)激增，極大地增加了后續(xù)偏移處理的計(jì)算量．

就實(shí)際中遇到的大量地質(zhì)目標(biāo)而言，盡管地表復(fù)雜、斷層等構(gòu)造斷裂發(fā)育，但除近地表以外的上覆地層的速度橫向變化并不是很劇烈，這就給疊前時(shí)間偏移的應(yīng)用提供了可能性．疊前時(shí)間偏移的主要優(yōu)點(diǎn)是只用一個(gè)等效速度（即均方根速度）描述地震波傳播，可對一類傾角、斷層較為復(fù)雜，但速度橫向變化不是很劇烈的構(gòu)造較好的成像［6］．獨(dú)立的、單一等效速度使得疊前時(shí)間偏移的速度建模變得簡單．本文進(jìn)一步發(fā)展疊前時(shí)間偏移方法，通過引入第二個(gè)等效參數(shù)描述近地表地震波傳播，發(fā)展可對炮、檢點(diǎn)不在同一水平面的三維地震資料直接進(jìn)行疊前時(shí)間偏移處理的成像方法，以此避免高速層出露和近地表速度建模這兩個(gè)困難．

現(xiàn)行的直接對起伏地表采集數(shù)據(jù)進(jìn)行偏移、避免使用靜校正技術(shù)的直接偏移方法主要是基于疊前深度偏移算法［7－8］．這一方法需要準(zhǔn)確的近地表速度和地下構(gòu)造的層速度模型，而這兩類速度又相互影響，與海上或水平地表情況下確定深度偏移的層速度模型相比，這類方法的速度建模更加困難．文獻(xiàn)［9］盡管在計(jì)算地震波走時(shí)時(shí)引入了平均速度和均方根速度，但這兩個(gè)速度是由層速度計(jì)算得到的，它們是成像點(diǎn)和炮點(diǎn)（或檢波點(diǎn)）的函數(shù)；由于在同一成像點(diǎn)有多個(gè)均方根速度存在（隨炮點(diǎn)或檢波點(diǎn)變化），這一速度是不能如疊前時(shí)間偏移那樣通過簡單的掃描方式得到的，只能通過先求取層速度模型，進(jìn)而解析計(jì)算得到．因此，文獻(xiàn)［9］仍屬于疊前深度偏移范疇，面臨著與疊前深度偏移相同的困難．董春暉等［10］發(fā)展了一個(gè)基于雙速度參數(shù)的二維直接疊前時(shí)間偏移方法，有效地將疊前時(shí)間偏移技術(shù)推廣應(yīng)用于起伏地表問題．此外，實(shí)際起伏地表問題總是三維的，炮、檢點(diǎn)不在同一方位角平面，導(dǎo)致走時(shí)、等效偏移距（浮動基準(zhǔn)面上的偏移距）的計(jì)算較二維情況的更復(fù)雜．本文則在確定兩個(gè)等效速度參數(shù)、處理三維起伏地表這兩方面提升了文獻(xiàn)［10］的方法，并且將這一方法應(yīng)用于三維起伏地表采集的實(shí)際地震資料．

本文的關(guān)鍵有兩點(diǎn)，一是從疊前深度偏移的單程波方程出發(fā)，解析給出了三維情況下由描述近地表和上覆層速度的兩個(gè)等效速度表達(dá)的地震波走時(shí)、幅值和等效偏移距，并發(fā)展了基于表插值的快速算法；二是發(fā)展了基于偏移道集同相軸的雙曲彎曲和局部微彎曲，分別確定兩個(gè)等效速度參數(shù)的方法和算法流程．文中也應(yīng)用了隨方位角和構(gòu)造傾角變化的三維時(shí)變偏移孔徑，這使得本文疊前時(shí)間偏移方法即有效地壓制了偏移噪音，也減少了偏移計(jì)算量．三維起伏地表的理論數(shù)據(jù)和實(shí)際地震資料證明了本文方法和處理流程的正確和有效性．

2 走時(shí)、振幅與等效速度參數(shù)

假設(shè)三維非均勻介質(zhì)可近似為層狀介質(zhì)，不失一般性，令炮點(diǎn)或檢波點(diǎn)在坐標(biāo)原點(diǎn)；在波數(shù)－頻率域，基于三維深度偏移的相移法，該炮點(diǎn)或檢波點(diǎn)的波場深度延拓可用時(shí)間深度表示為

若將第m層介質(zhì)界面作為基準(zhǔn)面，可近似式（1）中右端指數(shù)項(xiàng)中的相移量為

將式（2）代入（1）式并做空間傅里葉反變換，可得到空間－頻率域波場：

式（3）是一個(gè)二重震蕩積分，可利用穩(wěn)相點(diǎn)原理［11］求得漸進(jìn)解為

式（4）中，定義

其中px和py是射線參數(shù)的零點(diǎn)，可由下式解得：

引入新的變量pr，定義p0x＝prcosφ，p0y＝prsinφ，其中φ是坐標(biāo)點(diǎn)（x，y）的方位角，則（7）和（8）式可轉(zhuǎn)化為單變量的方程：

而代入式（6），則可進(jìn)一步得到地震波的幅值．

令a2＝0，即假設(shè)基準(zhǔn)面以上是明顯的低速帶，則式（11）簡化為

這就是應(yīng)用靜校正技術(shù)的走時(shí)公式．這表明式（11）涵蓋了現(xiàn)行的靜校正方法．若假定介質(zhì)為上、下層速度分別是V0和Vrms的雙層介質(zhì)，式（9）和（11）則可由Snell定律直接導(dǎo)出，這從另一個(gè)角度證明了式（9）和（11）的正確性；式（9）和（11）既擴(kuò)展了Snell定律的結(jié)果，也給出了非均勻介質(zhì)情況下等效速度參數(shù)的物理含義．

盡管可解析求解［12］式（9），但對三維偏移過程中大量的成像點(diǎn)、炮點(diǎn)、檢波點(diǎn)組合而言，直接計(jì)算涉及了太多的計(jì)算量．既然3個(gè)無量綱系數(shù)完全決定了式（9）的四次方程，我們發(fā)展了基于表插值的快速算法．3個(gè)無量綱系數(shù)均有明確的物理意義，就疊前時(shí)間偏移方法而言，其對大于70°角度傳播的高角度波場是不感興趣的，因此可令0≤a3≤tan70；即使高速層出露，近地表等效速度也不可能大于深層速度，因此可定義0≤a2≤1．0；疊前偏移成像更多的是關(guān)心中、深層，若假定從T2≥T1開始成像計(jì)算，則可定義0≤a1≤1．0．這樣，若對這3個(gè)無量綱系數(shù)等間距離散采樣，預(yù)先求解對應(yīng)的η并建立一個(gè)三維表，就可在偏移過程中通過計(jì)算3個(gè)無量綱系數(shù)，直接在表中拾取對應(yīng)的η．炮點(diǎn)至成像點(diǎn)和檢波點(diǎn)至成像點(diǎn)可用同一個(gè)表．這一基于表的算法可極大提高三維偏移的計(jì)算效率．

盡管（9）、（11）以及（6）式的走時(shí)和振幅值是基于（1）式的層狀介質(zhì)假設(shè)推出的，但通過允許兩個(gè)等效速度V0和Vrms以及基準(zhǔn)面橫向變化，式（9）、（11）以及（6）式的解可代表三維橫向弱非均勻介質(zhì)中地震波的走時(shí)和振幅值；對較強(qiáng)的速度橫向變化情況，只能采用基于層速度模型的深度域延拓方法；文獻(xiàn)［13］通過引入速度梯度等信息，對較強(qiáng)橫向變速情況下的走時(shí)計(jì)算進(jìn)行了討論．

等效速度參數(shù)Vrms和基準(zhǔn)面深度的橫向變化是由成像點(diǎn)的橫向坐標(biāo)（x，y）決定的．這樣，Vrms將與基準(zhǔn)面上的等效偏移距構(gòu)成常規(guī)疊前時(shí)間偏移的雙曲走時(shí)關(guān)系；如果令a2＝0，則本文算法可退化為浮動基準(zhǔn)面疊前時(shí)間偏移方法．另一個(gè)等效參數(shù)V0即可隨炮、檢點(diǎn)橫向變化，也可隨成像點(diǎn)橫向變化．

3 等效速度參數(shù)估計(jì)與等效偏移距計(jì)算

第2節(jié)的研究表明，可用兩個(gè)等效速度參數(shù)描述地震波在近地表和上覆層中的傳播；而如何準(zhǔn)確獲得這兩個(gè)等效參數(shù)，是本文方法實(shí)際應(yīng)用的關(guān)鍵問題．本文的兩個(gè)等效速度參數(shù)在每個(gè)成像點(diǎn)是惟一的，這兩個(gè)參數(shù)完全決定了該成像點(diǎn)處共反射點(diǎn)道集（CRP）的表現(xiàn)，因此可在偏移過程中根據(jù)該CRP道集的表現(xiàn)來決定它們；不同于文獻(xiàn)［9］的雙速度參數(shù)，它們是成像點(diǎn)和炮點(diǎn)（或檢波點(diǎn)）的函數(shù)，若干個(gè)不同的參數(shù)決定了該成像點(diǎn)處CRP道集的表現(xiàn)，因此不能由單個(gè)CRP道集來決定該速度參數(shù)．實(shí)際上，文獻(xiàn)［9］中的雙速度參數(shù)只是用于近似計(jì)算地震波的走時(shí)，它們是由已知的層速度模型解析計(jì)算得到的．

從偏移成像的角度出發(fā)，CRP道集中同相軸是否平直是確定兩個(gè)等效速度的標(biāo)準(zhǔn)，這也是本文確定等效速度的出發(fā)點(diǎn)．但兩個(gè)等效速度參數(shù)共同影響了同相軸的平直，只有區(qū)分這兩者對同相軸彎曲的獨(dú)特影響，才能有效地確定這兩個(gè)參數(shù)．

圖1 等效速度參數(shù)對走時(shí)影響分析示意圖Fig．1 Illustration of equivalent velocity parameter analysis using 2Dmodel including topography

圖2 等效速度參數(shù)對走時(shí)殘差影響Fig．2 Influence of equivalent velocity parameters on residual traveltime

本文用一個(gè)簡單的二維模型來考察兩個(gè)速度參數(shù)變化對走時(shí)的影響．如圖1所示，起伏地表按正弦函數(shù)變化，其波長為1000m，波鋒與波谷間的高程為200m；成像點(diǎn)在模型中心，對應(yīng)該成像點(diǎn)的基準(zhǔn)面如圖1中虛線所示，較地表高程曲線的波谷低0．1s（時(shí)間深度），而成像點(diǎn)在基準(zhǔn)面下1．5s處；基準(zhǔn)面上、下的介質(zhì)是均勻的，速度分別為1800m／s和2500m／s．圖2a給出了由（11）式求得的地面上炮點(diǎn)至成像點(diǎn)再反射到地表檢波點(diǎn)的地震波走時(shí)隨偏移距（水平距離）的變化曲線．若令V0＝1800m／s，Vrms＝2000，2500，3000m／s，利用式（11）做動校正，可得剩余走時(shí)曲線（如圖2b），而圖2c則給出了隨基準(zhǔn)面上的等效偏移距變化的剩余走時(shí)曲線；圖2b中的點(diǎn)線是最佳擬合剩余走時(shí)的雙曲線，而圖2c的點(diǎn)線與剩余走時(shí)曲線的實(shí)線重合．觀察圖2b和圖2c可發(fā)現(xiàn)，速度參數(shù)Vrms影響了同相軸的整體彎曲程度，而等效偏移距使得這種彎曲更好地展現(xiàn)了如常規(guī)時(shí)間偏移的雙曲特征．若令V0＝1600，1800，2000m／s，Vrms＝2500m／s，利用式（11）做動校正，可得剩余走時(shí)曲線（如圖2d）．觀察圖2d可發(fā)現(xiàn)，不準(zhǔn)確的速度參數(shù)V0在產(chǎn)生同相軸的小幅度彎曲的同時(shí)，導(dǎo)致了同相軸的局部繞曲．因此，CRP道集中同相軸的局部微彎曲是確定V0的指標(biāo)，而整體的雙曲彎曲可指示速度參數(shù)Vrms的正確與否．實(shí)際應(yīng)用中，同相軸的局部微彎曲可通過疊加道波形的寬窄來判斷，而雙曲彎曲使得可用常規(guī)的NMO方法拾取剩余速度．

上述討論表明，基準(zhǔn)面上的等效偏移距對確定速度參數(shù)是重要的；而三維起伏地表情況，特別是炮、檢點(diǎn)不在同一方位角平面上時(shí)，等效偏移距的計(jì)算較二維情況變的更復(fù)雜．圖3給出了三維情況下等效偏移距的示意圖，圖中S和G點(diǎn)分別代表起伏地表上的炮點(diǎn)和檢波點(diǎn)，AOB平面代表基準(zhǔn)面，I代表成像點(diǎn)．采用常規(guī)靜校正方法時(shí)，地震道的偏移距不變，如圖中的PQ；而本文方法考慮了地震波的實(shí)際傳播路徑，基準(zhǔn)面上的真實(shí)（等效）偏移距應(yīng)是MN．

已由式（9）或（10）求得由S傳至M再傳至I點(diǎn)的地震波的射線參數(shù)pr．對M傳至I點(diǎn)的地震波應(yīng)用式（10）（可認(rèn)為T1為零），地震波的射線參數(shù)pr不變，則簡單得到

式中（xs，ys）是炮點(diǎn)水平坐標(biāo)，（xg，yg）是檢波點(diǎn)水平坐標(biāo)，（x0，y0）是成像點(diǎn)水平坐標(biāo)，ηs和ηg分別是根據(jù)炮點(diǎn)和檢波點(diǎn)信息由表中拾取的η．式（14）表明利用基于表插值的快速算法可在得到走時(shí)和幅值的同時(shí)，得到等效偏移距；圖3中曲線SMI和GNI代表非均勻介質(zhì)中地震波的傳播路徑，因此式（14）可適用于非均勻介質(zhì)．

圖3 等效偏移距示意圖Fig．3 Illustration of equivalent offset on the datum

實(shí)際應(yīng)用中等效參數(shù)V0的估計(jì)有兩種情況：一是地表調(diào)查、微測井、折射或?qū)游鲮o校正等方法已給出了較準(zhǔn)確的近地表速度場；二是尚沒有做上述工作．對于第一種情況，令等效參數(shù)V0橫向通過較小的百分比掃描，根據(jù)各水平位置處剩余動校后疊加道波形的寬窄和道集的情況綜合確定各水平位置處的最佳百分比參數(shù)，對百分比參數(shù)橫向平滑，用平滑后的百分比參數(shù)修正近地表速度，即得到等效參數(shù)V0．對第二種情況，可以估計(jì)幾個(gè)近地表速度對最小偏移距數(shù)據(jù)進(jìn)行偏移，觀察最小偏移距的成像情況確定一個(gè)初始的等效參數(shù)V0，再將等效參數(shù)V0做速度百分比掃描，根據(jù)各水平位置處剩余動校后疊加道波形的寬窄和道集情況確定各水平位置處新的V0，橫向平滑后即是最終的近地表速度場．評估疊加道波形的寬窄時(shí)應(yīng)綜合考慮幾個(gè)不同深度的同相軸．

等效參數(shù)Vrms可根據(jù)求得的等效偏移距，在偏移形成的CRP道集上，利用常規(guī)的NMO技術(shù)拾取．詳見第5節(jié)的數(shù)據(jù)處理流程．

4 三維偏移與成像幅值權(quán)系數(shù)

基于式（10）的查表法僅給出炮點(diǎn)至成像點(diǎn)和檢波點(diǎn)至成像點(diǎn)的走時(shí)和振幅值；而偏移的目的是進(jìn)一步得到地下界面的反射系數(shù)圖像．若將單個(gè)地震道看作是僅有一個(gè)接收道的單炮記錄，則根據(jù)（4）式和基于（10）式的表插值算法，可分別得到炮域偏移的下行波場和反傳波場：

式中假設(shè)震源是一時(shí)間脈沖，f（ω）是接收信號的傅里葉變換．將式（15）代入波動方程疊前深度偏移的反褶積成像條件［14］，有成像結(jié)果

式中F′（t）是f（ω）對應(yīng)的時(shí)域函數(shù)的一階導(dǎo)數(shù)．式（16）的成像結(jié)果表明，對任一單道數(shù)據(jù)，對數(shù)據(jù)求一階導(dǎo)數(shù)；對成像區(qū)域的每一成像點(diǎn)，計(jì)算走時(shí)ts＋tg和成像權(quán)系數(shù)Ag／As；在數(shù)據(jù)的一階導(dǎo)數(shù)上拾取ts＋tg時(shí)刻的振幅值并乘上權(quán)系數(shù)，即完成該道數(shù)據(jù)的成像．

對全部地震道做上述操作，將成像結(jié)果累加，即完成了全部數(shù)據(jù)的成像．不同于忽略權(quán)系數(shù)的現(xiàn)行方法，式（16）的權(quán)系數(shù)Ag／As實(shí)現(xiàn)了正確補(bǔ)償?shù)卣鸩ǖ膸缀螖U(kuò)散效應(yīng)，使得深層成像更清晰．

5 偏移流程

本文偏移流程的關(guān)鍵是在偏移過程中確定兩個(gè)等效速度V0和Vrms，偏移流程如下：

1．綜合地表起伏情況、可能的地表信息和最小偏移距剖面，確定一個(gè)光滑的定浮動基準(zhǔn)面；

2．對地表調(diào)查、微測井、折射或?qū)游鲮o校正等方法已給出了較準(zhǔn)確的近地表速度場情況，使用該速度場作為初始V0；對無此類信息情況，根據(jù)經(jīng)驗(yàn)確定一均勻速度作為初始；

3．根據(jù)V0，在幾個(gè)典型CDP位置，將CMP道集靜校正到浮動基準(zhǔn)面，利用NMO確定初始；

4．對未經(jīng)靜校正處理的數(shù)據(jù)，依據(jù)初始V0和Vrms做式（16）的偏移計(jì)算，計(jì)算等效偏移距，形成等效偏移距下的共反射點(diǎn)道集；

5．利用初始Vrms，對共反射點(diǎn)道集做反動校，再做NMO動校拾取速度，對拾取的速度做空間平滑，作為新的Vrms；

6．利用初始V0和新的Vrms，對V0做百分比掃描，對較疏的CDP位置和主要的同相軸區(qū)域做局部偏移計(jì)算，利用第3節(jié)討論的方法確定新的V0；

7．在新的V0和Vrms下重新進(jìn)行偏移計(jì)算，形成等效偏移距下的共反射點(diǎn)道集．若剩余動校正量較小，進(jìn)行剩余動校正、切除和疊加，得到最終偏移剖面；若剩余動校正量較大，利用Vrms做反動校，再做動校，形成再次更新的Vrms，對V0和Vrms重復(fù)上述過程，得到最終偏移剖面；

8．由于浮動基準(zhǔn)面一般是非水平的，基于浮動基準(zhǔn)面的偏移成像將扭曲地下構(gòu)造圖像．計(jì)算浮動基準(zhǔn)面與選定的水平基準(zhǔn)面的垂直走時(shí)，用這一時(shí)間修正偏移圖像，可得到基于水平基準(zhǔn)面的偏移剖面．

6 數(shù)值算例

為了驗(yàn)證本文方法的正確性和工業(yè)應(yīng)用的可能性，分別應(yīng)用本文方法處理了三維起伏地表采集的理論模型數(shù)據(jù)和我國東部某工區(qū)的實(shí)際三維起伏地表采集資料，前者證明了本文理論和方法的正確性，后者表明了本文方法具有實(shí)際工業(yè)應(yīng)用價(jià)值．

6．1 三維起伏地表模型數(shù)據(jù)

三維起伏地表模型數(shù)據(jù)由東方地球物理公司提供．該數(shù)據(jù)近地表高程變化明顯，地表不存在明顯的降低速帶，圖4是模擬該數(shù)據(jù)集用到的速度模型在inline和crossline方向的切片，實(shí)線表示地表高程，最大高差達(dá)到590m，虛線是本文采用的浮動基準(zhǔn)面．該理論數(shù)據(jù)集共3940炮，每炮10條線接收，每條線201道，每炮共2010道，中間放炮，兩邊接收．inline方向道間距40m，crossline方向線距100m，最小偏移距50m，最大偏移距5100m，記錄總長度為5s，采樣率為4ms．圖5是該三維數(shù)據(jù)的典型單炮記錄，從中可以看到：炮記錄中的波場已被嚴(yán)重扭曲，同相軸不再是水平地表情況下的雙曲線．圖6給出了使用本文方法得到的最終偏移結(jié)果．5個(gè)界面的反射同相軸被很好的成像，成像結(jié)果與速度模型在形狀上不完全對應(yīng)是因?yàn)槌上窠Y(jié)果是時(shí)間剖面，這也是時(shí)間偏移的特點(diǎn)．由于利用了成像權(quán)系數(shù)，不用AGC深層構(gòu)造也得到了較好的成像；第5個(gè)同相軸較弱是該層的速度對比很小（見圖4速度模型），陡界面成像不好是因?yàn)榈卣鹩涗浀钠凭噍^小，導(dǎo)致不能接收到陡界面的反射信號所致．

圖4 三維理論數(shù)據(jù)速度模型，地表高程（實(shí)線），浮動基準(zhǔn)面（虛線）Fig．4 Velocity model，surface elevation（solid）and floating datum（dot）of 3Dsynthetic data

圖5 三維理論數(shù)據(jù)典型單炮記錄Fig．5 Typical shot record of 3Dsynthetic data

圖6 三維理論數(shù)據(jù)直接疊前時(shí)間偏移剖面Fig．6 Direct prestack time migration result of the 3Dsynthetic data

該理論模型盡管簡單，但由于地表不存在明顯的降低速帶，應(yīng)用常規(guī)的靜校正加浮動基準(zhǔn)面方法是不能正確成像的．

6．2 三維起伏地表實(shí)際數(shù)據(jù)

本節(jié)進(jìn)一步給出了應(yīng)用本文方法處理某工區(qū)三維實(shí)際數(shù)據(jù)的例子．圖7是該工區(qū)地表高程的等值線圖，可見地表起伏變化是較劇烈的．三維資料共1710炮，每炮16線接收，每線142道，道距40m，線距120m，中間放炮，記錄長度為6s，采樣率為4ms．該三維資料經(jīng)過壓制面波等去噪處理，但沒有進(jìn)行靜校正處理．圖8顯示了典型的單炮記錄，地表起伏導(dǎo)致的同相軸錯(cuò)動明顯可見．圖9給出了在兩個(gè)等效速度參數(shù)進(jìn)行更新過程中，典型共反射點(diǎn)道集中同相軸的變化情況，隨著速度參數(shù)的更新，同相軸變得更加平直和清晰．

圖7 三維實(shí)際數(shù)據(jù)地表高程Fig．7 Surface elevation of the 3Dfield data

我們沒有得到該工區(qū)地表速度的任何資料，因此憑經(jīng)驗(yàn)選取1800m／s作為初始．圖10是更新前后在一條inline線上的速度對比情況，圖11是最終得到的在一條inline線上的分布圖．圖12是在一條inline線上的偏移疊加剖面，同相軸連續(xù)性較好，斷層斷點(diǎn)清晰，原始資料中的同相軸錯(cuò)動被消除，這表明本文方法在偏移過程中正確地處理了起伏地表影響，實(shí)現(xiàn)了直接偏移．

本文偏移處理是在沒有任何近地表速度資料的情況下進(jìn)行的，這表明本文方法具有了自主建立近地表和上覆層速度的能力．

7 結(jié) 論

圖8 典型三維單炮記錄中單條測線的接收信號Fig．8 Seismic signature of one receiver line with respect to 3Dfield data

圖9 CRP道集隨著速度更新的變化情況Fig．9 Changes in CRP gather while updating the velocities

圖10 Inline方向速度參數(shù)V0更新情況Fig．10 Comparison between initial and updated V0parameter in inline

圖11 Inline方向速度參數(shù)Vrms更新結(jié)果Fig．11 Updated parameter Vrmsin inline

本文提出了一種可對起伏地表采集的三維地震資料直接進(jìn)行偏移成像的疊前時(shí)間偏移方法和流程．該方法可適應(yīng)高速層出露等不存在明顯低、降速帶的情況，能在偏移過程中自主建立（或修正）近地表速度模型和偏移使用的疊加速度模型．三維理論模型數(shù)據(jù)證明了方法的正確性和其處理近地表不存在明顯低、降速帶情況的能力．未經(jīng)靜校正處理的三維實(shí)際陸上資料的直接偏移結(jié)果表明，本文方法具有很好的工業(yè)化應(yīng)用能力．即使不考慮高速層出露，僅從避免了初至拾取、近地表速度建模等繁雜的靜校正流程而言，本文方法就有很好的實(shí)用價(jià)值．本文方法為我國陸上油氣勘探面臨的復(fù)雜地表問題指示了一個(gè)具有工業(yè)應(yīng)用價(jià)值的有效途徑．

圖12 實(shí)際數(shù)據(jù)直接疊前時(shí)間偏移疊加剖面的inline切片F(xiàn)ig．12 Stack section of direct prestack time migration in inline of 3Dfield data

（References）

［1］ Sun C W，Jing P G，Pu Y．Seismic imaging in the mountainous area of Southern China：Statics correction compared with migration from topography．SEG expanded abstract，2009．

［2］ 劉國峰，劉洪，李博等．山地地震資料疊前時(shí)間偏移方法及其GPU實(shí)現(xiàn)．地球物理學(xué)報(bào)，2009，52（12）：3101－3108．Liu G F，Liu H，Li B，et al．Method of prestack time migration of seismic data of mountainous regions and its GPU implementation．Chinese J Geophys．（in Chinese），2009，52（12）：3101－3108．

［3］ Cox M．Static Corrections for Seismic Reflection Surveys．Tulsa：Society of Exploration Geophysics，1999．

［4］ Bevc D．Flooding the topography：Wave－equation datuming of land data with rugged acquisition topography．Geophysics，1997，62（5）：1558－1569．

［5］ Kelamis P G，Erickson K E，Verschuur D J，Berkhout A J．Velocity－independent redatuming：A new approach to the near－surface problem in land seismic data processing．The Leading Edge，2002，21（8）：730－735．

［6］ Black J L，Brzostowski M A．Systematics of time－migration errors．Geophysics，1994，59（9）：1419－1434．

［7］ Reshef M．Depth migration from irregular surfaces with depth extrapolation methods．Geophysics，1991，56（1）：119－122．

［8］ Rajasekaran S，McMechan G A．Prestack processing of land data with complex topography．Geophysics，1995，60（6）：1875－1886．

［9］ Wiggins J W．Kirchhoff integral extrapolation and migration of nonplanar data．Geophysics，1984，49（8）：1239－1248．

［10］ 董春暉，張劍鋒．起伏地表下的直接疊前時(shí)間偏移．地球物理學(xué)報(bào)，2009，52（1）：239－244．Dong C H，Zhang J F．Prestack time migration including surface topography．Chinese J．Geophys．（in Chinese），2009，52（1）：239－244．

［11］ 李世雄．波動方程的高頻近似與辛幾何．北京：石油工業(yè)出版社，2001．Li S X．Wave Equation with High Frequency Approximation and Symplectic Geometry（in Chinese）．Beijing：Petroleum Industry Press，2001．

［12］ Zhang J F，Verschuur D J，Wapenaar C P A．Depth migration of shot records in heterogeneous，transversely isotropic media using optimum explicit operators．Geophysical Prospecting，2001，49（3）：287－299．

［13］ 張廉萍，劉洪，李幼銘．單程波李代數(shù)深度積分的精度分析和算法改進(jìn)．地球物理學(xué)報(bào)，2010，53（11）：2739－2746．Zhang L P，Liu H，Li Y M．The precision analysis and algorithm improvement in Lie algebra depth integral of oneway seismic wave．Chinese J．Geophys．（in Chinese），2010，53（11）：2739－2746．

［14］ 張宇．振幅保真的單程波方程偏移理論．地球物理學(xué)報(bào)，2006，49（5）：1410－1430．Zhang Y．The theory of true amplitude oneway wave equation migration．Chinese J．Geophys．（in Chinese），2006，49（5）：1410－1430．

3Dprestack time migration including surface topography

ZHANG Hao1，2，ZHANG Jian－Feng1

1 Key Laboratory of Petroleun Resources Research，Institute of Geology and Geophysics，Chinese Academy of Sciences，Beijing100029，China 2 Graduate University，Chinese Academy of Sciences，Beijing100049，China

We present a 3Dprestack time migration scheme and workflow that can image land seismic data without datum static corrections．The scheme mimics the influences of the nearsurface and overlay on wave propagation using two equivalent velocities．As a result，it is not necessary for wave path to be vertical incidence in the near－surface，as in conventional static corrections．This leads to the proposed scheme which can handle the more complex surface cases，such as the exposure of high velocity in the near－surface．The proposed scheme can update the two equivalent velocities with respect to the flats and consistency of the events．A synthetic 3D dataset and a field 3Ddataset without datum static corrections are used to demonstrate the scheme and the workflow．

Topography，3Dimaging，Near－surface velocity，Pre－stack time migration

P631收修定稿2010－12－16，2012－01－17收修定稿

國家自然科學(xué)基金重點(diǎn)項(xiàng)目（40930422）和國家科技重大專項(xiàng)（2011ZX05008－006）資助．

張浩，男，1983年生，中國科學(xué)院地質(zhì)與地球物理研究所博士研究生，主要從事疊前偏移成像研究．E－mail：zhanghao＠m(xù)ail．iggcas．a(chǎn)c．cn

張浩，張劍鋒．起伏地表采集數(shù)據(jù)的三維直接疊前時(shí)間偏移方法．地球物理學(xué)報(bào)，2012，55（4）：1335－1344，

10．6038／j．issn．0001－5733．2012．04．029．

Zhang H，Zhang J F．3Dprestack time migration including surface topography．Chinese J．Geophys．（in Chinese），2012，55（4）：1335－1344，doi：10．6038／j．issn．0001－5733．2012．04．029．

10．6038／j．issn．0001－5733．2012．04．029

（本文編輯 汪海英）