徐 浩，李 赫

（1.中國(guó)石化西北油田分公司勘探開發(fā)研究院，新疆 烏魯木齊830011；2.中國(guó)石油集團(tuán)東方地球物理勘探有限公司西南物探研究院，四川 成都610041）

地震沉積學(xué)的概念由ZENG等（1998）［1-2］首次提出，主要是應(yīng)用地震平面屬性特征來(lái)識(shí)別巖石及砂體宏觀特征、研究沉積體系發(fā)育演化，并彌補(bǔ)地震資料縱向分辨率的局限性，是目前薄層、薄互層砂體平面展布預(yù)測(cè)的重要方法之一，利用該方法對(duì)儲(chǔ)集體的內(nèi)幕結(jié)構(gòu)分析已經(jīng)成為巖性油氣藏勘探研究的主要突破口（陳旭，2010；劉化清等，2018）［3-4］。相位轉(zhuǎn)換技術(shù)、地層切片技術(shù)、地震數(shù)據(jù)分頻技術(shù)是地震沉積學(xué)的核心技術(shù) （董春梅等，2006； 曾洪流，2011）［5-6］。 近年來(lái)，國(guó)內(nèi)外學(xué)者在此經(jīng)典技術(shù)的基礎(chǔ)上輔以拓頻處理和地震屬性預(yù)測(cè)和反演等技術(shù)，在薄層砂體識(shí)別中得到了較好的運(yùn)用效果（常少英等，2015；黃林軍等，2012）［7-8］。

隨著勘探程度的日益提高， 傳統(tǒng)的反演手段已不能滿足儲(chǔ)層預(yù)測(cè)的需要， 基于雷克子波和楔形模型，構(gòu)建儲(chǔ)層厚度、振幅、頻率三者之間關(guān)系的AVF分頻反演技術(shù)應(yīng)運(yùn)而生，并得到了廣泛的運(yùn)用（季玉新等，2010；張廣智等，2015；姚振岸，2017；李偉等，2017）［9-12］。 此反演方法是將分頻信息和支持向量機(jī)的反演技術(shù)結(jié)合發(fā)展而成，優(yōu)勢(shì)在于可對(duì)地震數(shù)據(jù)進(jìn)行分頻分析，能充分運(yùn)用地震數(shù)據(jù)中的相對(duì)低頻和相對(duì)高頻成分， 通過(guò)向量機(jī)學(xué)習(xí)獲取其中有效信息，減少薄層砂體反演的不確定性，具備較高的分辨率（于建國(guó)等，2006；張延章等，2006；范洪軍等，2007）［13-15］。

托甫臺(tái)北區(qū)塊位于塔里木盆地塔河油田主體區(qū)，白堊系舒善河組為重要的儲(chǔ)油層系之一，該層系主要以辮狀河三角洲相沉積為主（王峻等，2012）［16］。在實(shí)際的勘探開發(fā)中，該層系面臨砂巖儲(chǔ)層厚度薄、橫向發(fā)育規(guī)模小且變化大、對(duì)比追蹤困難等難題；而舒善河組地震資料主頻低、信噪比和分辨率低、內(nèi)幕反射連續(xù)性差、低頻資料缺失較為嚴(yán)重，因此，薄層砂體地震反演及儲(chǔ)層預(yù)測(cè)難度較大（徐麗萍，2009）［17］。砂體和強(qiáng)振幅對(duì)應(yīng)關(guān)系差，使得傳統(tǒng)的振幅識(shí)別砂體方法在該層系無(wú)法應(yīng)用； 舒善河組內(nèi)幕砂泥巖從聲波時(shí)差（AC）曲線、密度（DEN）曲線以及縱波阻抗上無(wú)法區(qū)分開， 導(dǎo)致常規(guī)的反演無(wú)法在本區(qū)發(fā)揮作用，制約了薄砂體儲(chǔ)層的有效預(yù)測(cè)。本文以托甫臺(tái)北白堊系舒善河組為例，利用AVF分頻屬性反演方法，首先從巖石物理參數(shù)及基礎(chǔ)地震資料入手， 篩選目標(biāo)曲線， 對(duì)原始地震資料進(jìn)行拓頻處理； 然后通過(guò)BP神經(jīng)網(wǎng)絡(luò)算法獲得有效的高分辨率反演數(shù)據(jù)體，并在此基礎(chǔ)上識(shí)別和追蹤有利砂體。

1 概況

托甫臺(tái)北三維地震工區(qū)位于阿克庫(kù)勒凸起西南斜坡，哈拉哈塘凹陷東部。 區(qū)內(nèi)70余口井揭示白堊系地層，已發(fā)現(xiàn)的白堊系舒善河組油氣藏為受斷裂及構(gòu)造控制的巖性油氣藏（圖1）。 受印支-燕山構(gòu)造運(yùn)動(dòng)影響，中新生界地層發(fā)生反翹，塔河地區(qū)中新生界地層南東抬北西降。其壓應(yīng)力方向?yàn)镹E-SW向，持續(xù)性的穩(wěn)定區(qū)域擠壓導(dǎo)致局部發(fā)育張扭與壓扭作用以及工區(qū)南部下石炭統(tǒng)鹽丘的鹽拱作用，形成雁列型、地塹型、Y字型、階梯狀、階梯狀地塹等斷裂組合樣式，中生界斷裂成排成帶分布。 舒善河組地層厚約230 m，縱向上可劃分為下、中、上三段，出油層主要位于中段；沉積類型主要為辮狀河三角洲前緣亞相水下分流河道、席狀灘壩及湖泊相濱淺湖沉積，垂向上砂泥互層沉積，單砂層或砂層組橫向變化快，單砂體厚度薄，多數(shù)單砂體厚度集中在4～8 m之間；儲(chǔ)層巖性以細(xì)砂巖、極細(xì)砂巖為主，粉砂巖次之，儲(chǔ)層巖性細(xì)，但其儲(chǔ)集、滲濾能力較強(qiáng)，是巖性圈閉發(fā)育的有利層段。

圖1 托甫臺(tái)北工區(qū)位置

2 AVF分頻反演原理

不同厚度的地層，其調(diào)諧頻率不同。對(duì)于一個(gè)楔狀砂體模型， 用不同主頻的雷克子波與其褶積進(jìn)行正演模擬，得到一系列合成地震數(shù)據(jù)體，并可得到振幅與砂體厚度在不同頻率下的調(diào)諧曲線（圖2），建立不同時(shí)間厚度下振幅隨頻率變化（AVF）的關(guān)系（圖3）（于建國(guó)等，2006；李偉等，2017）［12-13］。 由AVF關(guān)系可知， 相同地層在不同的主頻頻率子波下會(huì)展現(xiàn)不同振幅特征的重要規(guī)律，但這種關(guān)系非常復(fù)雜，很難用一個(gè)顯示函數(shù)表示，需用支持向量機(jī)（SVM）非線性影射的方法在測(cè)井和地震子波分解剖面上找到這種關(guān)系，利用AVF信息進(jìn)行反演（于建國(guó)等，2006）［13］。

圖2 振幅與厚度在不同頻率時(shí)的調(diào)諧曲線

圖3 不同時(shí)間厚度下振幅隨頻率變化曲線

分頻反演首先要對(duì)地震資料進(jìn)行頻譜分析，掌握資料的主頻信息及有效頻帶范圍，同時(shí)根據(jù)井點(diǎn)數(shù)據(jù)統(tǒng)計(jì)，分析不同顯示級(jí)別井在其目的層的特征頻率。 根據(jù)特征頻率、主頻及頻帶寬度范圍設(shè)計(jì)合適的尺度進(jìn)行分頻，產(chǎn)生不同頻段的數(shù)據(jù)體，從而達(dá)到分頻的目的。 對(duì)于分頻后的數(shù)據(jù)體，利用支持向量機(jī)（SVM）的方法計(jì)算出不同厚度下振幅與頻率（AVF）之間的關(guān)系，將AVF關(guān)系引入反演，從而建立起測(cè)井波阻抗曲線與地震波形間的非線性映射關(guān)系， 并基于Geoscope軟件得到反演結(jié)果（于建國(guó)等，2006）［13］。

3 AVF分頻反演實(shí)現(xiàn)

3.1 拓頻處理

拓頻高分辨率處理技術(shù)是一項(xiàng)全新的地震數(shù)據(jù)處理方法，主要將低頻子波形成的地震數(shù)據(jù)轉(zhuǎn)換為由高頻子波形成的地震數(shù)據(jù)，在保證地震數(shù)據(jù)的時(shí)頻特征及波組特征不發(fā)生變化的情況下拓寬頻帶，提高相應(yīng)的地震分辨率，從而達(dá)到更精確地識(shí)別薄層沉積砂體的目的 （袁紅軍等，2008； 孫夕平等，2010；常少英等，2015）［7，18-19］。 本區(qū)頻譜寬度集中在20～60 Hz之間，主頻28 Hz，原始地震數(shù)據(jù)分辨率較低，地震標(biāo)定原始剖面與合成地震道吻合程度也低。因此，需要對(duì)原始地震資料進(jìn)行拓頻處理。拓頻方法主要有三種：（1）合理保留低頻；（2）相對(duì)增強(qiáng)調(diào)諧能量；（3）增強(qiáng)高頻弱信號(hào)處理。 根據(jù)實(shí)際資料分析，本文主要運(yùn)用調(diào)諧能力增強(qiáng)法進(jìn)行拓頻。

拓頻后地震資料的主頻為35 Hz， 頻帶寬度為20～80 Hz。 從拓頻處理前后地震道與合成地震記錄的對(duì)比情況可以看出（圖4）：拓頻后的資料與合成記錄對(duì)比特征較好。 以下工作均是基于此拓頻體的基礎(chǔ)上進(jìn)行的。

圖4 拓頻前后地震道與合成地震記錄對(duì)比

3.2 特征頻率及曲線優(yōu)選

主要從兩個(gè)方面優(yōu)選特征頻率：（1） 有油氣顯示井和無(wú)油氣顯示井的頻率區(qū)分；（2） 拓頻后的資料主頻率。 兩者相結(jié)合選擇有效頻率刻畫砂體及含油氣性特征。 通過(guò)對(duì)研究區(qū)白堊系舒善河組46口鉆井的頻率進(jìn)行統(tǒng)計(jì)分析，結(jié)果表明油氣顯示井頻率主要集中在33～38 Hz，平均頻率35 Hz，而27口無(wú)油氣顯示井平均頻率為48 Hz左右，兩者差異較大（表1）；拓頻后地震資料的主頻為35 Hz。 綜上分析，本次分頻反演的特征頻率確定為35 Hz。

從工區(qū)內(nèi)目的層段砂泥巖的巖石物理參數(shù)分析可以看出（圖5），砂泥速度、密度和波阻抗差別均較小，直接運(yùn)用波阻抗進(jìn)行反演，取得的效果不能達(dá)到預(yù)期目標(biāo)。 從該工區(qū)內(nèi)70口井的測(cè)井資料綜合分析，自然電位（SP）曲線對(duì)砂泥區(qū)分響應(yīng)特征較好。

表1 研究區(qū)目的層不同顯示級(jí)別與頻率對(duì)應(yīng)關(guān)系表

圖5 目的層砂泥巖巖石物理參數(shù)分析

3.3 反演實(shí)現(xiàn)過(guò)程

首先對(duì)工區(qū)內(nèi)參與反演的67口井的SP曲線進(jìn)行歸一化處理，在保證一定區(qū)域內(nèi)砂泥巖反射特征相似的前提下，進(jìn)行合成記錄的制作和調(diào)整，得到合理的目標(biāo)曲線。 另外通過(guò)對(duì)地震資料拓頻得到與合成記錄對(duì)應(yīng)關(guān)系較好的拓頻數(shù)據(jù)體，并對(duì)拓頻數(shù)據(jù)體作外部屬性運(yùn)算，如瞬時(shí)屬性、分頻屬性、頻率屬性等，也可建立模型如相控模型或者井控模型作為外部屬性對(duì)反演過(guò)程進(jìn)行約束。 然后選擇合適的學(xué)習(xí)方法（BP神經(jīng)網(wǎng)絡(luò)、支持向量機(jī)和粒子群算法）分頻段建立非線性映射關(guān)系，最后得到反演結(jié)果。

本次參與反演的數(shù)據(jù)體為拓頻數(shù)據(jù)體，67口測(cè)井?dāng)?shù)據(jù)，外部屬性為井控模型。 選取工區(qū)內(nèi)25口井建立三條連井剖面， 通過(guò)聯(lián)井小層對(duì)比的方法，建立砂體識(shí)別模式，并將其作為井控模型約束整個(gè)數(shù)據(jù)體進(jìn)行反演。 AVF分頻反演選取的方法為BP神經(jīng)網(wǎng)絡(luò)，其優(yōu)勢(shì)是能夠通過(guò)學(xué)習(xí)獲得多個(gè)輸入和輸出數(shù)據(jù)的映射關(guān)系，并利用這種網(wǎng)絡(luò)關(guān)系得到對(duì)應(yīng)的輸出結(jié)果。 此種方法運(yùn)用過(guò)程中，調(diào)節(jié)算法和參數(shù)，通過(guò)檢查學(xué)習(xí)曲線和目標(biāo)曲線的相關(guān)度，決定學(xué)習(xí)效果好壞，從而得到滿意的反演結(jié)果。

4 應(yīng)用效果分析

在上述工作的基礎(chǔ)上，完成了工區(qū)內(nèi)67口井白堊系舒善河組的分頻反演，得到了反演數(shù)據(jù)體。 后期的效果分析工作主要從兩個(gè)方面進(jìn)行：（1） 參與反演的井與反演體的吻合度分析；（2） 泥巖百分比含量曲線讀取的砂巖百分比與反演體得到的砂體厚度之間的擬合關(guān)系。

從反演結(jié)果可以看出（圖6），反演結(jié)果與SP曲線吻合度較高。 反演剖面的縱橫向分辨率高，砂體橫向變化合理，尤其是在單砂體的刻畫中，達(dá)到了預(yù)期目標(biāo)。 從選取的兩口未參與反演的井驗(yàn)證來(lái)看，曲線與反演數(shù)據(jù)吻合度較高，主要砂體的識(shí)別程度較好，證實(shí)了反演的可靠性。

圖6 過(guò)驗(yàn)證井的反演剖面

此外， 拾取單井目的層段泥質(zhì)含量曲線數(shù)值，與反演數(shù)據(jù)體上拾取的砂巖厚度數(shù)據(jù)進(jìn)行擬合分析，結(jié)果如圖7所示，參與反演井及后期驗(yàn)證井的泥質(zhì)含量與AVF統(tǒng)計(jì)砂體厚度擬合整體呈線性特征，進(jìn)一步證實(shí)了反演方法的有效性。

圖7 泥質(zhì)含量與AVF反演統(tǒng)計(jì)砂體厚度擬合

在AVF反演體的基礎(chǔ)上， 通過(guò)帶通濾波運(yùn)算，得到常規(guī)振幅地震數(shù)據(jù)體，再運(yùn)用出油氣井對(duì)其進(jìn)行標(biāo)定追蹤，在本區(qū)舒善河組中段共計(jì)標(biāo)定出三套儲(chǔ)蓋組合。 其中Oil1圈閉和Oil3圈閉為主要出油氣層段，通過(guò)追蹤兩者的層面，刻畫出有利單砂體的分布特征，將其疊合在一起（圖8），可為下一步勘探提供重要依據(jù)。 結(jié)合波形指示反演以及“甜點(diǎn)”預(yù)測(cè)結(jié)果，綜合優(yōu)選H點(diǎn)為工區(qū)內(nèi)最為有利的部署位置。

圖8 斷裂及各層有利砂體疊合

5 結(jié)論

（1）AVF分頻屬性反演作為一項(xiàng)全新的技術(shù)手段，目前已被廣泛運(yùn)用于塔北薄砂體儲(chǔ)層的預(yù)測(cè)和識(shí)別中。 通過(guò)地震資料拓頻處理、特征頻率和曲線優(yōu)選建立最優(yōu)井控模型， 應(yīng)用模型約束AVF分頻屬性反演，采用BP神經(jīng)網(wǎng)絡(luò)法對(duì)學(xué)習(xí)曲線和原始曲線進(jìn)行監(jiān)控，從而得到有效的反演結(jié)果。

（2）模型控制AVF分頻屬性反演在托甫臺(tái)北地區(qū)舒善河組薄砂體識(shí)別中的應(yīng)用效果較好。 反演所得的Oil1和Oil3砂體展布特征得到了多口鉆井驗(yàn)證。 為下一步鉆井提供有力的技術(shù)支撐。