一種基于模型統(tǒng)計(jì)拾取的深水obc資料多次波壓制方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明設(shè)及石油天然氣地震勘探領(lǐng)域，特別是關(guān)于一種基于模型統(tǒng)計(jì)拾取的深水 0BC資料多次波壓制方法。
【背景技術(shù)】
[0002] 多次波是海上地震勘探中最常見的一種規(guī)則干擾波，它源于震源激發(fā)的地震波在 幾個(gè)強(qiáng)反射界面之間的多次反射。多次波容易導(dǎo)致在接收的地震剖面上，來自地下深部目 的層的有效反射信號(hào)被多次波強(qiáng)能量所掩蓋，進(jìn)而使疊加速度場(chǎng)求取不準(zhǔn)確，在偏移剖面 上形成很多假象。為此，行業(yè)內(nèi)先后發(fā)展了預(yù)測(cè)反權(quán)積、二維濾波、SRME(Surface-Related MultipleElimination,表面相關(guān)多次波去除）等一系列壓制多次波的技術(shù)方法。該些方 法都有各自的基本假設(shè)和適用條件，在一定范圍內(nèi)取得了良好的應(yīng)用效果。然而，在對(duì)深水 0BC(0ceanBottomC油le，海底電纜）采集的多波多分量地震資料壓制多次波時(shí)，上述該些 方法卻不能得到滿意的多次波壓制效果。例如預(yù)測(cè)反權(quán)積壓制鳴震，通過設(shè)計(jì)一個(gè)反濾波 算子壓縮地震子波，對(duì)多次波的周期性要求很高，且只能適用于淺水近道的情況。二維濾波 根據(jù)動(dòng)校正后的剩余時(shí)差分離一次波和多次波，先要拾取均方根速度，再到變換域中濾除 多次波，可是0BC資料的放炮點(diǎn)和接收點(diǎn)不在同一個(gè)基準(zhǔn)面上，當(dāng)海水較深時(shí)，靜校正誤差 很大，如圖1所示，從（d)中可W看出0BC資料與拖纜采集資料的時(shí)差，既隨偏移距變化，又 隨反射層深度變化，埋深越淺，變化越大，最大時(shí)差為海水層的*。/2 (t。是海水層的自激自 收時(shí)間）。所從通過簡(jiǎn)單的高程靜校正，將0BC資料的接收點(diǎn)校正到海面，會(huì)產(chǎn)生很大的誤 差，且海水越深，誤差越大。該會(huì)導(dǎo)致均方根速度拾取不準(zhǔn)確，從而影響對(duì)動(dòng)校剩余時(shí)差的 區(qū)分，而且無論F-K變換還是Radon變換都會(huì)引入假頻，降低資料的信噪比。且由于近道的 一次波與多次波剩余時(shí)差很小，所W該方法對(duì)近道效果不好，容易傷害到有效波。此外，對(duì) 0BC橫波分量而言，縱波激發(fā)橫波接收，其均方根速度拾取本身也是十分困難的。SRME是目 前海上拖纜采集資料去多次波效果最好，應(yīng)用最廣的一項(xiàng)技術(shù)，它最大的特點(diǎn)是無需任何 地下信息，僅利用地震數(shù)據(jù)自身進(jìn)行時(shí)空權(quán)積來預(yù)測(cè)多次波。實(shí)際應(yīng)用中，要注意數(shù)據(jù)規(guī)則 化和提供近道數(shù)據(jù)。可是，將SRME用于深水0BC資料同樣得不到滿意的多次波壓制效果。 該是因?yàn)?，SRME將多次波看作是由兩段在"地表反射點(diǎn)"相連的一次反射路徑組成的（"地 表反射點(diǎn)"位置射線傳播遵循Snell定律）。而0BC數(shù)據(jù)是在海面放炮，海底接收，缺少了從 接收點(diǎn)到炮點(diǎn)的射線路徑。雖然可W通過靜校正將海底電纜校正到炮點(diǎn)所在的基準(zhǔn)面上， 但當(dāng)海水較深時(shí)，靜校正誤差很大，且不滿足地震波傳播規(guī)律，該使得SRME對(duì)0BC資料的中 遠(yuǎn)道多次波壓制效果很差，如圖2所示，隨偏移距增加，預(yù)測(cè)的多次波誤差迅速增大。
[0003] 由此可見，目前生產(chǎn)中還缺乏有效的針對(duì)0BC資料的多次波壓制技術(shù)。將常規(guī)多 次波壓制技術(shù)應(yīng)用到0BC資料的主要困難是；①炮點(diǎn)和接收點(diǎn)不在同一個(gè)基準(zhǔn)面上；②靜 校正量時(shí)變空變；⑨射線路徑不對(duì)稱；④轉(zhuǎn)換波速度拾取困難；⑥多次波能量越強(qiáng)壓制越 困難。現(xiàn)今主流的多次波壓制技術(shù)包都含多次波預(yù)測(cè)和多次波相減兩個(gè)步驟。而其中的多 次波預(yù)測(cè)的準(zhǔn)確性常常是方法成功的關(guān)鍵。

【發(fā)明內(nèi)容】

[0004] 針對(duì)上述問題，本發(fā)明的目的是提供一種能夠有效壓制深水OBC地震資料中多次 波干擾的基于模型統(tǒng)計(jì)拾取的深水OBC資料多次波壓制方法。
[0005] 為實(shí)現(xiàn)上述目的，本發(fā)明采取W下技術(shù)方案；一種基于模型統(tǒng)計(jì)拾取的深水OBC 資料多次波壓制方法，包括W下步驟；1)對(duì)原始地震資料進(jìn)行預(yù)處理；2)對(duì)預(yù)處理后的地 震資料進(jìn)行道集分選，使得多次波在分選后的地震道集上呈現(xiàn)周期性特征或在速度譜上呈 現(xiàn)低速特征，其中，對(duì)于2-D地震資料，按帶符號(hào)的共偏移距剖面進(jìn)行道集分選；對(duì)于3-D地 震資料，按CDP進(jìn)行道集分選；3)通過2-D共偏移距道集剖面上反射同相軸的周期性特征， 或通過3-D共CDP道集剖面結(jié)合速度譜能量團(tuán)信息識(shí)別出低速的多次波強(qiáng)能量，確定地震 資料中所有需要壓制的多次波同相軸；4)讀取其中一個(gè)道集，在該道集上拾取出需要壓制 的多次波同相軸；5)平滑拾取的多次波同相軸，濾除多次波同相軸中偏離的局部極值；6) W平滑后的多次波同相軸為中線，對(duì)共偏移距道集上的各道做靜態(tài)時(shí)移，將平滑后的多次 波同相軸拉平，再沿多次波同相軸分時(shí)窗統(tǒng)計(jì)時(shí)窗內(nèi)各道多次波的振幅特征，建立多次波 模型；7)從拉平的地震資料中減去步驟6)所建立的多次波模型，得到壓制多次波后的一次 波數(shù)據(jù)；8)將步驟7)得到的一次波數(shù)據(jù)沿步驟6)中的多次波同相軸中線反向靜態(tài)時(shí)移， 校正回原來的位置，得到壓制多次波后的地震道集；9)W步驟5)得到的平滑后的多次波同 相軸為參考，在下一個(gè)道集上重復(fù)步驟4)~步驟9)，直至完成地震資料全部道集的多次波 壓制處理，得到壓制多次波后的地震資料。
[0006] 所述步驟1)中的預(yù)處理包括去除壞炮、壞道，地震道極性反轉(zhuǎn)校正，去除隨機(jī)噪 音，球面擴(kuò)散補(bǔ)償，炮集間的能量均衡。
[0007] 所述步驟4)在該道集上拾取出需要壓制的多次波同相軸，根據(jù)多次波同相軸特 點(diǎn)分為W下兩種情況；①針對(duì)海底一次反射波能量強(qiáng)、同相軸連續(xù)的海底多次波，首先，使 用計(jì)算機(jī)自動(dòng)拾取海底一次反射波；然后，根據(jù)海底一次反射波的到達(dá)時(shí)間T。和當(dāng)前地震 道集的偏移距ofs計(jì)算出海底一階、二階多次波的達(dá)到時(shí)間Ti和T2，
[0010] 其中，V代表聲波在海水中的傳播速度；最后，分別沿Ti和T2在上下各1/4子波波 長(zhǎng)范圍內(nèi)根據(jù)多次波同相軸的子波極性和振幅相位特征，自動(dòng)拾取海底一階、二階多次波 同相軸；②針對(duì)局部強(qiáng)能量多次波，首先，在多次波同相軸上按相等的水平間隔人工拾取若 干參考點(diǎn)，W及人工拾取多次波同相軸的左右端點(diǎn)，人工拾取參考點(diǎn)時(shí)，在其它強(qiáng)能量反射 軸斜交所拾取的多次波同相軸的交點(diǎn)附近適當(dāng)添加參考點(diǎn)；然后，采用樣條插值計(jì)算參考 基線，并根據(jù)多次波同相軸的子波極性、振幅相位特征或時(shí)窗內(nèi)相鄰道的相關(guān)分析信息，沿 著參考基線逐點(diǎn)追蹤出要拾取的多次波同相軸。
[0011] 所述步驟5)中平滑拾取的多次波同相軸時(shí)，采用中值濾波。
[0012] 所述步驟6)中統(tǒng)計(jì)時(shí)窗內(nèi)各道多次波的振幅特征采用平均值或中值疊加方法。
[0013] 本發(fā)明由于采取W上技術(shù)方案，其具有W下優(yōu)點(diǎn)；1、本發(fā)明包括多次波模型建立 和相減，即首先根據(jù)道集上多次波的走時(shí)特征識(shí)別并解釋多次波，再沿同相軸統(tǒng)計(jì)多次波 振幅特征，采用中值疊加或統(tǒng)計(jì)方法建立多次波模型，然后從原始地震資料中減去多次波， 從而達(dá)到多次波壓制的目的，可W處理放炮點(diǎn)位置與接收點(diǎn)位置不在同一個(gè)水平面上的情 況，避開了二維濾波和SRME等常規(guī)方法處理深水OBC資料時(shí)需要面對(duì)的靜校正和速度拾 取兩個(gè)難題，不僅適用于深水OBC地震資料，對(duì)常規(guī)海上拖纜采集地震資料也同樣有效。如 圖3所示，對(duì)Marmousi-II模型OBC縱波資料壓制海底多次波處理，該模型海水水深450m， 海底水平，成像剖面中的海底多次波得到了有效壓制，此時(shí)，諸如預(yù)測(cè)反權(quán)積、二維濾波和 SRME等常用方法都不滿足假設(shè)條件，不能取得很好的多次波壓制效果。2、本發(fā)明采用不同 的方式對(duì)地震資料進(jìn)行分選，對(duì)于