徐錦承， 張劍鋒

1 中國科學(xué)院地質(zhì)與地球物理研究所, 北京　100029 2 中國科學(xué)院大學(xué), 北京　100049

?

基于偏移成像道集的剩余靜校正方法

徐錦承1,2， 張劍鋒1

1 中國科學(xué)院地質(zhì)與地球物理研究所, 北京100029 2 中國科學(xué)院大學(xué), 北京100049

摘要針對(duì)陸上地震資料處理的靜校正問題，提出了一種基于偏移成像道集的剩余靜校正方法.與傳統(tǒng)的由動(dòng)校正后的CMP道集中拾取剩余時(shí)差不同，本文基于偏移成像道集求取剩余時(shí)差，避免了復(fù)雜情況下同相軸歸位不準(zhǔn)確導(dǎo)致的剩余時(shí)差拾取誤差.通過生成隨炮點(diǎn)和檢波點(diǎn)位置變化的偏移道集，實(shí)現(xiàn)了由偏移道集中直接拾取炮、檢點(diǎn)的地表一致性剩余時(shí)差；該炮、檢點(diǎn)偏移道集只在指定的局部時(shí)窗生成，并不增加大的計(jì)算量.二維和三維實(shí)際數(shù)據(jù)測試表明了該方法的有效性和實(shí)用性.

關(guān)鍵詞剩余靜校正；疊前偏移；偏移成像道集

The method of residual static corrections is used to improve the resolution of pre-stack migration. We obtain the residual time shifts by migrated gathers rather than CMP gathers in the conventional method, which can reduce the artificial picking error by incorrect migration or indistinct events. Firstly, we generate the shot-domain and receiver-domain migrated gathers in the local time window. The de-noised migrated result is used as the reference gathers section. Then, the surface consistent residual time shifts are obtained by the shot-domain and receiver-domain migrated gathers. We need to calculate the residual statics according to shot and receiver positions separately. Same as conventional residual static correction, we pick up the residual time shifts in the migrated gathers by correlation with the reference gathers. After this, we remove the outliers and compute the arithmetic mean of the residual time shifts in each shot and receiver. This whole process can be carried out only in the selected time windows, which can save the calculation obviously. Finally, we put those residual statics of every shot and receiver into data in the migration.

We test this method on two field data set examples. One is 2.5D real seismic data recorded with a single receiver line. Another is the real data of 3D acquisition. The proposed method can be used both in 2D and 3D field data sets, which can be easily combined to conventional seismic data processing. The CRP gathers comparison show that this approach improves the quality greatly. The events become more and more continuous and clear. The comparisons of migrated stack section clearly show that by imposing the proposed residual static corrections, we can image clearer faults and more continuous events.

After theoretical investigation and real data application, we have the following conclusions: (1) The proposed method is suitable and effective for 2D and 3D real land seismic data. (2) This method can be easily combined with the existing seismic data processing flow. (3) The proposed process flow does not increase amount of computation. (4) The results show the residual static correction based on migrated gathers can obviously improve the S/N and resolution of pre-stack time migration.

1引言

靜校正是陸上地震資料處理的一個(gè)關(guān)鍵問題.靜校正通常包含野外靜校正和剩余靜校正.野外靜校正主要解決的是長波長時(shí)差(低頻部分)，而剩余靜校正主要解決短波長時(shí)差(高頻部分)問題.

傳統(tǒng)的野外靜校正假定炮點(diǎn)、檢波點(diǎn)的出射和入射波是垂直傳播的，這一假設(shè)對(duì)近地表存在明顯的低速層是成立的.由于應(yīng)用簡單，已在實(shí)際勘探中得到廣泛的應(yīng)用.對(duì)于高速層出露、不存在明顯低速帶的一類地質(zhì)目標(biāo)，就不能采用這類方法，必需考慮地震波在近地表地層中的實(shí)際傳播路徑；為此發(fā)展了基于波場延拓的基準(zhǔn)面靜校正方法和考慮實(shí)際采集表面的真地表(直接)疊前偏移方法(Reshef, 1991; Rajasekaran and McMechan, 1995; 董春暉和張劍鋒, 2009).但這兩類方法都不能消除由于近地表速度快速橫向變化、測量和采集誤差導(dǎo)致的短波長扭曲，必需進(jìn)一步應(yīng)用剩余靜校正.

剩余靜校正的引入既可以消除測量和觀測誤差，也可以補(bǔ)償近地表速度模型和靜校正方法本身的誤差，對(duì)提高成像分辨率有明顯應(yīng)用效果.目前大多數(shù)剩余靜校正方法都是通過反射波信息采用統(tǒng)計(jì)相關(guān)技術(shù)來提高疊加質(zhì)量(Marsden, 1993).其求解方法主要有線性反演和非線性反演兩類.其中線性反演法是對(duì)旅行時(shí)分解計(jì)算炮點(diǎn)和檢波點(diǎn)的時(shí)移量，主要代表有廣義線性反演法(Wiggins et al., 1976)和最大能量法(Ronen and Claerbout, 1985)；非線性反演方法主要代表有遺傳算法和模擬退火法(Rothman, 1985; Stork and Kusuma, 1992)以及將二者混合的全局優(yōu)化方法(井西利等, 2002；茍量等, 2005).

現(xiàn)行的剩余靜校正方法都是在基準(zhǔn)面靜校正(野外靜校正)和動(dòng)校正后的CMP道集上拾取剩余時(shí)差.復(fù)雜反射情況下，并不是所有反射同相軸都在CMP道集上呈現(xiàn)雙曲走時(shí)，因此動(dòng)校正并不能拉平所有同相軸；而且低信噪比情況下，也很難準(zhǔn)確判定弱幅值的同相軸.這些因素導(dǎo)致直接在(動(dòng)校正后的)CMP道集上拾取剩余時(shí)差易于引入誤差，也存在著使同相軸產(chǎn)生錯(cuò)位的風(fēng)險(xiǎn).由于偏移考慮了反射界面傾斜，斷點(diǎn)繞射等復(fù)雜情況，可使反射波正確歸位；所有反射同相軸在偏移后道集上均應(yīng)呈現(xiàn)平直狀態(tài)，而部分疊加也提高了偏移道集的信噪比，因此基于偏移后的道集拾取剩余時(shí)差是一個(gè)更好的選擇.

疊前時(shí)間偏移目前已成為工業(yè)界廣泛應(yīng)用的關(guān)鍵技術(shù)，該方法可對(duì)一類傾角、斷層較為復(fù)雜但速度橫向變化不是很劇烈的構(gòu)造較好成像.與疊前深度偏移相比，除具有較高的計(jì)算效率外，其主要的優(yōu)點(diǎn)是只需使用疊加速度；這樣可簡單地通過速度掃描等方式得到恰當(dāng)?shù)乃俣饶Ｐ?，回避了使用疊前深度偏移方法面臨的一個(gè)主要困難.因此本文以波動(dòng)方程架構(gòu)的積分疊前時(shí)間偏移為基礎(chǔ)(趙偉等, 2010)，發(fā)展基于偏移后道集拾取剩余時(shí)差的剩余靜校正方法；但這一思路可同樣應(yīng)用于基于KKF積分的疊前深度偏移方法.若進(jìn)一步應(yīng)用反偏移方法，基于偏移道集的剩余靜校正方法也可用于生成靜校正的疊前資料.

常規(guī)疊前時(shí)間偏移生成的偏移道集是按偏移距大小構(gòu)成的，不同炮點(diǎn)和檢波點(diǎn)的剩余時(shí)差在這一道集中已經(jīng)被混合，難以區(qū)分；為基于疊前時(shí)間偏移的偏移道集拾取剩余時(shí)差，必須形成新的、按炮點(diǎn)和檢波點(diǎn)組成的偏移道集.現(xiàn)行的剩余靜校正方法一般假定地表一致性，即要求同一位置的炮點(diǎn)或檢波點(diǎn)有相同的剩余時(shí)差，這一假定對(duì)保證剩余靜校正的穩(wěn)定性，避免引起同相軸的扭曲和錯(cuò)位有重要作用.為在基于偏移道集的剩余靜校正方法中實(shí)現(xiàn)地表一致性要求，本文提出了分別基于形成炮點(diǎn)和檢波點(diǎn)偏移道集拾取剩余時(shí)差的剩余靜校正方法.

文中以實(shí)際的起伏地表采集資料為例討論所發(fā)展的基于偏移道集的剩余靜校正方法和使用流程.實(shí)際資料的成像效果證明了本文方法的有效性.該方法易于結(jié)合到現(xiàn)行的地震數(shù)據(jù)處理流程；由于炮、檢點(diǎn)偏移道集只在指定的局部時(shí)窗生成，因此并不增加大的計(jì)算量.

2方法原理

2.1炮、檢點(diǎn)偏移道集

T≥T0(x,y,xm,ym,h)，

(1)

在式(1)的共反射點(diǎn)道集中，具有不同炮點(diǎn)和檢波點(diǎn)的偏移道已經(jīng)疊加(只要它們有相同的偏移距)；因此不同炮點(diǎn)和檢波點(diǎn)的剩余時(shí)差在這一道集中已經(jīng)被混合，難以區(qū)分.為拾取單個(gè)炮點(diǎn)和檢波點(diǎn)的剩余時(shí)差，必須形成新的按炮點(diǎn)和檢波點(diǎn)組成的偏移道集，表達(dá)式為

T≥T0(x,y,xm,ym,h)，

(2)

T≥T0(x,y,xm,ym,h)，

(3)

其中(xs,ys)和(xg,yg)分別是地震道的炮點(diǎn)和檢波點(diǎn)坐標(biāo).在式(2)和(3)中，具有相同炮點(diǎn)或相同檢波點(diǎn)的偏移道集將被疊加.由于式(2)和(3)的炮點(diǎn)和檢波點(diǎn)偏移道集不能保留準(zhǔn)確的剩余動(dòng)校正，也不能區(qū)分長偏移距的拉伸，在形成炮點(diǎn)和檢波點(diǎn)偏移道集時(shí)必須同時(shí)完成剩余動(dòng)校正和拉伸切除.理論上，在成像剖面的每個(gè)水平位置(x,y)，都有關(guān)于全部炮點(diǎn)或檢波點(diǎn)的偏移道集；實(shí)際上，由于偏移孔徑以及拉伸切除的影響，在每個(gè)水平位置僅有關(guān)于部分炮點(diǎn)或檢波點(diǎn)的有效偏移道集.當(dāng)計(jì)算用于拾取剩余時(shí)差的炮、檢點(diǎn)偏移道集時(shí)，可選取更小的偏移孔徑.

與常規(guī)剩余靜校正方法一樣，基于偏移道集的剩余靜校正也是在局部時(shí)窗內(nèi)進(jìn)行的，因此式(2)和(3)的偏移僅需對(duì)給定的時(shí)窗成像.一般選擇具有清晰的、連續(xù)或部分連續(xù)的同相軸的區(qū)域?yàn)橹付〞r(shí)窗.時(shí)窗的時(shí)間深度可隨CDP位置變化.根據(jù)資料品質(zhì)的不同，在同一CDP位置可選取幾個(gè)不同時(shí)間深度的時(shí)窗，這樣可將不同窗口中求出的剩余靜校正量平均，以取得更穩(wěn)定的結(jié)果.

2.2地表一致性的炮、檢點(diǎn)剩余時(shí)差

地表一致性假定要求同一位置的炮點(diǎn)或檢波點(diǎn)有相同的剩余時(shí)差，這一假定對(duì)壓制噪聲影響，保證剩余靜校正的穩(wěn)定性，避免引起同相軸的扭曲和錯(cuò)位有重要作用.現(xiàn)行的剩余靜校正方法通過聯(lián)立求解超定的時(shí)差分解方程(Taner et al., 1974; Larner et al., 1979)，獲得滿足地表一致性假設(shè)的炮點(diǎn)和檢波點(diǎn)剩余時(shí)差.不同于動(dòng)校后的CMP道集，炮點(diǎn)偏移道集中每道的剩余時(shí)差分解為

(4)

其中下標(biāo)j代表炮點(diǎn)，i代表檢波點(diǎn)，k代表成像的水平位置(CDP點(diǎn))，Sj和Ri分別是對(duì)應(yīng)炮點(diǎn)和檢波點(diǎn)位置的、滿足地表一致性假設(shè)的炮點(diǎn)和檢波點(diǎn)剩余時(shí)差，Nk為噪聲，Ωjk為對(duì)應(yīng)炮和成像點(diǎn)的炮偏移的菲涅耳帶，Nf為Ωjk內(nèi)包含檢波點(diǎn)的個(gè)數(shù).準(zhǔn)確確定式(4)中與菲涅耳帶有關(guān)的剩余時(shí)差貢獻(xiàn)是困難的，因此不能沿用常規(guī)剩余靜校正的方法.

對(duì)同一炮點(diǎn)在不同CDP點(diǎn)的偏移道集中的剩余時(shí)差做算術(shù)平均，公式為

(5)

其中M代表成像位置的個(gè)數(shù)，L代表檢波點(diǎn)的個(gè)數(shù).假設(shè)噪聲和剩余時(shí)差滿足高斯分布，式(5)中右端的第2和第3項(xiàng)將趨于零；由此可得滿足地表一致性假設(shè)的炮點(diǎn)剩余時(shí)差為

(6)

同理可得滿足地表一致性假設(shè)的檢波點(diǎn)剩余時(shí)差為

(7)

其中tik是對(duì)應(yīng)第k個(gè)CDP點(diǎn)(成像水平位置)的檢波點(diǎn)偏移道集中對(duì)應(yīng)第i個(gè)檢波點(diǎn)的剩余時(shí)差.

實(shí)際應(yīng)用中，由于偏移孔徑可切除拉伸作用，每個(gè)炮點(diǎn)或檢波點(diǎn)對(duì)應(yīng)的有效M是有限的；此外，在做算術(shù)平均前，應(yīng)采用中值濾波方法剔出異常的剩余時(shí)差.

2.3剩余時(shí)差拾取

與常規(guī)剩余靜校正方法相同，在基于炮、檢點(diǎn)偏移道集的剩余時(shí)差拾取中(即決定tjk和tik時(shí))，將利用偏移道與參考道在指定時(shí)窗中的互相關(guān)來估計(jì)剩余時(shí)差.本文取偏移疊加剖面中的成像道作為參考道集.為提高參考道的信噪比，將對(duì)偏移疊加剖面應(yīng)用FX和FK濾波等方法衰減隨機(jī)和相干噪聲；這樣得到的參考道較常規(guī)剩余靜校正方法使用的參考道有更高的信噪比.

對(duì)每個(gè)水平位置k，從該點(diǎn)的炮點(diǎn)偏移道集中選取對(duì)應(yīng)的偏移道gj(T)，從噪聲衰減后的偏移疊加剖面中選取對(duì)應(yīng)的成像道pk(T)，可求得互相關(guān)結(jié)果為

(9)

用檢波點(diǎn)偏移道集的偏移道qi(T)和成像道pk(T)，同理可決定檢波點(diǎn)的剩余時(shí)差tik.這樣就可由式(6)和(7)求得地表一致性假設(shè)的炮點(diǎn)和檢波點(diǎn)剩余時(shí)差.

在剩余時(shí)差拾取過程中，有兩個(gè)因素決定著剩余時(shí)差的拾取精度，即時(shí)窗的選取和互相關(guān)函數(shù)峰值點(diǎn)的判斷.在確定時(shí)窗時(shí)，應(yīng)選取反射波組能量較強(qiáng)，特征明顯，信噪比高的層位；而時(shí)窗的寬度則和數(shù)據(jù)的信噪比有關(guān)，信噪比較低的情況下應(yīng)選擇較窄的時(shí)窗，以提高相關(guān)計(jì)算的準(zhǔn)確度，但也要注意過窄的時(shí)窗可能會(huì)造成假象.

2.4工作流程

本文剩余靜校正量的求取在疊前時(shí)間偏移的基礎(chǔ)上進(jìn)行，計(jì)算流程如圖1 所示，主要有以下幾個(gè)步驟：

(1)產(chǎn)生參考道.對(duì)疊前地震數(shù)據(jù)應(yīng)用疊前時(shí)間偏移，對(duì)生成的偏移疊加剖面應(yīng)用FX和FK濾波等方法衰減隨機(jī)和相干噪聲，偏移疊加剖面上的成像可作為參考道.

(2)確定計(jì)算時(shí)窗.在偏移疊加剖面上選定清晰的同相軸，據(jù)此在各個(gè)CDP處確定互相關(guān)計(jì)算的時(shí)窗.

(3)形成局部炮、檢點(diǎn)偏移道集.根據(jù)式(2)和(3)進(jìn)行指定時(shí)間窗口的局部偏移，形成各CDP處隨炮點(diǎn)坐標(biāo)變化的局部炮點(diǎn)偏移道集和隨檢波點(diǎn)坐標(biāo)變化的局部檢波點(diǎn)偏移道集.

圖1　剩余靜校正流程圖Fig.1　Flow chart of residual static correction

(4)計(jì)算剩余時(shí)差.將參考道與對(duì)應(yīng)的炮、檢點(diǎn)偏移道集利用式(8)和(9)計(jì)算得到剩余時(shí)差.

(5)求取剩余靜校正量.由式(6)和(7)將各組剩余時(shí)差剔除異常值后做算術(shù)平均即可得到各炮、檢點(diǎn)的剩余靜校正量.

實(shí)際應(yīng)用中炮點(diǎn)和檢波點(diǎn)可按區(qū)域分組計(jì)算，通過將相同組的偏移道進(jìn)行疊加，可以有效壓制噪聲提高信噪比；對(duì)于地表情況比較復(fù)雜的數(shù)據(jù)則不宜采用大范圍分組，因?yàn)橄噜徟诨驒z波點(diǎn)可能存在比較大的差別從而改變了成像道集的波形，不利于互相關(guān)準(zhǔn)確求取剩余時(shí)差.

3實(shí)際數(shù)據(jù)測試

為了驗(yàn)證本文方法及流程的實(shí)用性，本文給出兩個(gè)實(shí)際數(shù)據(jù)算例，一個(gè)是采單線接收的2.5D實(shí)際地震數(shù)據(jù)，另一個(gè)是三維采集的實(shí)際數(shù)據(jù).

3.1單線接收數(shù)據(jù)算例

該數(shù)據(jù)共560炮，每炮240道接收，記錄長度為6 s，采樣間隔為4 ms.圖2是疊前時(shí)間偏移得到的偏移疊加剖面.該剖面在成像區(qū)域內(nèi)連續(xù)性較好，波形穩(wěn)定，對(duì)其進(jìn)行去噪后將它作為剩余靜校正處理的參考道集.圖中用紅線標(biāo)注的是選取的計(jì)算時(shí)窗的中心位置，黃色的長方形方框是用于放大對(duì)比的主要區(qū)域.

圖3展示了隨偏移距變化的共反射點(diǎn)道集、隨炮點(diǎn)組變化的成像道集和隨檢波點(diǎn)組變化的成像道集.從成像道集上可以觀察出同相軸已比較平直，但局部還存在剩余靜校正量.圖4是通過互相關(guān)計(jì)算出的每個(gè)炮點(diǎn)組在不同CDP點(diǎn)位置處的剩余時(shí)差.可以看出，對(duì)于一個(gè)炮點(diǎn)組或檢波點(diǎn)組來說，相鄰CDP點(diǎn)處的剩余時(shí)差比較接近；相同炮點(diǎn)組和檢波點(diǎn)組的剩余時(shí)差在不同CDP位置處也比較相近，這可以當(dāng)作剩余時(shí)差求取是否正確及異常值剔除是否合理的一個(gè)參考標(biāo)準(zhǔn).圖5是按照不同CDP位置做算術(shù)平均得到關(guān)于炮點(diǎn)組位置和檢波點(diǎn)組位置的剩余靜校正量.圖中炮點(diǎn)組和檢波點(diǎn)組的最大剩余靜校正量都約為8ms，且整體算數(shù)平均值都接近于零，表明該數(shù)據(jù)本身剩余靜校正量比較小，計(jì)算得到的剩余靜校正量也比較合理.

圖2　剩余靜校正參考道集剖面Fig.2　Reference gather section for residual static correction

圖3　基于不同域疊加的成像道集(a) 共反射點(diǎn)道集； (b) 共炮點(diǎn)組成像道集； (c) 共檢波點(diǎn)組成像道集.Fig.3　Imaging gathers based on stacked various domains(a) CRP gather； (b) Common shot gather； (c) Common receiver gather.

圖4　炮點(diǎn)組在不同CDP點(diǎn)位置處的剩余靜校正量Fig.4　Residual shot statics located at different CDP

圖5　關(guān)于炮點(diǎn)組位置和檢波點(diǎn)組位置的剩余靜校正量曲線圖(a) 關(guān)于炮點(diǎn)組位置； (b) 關(guān)于檢波點(diǎn)組位置.Fig.5　Curves of residual statics(a) On shot position； (b) On receiver position.

圖6　剩余靜校正施加前后成像道集的對(duì)比圖(a) 沒做剩余靜校正； (b) 做了剩余靜校正.Fig.6　Comparison of imaging gather before (a) and after (b) residual static correction

圖7　剩余靜校正后再次偏移成像結(jié)果Fig.7　Migration imaging result after residual static correction

圖8　剩余靜校正前后偏移成像結(jié)果對(duì)比(a) 沒做剩余靜校正； (b) 做了剩余靜校正.顯示的是圖2和7長方形方框的放大結(jié)果.Fig.8　Comparison of migration imaging results before and after residual static correction(a) Before correction； (b) After correction. Enlarged details of rectangles in Fig.2 and 7.

圖10　施加剩余靜校正后測線方向和垂直測線方向的疊前時(shí)間偏移結(jié)果Fig.10　Inline and cross-line imaging results by pre-stack time migration with residual static correction

我們將上述剩余靜校正結(jié)果施加到數(shù)據(jù)中再次進(jìn)行疊前時(shí)間偏移.圖6是本文剩余靜校正施加前后成像道集的對(duì)比圖，可以看出道集質(zhì)量有明顯提高.圖7是本文剩余靜校正施加后的偏移疊加剖面，可以看出信噪比和分辨率都得到了明顯提高，構(gòu)造細(xì)節(jié)更加豐富.圖8是圖7和圖2中方框所在位置的放大對(duì)比圖，可以看出主要構(gòu)造部位有了更清晰地刻畫.

3.2三維實(shí)際數(shù)據(jù)算例

為了進(jìn)一步驗(yàn)證方法的實(shí)用性，我們采用本文方法處理了一個(gè)我國東北部的三維實(shí)際地震資料.三維數(shù)據(jù)的處理流程跟2.5D數(shù)據(jù)類似.圖9是沒有做偏移后剩余靜校正時(shí)沿測線方向和垂直測線方向的疊前時(shí)間偏移結(jié)果；圖10是施加偏移后剩余靜校正后沿測線方向和垂直測線方向的疊前時(shí)間偏移結(jié)果.對(duì)比兩圖可以看出成像效果有了明顯提高.經(jīng)過應(yīng)用本文的偏移后剩余靜校正，得到了更清晰的斷層(橢圓指示區(qū)域)和連續(xù)的同相軸(方框指示區(qū)域).

4結(jié)論

本文將剩余靜校正量求取由疊前的CMP道集推進(jìn)到偏移后的偏移道集，提出了一種基于偏移道集的地表一致性剩余靜校正方法.文中結(jié)合疊前時(shí)間偏移給出了基于偏移道集的剩余靜校正方法.為獲得地表一致性的炮、檢點(diǎn)剩余時(shí)差，發(fā)展了疊前時(shí)間偏移生成炮、檢點(diǎn)偏移道集，剩余時(shí)差拾取等系列方法.炮、檢點(diǎn)偏移道集和剩余時(shí)差拾取僅涉及小時(shí)窗內(nèi)的偏移和互相關(guān)計(jì)算，并不增加大的計(jì)算量.文中給出了詳細(xì)的資料處理流程，也給出二維、三維起伏地表采集資料的實(shí)際應(yīng)用例子.實(shí)際資料處理結(jié)果表明，本文方法明顯改善了疊前時(shí)間偏移成像的分辨率，斷點(diǎn)更加清晰.該方法可容易結(jié)合到現(xiàn)行的地震資料處理流程中.

References

Dong C H, Zhang J F. 2009. Prestack time migration including surface topography.ChineseJ.Geophys. (in Chinese), 52(1): 239-244.Gou L, Peng Z M, Wang K B. 2005. Study of method of hybrid optimum automatic residual static corrections.OGP(in Chinese), 40(4): 454-458. Jing X L, Yang C C, Li Y M, et al. 2002. A global optimized algorithm for seismic residual statics corrections.ChineseJ.Geophys. (in Chinese), 45(5): 707-713.

Larner K L, Gibson B R, Chambers R, et al. 1979. Simultaneous estimation of residual static and crossdip corrections.Geophysics, 44(7): 1175-1192.Marsden D. 1993. Static corrections-a review(part Ⅲ).TheLeadingEdge, 12(3): 210-216.

Rajasekaran S, McMechan G A. 1995. Prestack processing of land data with complex topography.Geophysics, 60(6): 1875-1886.

Reshef M. 1991. Depth migration from irregular surfaces with depth extrapolation methods.Geophysics, 56(1): 119-122.

Ronen J, Claerbout J F. 1985. Surface-consistent residual statics estimation by stack-power maximization.Geophysics, 50(12): 2759-2767.

Rothman D H. 1985. Nonlinear inversion, statistical mechanics, and residual statics estimation.Geophysics, 50(12): 2784-2796. Stork C, Kusuma T. 1992. Hybrid genetic autostatics: new approach for large-amplitude statics with noisy data.∥ Expanded Abstracts of 62nd SEG Meeting, 1127-1131. Taner M T, Koehler F, Alhilali K A. 1974. Estimation and correction of near-surface time anomalies.Geophysics, 39(4): 441-463.

Wiggins R A, Larner K L, Wisecup R D. 1976. Residual statics analysis as a general linear inverse problem.Geophysics, 41(5): 922-938.

Zhang J F, Zhang J J, Zhang H, et al. 2010. Three dimensional pre-stack time migration for anisotropic media. CN: ZL20101059716.0 (in Chinese).Zhao W, Zhang J F, Hao Z J, et al. 2010. Three dimensional true amplitude pre-stack time migration. CN: ZL201010288431.4 (in Chinese).

附中文參考文獻(xiàn)

董春暉, 張劍鋒. 2009. 起伏地表下的直接疊前時(shí)間偏移. 地球物理學(xué)報(bào), 52(1): 239-244.

茍量, 彭真明, 王克斌. 2005. 混合優(yōu)化自動(dòng)剩余靜校正方法研究. 石油地球物理勘探, 40(4): 454-458.

井西利, 楊長春, 李幼銘等. 2002. 地震靜校正全局最優(yōu)化問題的求解. 地球物理學(xué)報(bào), 45(5): 707-713.

張劍鋒, 張江杰, 張浩等. 2010. 各向異性三維疊前時(shí)間偏移方法. 中國專利: ZL201010597160.0.

趙偉, 張劍鋒, 郝振江等. 2010. 三維保幅疊前時(shí)間偏移方法. 中國專利: ZL201010288431.4.

(本文編輯張正峰)

基金項(xiàng)目國家自然科學(xué)基金面上項(xiàng)目(41174112)和國家自然科學(xué)基金重點(diǎn)項(xiàng)目(41330316)聯(lián)合資助.

作者簡介徐錦承,男,1987年生,中國科學(xué)院地質(zhì)與地球物理研究所博士研究生,主要從事疊前偏移成像研究. E-mail:xujincheng@mail.iggcas.ac.cn

doi:10.6038/cjg20160222 中圖分類號(hào)P631

收稿日期2014-03-14，2015-09-07收修定稿

Residual static correction based on migrated gathers

XU Jin-Cheng1,2， ZHANG Jian-Feng1

1InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China2UniversityofChineseAcademyofSciences,Beijing100049,China

AbstractStatic correction is a very important topic in land seismic exploration. Field and datum static correction can eliminate time distortions caused by variable elevations of the source and receiver. However, there still exist small time distortions, which may affect the continuity and resolution of events in migration results. So residual static correction is necessary in the conventional static correction processing. Using migrated gathers to pick up the residual time shifts of the shot and receiver, a new residual static corrections method is presented.

KeywordsResidual static corrections； Pre-stack migration； Migrated gathers

徐錦承， 張劍鋒. 2016. 基于偏移成像道集的剩余靜校正方法.地球物理學(xué)報(bào)，59(2):655-663,doi:10.6038/cjg20160222.

Xu J C, Zhang J F. 2016. Residual static correction based on migrated gathers.ChineseJ.Geophys. (in Chinese),59(2):655-663,doi:10.6038/cjg20160222.