楊小龍

(中國(guó)石化華東油氣分公司勘探開(kāi)發(fā)研究院，江蘇 南京 210007)

0 引 言

四川桂花三維頁(yè)巖氣地震勘探研究區(qū)地處四川東南部瀘州境內(nèi)，面積200 km2。該區(qū)位于四川盆地南緣，地表起伏較大，最低海拔521 m，最高海拔1 816 m，相對(duì)高差最大1 300 m，山地靜校正問(wèn)題突出；表層巖性橫向變化快，地表出露巖性多樣，由北向南依次發(fā)育白堊紀(jì)砂泥巖、侏羅紀(jì)砂巖、三疊紀(jì)灰?guī)r、二疊紀(jì)灰?guī)r及志留紀(jì)頁(yè)巖。南部灰?guī)r占比大，采集資料信噪比較低；地下局部地層傾角大，速度場(chǎng)呈切割條帶狀變化，陡傾構(gòu)造精確成像難度較大。

區(qū)內(nèi)地震勘探程度較低，以前未進(jìn)行三維地震勘探，周邊只與少數(shù)二維側(cè)線相交。二維地震數(shù)據(jù)由于采集方法和處理技術(shù)的限制，處理成果斷點(diǎn)不清，難以實(shí)現(xiàn)精細(xì)解釋。

針對(duì)區(qū)內(nèi)灰?guī)r區(qū)低信噪比、高陡構(gòu)造等特點(diǎn)，借鑒南方山地資料近幾年的處理方法研究成果，重點(diǎn)采用山地靜校正、弱信號(hào)保護(hù)技術(shù)、五維數(shù)據(jù)重構(gòu)技術(shù)及基于構(gòu)造約束網(wǎng)格層析建模的疊前深度偏移成像技術(shù)，提高資料處理品質(zhì)，為研究區(qū)地震資料的深度解釋和地震勘探提供可靠的基礎(chǔ)資料。

1 處理方法

1.1 技術(shù)路線

針對(duì)研究區(qū)采集密度低、南部灰?guī)r資料品質(zhì)低、地層傾角大等特點(diǎn)，通過(guò)大量測(cè)試和分析，提出適用于該區(qū)資料的處理方法和思路(圖1)。

圖1 桂花地區(qū)三維地震資料處理關(guān)鍵技術(shù)Fig.1 Key technology of 3D seismic data processing in Guihua area

1.2 山地靜校正

在靜校正方面，由于炮道密度低、研究區(qū)山地地形起伏大、地表出露巖性多變、灰?guī)r區(qū)原始地震資料信噪比低和一致性差等原因，影響了區(qū)內(nèi)的初至統(tǒng)計(jì)效應(yīng)，進(jìn)而影響靜校正的計(jì)算精度，因此需進(jìn)行精細(xì)初至拾取及靜校正參數(shù)的全面測(cè)試，保證靜校正量的準(zhǔn)確度。通過(guò)大量測(cè)試發(fā)現(xiàn)，桂花地區(qū)靜校正對(duì)替換速度較為敏感，4 000 m/s時(shí)南部灰?guī)r區(qū)成像效果最佳；不同軟件差異很大，采用ToModel軟件層析靜校正灰?guī)r區(qū)高陡構(gòu)造的成像效果最好(圖2)。

圖2 桂花地區(qū)靜校正前(a)后(b)剖面對(duì)比圖Fig.2 Profile comparison before (a)and after (b) static correction in Guihua area

1.3 灰?guī)r區(qū)弱信號(hào)保護(hù)處理

研究區(qū)南部出露巖性大部分為三疊紀(jì)及二疊紀(jì)灰?guī)r，激發(fā)接收條件較差，有效反射信號(hào)能量很弱、干擾嚴(yán)重，地震資料的信噪比較低，資料品質(zhì)較差，記錄的有效信息具有弱振幅、弱連續(xù)性等特點(diǎn)(曾梅，2003；李蘇光等，2013；鄧攻等，2015)。采用常規(guī)處理方法時(shí)，有效信號(hào)在數(shù)據(jù)處理中會(huì)隨干擾波一同被壓制，故采用子波極性重構(gòu)補(bǔ)償弱信號(hào)處理技術(shù)(圖3)，建立信噪比模型，保護(hù)弱小信號(hào)。通過(guò)對(duì)灰?guī)r區(qū)低信噪比單炮進(jìn)行弱信號(hào)保護(hù)處理，噪音得到有效壓制，有效地震反射波組清晰可見(jiàn)，資料信噪比得到有效提升(圖4)。

圖3 灰?guī)r地區(qū)弱信號(hào)保護(hù)處理流程Fig.3 Weak signal protection process in limestone area

圖4 灰?guī)r地區(qū)弱信號(hào)保護(hù)處理前(a)后(b)效果圖Fig.4 Effect comparison before (a)and after (b) weak signal protection treatment in limestone area

1.4 五維數(shù)據(jù)重構(gòu)

桂花地區(qū)三維地震采集為變覆蓋次數(shù)觀測(cè)系統(tǒng)，覆蓋次數(shù)從72次至172次變化分布；北部與南部覆蓋次數(shù)存在差異，南部資料信噪比較低。為解決地震資料成像質(zhì)量問(wèn)題，采用MPFI五維數(shù)據(jù)重構(gòu)等處理方法，提高南部灰?guī)r資料信噪比，確保成像效果。五維數(shù)據(jù)規(guī)則化通過(guò)調(diào)整橫向采樣間隔，提高地下復(fù)雜繞射波場(chǎng)及高陡傾角構(gòu)造的采樣密度，降低偏移成像產(chǎn)生的假頻信息和偏移噪聲。通過(guò)傅立葉變換將數(shù)據(jù)從時(shí)間空間域轉(zhuǎn)到空間頻率域，選擇不同維度信息進(jìn)行數(shù)據(jù)插值重構(gòu)。五維插值充分利用主測(cè)線、聯(lián)絡(luò)線、時(shí)間、偏移距和方位角5個(gè)維度信息(崔永福等，2016；齊鵬等，2018)。時(shí)間軸上由于對(duì)地震資料采用基本均勻采樣，因此無(wú)需進(jìn)行指定測(cè)試；主測(cè)線與聯(lián)絡(luò)線窗口參數(shù)分別定義主測(cè)線和聯(lián)絡(luò)線方向的炮點(diǎn)—檢波點(diǎn)長(zhǎng)度，會(huì)影響每組迭代過(guò)程中映射到傅里葉域的系數(shù)數(shù)量和相應(yīng)系數(shù)的計(jì)算結(jié)果，因此需進(jìn)行測(cè)試。采用31×31窗口。

采用五維數(shù)據(jù)重構(gòu)方法，重點(diǎn)提高砂巖區(qū)資料信噪比一致性，處理后的剖面(圖5)顯示，各波場(chǎng)信息更豐富，能量、信噪比一致性更好,該法適用于研究區(qū)不規(guī)則觀測(cè)系統(tǒng)、構(gòu)造復(fù)雜陡傾角的數(shù)據(jù)成像。

圖5 五維數(shù)據(jù)重構(gòu)處理前(a)后(b)剖面對(duì)比圖Fig.5 Profile comparison before (a)and after (b) five-dimensional data reconstruction processing

1.5 基于構(gòu)造約束網(wǎng)格層析建摸

研究區(qū)具有地下構(gòu)造簡(jiǎn)單、斷裂不發(fā)育的特點(diǎn)，但傾角較陡；淺層資料信噪比及有效覆蓋次數(shù)較低，層析建模求解方程中多解性大，迭代速度較慢，獲取的速度模型不夠準(zhǔn)確。以縱向速度變化明顯的地震標(biāo)志層位為約束，由淺至深逐層進(jìn)行速度反演，可以建立該區(qū)更加準(zhǔn)確的速度模型(任俊興等，2020)。

1.5.1 初始速度建模 按照速度建模流程(圖6)，首先結(jié)合地質(zhì)認(rèn)識(shí)和構(gòu)造解釋對(duì)時(shí)間域地震剖面進(jìn)行層位標(biāo)定，以區(qū)內(nèi)龍馬溪組、梁家店組等大套地層為框架建立構(gòu)造模型，同時(shí)進(jìn)行平滑處理消除異常點(diǎn)，然后對(duì)速度進(jìn)行時(shí)深轉(zhuǎn)化，形成深度域的速度場(chǎng)，結(jié)合構(gòu)造模型的約束生成初始深度域速度模型(任俊興等，2020)。圖7顯示，速度模型尚不能真實(shí)反映地下的速度場(chǎng)，因此需要后續(xù)速度迭代和優(yōu)化處理。

圖6 疊前深度偏移初始速度模型建立流程Fig.6 Flow chart of initial velocity model of pre-stack depth migration

圖7 構(gòu)造層位(a)和初始速度模型(b)(深度域)Fig.7 Structural horizon (a)and initial velocity model (b)(depth domain)

1.5.2 剩余曲率拾取 拾取反射波同相軸剩余曲率之前，在不改變反射同相軸剩余曲率的前提下對(duì)CRP道集進(jìn)行凈化處理，以提高道集信噪比。分構(gòu)造區(qū)域?qū)Σ煌氖Ｓ嗲适叭?shù)進(jìn)行優(yōu)化，例如根據(jù)不同的構(gòu)造類型設(shè)置不同的剩余曲率拾取范圍、地層傾角范圍、相似值等，以適應(yīng)剩余曲率拾取的需求(袁剛等，2013；任俊興等，2020)。利用地質(zhì)層位和過(guò)井控制，選出適用于研究區(qū)的剩余曲率值(圖8)進(jìn)行網(wǎng)格層析反演。

圖8 道集處理前(a)后(b)剩余曲率拾取對(duì)比圖Fig.8 Comparison of residual curvature picking before (a)and after (b)gather processing

1.5.3 速度模型優(yōu)化 首先采用無(wú)構(gòu)造約束的網(wǎng)格層析反演進(jìn)行速度迭代，形成橫向速度變化較準(zhǔn)確的速度模型；再進(jìn)行構(gòu)造約束的網(wǎng)格層析反演，提高速度模型的縱向反演精度(李來(lái)林等，2004；胡英等，2006)；以此速度為疊前深度偏移的最終速度模型，進(jìn)行最終疊前深度偏移成像處理。

研究區(qū)屬南方山地，區(qū)內(nèi)普遍發(fā)育中古生界地層，地質(zhì)原因?qū)е碌卣鸩ǚ瓷渌俣瓤v向變化快且產(chǎn)生速度倒轉(zhuǎn)，因此在速度反演過(guò)程中，首先以龍?zhí)督M底為界面，先反演上覆地層速度，直至以上地層的CRP道集同相軸拉平，然后以梁家組底為界面進(jìn)行第二輪反演。按照逐層反演的思路，對(duì)4套標(biāo)志層進(jìn)行逐層反演，直至所有層的同相軸均被拉平，得到的速度模型會(huì)更準(zhǔn)確(袁剛等，2013；任俊興等，2020)。構(gòu)造約束速度建模前后剖面對(duì)比圖(圖9)顯示，構(gòu)造約束速度剖面與地質(zhì)構(gòu)造吻合度更高，更能反映地下真實(shí)速度場(chǎng)，利于資料的精細(xì)解釋。

圖9 構(gòu)造約束前(a)后(b)速度剖面對(duì)比圖Fig.9 Velocity profile comparison before (a) and after (b)structural constraints

2 應(yīng)用效果

通過(guò)山地靜校正、弱信號(hào)保護(hù)處理及五維數(shù)據(jù)重構(gòu)，有效改善了地震剖面有效頻段信噪比——從2提升到3(圖5、圖10)。采用基于構(gòu)造約束網(wǎng)格層析疊前深度偏移方法解析后，主要目的層龍馬溪組等的波組特征清晰連續(xù)、構(gòu)造完整、保真度高，能滿足區(qū)內(nèi)構(gòu)造解釋的需要。通過(guò)仁頁(yè)X井層位標(biāo)定發(fā)現(xiàn)，目的層波組與實(shí)際地層對(duì)應(yīng)良好(圖11)，實(shí)鉆與地震層位相差30 m(在2%誤差以內(nèi))，符合預(yù)期，可以支撐后續(xù)井位部署。茅口組底面地震異常反射波具上拱、下拉、扭動(dòng)等特征，均方根振幅值中等—低(圖12)。結(jié)合振幅屬性及地震反射波形，初步圈定異常區(qū)域若干個(gè)，具有一定的勘探意義。

圖10 處理前后剖面有效頻段信噪比對(duì)比圖Fig.10 SNR comparison of effective frequency band profile before (red line)and after (blue line)processing

圖11 構(gòu)造約束疊前深度偏移剖面井位對(duì)比圖Fig.11 Well location comparison of pre-stack depth migration profile with structural constraints

圖12 三維地震反射異常點(diǎn)均方根振幅屬性圖Fig.12 Root mean square amplitude attribute diagram of three dimensional seismic reflection anomaly

3 結(jié) 論

(1)通過(guò)大量參數(shù)測(cè)試和軟件對(duì)比發(fā)現(xiàn)，ToModel軟件層析靜校正更適用于解析桂花地區(qū)山地資料，替換速度為4 000 m/s時(shí)灰?guī)r區(qū)成像效果最好。

(2)采用弱信號(hào)保護(hù)技術(shù)處理灰?guī)r區(qū)資料，更大程度地保留了原始地震有效反射波組信息，提高了資料的利用率；采用五維數(shù)據(jù)重構(gòu)處理方法解決因覆蓋次數(shù)較低和偏移距分配不均等引起的研究區(qū)南北資料品質(zhì)差異，信噪比及資料一致性得到改善，信噪比提高了1個(gè)點(diǎn)，避免了偏移剖面的畫弧現(xiàn)象。

(3)采用構(gòu)造約束速度建模和疊前深度偏移成像技術(shù)能夠保證成像效果，查明了桂花區(qū)塊的構(gòu)造細(xì)節(jié)、斷層分布及地層展布特征；獲得了區(qū)內(nèi)較連續(xù)的五峰組—龍馬溪組地震反射波組，滿足灰?guī)r高陡構(gòu)造準(zhǔn)確成像的要求，實(shí)鉆和地震剖面深度誤差控制在2%內(nèi)。此次三維處理成果可進(jìn)行構(gòu)造解釋和精細(xì)評(píng)價(jià)，為井位部署提供支撐。