杜 偉

1 研究區(qū)概況

K1地區(qū)為沙丘、沙梁地貌，目的層深4 000～6 000 m，地表巨厚的沙丘對(duì)地震信號(hào)具有強(qiáng)烈的吸收衰減作用，尤其對(duì)地震波的高頻成分吸收作用更為強(qiáng)烈，造成地震資料中深層分辨率、信噪比低，地震資料解釋困難。由于地震部署時(shí)間的差異和地震采集參數(shù)的不同，在連片處理后，橫向上，不同區(qū)塊三維地震數(shù)據(jù)交界處出現(xiàn)一些振幅能量變化不均勻現(xiàn)象，存在明顯的能量分界線。縱向上，目的層段內(nèi)深淺層地震反射能量差異大，局部出現(xiàn)強(qiáng)能量及能量空間方向突變等現(xiàn)象。反射能量的不均衡不僅造成地震數(shù)據(jù)分辨率及信噪比低，同時(shí)也使地震屬性分析、反演結(jié)果失真，嚴(yán)重制約了后續(xù)地震資料的構(gòu)造解釋和儲(chǔ)層預(yù)測(cè)研究。

針對(duì)地震數(shù)據(jù)體存在的這些問(wèn)題，制定了能量均衡處理、提高分辨率處理和提高信噪比處理“三步走”的處理流程，其實(shí)現(xiàn)方法分別為層位控制下的振幅增益補(bǔ)償計(jì)算、混合相位子波反褶積和傾角導(dǎo)向控制下的中值濾波。實(shí)際地震資料處理結(jié)果表明，在保證原始數(shù)據(jù)信息的同時(shí)，該處理流程大幅度提高了三維地震資料的品質(zhì)，為后期構(gòu)造解釋及儲(chǔ)層預(yù)測(cè)奠定了基礎(chǔ)。

2 地震數(shù)據(jù)目標(biāo)處理技術(shù)

2.1 能量均衡處理

能量均衡處理主要是為了改善縱橫方向上存在的非地質(zhì)因素引起的地震振幅的強(qiáng)弱差異，在保留振幅相對(duì)強(qiáng)弱關(guān)系的前提下使地震剖面整體振幅趨于一定范圍，從而滿足地震解釋的需要[1]。從理論上來(lái)說(shuō)，要區(qū)分地震振幅的變化是地質(zhì)因素引起還是非地質(zhì)變化引起的是基本不可能的，但根據(jù)對(duì)地震資料采集、處理過(guò)程的了解和地下地質(zhì)情況的掌握，是能夠根據(jù)經(jīng)驗(yàn)對(duì)地震振幅變化的合理性進(jìn)行判斷的。另外，根據(jù)解釋目的的需要，往往需要對(duì)目的層地震振幅進(jìn)行適當(dāng)?shù)哪芰炕謴?fù)及補(bǔ)償。

在疊后地震數(shù)據(jù)處理中，縱向振幅增益補(bǔ)償計(jì)算可有效地均衡地震道淺、中、深層能量的差異，其實(shí)現(xiàn)方法是縱向振幅增益補(bǔ)償計(jì)算。首先，選取目的層頂?shù)追秶?，利用地震解釋層位及?nèi)插層位將目的層劃分為多個(gè)近等厚的小時(shí)窗，選取目的層頂部的小時(shí)窗作為參考時(shí)窗，計(jì)算其振幅值，如平均絕對(duì)振幅、均方根振幅等。其次，計(jì)算其它各小時(shí)窗內(nèi)振幅值，將參考時(shí)窗內(nèi)振幅值分別與其它小時(shí)窗內(nèi)振幅值相除求取比值，也即比例因子。最后，利用比例因子對(duì)其它各時(shí)窗內(nèi)振幅數(shù)據(jù)進(jìn)行相乘，即可將各個(gè)小時(shí)窗內(nèi)的振幅值計(jì)算為統(tǒng)一水平，從而實(shí)現(xiàn)縱向上的能量均衡。

圖1中目的層段厚度大且地層起伏較大，目的層內(nèi)垂向上能量差異較大，利用地震解釋層位作為控制時(shí)窗，進(jìn)行縱向能量均衡計(jì)算，較好地實(shí)現(xiàn)了目的層段內(nèi)整體能量的均衡。

圖1 縱向能量均衡處理前（左）后（右）地震剖面對(duì)比

橫向上的能量不均衡性較為復(fù)雜，除去地質(zhì)因素外，一般體現(xiàn)在兩個(gè)方面：一是三維數(shù)據(jù)體的邊緣部分，相對(duì)于探區(qū)滿覆蓋地區(qū)而言，探區(qū)邊緣因覆蓋次數(shù)不足，容易造成疊加后地震反射能量整體偏弱；二是分塊采集的三維地震數(shù)據(jù)體在拼接過(guò)程中，拼接處容易留下能量差異的痕跡。針對(duì)前者，時(shí)不變道均衡處理可改善數(shù)據(jù)體邊緣能量弱等問(wèn)題；針對(duì)后者，通過(guò)計(jì)算拼接線兩邊振幅比值，利用比例因子法進(jìn)行反計(jì)算，即可消除拼接處的能量不均衡性。圖2中左圖存在明顯的數(shù)據(jù)拼接痕跡，兩邊能量差異明顯，對(duì)兩邊振幅值進(jìn)行比例因子計(jì)算后，兩邊反射能量校正到了同一水平。

圖2 拼接數(shù)據(jù)體能量均衡處理前（左）后（右）地震剖面對(duì)比

2.2 提高分辨率處理

近年來(lái)，進(jìn)行疊后地震數(shù)據(jù)拓頻處理得到了廣泛的應(yīng)用[2-5]，其理論基礎(chǔ)較為完善，技術(shù)方法也逐漸成熟。相對(duì)于常規(guī)處理，疊后地震數(shù)據(jù)拓頻處理具有成本低、見(jiàn)效快等優(yōu)勢(shì)，無(wú)疑是一種提高分辨率的快速補(bǔ)救方法。

本次提高分辨率處理采用的是混合相位子波反褶積技術(shù)，全部過(guò)程由地質(zhì)放大鏡軟件完成。處理的關(guān)鍵是在井資料的控制下，提取理想的混合相位子波。首先，在忽略噪聲影響的情況下，將井旁地震道的頻譜取對(duì)數(shù)轉(zhuǎn)化為線性系統(tǒng)，通過(guò)設(shè)計(jì)低通濾波器，子波和反射系數(shù)在復(fù)賽譜域可以得到分離，從而達(dá)到提取子波的目的[6]。然后，基于譜模擬的思想，根據(jù)地震數(shù)據(jù)的頻譜分布，自定義一個(gè)期望輸出的子波，用雷克子波、余氏子波或人工定義的其它子波，確保信噪比水平在合理范圍[7]。最后，利用期望輸出的子波與提取子波的振幅譜相除，得到反射子波；反射子波與地震道進(jìn)行褶積，得到拓頻后的地震記錄。

圖3是拓頻處理前后地震剖面的對(duì)比，可以看出，拓頻后目的層主頻從25 Hz提高到了35 Hz，頻帶寬度也由8～40 Hz拓寬到了8～50 Hz，從地震剖面對(duì)比來(lái)看，在保持強(qiáng)反射界面不變的同時(shí)，拓頻后地震剖面同相軸更細(xì)，數(shù)量更多，反映的地質(zhì)細(xì)節(jié)更豐富、清晰，縱向上的分辨率得到了提高。

以上過(guò)程是在復(fù)賽譜中交互完成，提取的子波振幅譜準(zhǔn)確客觀，在實(shí)現(xiàn)頻譜拓寬的同時(shí)，還可以使剖面零相位化。由于地震資料分辨率與信噪比存在天然的矛盾性，過(guò)高的拓頻容易帶來(lái)噪聲和產(chǎn)生假的同相軸，故對(duì)疊后資料的拓頻處理通常需要謹(jǐn)慎對(duì)待。提高分辨率處理過(guò)程中尤其要注意質(zhì)量控制，可利用處理前后信噪比譜作參考，利用測(cè)井曲線、合成地震記錄等作監(jiān)控，防止虛假地震同相軸的產(chǎn)生，確保地震資料能夠最大限度地真實(shí)反映地質(zhì)現(xiàn)象，從而找到分辨率和保真度之間的最佳平衡點(diǎn)。

圖3 拓頻處理前（左）后（右）地震剖面對(duì)比

2.3 提高信噪比處理

由于地表及采集因素的影響，原始地震數(shù)據(jù)本身信噪比較低，通過(guò)前述能量均衡處理和提高分辨率處理后，不可避免地會(huì)放大原始噪聲并產(chǎn)生新的噪聲，因此，疊后去噪處理則顯得尤為重要。去噪處理主要基于地震數(shù)據(jù)體的濾波計(jì)算，數(shù)據(jù)濾波有很多種方法[8-10]，各種濾波方法都有一定的優(yōu)點(diǎn)和適用性，考慮到構(gòu)造解釋的需要，選擇基于傾角控制下的構(gòu)造導(dǎo)向中值濾波方法[11]。該方法的特點(diǎn)是沿地層反射傾角方向進(jìn)行追蹤對(duì)比或掃描計(jì)算，能獲得更加準(zhǔn)確的地質(zhì)目標(biāo)信息，在增強(qiáng)地震同相軸連續(xù)性的同時(shí)，對(duì)斷點(diǎn)的保持也得到了加強(qiáng)，更有利于層位和斷層的解釋[12-13]。

圖4是濾波前后地震剖面的對(duì)比，可以看出，濾波前地震剖面上噪聲較多，同相軸連續(xù)性差；濾波處理后，地震剖面上噪音得到明顯壓制，地震剖面更為“干凈”，地震同相軸的連續(xù)性也得到加強(qiáng)，而且地質(zhì)體特征和斷層邊界信息得到了很好的保持。

圖4 濾波處理前（左）后（右）地震剖面對(duì)比

3 效果分析

K1地區(qū)古生界碎屑巖以砂泥沉積為主，單砂體厚度薄，目的層埋深為4 000～6 000 m。地震數(shù)據(jù)體在目的層段主頻較低，約25 Hz，頻帶寬度為8～40 Hz。從原始地震資料品質(zhì)來(lái)看，其在反射能量、縱向分辨率上均無(wú)法滿足構(gòu)造解釋和儲(chǔ)層預(yù)測(cè)研究的需要（圖5左）。經(jīng)過(guò)前述三個(gè)步驟的處理后（圖5右），地震數(shù)據(jù)體在能量均衡性方面有很大改善，分辨率和信噪比均得到了提高，處理后的地震剖面更加協(xié)調(diào)，更有利于構(gòu)造解釋及儲(chǔ)層預(yù)測(cè)。

4 結(jié)論

（1）針對(duì)復(fù)雜三維連片探區(qū)，根據(jù)不同解釋任務(wù)的需要，疊后三維地震數(shù)據(jù)通常需要進(jìn)行目標(biāo)處理。本文中應(yīng)用到的能量均衡處理、拓頻處理及去噪處理都是針對(duì)性的地震資料處理技術(shù)，實(shí)際資料處理結(jié)果表明，將這三者結(jié)合起來(lái)作為一個(gè)完整的處理流程通常能更加顯著地改善原始地震數(shù)據(jù)的品質(zhì)。

圖5 原始地震剖面（左）與目標(biāo)處理后地震剖面（右）對(duì)比

（2）處理過(guò)程中，計(jì)算時(shí)窗的選取、拓頻參數(shù)的定義及去噪方法和參數(shù)的選取都是處理過(guò)程中的關(guān)鍵環(huán)節(jié)，且每一個(gè)處理環(huán)節(jié)都容易抹去有效信息，同時(shí)帶來(lái)虛假信息。如何平衡分辨率與信噪比兩方面的矛盾，不僅涉及到大量參數(shù)試驗(yàn)，還需要很多過(guò)程的質(zhì)量監(jiān)控，甚至包括處理解釋人員的經(jīng)驗(yàn)及對(duì)該區(qū)地質(zhì)特征的認(rèn)識(shí)程度。

（3）相對(duì)于常規(guī)疊前處理而言，疊后地震數(shù)據(jù)體目標(biāo)處理過(guò)程相對(duì)簡(jiǎn)單，而且也有較多成熟商業(yè)軟件支持，具有投入成本低、見(jiàn)效快等優(yōu)勢(shì)，在目前低油價(jià)、追求效益勘探的大背景下，值得推廣應(yīng)用。

參考文獻(xiàn)

[1] 趙志萍，徐輝. 三維地震資料連片的一致性處理技術(shù)[J].油氣地球物理，2003，1（2）：17–20.

[2] 王有濤，陳學(xué)國(guó)，郝志偉，等. 高分辨率反演在春風(fēng)油田儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J]. 斷塊油氣田，2013，20（3）：293–295.

[3] 張會(huì)卿，聶國(guó)振，燕云. 基于拓頻的地質(zhì)統(tǒng)計(jì)反演技術(shù)在馬東東地區(qū)的應(yīng)用[J]. 石油地質(zhì)與工程，2014，28（2）：52–54.

[4] 周宗良，吳剛，肖建玲，等. 地震拓頻處理技術(shù)在歧口深層氣藏儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J]. 工程地球物理學(xué)報(bào)，2015，12（5）：604–609.

[5] 楊瑞召，趙爭(zhēng)光，馬彥龍，等. 利用譜藍(lán)化和有色反演分辨薄煤層[J]. 天然氣地球科學(xué)，2013，24（1）：156–160.

[6] 李國(guó)發(fā)，彭蘇萍，高日勝，等. 復(fù)賽譜域提取混合相位子波的方法[J]. 天然氣工業(yè)，2005，25（1）：85–87.

[7] 姜平，于興河，黃月銀，等. 儲(chǔ)層精細(xì)描述在東方1－1氣田中的應(yīng)用[J]. 地學(xué)前緣，2012，19（2）：87–94.

[8] 袁恩輝，陳鳳來(lái)，齊志斌，等. 奇異值分解在去除特殊干擾波中的研究與應(yīng)用[J]. 非常規(guī)油氣，2015，2（5）：20–25.

[9] 董志華，周家雄，隋波，等. 基于擴(kuò)散濾波方法的小斷層識(shí)別技術(shù)[J]. 石油地球物理勘探，2013，48（2）：758–762.

[10] 仲旭艷，韓其言. 中值濾波及均值濾波在地震異常波提取中的應(yīng)用[J]. 地下水，2017，39（2）：75–88.

[11] 趙麗婭，傅群，潘海濱. OpendTect軟件傾角控制模塊在地震解釋中的應(yīng)用[J]. 海洋地質(zhì)動(dòng)態(tài)，2008，24（4）：37–41.

[12] 尹川，杜向東，趙汝敏，等. 基于傾角控制的構(gòu)造導(dǎo)向?yàn)V波及其應(yīng)用[J]. 地球物理學(xué)進(jìn)展，2014，29（6）：2 818–2 822.

[13] 劉洋，王典，劉財(cái)，等. 局部相關(guān)加權(quán)中值濾波技術(shù)及其在疊后隨機(jī)噪聲衰減中的應(yīng)用[J]. 地球物理學(xué)報(bào)，2011，54（2）：358–366.