屈 元, 刁永波, 李睿寧,2, 張明柱, 宋繼勝
(1.東方地球物理勘探有限公司 西南物探研究院,成都 610213;2.東方地球物理勘探有限公司 西南物探分公司,成都 610213)
近年,隨著經(jīng)濟(jì)形勢(shì)下行壓力增大,石油行業(yè)上游勘探投資日趨縮緊。物探市場(chǎng)競(jìng)爭(zhēng)日益激烈且項(xiàng)目要求越來(lái)越高,石油物探是集技術(shù)密集型與勞動(dòng)密集型于一體的勘探上游行業(yè)。降本增效精細(xì)勘探成為地震勘探新常態(tài)。
地震勘探檢波器組合法是利用干擾波與有效波在傳播方向上的差別提出。檢波器組合作為組合的一種重要形式,在壓制干擾波提高資料信噪比起著重要作用,但是由于檢波器組合往往需要一定組合基距,且要求一定組合個(gè)數(shù)。這就使得在地表復(fù)雜工區(qū)中,由于地表高程、地表巖性不一致造成同一組檢波器每一只接收的地震起跳時(shí)間有差異,每只檢波器的耦合性也不同,加之受不同檢波器的靈敏度的誤差等原因,檢波器組合方式體現(xiàn)出局限性[2-3]。靜校正誤差的各個(gè)時(shí)差值實(shí)際上造成了一種濾波作用,并造成各疊加道子波的相位特性的不一致(子波的波形不一致),使水平疊加中不能實(shí)現(xiàn)“同相疊加”。單只高精度檢波器則規(guī)避了上述局限,在復(fù)雜地表山地地區(qū)適用性更強(qiáng),但噪聲發(fā)育更嚴(yán)重。
研究區(qū)域位于塔里木盆地KP地區(qū),為首塊在塔里木盆地山地地表開(kāi)展的單檢接收的地震采集區(qū)。工區(qū)地表高程范圍變化大,山體區(qū)地形切割大,局部高差最大達(dá)840 m(圖1)。
圖1 KP地區(qū)高程截面圖Fig.1 The height of cross section figure of KP area
由于單檢具有噪聲發(fā)育更重且時(shí)常與復(fù)雜地表共生的施工放樣特點(diǎn)。為此,如何針對(duì)復(fù)雜地表?xiàng)l件下的單點(diǎn)檢波器接受的室內(nèi)資料處理是地震資料處理人員面臨的新問(wèn)題。
針對(duì)性處理措施包含以下3個(gè)技術(shù)手段:
靜校正處理在地震資料成像處理中扮演著重要角色,是做好地震資料成像處理的基礎(chǔ)。
勘探實(shí)踐表明層析靜校正對(duì)地形復(fù)雜,地表巖性變化劇烈的工區(qū)適應(yīng)能力更強(qiáng),在塔西南高難山體區(qū)常采用層析反演靜校正。
層析反演中網(wǎng)格大小決定反演模型分辨率的高低,一定限度內(nèi),模型分辨率越高,能夠解決的資料高頻靜校正問(wèn)題就越強(qiáng)。為了滿足覆蓋次數(shù)及分辨率需要,野外采集高精度單點(diǎn)檢波器與小道距接收方式相伴而生。
在室內(nèi)靜校正處理中,通過(guò)對(duì)不同網(wǎng)格層析反演對(duì)比,小網(wǎng)格反演在刻畫(huà)低降速層上更準(zhǔn),利于資料處理人員精細(xì)更加精細(xì)地表層層析模型建立。
在層析模型的建立中,通過(guò)對(duì)不同網(wǎng)格大小層析反演的計(jì)算和應(yīng)用對(duì)比(圖2)。由圖2可以看出,較小網(wǎng)格反演層析模型對(duì)資料低降速層刻畫(huà)更清晰。資料成像效果更佳(圖3)。
在單檢小道距采集的室內(nèi)成像處理中,應(yīng)在確保迭代收斂和反演結(jié)果合理的前提下,盡量選擇小的網(wǎng)格。
山地地震資料由于地表起伏劇烈,單炮資料噪聲線性規(guī)律受到破壞,且起伏地表落差大,丟道現(xiàn)象普遍,造成觀測(cè)系統(tǒng)不規(guī)則。在高難山地地震資料去噪處理中,首先應(yīng)對(duì)資料的規(guī)則噪聲進(jìn)行線性規(guī)律恢復(fù),在噪聲壓制前,對(duì)資料進(jìn)行靜校正處理,如圖4所示。對(duì)比圖4可知,單炮靜校正處理后FK譜噪聲特性得到較好恢復(fù),利于后期線性噪聲壓制。
圖2 不同網(wǎng)格層析模型對(duì)比Fig.2 The contrast of tomographic model in different grid(a)30m*30m;(b)7.5m*7.5m
圖3 不同網(wǎng)格層析疊加對(duì)比Fig.3 The contrast of tomographic stack in different grid(a)30m*30m;(b)7.5m*7.5m
圖4 靜校正應(yīng)用前后單炮及FK譜對(duì)比Fig.4 Static correction of shot and FK spectrum comparison before and after application(a)靜校正應(yīng)用前單炮;(b)靜校正應(yīng)用后單炮;(c)靜校正應(yīng)用前FK譜;(d)靜校正應(yīng)用后FK譜
圖5 NUCNS去噪與FXCN噪聲壓制效果對(duì)比Fig.5 The noise suppression effect compare between NUCNS and FXCN(a)FXCN噪后單炮;(b)NUCNS噪后單炮
在山地單點(diǎn)檢波器采集資料的去噪處理中,采用omega新去噪模塊NUCNS,它是趨于真實(shí)坐標(biāo)信息驅(qū)動(dòng),對(duì)不規(guī)則采集系統(tǒng)在防假頻處理中效果更佳。常規(guī)FXCN去噪高陡山地資料的非規(guī)則采樣單炮進(jìn)行噪聲壓制易產(chǎn)生空間假頻率。NUCNS模塊在噪聲壓制中能夠較好地控制假頻的產(chǎn)生,特別是在丟道較多的區(qū)域,效果更甚(圖5)。
圖6 反向?yàn)V波去噪處理路程框圖Fig.6 The inverse filtering denoising block diagram
除此之外,由于單檢在組合壓噪方面的存在缺失,資料噪聲普遍較“串”接收方式更重。資料處理人員在低信噪?yún)^(qū)地震資料去噪處理中,往往采用較為強(qiáng)化的噪聲壓制參數(shù)進(jìn)行噪聲壓制。但是,這種強(qiáng)化的壓噪方法往往會(huì)造成對(duì)資料有效信號(hào)的更大傷害,如何進(jìn)行有效地信噪分離是一個(gè)難點(diǎn)。在處理此類(lèi)區(qū)域噪聲壓制中,采用“反向?yàn)V波”的思路(圖6),對(duì)噪聲資料輸出一側(cè)反向通過(guò)DIP濾掉低角度信息,再利用不含有低角度的噪聲信息與原始數(shù)據(jù)進(jìn)行取差,從而達(dá)到對(duì)有效信息的保護(hù)(圖7)。
圖7 反向?yàn)V波去噪處理單炮對(duì)比Fig.7 Inverse filtering denoising shot comparison(a)noise1;(b)noise2;(c)差值
在復(fù)雜起伏地表工區(qū),基于地震波垂直上下傳播的條件不成立 ,其共深度點(diǎn)時(shí)距曲線也不可以近似地認(rèn)為是雙曲線,為此引入浮動(dòng)基準(zhǔn)面就是為了滿足在地表復(fù)雜工區(qū),靜校正的計(jì)算[7-8],浮動(dòng)基準(zhǔn)面的選取是在一段或相臨幾個(gè)CMP道集的炮點(diǎn)和接收點(diǎn)所涉及的范圍內(nèi),確定一個(gè)時(shí)間域的地形平均面。在這個(gè)面上,時(shí)距曲線可以近似地認(rèn)為是雙曲線,從而進(jìn)行速度分析、動(dòng)校正與CMP疊加。
單點(diǎn)檢波器的應(yīng)用是為了配合野外高難度探區(qū)的采集施工放樣,可以說(shuō)單檢的應(yīng)用與復(fù)雜地表是共生的。在山地地震資料處理中,由于地形起伏大且丟道現(xiàn)象頻發(fā),因此在地形起伏大或炮檢距較大的情況下,試驗(yàn)合理的平滑參數(shù),選擇合適的浮動(dòng)平滑時(shí)間地形面(CMP疊加參考面)顯得十分必要。常規(guī)圓滑浮動(dòng)面的計(jì)算方法,是對(duì)基于同一CMP道集內(nèi)所有炮檢點(diǎn)信息均參與平均、圓滑。然而,由于山地探區(qū)地表復(fù)雜,造成資料丟道多,特別是在構(gòu)造主體,地形起伏大,野外施工條件惡劣,地震資料近道有效信息少,遠(yuǎn)偏移距成像,基于OMEGA系統(tǒng)自動(dòng)算取的CMP面存在異常,易造成成像虛假。為此采用人工分解方法從新計(jì)算CMP浮動(dòng)面,在計(jì)算CMP浮動(dòng)面中,按一定搜索范圍統(tǒng)計(jì)相關(guān)炮、檢點(diǎn)高程信息進(jìn)行平滑,避免了軟件自動(dòng)按共中心點(diǎn)道集搜取信息在丟道較多區(qū)域造成浮動(dòng)面與地形時(shí)間線時(shí)差過(guò)大,造成時(shí)間域構(gòu)造假象的問(wèn)題。通過(guò)對(duì)比,人工計(jì)算CMP浮動(dòng)面后的地震資料超道集成像質(zhì)量得到進(jìn)一步改善(圖8),疊加剖面時(shí)間域成像可靠性更高(圖9)。值得注意的是人工計(jì)算CMP面搜索半徑的確定,既要一定范圍的圓滑,又要盡量保持與地形線接近。通過(guò)新老資料對(duì)比(圖10),新資料在信噪比和分辨率上都較老資料有較大提升,特別是在構(gòu)造主體部位新處理資料同相軸連續(xù)性更好,成像可靠性更高。
圖8 基于人工分解CMP面與自動(dòng)分解CMP面對(duì)資料速度拾取道集的影響Fig.8 The influence of velocity pick gathers between automatic decomposition and artificial decomposition of CMP surface(a)自動(dòng)計(jì)算;(b)人工計(jì)算
圖9 基于人工分解CMP面與自動(dòng)分解CMP面成像效果對(duì)比Fig.9 The imaging contrast effect between automatic decomposition and artificial decomposition of CMP surface(a)自動(dòng)分解;(b)人工分解
圖10 KP地區(qū)地震資料新老剖面對(duì)比Fig.10 The contrast of new and old section of seismic data in KP area(a)老資料;(b)新資料
通過(guò)本次研究得到以下認(rèn)識(shí):
1)針對(duì)單點(diǎn)檢波器往往具有小面元高密度放樣特點(diǎn),采用小面元層析反演靜校正計(jì)算能夠較好解決資料低頻成像問(wèn)題。
2)針對(duì)單點(diǎn)檢波器具有噪聲發(fā)育的特點(diǎn),采用分步、分域的反向?yàn)V波非規(guī)則噪聲壓制保幅去噪手段。
3)針對(duì)山地單檢與復(fù)雜地表共生的特點(diǎn)采用了基于人工分解的CMP小圓滑面成像技術(shù)等室內(nèi)處理針對(duì)性措施,能夠較好地恢復(fù)地層真實(shí)信息,為地震資料處理后期基于起伏地表偏移打下好的基礎(chǔ)。