朱寶衡，謝春雨

（中國(guó)石化上海海洋油氣分公司勘探開發(fā)研究院，上海 200120）

河流相儲(chǔ)層是我國(guó)東部老油田的主要儲(chǔ)集類型，內(nèi)部構(gòu)型復(fù)雜、儲(chǔ)層厚度薄、巖相變化快、橫向變化劇烈，預(yù)測(cè)識(shí)別難度大。隨著勘探程度的加深，復(fù)雜的巖性油氣藏已經(jīng)成為目前主要勘探目標(biāo)，而對(duì)薄儲(chǔ)層的準(zhǔn)確識(shí)別是巖性油氣藏的一個(gè)技術(shù)難題。傳統(tǒng)的儲(chǔ)層預(yù)測(cè)方法（如常規(guī)地震屬性、時(shí)頻分析等）往往都因?yàn)槠溆邢薜姆直媛?，在識(shí)別薄儲(chǔ)層中發(fā)揮的作用極為有限［1］。地震反演是當(dāng)今儲(chǔ)層預(yù)測(cè)的主要手段。比較成熟的廣義線性反演、遞歸反演、譜反演和約束稀疏脈沖反演可以在一定程度上描述巖性油氣藏的橫向變化，然而反演結(jié)果依然具有很強(qiáng)的多解性，在識(shí)別薄層方面發(fā)揮的作用依然十分有限［2-7］。

為了減少反演結(jié)果的多解性和不穩(wěn)定性，層序界面約束和貝葉斯理論被應(yīng)用于反演［8］。層序界面對(duì)反演的約束作用主要體現(xiàn)在兩個(gè)方面：①反演初始低頻模型的建立；②反演迭代過程的控制。Niepsuj［9］等探究了不同初始模型對(duì)于最終反演結(jié)果的影響；Kemper［10］在地震反演中將不同的層序界面劃分出的地質(zhì)單元賦予不同的初始模型，強(qiáng)調(diào)了層序界面對(duì)于反演初始模型建立的重要性。自從Tarantola［11-12］提出貝葉斯反演以來(lái)，貝葉斯理論被廣泛的應(yīng)用于地球物理反演的目標(biāo)函數(shù)中。Buland［13］、Eidsvik［14］等將貝葉斯理論應(yīng)用到疊前流體因子的反演中。Downtown［15］等通過比較不同的先驗(yàn)分布函數(shù)，發(fā)現(xiàn)與高斯先驗(yàn)分布相比，單變量的柯西分布能夠更好地約束反演過程。Wubshet［16］等利用多參數(shù)柯西分布約束疊前三參數(shù)反演，取得了良好的效果。另外，為了提高識(shí)別薄層的分辨率，地震分頻技術(shù)被廣泛地應(yīng)用于地震資料解釋與儲(chǔ)層預(yù)測(cè)中。龔洪林［17］等利用分頻屬性成功的識(shí)別出了碳酸鹽巖地層中的薄儲(chǔ)層。馬志霞［18］等利用gabor-morlet變換對(duì)地震數(shù)據(jù)進(jìn)行分頻并進(jìn)行儲(chǔ)層預(yù)測(cè)。此外，在反演過程中引入多尺度目標(biāo)函數(shù)并將目標(biāo)函數(shù)的求解放入一系列嵌套的問題中，利用不同尺度目標(biāo)函數(shù)的組合提升尋優(yōu)能力［19-20］。李坤［21］等將平滑模型約束引入貝葉斯多尺度反演過程中，降低了反演的多解性；陶倩倩［22］等利用分頻反演對(duì)潿西南凹陷的烴源巖總有機(jī)碳TOC含量進(jìn)行了預(yù)測(cè)。

本文提出頻率域分頻多尺度反演方法，首先通過對(duì)研究區(qū)進(jìn)行巖石物理統(tǒng)計(jì)，分析目的層的波阻抗特征，以便利用阻抗對(duì)儲(chǔ)層和非儲(chǔ)層進(jìn)行有效的區(qū)分，證明了用反演進(jìn)行儲(chǔ)層預(yù)測(cè)的可行性。其次結(jié)合地震資料頻譜特征的分析結(jié)果，將地震資料分為低、中、高三個(gè)頻段依次進(jìn)行三個(gè)階段的反演，前一階段的反演結(jié)果作為下一階段的初始模型并利用層序界面對(duì)反演過程進(jìn)行約束，通過不同頻率的逐級(jí)迭代方式提升反演結(jié)果的分辨率。

在實(shí)際資料的反演過程中，通過對(duì)勝利油田三合村凹陷的三維地震數(shù)據(jù)進(jìn)行頻譜分析確定主頻并劃定相應(yīng)高、中、低反演頻帶的頻帶范圍。最后，在層序界面的控制下，通過引入貝葉斯框架，假設(shè)反射系概率似然函數(shù)滿足高斯分布，地下反射系數(shù)序列滿足柯西分布，從而降低反演的多解性。對(duì)勝利油田三合村凹陷的三維地震數(shù)據(jù)進(jìn)行反演，對(duì)薄儲(chǔ)層具有好的識(shí)別效果。

1 層序界面約束的分頻多尺度反演方法

在分頻多尺度反演過程中，當(dāng)似然函數(shù)滿足高斯分布，可引入貝葉斯理論來(lái)降低反演的多解性。另外，在實(shí)際三維地震資料反演過程中通過層序約束，然后逐道外推，前一階段的反演結(jié)果作為下一階段的初始模型，以此逐級(jí)迭代來(lái)提高反演的分辨率。

1.1 頻率域反演原理

忽略投射損失、球面擴(kuò)散以及多次波，地震道可以表示為：

其中，ω(t)為地震子波，r(t)時(shí)間域?yàn)榉瓷湎禂?shù)序列，n(t)為非反射噪聲。對(duì)式（1）兩端做傅里葉變換可得：

其中，ω為地震道頻率，Hz；S(ω)為地震道頻譜，R(ω)為反射系數(shù)序列的頻譜，N(ω)為噪聲頻譜。假設(shè)地震資料為平穩(wěn)信號(hào)(Q=∞)，可將式（2）的前一項(xiàng)寫成矩陣形式：

其中，fi（i=1，2，L，M）為地震資料頻率，M為分頻反演時(shí)參與反演的數(shù)據(jù)體個(gè)數(shù)，個(gè)。因此，式（2）可以寫為：

其中，m為待反演參數(shù)向量，S為地震道傅里葉變換序列。

1.2 貝葉斯反演理論

根據(jù)貝葉斯反演理論，聯(lián)合式（5）和（6）可以得到目標(biāo)函數(shù)：

對(duì)式（7）作進(jìn)一步的優(yōu)化處理，可以得到最終的反演方程為：

其中m為待反演阻抗向量，g/cm3·m/s；S為地震道傅里葉變換序列：

分頻反演的具體過程如圖1，將地震數(shù)據(jù)分為7～20 Hz、7～30 Hz、7～45 Hz三個(gè)數(shù)據(jù)體分別進(jìn)行反演。

圖1 分頻反演流程

1.3 層序界面對(duì)反演外推的控制作用

根據(jù)層序界面建立用于反演約束的初始低頻模型，在對(duì)每個(gè)地震數(shù)據(jù)體的反演過程中，采用逐道外推的方式進(jìn)行反演，外推的過程從井旁道開始并逐步覆蓋整個(gè)研究區(qū)。對(duì)于井旁道反演的初始模型由測(cè)井曲線的低頻分量獲得；對(duì)于非井旁道，前一道的反演結(jié)果作為該道反演的初始模型，整個(gè)外推的過程以三維面元的形式進(jìn)行。外推反演的優(yōu)勢(shì)在于融合了測(cè)井、地質(zhì)以及地震波形變化的信息，減少了反演結(jié)果的多解性。其次層序界面在反演過程中根據(jù)低頻背景約束地震道外推反演中每個(gè)采樣點(diǎn)的起時(shí)位置，外推過程中層序界面還約束反演時(shí)窗的范圍。

2 實(shí)例應(yīng)用

2.1 三合村凹陷地質(zhì)背景

研究區(qū)位于勝利油田三合村凹陷，是渤海灣盆地沾化凹陷的一個(gè)次級(jí)單元，緊鄰渤南凹陷、孤南凹陷和富林凹陷（圖2）。渤海灣盆地在構(gòu)造形成過程中經(jīng)歷過抬升、下降、“斷-拗”沉積發(fā)育階段，在郯廬斷裂拉張和盆地本身南北向伸展的共同作用下最終形成“兩坳一隆”的基本構(gòu)造格局，在形成時(shí)期受強(qiáng)烈的擠壓作用以及塊斷隆升活動(dòng)的影響，導(dǎo)致擠壓褶皺、逆沖斷層和負(fù)反轉(zhuǎn)構(gòu)造發(fā)育，其中北部發(fā)育一系列階梯狀的具有東西走向的斷層，對(duì)古近系、新近系儲(chǔ)層的沉積具有控制作用。

圖2 三合村凹陷局部構(gòu)造

本次研究的目的層位于三合村凹陷中新世館陶組館上段，埋藏深度在2 800～3 000 m，地層起伏平緩，地層傾角小于5°，主要發(fā)育曲流河相沉積，河道砂體橫向上變化較快，縱向上相互疊合，斷層為油氣的運(yùn)移提供通道，油氣通過主要的油源斷層孤南斷裂由北向南運(yùn)移，主要物源由工區(qū)北部向南部推進(jìn)，受斷層活動(dòng)的影響，砂體呈塊狀、條狀分布［26］。

2.2 研究區(qū)成藏特征概述

研究區(qū)內(nèi)目的層本身埋深較淺，本地?zé)N源巖發(fā)育但不成熟，油源來(lái)自渤南凹陷沙三段三角洲及沙四段湖相沉積體系中發(fā)育的烴源巖。明化鎮(zhèn)組底部的泥巖為儲(chǔ)集體提供良好的蓋層。東營(yíng)組和沙一段之間的不整合面以及研究區(qū)北部發(fā)育的斷層為油源提供了運(yùn)移的通道，油氣沿不整合面和北部的斷層運(yùn)移至館上段。由于埋藏深度較淺，儲(chǔ)層受成巖作用影響較弱，砂體的孔滲條件較好。館陶組平均孔隙度為31.3%，平均滲透率為400×10-3μm2，屬于高孔隙度-高滲透性油藏。但是館上段儲(chǔ)層油水關(guān)系復(fù)雜，砂體連通性較差，油藏大多受砂體限制，延伸范圍較小。

利用三合村地區(qū)采集的地震資料，對(duì)主要目的層段進(jìn)行頻譜分析（圖3），可以看到地震主頻為20 Hz，頻帶寬度為5～45 Hz，如此我們可以將地震數(shù)據(jù)分為7～20 Hz、7～30 Hz、7～45 Hz三個(gè)數(shù)據(jù)體分別進(jìn)行反演。

圖3 地震資料頻譜特征

圖4為儲(chǔ)層、非儲(chǔ)層縱波阻抗-密度交匯圖，可以看出儲(chǔ)層、非儲(chǔ)層的阻抗具有較好的區(qū)分性，儲(chǔ)層具有高阻抗特性。

2.3 地震層序界面解釋

研究區(qū)目的層處在盆地發(fā)育晚期，基底的沉降速率明顯減弱，可容性空間降低，河道疊置程度加大導(dǎo)致砂體側(cè)向連續(xù)性變差，在地震剖面上會(huì)產(chǎn)生不連續(xù)的反射。通過對(duì)館上段地震反射結(jié)構(gòu)、地震反射構(gòu)型以及地震反射外部形態(tài)的分析，劃分出兩個(gè)三級(jí)層序界面：館上段頂、底界面。層序界面解釋的精細(xì)程度會(huì)在一定程度上影響反演結(jié)果的分辨率，因此在兩個(gè)三級(jí)層序界面之間劃分出四級(jí)層序界面作為反演的約束界面（圖5），增加層序界面約束可靠性，提高反演分辨率。

圖5 Line156過井層序界面分布情況

3 分頻反演結(jié)果分析

在分頻過程中，基于研究區(qū)地震資料頻譜分析結(jié)果，將地震數(shù)據(jù)分頻得到7～20 Hz、7～30 Hz和7～45 Hz三個(gè)頻帶的數(shù)據(jù)進(jìn)行反演。圖6a、6b、6c為不同頻帶的地震數(shù)據(jù)依次反演結(jié)果的過井剖面，可以看出低頻分量反演結(jié)果描繪了地層大致的起伏特征，儲(chǔ)層段的大致分布位置，但是對(duì)于薄層刻畫不明顯（圖6a）。隨著參與反演的地震資料高頻分量的增加，反演結(jié)果對(duì)于薄儲(chǔ)層的刻畫能力明顯增強(qiáng)（圖6b、6c）。

圖6 分頻反演連井剖面

從連井反演剖面中可以看出，目的層段內(nèi)主要發(fā)育砂泥巖薄互層，受河道砂體疊置的影響，砂體的橫向連續(xù)性較差，儲(chǔ)層橫向延伸范圍較小。反演剖面清楚地描述出了井間砂體的連通情況以及尖滅特征。

基于分頻反演成果，通過提取主要目的層段的屬性切片，獲得砂體在平面的形態(tài)展布，并將提取的屬性切片與研究區(qū)構(gòu)造圖疊合（圖7），可以看出砂體形態(tài)在平面上呈細(xì)小條帶狀，發(fā)育三角洲前緣水下分支河道，西北物源，河道向東南方延展，砂體厚度自北向南逐漸變薄。

圖7 分頻反演砂體平面屬性與構(gòu)造疊合

利用分頻反演的結(jié)果對(duì)三合村地區(qū)的河道形態(tài)輪廓進(jìn)行刻畫（圖8），反演結(jié)果與已鉆井結(jié)果吻合（K71-48、K71-28井處于分支河道主體部位，K71-3井位于河道間）。綜合反演剖面和儲(chǔ)層河道砂體在平面上的展布形態(tài)、規(guī)模大小等可為下一步高效勘探提供依據(jù)。

圖8 利用反演結(jié)果對(duì)墾西地區(qū)三合村工區(qū)河道刻畫

4 結(jié)論

（1）頻率域分頻多尺度反演方法，保證了低頻信息準(zhǔn)確性的基礎(chǔ)上進(jìn)一步提高反演分辨率。

（2）針對(duì)墾西地區(qū)三合村凹陷河流相儲(chǔ)層橫向變化快，砂體厚度薄的特點(diǎn)，利用頻率域分頻多尺度反演方法進(jìn)行儲(chǔ)層預(yù)測(cè)，取得了比較好的應(yīng)用效果。

（3）該方法能對(duì)橫向變化快的薄互層砂體有效地識(shí)別，提升了疊后反演對(duì)薄儲(chǔ)層的預(yù)測(cè)能力，清晰刻畫出河道砂體邊界，預(yù)測(cè)砂體平面展布，對(duì)于同樣地質(zhì)背景探區(qū)具有一定的參考意義。